backsolve... --> class()es for special matrices ?

Martin Maechler Martin Maechler <maechler@stat.math.ethz.ch>
Fri, 22 Jan 1999 09:08:50 +0100


>>>>> "JonR" == Jonathan Rougier <J.C.Rougier@durham.ac.uk> writes:

    JonR> ...  By the way, I have `solve'
    JonR> methods for triangular matrices and variance matrices -- would
    JonR> you be interested?

	{ Jonathan, I hope it's okay if I CC this to R-devel;  
	  this must be of a wider interest }

Ye.e..s;  
for triangular ones, there's the 4 different ways of 
backsolve()/forwardsolve()..

Well, we (R core) have  discussed matrix classes months and years ago..
and haven't really started doing anything about it.

I think we then agreed to follow Splus'  Matrix library
and use class "Matrix"
(and classes inheriting from "Matrix", i.e.
	class == c("specialMatrix", "Matrix")
).

If we really go into this, then 
we'd probably want to use the SAME class names that S-plus' library(Matrix)
is using. 
Also, I remember that people (on the old R-testers list?) 
very much desired to have 1-column or 1-row matrices instead of vectors
when working within the "Matrix" class  [= behavior of  Octave/Matlab ?]

Other opinions / (code) contributions ?

Martin Maechler <maechler@stat.math.ethz.ch>	http://stat.ethz.ch/~maechler/
Seminar fuer Statistik, ETH-Zentrum SOL G1;	Sonneggstr.33
ETH (Federal Inst. Technology)	8092 Zurich	SWITZERLAND
phone: x-41-1-632-3408		fax: ...-1086			<><
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._