Date: Thu, 20 Feb 1997 10:23:53 -0600 (CST)
From: "Jesse W. Bennett" <jesse AT lenny DOT dseg DOT ti DOT com>
Reply-To: Jesse Bennett <jbennett AT ti DOT com>
To: Eli Zaretskii <eliz AT is DOT elta DOT co DOT il>
cc: Jesse Bennett <jbennett AT ti DOT com>, kagel AT dg1 DOT bloomberg DOT com,
        djgpp AT delorie DOT com
Subject: Re: Netlib code [was Re: flops...]
In-Reply-To: <Pine.SUN.3.91.970220091557.25071N-100000@is>
Message-ID: <Pine.LNX.3.91.970220101230.25706A-100000@lenny.dseg.ti.com>
MIME-Version: 1.0
Content-Type: TEXT/PLAIN; charset=US-ASCII

On Thu, 20 Feb 1997, Eli Zaretskii wrote: 
> On Wed, 19 Feb 1997, Jesse W. Bennett wrote:
> 
> > to solve X' = [inv(A)B]' = B'inv(A').  I see no easy way to do this
> > without either transposing the C arrays or rewriting a large amount of
> > code (actually I could solve X = inv(A')B' and transpose the X array to
> 
> Can't you edit the source and switch the indices where it accesses the 
> arrays?

Yes, but this occurs in *many* places since this is a very large and
complex piece of code.  The cited example was representative; the problem
is not as localized as the example suggests.  Many of the affected
routines are general library functions I have written over the years and
changing them will break other code that is otherwise unrelated.  If the
problem was isolated to a single source file (or even just a few) this
would be the obvious solution. 

Best Regards,
Jesse