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
53 3 0 0
1 1 -1 1
2 1 1 1
3 -1 1 1
4 -1 -1 1
5 1 -1 -1
6 1 1 -1
7 -1 1 -1
8 -1 -1 -1
9 -0.214 0.214 -0.214
10 -0.214 -0.214 -0.214
11 0.214 -0.214 -0.214
12 0.214 0.214 -0.214
13 -0.214 0.214 0.214
14 -0.214 -0.214 0.214
15 0.214 -0.214 0.214
16 0.214 0.214 0.214
17 -0.037350531086211003 -0.037350531086211003 -1
18 0 0 0.214
19 0 0.214 0
20 0.16391322700620581 -1 -0.16391322700620581
21 0 0 1
22 -1 0 0
23 -1 0 -1
24 1 -1 0
25 0 -1 -1
26 0 1 0
27 1 0 0
28 0 0.214 -0.214
29 0 0 -0.214
30 0.214 0 0.214
31 0.214 0 0
32 0.214 0.214 0
33 0 0.214 0.214
34 -0.214 0.214 0
35 0 -0.214 0.214
36 0 -0.214 0
37 -0.214 0 0.214
38 -0.214 0 0
39 -1 0 1
40 0 -1 1
41 0 1 1
42 1 0 1
43 -1 1 0
44 -0.214 0 -0.214
45 0 -0.214 -0.214
46 0.214 0 -0.214
47 0.214 -0.214 0
48 -1 -1 0
49 1 0 -1
50 1 1 0
51 0 1 -1
52 -0.214 -0.214 0
53 0 0 0
# Generated by tetgen -pgaA trycube.poly
192 4 1
1 10 23 17 8 2
2 47 24 15 20 2
3 51 28 9 17 2
4 32 26 12 19 2
5 10 23 8 22 2
6 30 21 42 16 2
7 16 27 42 30 2
8 47 27 15 24 2
9 35 18 14 21 2
10 32 27 31 12 2
11 14 39 4 21 2
12 28 29 17 12 2
13 1 27 24 15 2
14 12 28 19 26 2
15 16 27 32 50 2
16 16 33 26 19 2
17 19 34 13 26 2
18 25 45 10 20 2
19 20 24 15 40 2
20 16 26 33 41 2
21 45 25 10 17 2
22 40 14 4 21 2
23 30 16 18 21 2
24 34 22 38 13 2
25 23 44 9 22 2
26 12 49 6 27 2
27 44 23 9 17 2
28 15 27 30 42 2
29 14 48 4 22 2
30 22 9 38 44 2
31 5 11 20 24 2
32 9 43 26 7 2
33 14 22 39 37 2
34 13 22 37 39 2
35 39 37 13 21 2
36 10 29 45 17 2
37 48 14 52 22 2
38 13 26 41 33 2
39 14 48 20 4 2
40 28 51 9 26 2
41 30 27 31 16 2
42 29 12 46 17 2
43 49 46 11 27 2
44 12 27 50 32 2
45 16 33 18 21 2
46 18 13 37 21 2
47 36 10 45 20 2
48 14 39 22 4 2
49 28 51 26 12 2
50 11 25 5 20 2
51 9 26 19 28 2
52 16 50 2 27 2
53 40 14 20 4 2
54 9 17 23 51 2
55 40 1 15 21 2
56 9 51 7 26 2
57 32 26 50 12 2
58 34 13 43 22 2
59 13 41 26 3 2
60 35 14 36 20 2
61 9 23 22 7 2
62 46 49 17 12 2
63 11 5 27 24 2
64 10 25 8 17 2
65 11 5 25 49 2
66 10 8 25 48 2
67 13 26 34 43 2
68 9 7 51 23 2
69 9 26 43 34 2
70 41 33 16 21 2
71 12 50 27 6 2
72 11 49 27 5 2
73 16 42 27 2 2
74 9 17 28 29 2
75 9 44 17 29 2
76 30 27 15 31 2
77 15 42 30 21 2
78 15 18 21 30 2
79 46 27 12 31 2
80 15 31 27 47 2
81 16 27 31 32 2
82 16 50 32 26 2
83 16 19 26 32 2
84 13 21 18 33 2
85 13 33 41 21 2
86 13 19 26 33 2
87 13 22 38 37 2
88 9 34 43 22 2
89 9 19 26 34 2
90 15 21 18 35 2
91 15 40 35 20 2
92 52 14 20 36 2
93 24 1 15 40 2
94 15 40 21 35 2
95 15 20 35 36 2
96 15 47 20 36 2
97 14 21 18 37 2
98 14 37 39 21 2
99 34 38 22 9 2
100 13 39 3 22 2
101 14 22 37 38 2
102 14 52 22 38 2
103 13 39 21 3 2
104 40 14 21 35 2
105 40 20 14 35 2
106 16 41 2 26 2
107 13 41 3 21 2
108 16 41 21 2 2
109 15 42 1 27 2
110 16 42 2 21 2
111 15 42 21 1 2
112 13 43 3 26 2
113 9 43 7 22 2
114 13 43 22 3 2
115 10 22 38 44 2
116 10 23 22 44 2
117 10 17 23 44 2
118 10 29 17 44 2
119 11 20 36 45 2
120 11 25 20 45 2
121 11 17 25 45 2
122 11 29 17 45 2
123 11 17 29 46 2
124 46 49 11 17 2
125 11 31 27 46 2
126 11 27 31 47 2
127 11 24 27 47 2
128 11 20 24 47 2
129 11 36 20 47 2
130 10 20 48 25 2
131 14 48 52 20 2
132 10 48 22 8 2
133 11 17 49 25 2
134 12 49 17 6 2
135 46 49 12 27 2
136 12 50 6 26 2
137 16 50 26 2 2
138 12 51 6 17 2
139 12 51 26 6 2
140 28 51 12 17 2
141 10 20 36 52 2
142 10 48 20 52 2
143 10 22 48 52 2
144 10 38 22 52 2
145 36 35 15 53 1
146 15 35 18 53 1
147 30 15 18 53 1
148 30 31 15 53 1
149 15 31 47 53 1
150 47 36 15 53 1
151 11 45 36 53 1
152 29 45 11 53 1
153 11 47 31 53 1
154 36 47 11 53 1
155 11 46 29 53 1
156 31 46 11 53 1
157 36 14 35 53 1
158 14 18 35 53 1
159 31 30 16 53 1
160 16 30 18 53 1
161 38 37 14 53 1
162 14 37 18 53 1
163 14 36 52 53 1
164 52 38 14 53 1
165 10 44 38 53 1
166 29 44 10 53 1
167 10 52 36 53 1
168 38 52 10 53 1
169 36 45 10 53 1
170 10 45 29 53 1
171 16 18 33 53 1
172 33 19 16 53 1
173 16 32 31 53 1
174 19 32 16 53 1
175 31 12 46 53 1
176 12 29 46 53 1
177 31 32 12 53 1
178 32 19 12 53 1
179 28 12 19 53 1
180 29 12 28 53 1
181 19 34 9 53 1
182 9 34 38 53 1
183 38 44 9 53 1
184 9 44 29 53 1
185 9 28 19 53 1
186 9 29 28 53 1
187 38 34 13 53 1
188 13 34 19 53 1
189 37 38 13 53 1
190 13 18 37 53 1
191 13 33 18 53 1
192 19 33 13 53 1
# Generated by tetgen -pgaA trycube.poly
------------------------------------------------------------------------------
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
