Hi R,
I tried debugging sfepy/terms/extmods/termsNavierStokes.c but could not locate the error.
Can you tell me when does sfepy raises a ccore error (I saw it in some .pyx and .c files)
Can you kindly send me the .vtk file generated by the script so that meanwhile I can run an analysis and add it to my mid term report.
Regards,
Ankit
On Friday, 26 July 2013 20:29:27 UTC+5:30, Robert Cimrman wrote:
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.
>>>
>>>
>