Re: [Getfem-users] please help with VTK deformed mesh export

[julien pommier]

Hi Andriy,
this was a non trivial bug. I just commited revision 2508 to fix it:
http://svn.gna.org/viewcvs/getfem?rev=2508&view=rev

however if you want a quick and ugly fix, just look for the clear() 
method of fem_precomp_pool in getfem_fem.h and remove the call to 
"del_stored_object"

--
Julien

Andriy Andreykiv wrote:
> Hi,
>
>         I'm trying to export a deformed mesh in a VTK format.
> For the start I was using elastostatic.cc from the tests directory of Getfem.
> I put the the following code in the main() of elastostatic.cc:
>
>       getfem::vtk_export exp("elastostatic.vtk", true);
>       getfem::mesh_slice_cv_dof_data<plain_vector> base(p.mf_u, U);
>       getfem::slicer_apply_deformation a1(base);
>       getfem::stored_mesh_slice sl;
>       sl.build(p.mf_u.linked_mesh(),a1,1);
>       exp.exporting(sl);
>       exp.write_point_data(p.mf_u, U, "elastostatic_displacement");
>
> The build goes fine, but the execution gives an error.
> The error I get is: Error in dal_static_stored_objects.cc, line 199 :
> This object is not stored : 0x9dcdef8
> Obviously, sl.build doesn't work in this case, but I don't understand why.
>
> I was also trying the first method from the two, explained in the manual, by 
> using getfem::slicer objects (see the code bellow), but the result is the 
> same.
>
>       getfem::vtk_export exp("elastostatic.vtk", true);
>       getfem::mesh_slice_cv_dof_data<plain_vector> base(p.mf_u, U);
>       getfem::slicer_apply_deformation a1(base);
>       getfem::stored_mesh_slice sl;
>       getfem::slicer_build_stored_mesh_slice a2(sl);
>       getfem::mesh_slicer slicer(p.mf_u.linked_mesh( ));
>       slicer.push_back_action(a1);
>       slicer.push_back_action(a2);
>       slicer.exec(1);
>       exp.exporting(sl);
>       exp.write_point_data(p.mf_u, U, "elastostatic_displacement");     
>  
> If you have any guess how to fix this I would really appreciate your help.
>
> Best regards,
>                                 Andriy
>
> P.S All the attempts to export the results without the slice_apply_deformation 
> class are successfull.
>       
>
> __________________________________________
> Andriy Andreykiv
> Delft University of Technology
> Faculty of Mechanical Engineering 
> Fundamentals of Microsystems
> Mekelweg 2
> 2628 CD  Delft
> The Netherlands
>
> E-mail  :  address@hidden
> Tel.    :   (31) 15 2786818
> Fax.    :   (31) 15 2789475
> www     :   http://www-tm.wbmt.tudelft.nl/~andrico
> private:  (31) 6 47376804
>
> _______________________________________________
> Getfem-users mailing list
> address@hidden
> https://mail.gna.org/listinfo/getfem-users
>   


--
Julien Pommier, bureau 111
GMM, INSA Toulouse, tél:05 61 55 93 42