[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
solving Ax=b
From: |
Sven Khatri |
Subject: |
solving Ax=b |
Date: |
Sat, 4 Jan 2003 01:00:26 -0800 |
Hi All,
I'm looking to solve a problem that is equivalent to solving for x in
the system of linear equations:
Ax = b
where A \in R^{n \times n}, x \in R^n, b \in R^n. But the catch is
that n is too large to explicitly construct A (even in Matlab's sparse
implementation) but I can construct a function f so that f(x) = Ax. Is
there an easy way (within octave) to solve for x? I can solve the problem
by introducing \tilde{A} = I - A, so that the solution to
x = \tilde{A}x + b
is a solution to the above problem and \tilde{A} is a contraction and
solving the problem iteratively but this computation is VERY slow.
thanx...
Sven
PS I hope you all are comfortable with the above latex notation
-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.
Octave's home on the web: http://www.octave.org
How to fund new projects: http://www.octave.org/funding.html
Subscription information: http://www.octave.org/archive.html
-------------------------------------------------------------
- solving Ax=b,
Sven Khatri <=