Implement MatSparse<GPU2>::multiply(MatDense<GPU2>), Matdense<GPU2>::multiplyLeft(MatSparse<GPU2>).

src/faust_linear_operator/GPU2/faust_MatDense_gpu.cpp.in

@@ -481,6 +481,12 @@ namespace Faust
 			}
 		}
 
+		template<>
+		void Faust::MatDense<@FAUST_SCALAR_FOR_GM@,GPU2>::multiplyLeft(MatSparse<@FAUST_SCALAR_FOR_GM@,GPU2> const& S, const char transS)
+		{
+			S.multiply(*this);
+		}
+
 	template<>
 		void Faust::MatDense<@FAUST_SCALAR_FOR_GM@, GPU2>::operator=(const MatDense<@FAUST_SCALAR_FOR_GM@, Cpu>& A)
 		{

src/faust_linear_operator/GPU2/faust_MatDense_gpu.h

@@ -79,6 +79,7 @@ namespace Faust
 				void multiplyRight(const MatDense<FPP, GPU2>& A);
 				//! \brief Replace (this) by S * (this)
 				void multiplyLeft(const MatSparse<FPP, Cpu>& S, const char transS='N');
+				void multiplyLeft(const MatSparse<FPP, GPU2>& S, const char transS='N');
 				void multiply(const Vect<FPP, GPU2>& vec, Vect<FPP, GPU2>& out_vec) const;
 				//! \brief compute MatDense-vector multiplication
 				//! \param vec : the vector

src/faust_linear_operator/GPU2/faust_MatSparse_gpu.cpp.in

@@ -486,6 +486,18 @@ namespace Faust
 			vec.gpu_mat = gpu_dmat;
 		}
 
+	template<>
+		void MatSparse<FSFG, GPU2>::multiply(MatDense<FSFG,GPU2>& mat, char opThis/*='N'*/) const
+		{
+			gm_Op gop_this;
+			char2gm_Op(opThis, gop_this);
+			auto spm_funcs = GPUModHandler::get_singleton()->spm_funcs((@FAUST_SCALAR_FOR_GM@)(0));
+			auto dsm_funcs = GPUModHandler::get_singleton()->dsm_funcs((@FAUST_SCALAR_FOR_GM@)(0));
+			auto gpu_dmat = spm_funcs->mul_gpu_dsm_ext(gpu_mat, mat.gpu_mat, nullptr, gop_this, OP_NOTRANSP);
+			dsm_funcs->free(mat.gpu_mat);
+			mat.gpu_mat = gpu_dmat;
+		}
+
 	template<>
 		size_t Faust::MatSparse<FSFG,GPU2>::getNBytes() const
 		{

src/faust_linear_operator/GPU2/faust_MatSparse_gpu.h

@@ -92,6 +92,7 @@ namespace Faust
 				std::string to_string(const bool transpose=false, const bool displaying_small_mat_elts=false) const;
 				MatType getType() const;
 				void multiply(Vect<FPP,GPU2>& vec, char opThis='N') const;
+				void multiply(MatDense<FPP,GPU2>& vec, char opThis='N') const;
 				~MatSparse();
 
 			MatSparse<FPP,GPU2>* get_rows(faust_unsigned_int row_id_start, faust_unsigned_int num_rows) const;