Path: utzoo!utgpu!news-server.csri.toronto.edu!bonnie.concordia.ca!ccu.umanitoba.ca!herald.usask.ca!alberta!brazeau.ucs.ualberta.ca!unixg.ubc.ca!ubc-cs!uw-beaver!cornell!rochester!pt.cs.cmu.edu!o.gp.cs.cmu.edu!fs7.ece.cmu.edu!faraday.ece.cmu.edu!winstead
From: winstead@faraday.ece.cmu.edu (Charles Holden Winstead)
Newsgroups: comp.lang.fortran
Subject: Re: Need complex matrix solver for CM
Message-ID: <1991Jun28.061454.29461@fs7.ece.cmu.edu>
Date: 28 Jun 91 06:14:54 GMT
References: <1991Jun27.161918.10218@leland.Stanford.EDU> <1991Jun27.190724.26161@ariel.unm.edu> <7887@mace.cc.purdue.edu>
Sender: news@fs7.ece.cmu.edu (USENET News System)
Organization: Electrical and Computer Engineering, Carnegie Mellon
Lines: 21

>In article <1991Jun27.190724.26161@ariel.unm.edu> john@spectre.unm.edu (John Prentice) writes:
>linear system twice as large. If you need a High School math reminder:
                                             ^^^^^^^^^^^^^^^^^^^^^^^^^^

Yo dude, grow up!  I agree that it's a simple task of changing a complex 
linear system into a real one, but with a language that deals with complex
numbers as nicely as fortran does, let's agree that there are some benefits
to not switching from complex to real and back again repeatedly.  Why not
just look through the code that works for real matrices and change the 
corresponding variable declaration from real to complex.

>
>Ax=b ==>
>
>+ A(r)  -A(i) +  x(r)          b(r)
>|             |           =    
>+ A(i)   A(r) +  x(i)          b(i)

Gee thanks.

-Chuck