Path: utzoo!utgpu!news-server.csri.toronto.edu!cs.utexas.edu!qt.cs.utexas.edu!yale.edu!think.com!spool.mu.edu!uunet!stanford.edu!leland.Stanford.EDU!hanauma!francis
From: francis@hanauma.stanford.edu (Francis Muir)
Newsgroups: comp.lang.fortran
Subject: Re: Need complex matrix solver for CM
Message-ID: <1991Jun27.201955.4245@leland.Stanford.EDU>
Date: 27 Jun 91 20:19:55 GMT
Sender: news@leland.Stanford.EDU (Mr News)
Organization: Stanford University, Department of Geophysics
Lines: 31

John Prentice writes:

  Francis Muir writes:

    John Prentice writes:

      I need a routine to solve a dense complex valued linear system on the
      CM2, probably using a LU decomposition and Gaussian elimination.  From
      what I know of it, CMSSL only has software like this for real valued
      matrices (please correct me if I am wrong), while I need complex valued
      ones.  Any help would be most appreciated.

    Re-write the code in REAL; it will probably illuminate.
						
  This is an amazingly presumptuous, not to mention useless, comment.  If I
  have to write it in REAL, I will.  But to do so completely obscures the
  elegant structure of the matrices I am solving and requires both artificial
  and inelegant constructs.  I would vastly prefer to have software to handle
  the complex matrices as they are and keep my code simple.

?

If you were to look at:

	CMSSL for CM FORTRAN, Version 2.1, March 1991

you would see in Chapter 2, Linear Algebra Routines, that all routines handle
COMPLEX as well as REAL. For dense systems CMSSL offers both Gauss-Jordan and
QR.

							Fido