On 26 May 2010, at 09:10, Jaroslav Hajek wrote: > > OK. Please note that your C++ version only works in real doubles, > whereas a properly written m-file function will probably work well in > both single and double (that is typically automatic, one only needs to > be careful about eps() and tolerances) and both real and complex (that > is sometimes harder, requires using ' and .' properly, sometimes it is > not possible).
Jaroslav, I just commited an pgmres.m and mgorth.[m,oct] to linear-algebra. pgmres.m seems to be almost as fast as pgmres.oct, as long as mgorth.oct is used but is slower if mgorth.m is used. On the other hand the oct file version of mgorth suffers from the problems you noted above. What do you suggest? c. ------------------------------------------------------------------------------ _______________________________________________ Octave-dev mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/octave-dev
