bug (resolu) transpose matlab mais (A'*x) lent par rapport a (trans_A=A';trans_A*x

misc/demo/tools/greed_omp_chol.m

@@ -217,7 +217,7 @@ end
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
 if          isa(A,'float')      P =@(z) A*z;  Pt =@(z) A'*z;
-elseif      isobject(A)         P =@(z) A*z;  Pt =@(z) A'*z;
+elseif      isobject(A)         P =@(z) A*z; trans_A=A'; Pt =@(z) trans_A*z;%% modification nicolas bellot
 elseif      isa(A,'function_handle') 
     try
         if          isa(Pt,'function_handle'); P=A;

wrapper/matlab/src/mexFaust.cpp

@@ -52,9 +52,7 @@ void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
         
 		std::vector<Faust::MatSparse<FFPP,Cpu> > vec_spmat;
 		mwSize nb_element = mxGetNumberOfElements(prhs[1]);
-		if (nb_element == 0)
-			mexWarnMsgTxt("Empty cell array.");
-    		else
+		if (nb_element != 0)
 		{
 			mxArray * mxMat;
 			for (int i=0;i<nb_element;i++)
@@ -132,13 +130,9 @@ void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
 	{
 
 
-		// if (nlhs != 1)
-		// {
-			// (*core_ptr).transpose();
-		// }
-		// return;
-		if (nlhs != 1)
-			(*core_ptr).transpose();
+
+		if (nlhs > 1)
+			mexErrMsgTxt("transpose : too much output argument");
 		else
 		{
 			Faust::Transform<FFPP,Cpu>* F = new Faust::Transform<FFPP,Cpu>((*core_ptr));