2009/7/20 Robert Cimrman cimr...@ntc.zcu.cz:
Robert Cimrman wrote:
Antonino Ingargiola wrote:
Hi Robert,
2009/7/20 Robert Cimrman cimr...@ntc.zcu.cz:
Hi Antonio,
I am just preparing a new release and think that some of the umfpack related issues have been fixed already in the Git version - could you check that out? You use 2009.2, right?
You do not need Git to try it, just grab [1].
[1] http://git.sympy.org/?p=sfepy.git%3Ba=snapshot%3Bh=HEAD%3Bsf=tgz
I've tried the latest version. Now I don't have the problem on test_quadratures.py anymore. (I've configured site_cfg_template.py to use the local numpy 1.3.0 installation).
However now runTest.py freezes the terminal during the test_input_navier_stokes.py test. The output is appended.
Thanks for your quick answer.
This might be related to [1] (at the bottom of the message) - do you have umfpack installed?
Yes, as mentioned I didn't find UMFPACK anymore in the latest ubuntu. However searching better I found that the library is now included in the libsuitesparse package. I have installed the development package libsuiteparse-dev, then installed your umfpack scikit from svn.
Recompiling sfepy now I have only one error in tests/test_input_acoustics.py. Here it is the output:
$ ./runTests.py --debug tests/test_input_acoustics.py <<< tests/test_input_acoustics.py sfepy: left over: ['__builtins__', '_filename', '__file__', '__package__', 'TestInput', 'input_name', '__name__', 'output_name', '__doc__'] sfepy: left over: ['w', '__builtins__', 'rho', '__name__', 'v_n', 'c', '__doc__', '__package__', '_filename', '__file__']
test instance prepared (1 test(s))
... solving ../input/acoustics.py... sfepy: reading mesh (database/acoustics.mesh)... sfepy: ...done in 0.00 s sfepy: setting up domain edges... sfepy: ...done in 0.00 s sfepy: creating regions... sfepy: leaf Gamma_out region_Gamma_out__0 sfepy: leaf Omega region_Omega__1 sfepy: leaf Gamma_in region_Gamma_in__2 sfepy: ...done in 0.01 s sfepy: equation "Acoustic pressure": sfepy: 117649.0 * dw_laplace.ivol.Omega( one.one, q, p ) - 1000000.0 * dw_mass_scalar.ivol.Omega( q, p ) - 343000j * dw_surface_mass_scalar.isurf.Gamma_out( q, p ) = 182355950j * dw_surface_integrate.isurf.Gamma_in( q ) sfepy: setting up dof connectivities... sfepy: ...done in 0.00 s sfepy: describing geometries... sfepy: ...done in 0.01 s sfepy: using solvers: nls: newton ls: ls sfepy: matrix shape: (61, 61) sfepy: assembling matrix graph... sfepy: ...done in 0.00 s sfepy: matrix structural nonzeros: 445 (1.20e-01% fill) sfepy: updating materials... sfepy: one sfepy: ...done in 0.00 s /home/anto/src/python/lib/python2.6/site-packages/scipy/linsolve/__init__.py:4: DeprecationWarning: scipy.linsolve has moved to scipy.sparse.linalg.dsolve warn('scipy.linsolve has moved to scipy.sparse.linalg.dsolve', DeprecationWarning) sfepy: nls: iter: 0, residual: 1.116698e+07 (rel: 1.000000e+00) warning: singular matrix zeroing nan and inf entries... warning: (almost) singular matrix! (estimated cond. number: inf) sfepy: rezidual: 0.00 [s] sfepy: solve: 0.00 [s] sfepy: matrix: 0.00 [s] sfepy: linear system not solved! (err = 1.116698e+07) sfepy: nls: iter: 1, residual: 1.116698e+07 (rel: 1.000000e+00) ... ../input/acoustics.py solved --- test_input: failed! !!! 1 test failed 1 test file executed in 0.51 s, 1 failure(s) of 1 test(s)
Many thanks for your help,
~ Antonio