[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

sparse matrix improvements in 2.9.x

From: Shai Ayal
Subject: sparse matrix improvements in 2.9.x
Date: Sat, 10 Jun 2006 08:04:11 +0300
User-agent: Thunderbird (Windows/20060516)

Hi all,

I have to invert a sparse matrix using the "\" operator. This matrix is the result of a discretization of a Poisson equation using finite differencing. When I increase the resolution of the discretization above a certain threshold the inversion stops to work. I can currently invert a ~200x200 matrix, but not 300x300 derived from the same physical problem (i.e. I just decrease the cell size). I am almost sure this is not a memory problem (I have 1GB ram). So this is probably a numerical problem. It might be something in my discretization code, or in the inversion code.

I am currently using 2.1.72 with octave-forge. Should I upgrade to 2.9.x? Is the sparse inversion there "better" in that it tries to do some clever decompositions (like matlab)?


reply via email to

[Prev in Thread] Current Thread [Next in Thread]