
Antonino Ingargiola wrote:
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;a=snapshot;h=HEAD;sf=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.
Whoops, I missed that in your e-mail. :)
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.
Good!
Recompiling sfepy now I have only one error in tests/test_input_acoustics.py. Here it is the output:
I guess you needed no recompilation in this case - sfepy uses umfpack only as an external package with no link-time dependence.
$ ./runTests.py --debug tests/test_input_acoustics.py <<< tests/test_input_acoustics.py 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)
This is strange, look at what I get:
sfepy: nls: iter: 0, residual: 1.116698e+07 (rel: 1.000000e+00) sfepy: rezidual: 0.00 [s] sfepy: solve: 0.00 [s] sfepy: matrix: 0.00 [s] sfepy: nls: iter: 1, residual: 8.460917e-08 (rel: 7.576731e-15) ... ../input/acoustics.py solved
- no singular matrix warning. I have tested with numpy 1.2.1, scipy 0.7.0, umfpack 5.2.0 on gentoo, latest git version of sfepy. But I do not use the scikit (good old deprecated scipy.sparse.linalg.dsolve.umfpack). Let me check with the latest scikit
- yes, it works too.
What is you swig version?
Many thanks for your help,
Glad to help, r.