El 18/07/2012, a las 15:05, Jed Brown escribió: > MatMultHermitianTranspose could have a default implementation in terms of the others. Then it would work in all cases, right? > Do you mean VecConjugate(x) plus MatMultTranspose(A,x,y) plus VecConjugate(y) ?? I guess this would work. Jose