I found the problem. GMVIO changes the order of the tetrahedra
vertices:
*.ele file:
192 4 1
1 10 23 17 8 2
2 47 24 15 20 2
3 51 28 9 17 2
*.gmv file:
cells 192
tet 4
10 17 23 8
tet 4
47 15 24 20
tet 4
51 9 28 17
In contrast to this, xdr_io works fine.
Robert
Am Montag, den 08.08.2011, 17:29 +0200 schrieb robert:
> Hello,
>
> I have seen that there has already been a topic like this one some years
> ago. However, I couldn't find a solution for the problem.
>
> At the end of a calculation I want to evaluate the results by using
> system.point_value(int,Point). Everything works fine if I do it
> immediately. However, if I write the equation system and the mesh and
> read it afterwards:
>
> std::cout << "Calculate point value: ";
>
> equation_systems.write("equation_systems");
> mesh.write("mesh.gmv");
> Mesh old_mesh;
> EquationSystems old_equation_systems(old_mesh);
> old_mesh.read("mesh.gmv");
> old_equation_systems.read("equation_systems");
> TransientLinearImplicitSystem & old_system =
> old_equation_systems.get_system<TransientLinearImplicitSystem>
> ("Convection-Diffusion");
> Point point(0.,0.,0.);
> Real value = old_system.point_value(0,point);
>
> std::cout << "\t\t\t\t"<<value << " done" << std::endl;
>
>
> I get the following error message:
>
> ERROR: negative Jacobian: -0.00980034 in element 144
> [0] src/fe/fe_map.C, line 362, compiled Aug 4 2011 at 17:19:33
> terminate called after throwing an instance of 'libMesh::LogicError'
> what(): Error in libMesh internal logic.
>
> I think there might be a problem with accuracy when writing and reading
> the equation system. Is there any better way to evaluate the data?
>
> I have attached a sample mesh and a slightly modified version of
> example 9 which gives me the above error.
>
> Thank you in advance,
> Robert
>
> ------------------------------------------------------------------------------
> BlackBerry® DevCon Americas, Oct. 18-20, San Francisco, CA
> The must-attend event for mobile developers. Connect with experts.
> Get tools for creating Super Apps. See the latest technologies.
> Sessions, hands-on labs, demos & much more. Register early & save!
> http://p.sf.net/sfu/rim-blackberry-1
> _______________________________________________ Libmesh-users mailing list
> [email protected]
> https://lists.sourceforge.net/lists/listinfo/libmesh-users
------------------------------------------------------------------------------
BlackBerry® DevCon Americas, Oct. 18-20, San Francisco, CA
The must-attend event for mobile developers. Connect with experts.
Get tools for creating Super Apps. See the latest technologies.
Sessions, hands-on labs, demos & much more. Register early & save!
http://p.sf.net/sfu/rim-blackberry-1
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
