Dear libmesh users and developers,
I'm having an issue with the contains_point function in the Elem class. I'm
some inconsistencies. Namely, For the same mesh and same points (not nodes)
to test whether the points are in the mesh, I'm getting true sometimes and
false some other time. To provide an example of my results:
(inMesh? means in any of the elements? The .ele file is provided after this
file and the mesh is a sphere of radius 0.4)
PointNum x y z inMesh?
1 0 0 0 1
2 0.01 0.01 0.01 1
3 0.01 -0.01 0.01 0
4 0.01 0.01 -0.01 0
5 0.01 0.01 0.01 0
6 0.01 0.05 0.01 0
7 0.01 -0.05 0.01 0
8 0.02 -0.03 0.04 0
9 0.01 0.08 0.05 0
10 1 1 1 0
Point 2 and 5 are identical, but are giving different results! I'm not using
parallelization here.
The mesh (in tetgen .node and .ele format, should be self-explanatory) is
provided below:
test.1.node
19 3 2 0
1 0.40000000000000002 0 0 1 0
2 0 0.40000000000000002 0 1 0
3 -0.40000000000000002 0 0 1 0
4 0 -0.40000000000000002 0 1 0
5 0 0 -0.40000000000000002 1 0
6 0 0 0.40000000000000002 1 0
7 0.28284271247484399 0.28284271247439408 0 1 0
8 -0.28284271247439408 0.28284271247484399 0 1 0
9 -0.28284271247484399 -0.28284271247439408 0 1 0
10 0.28284271247439408 -0.28284271247484399 0 1 0
11 0 0.28284271247484399 -0.28284271247439408 1 0
12 0 -0.28284271247439408 -0.28284271247484399 1 0
13 0 -0.28284271247484399 0.28284271247439408 1 0
14 0 0.28284271247439408 0.28284271247484399 1 0
15 0.28284271247484399 0 0.28284271247439408 1 0
16 -0.28284271247439408 0 0.28284271247484399 1 0
17 -0.28284271247484399 0 -0.28284271247439408 1 0
18 0.28284271247439408 0 -0.28284271247484399 1 0
19 0 0 0 1 0
# Generated by tetgen test.node
test.1.ele:
32 4 0
1 9 3 19 17
2 9 19 12 17
3 8 11 19 17
4 18 19 12 10
5 10 1 19 15
6 2 19 8 11
7 12 17 19 5
8 16 19 9 3
9 12 19 4 10
10 8 19 16 3
11 18 19 7 11
12 4 12 9 19
13 19 10 18 1
14 13 10 19 15
15 7 19 2 11
16 19 17 11 5
17 19 3 8 17
18 12 19 18 5
19 19 11 18 5
20 7 19 18 1
21 2 7 14 19
22 4 19 13 10
23 16 6 13 19
24 14 6 16 19
25 4 9 13 19
26 13 9 16 19
27 2 14 8 19
28 14 16 8 19
29 13 19 6 15
30 6 19 14 15
31 14 19 7 15
32 19 1 7 15
# Generated by tetgen test.node
Any help will be greatly appreciated.
Thanks,
Karen
------------------------------------------------------------------------------
This SF.net Dev2Dev email is sponsored by:
Show off your parallel programming skills.
Enter the Intel(R) Threading Challenge 2010.
http://p.sf.net/sfu/intel-thread-sfd
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
