On Jul 10, 2008, at 9:58 AM, Rashmi Ranjan Mohanty wrote:
I am running a Fipy code with a system size of 400x400 mesh points
on a Mac Pro 3GHz Dual-Intel Xeon processor with 2GB memory. After
running few steps it stops with the following error message:
python(1013) malloc: *** vm_allocate(size=7921664) failed (error
code=3)
python(1013) malloc: *** error: can't allocate region
python(1013) malloc: *** set a breakpoint in szone_error to debug
malloc fails for local dworkptr[].Traceback (most recent call last):
File "isothermal-new2D.py", line 355, in <module>
compRes = compEq.sweep(var = comp, boundaryConditions = BCsC,
dt=time, solver=solver)
File "/sw/lib/python2.5/site-packages/fipy/terms/term.py", line
153, in sweep
self._solveLinearSystem(var, solver, matrix, RHSvector)
File "/sw/lib/python2.5/site-packages/fipy/terms/term.py", line
102, in _solveLinearSystem
solver._solve(matrix, array, RHSvector)
File "/sw/lib/python2.5/site-packages/fipy/solvers/
linearLUSolver.py", line 92, in _solve
LU = superlu.factorize(L._getMatrix().to_csr())
MemoryError
I guess it has something to do with memory, but could not figure out
a solution.
I would appreciate your help.
My suspicion would lie on your use of the LU solver. The PySparse LU
solver is robust, and often faster than the iterative solvers, but it
does so by using a great deal of memory.
Have you tried other solvers?
