[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
RE: LU decomposition/backsubstitution
From: |
Ted Harding |
Subject: |
RE: LU decomposition/backsubstitution |
Date: |
Mon, 06 Dec 1999 11:24:28 -0000 (GMT) |
On 06-Dec-99 Thomas Hoffmann wrote:
> To solve a linear system A x = b, I use "left division" x=A\b .
>
> But given I have a system with fixed A and several (succesively
> computed) b's, I would like to use something like a precomputed
> decomposition of A which I can apply my b's to, to get the solutions
> for every b with minimal computational effort.
Surely setting
B = A\U
where U is an identity matrix of appropriate dimensions, would give
you what you are looking for (to within rounding errors, anyway):
B*b = A\U*b = A\b
Ted.
--------------------------------------------------------------------
E-Mail: (Ted Harding) <address@hidden>
Date: 06-Dec-99 Time: 11:24:28
2205332 seconds left
------------------------------ XFMail ------------------------------
-----------------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.
Octave's home on the web: http://www.che.wisc.edu/octave/octave.html
How to fund new projects: http://www.che.wisc.edu/octave/funding.html
Subscription information: http://www.che.wisc.edu/octave/archive.html
-----------------------------------------------------------------------