Hi Alec,
On 09/25/2012 08:23 PM, Alec Kalinin wrote:
Hello SfePy,
I have found the reason of error! The mesh file was generated in Gmsh tool and contains the nodes that were not a part of any tetra. So these standalone nodes were the reason of the singular matrix.
It's good you found the reason.
I think, there are two ways to avoid this error:
- Write a simple python script that preproceses mesh file and remove standalone nodes.
- Write a simple python script that preproceses mesh file and add proper group number to the nodes. Then in problem description file we can use the node selection code like 'nodes of group 1'.
It could be made so that 1. is done in sfepy sutomatically, also giving a warning when such a situation occurs. I will try to fix this soon.
Thanks. r.
Sincerely, Alec.
On Tuesday, September 25, 2012 6:11:51 PM UTC+4, Alec Kalinin wrote:
Hello SfePy,
During the problem solution I got the error:
sfepy: nls: iter: 0, residual: 5.088014e+00 (rel: 1.000000e+00) warning: singular matrix zeroing nan and inf entries... /usr/lib/python2.7/dist-packages/scipy/sparse/linalg/dsolve/umfpack/umfpack.py:556: RuntimeWarning: divide by zero encountered in double_scalars econd = 1.0 / self.info[UMFPACK_RCOND] warning: (almost) singular matrix! (estimated cond. number: inf) sfepy: rezidual: 0.01 [s]
I attached the problem description and the mesh for this problem. What is the reason of the error? May be the mesh is not very detailed?
Sincerely, Alexander Kalinin