Hi Ankit,
the failing function is in sfepy/terms/extmods/termsNavierStokes.c - try debugging it (printf...). I cannot reproduce it, sorry :]
r.
On 07/26/2013 03:13 PM, Ankit Mahato wrote:
On my other system convect_build_vtg() is again causing the error. But this time the Traceback is a bit different with mem_free_mem(): error exit!
convdiff: left over: ['os', '__builtins__', '__doc__', '__name__', '__package__', 'verbose', 'sys', 'data_dir', 'cwd', '_filename', 'main', '__file__'] convdiff: reading mesh (/usr/local/lib/python2.7/dist-packages/sfepy/meshes/2d/rectangle_fine_quad.mesh)... convdiff: ...done in 0.03 s convdiff: creating regions... convdiff: Wall convdiff: Top convdiff: Surface convdiff: Driven convdiff: Entry convdiff: Omega convdiff: ...done in 0.05 s convdiff: using solvers: ts: no ts nls: newton ls: ls convdiff: equation "balance": convdiff: + dw_div_grad.5.Omega(m.viscosity, v, u) + dw_convect.5.Omega(v, u) - dw_stokes.5.Omega(v, p) = 0 convdiff: equation "incompressibility": convdiff: dw_stokes.5.Omega(u, q) = 0 convdiff: setting up dof connectivities... convdiff: ...done in 0.00 s convdiff: updating variables... convdiff: ...done convdiff: matrix shape: (44949, 44949) convdiff: assembling matrix graph... convdiff: ...done in 0.24 s convdiff: matrix structural nonzeros: 1747752 (8.65e-04% fill) convdiff: updating materials... convdiff: m convdiff: ...done in 0.02 s convdiff: nls: iter: 0, residual: 2.002082e-02 (rel: 1.000000e+00) convect_build_vtg(): ERR_Switch mem_free_mem(): error exit! mem_free_mem(): error exit! mem_free_mem(): error exit! mem_free_mem(): error exit! mem_free_mem(): error exit! mem_free_mem(): error exit! mem_free_mem(): error exit! mem_free_mem(): error exit! Traceback (most recent call last): File "examples/standalone/convdiff/convective_diffusive.py", line 131, in <module> main() File "examples/standalone/convdiff/convective_diffusive.py", line 113, in main flow = problem.solve() File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/problemDef.py", line 933, in solve vec = solvers.nls(vec0) File "/usr/local/lib/python2.7/dist-packages/sfepy/solvers/nls.py", line 345, in __call__ mtx_a = fun_grad(vec_x) File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/evaluate.py", line 66, in eval_tangent_matrix mtx = pb.equations.eval_tangent_matrices(vec, mtx) File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/equations.py", line 640, in eval_tangent_matrices self.evaluate(mode='weak', dw_mode='matrix', asm_obj=tangent_matrix) File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/equations.py", line 526, in evaluate asm_obj=asm_obj) File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/equations.py", line 766, in evaluate ret_status=True) File "/usr/local/lib/python2.7/dist-packages/sfepy/terms/terms.py", line 1473, in evaluate diff_var, **kwargs) File "/usr/local/lib/python2.7/dist-packages/sfepy/terms/terms.py", line 1311, in eval_real status = self.call_function(out, fargs) File "/usr/local/lib/python2.7/dist-packages/sfepy/terms/terms.py", line 1296, in call_function raise ValueError('term evaluation failed! (%s)' % self.name) ValueError: term evaluation failed! (dw_convect)
On Friday, 26 July 2013 17:42:57 UTC+5:30, Ankit Mahato wrote:
On Friday, 26 July 2013 17:39:52 UTC+5:30, Robert Cimrman wrote:
On 07/26/2013 01:55 PM, Ankit Mahato wrote:
On Friday, 26 July 2013 17:17:09 UTC+5:30, Robert Cimrman wrote:
On 07/26/2013 01:40 PM, Ankit Mahato wrote:
On Friday, 26 July 2013 16:54:45 UTC+5:30, Robert Cimrman wrote: > > Do you have the latest code with Navier-Stokes terms updated for 2D? > > Yes. Also the navier_stokes2D.py is working perfectly fine.
Strange. It works for me. Try "make clean; make"...
yes R. Did it. still throwing the same Runtime error. Other standalone codes are working fine. Is it working on your 32-bit box.
Yes, it works also on the 32 bit box.
I will try it on another system and check it.
r.