Dear Nek Users,
I am trying to get periodicity to work in the lastest version of Nek. I've made an extremely simple mesh within prenek (prex) in polar coordinates and cannot seem to run this case. When I run genmap, I get the error below. Did I do something wrong in setting up this mesh?
>> genmap
reading a.rea
Input mesh tolerance (default 0.2):
NOTE: smaller is better, but generous is more forgiving for bad meshes.
0.01
reading mesh data ...
start locglob_lexico: 4 2 8 1.0000000000000000E-002
locglob: 1 1 8
locglob: 2 5 8
locglob: 1 6 8
locglob: 2 6 8
done locglob_lexico: 6 6 8 2
start periodic vtx: 2 6
1 4 4 Matrix: x0
1 x0 1.0000E+00 1.0000E+00 5.4347-323 4.9407-324
2 x0 -3.5355E-01 3.5355E-01 6.9423-310 6.9423-310
3 x0 3.5355E-01 -3.5355E-01 6.9529-310 6.9529-310
4 x0 6.9423-310 6.9529-310 2.3477-313 3.1632-317
2 4 4 Matrix: x1
1 x1 1.0000E+00 1.0000E+00 6.3240-322 0.0000E+00
2 x1 -3.5355E-01 3.5355E-01 0.0000E+00 6.9423-310
3 x1 -3.5355E-01 3.5355E-01 6.9423-310 6.9529-310
4 x1 6.9423-310 6.9423-310 6.9423-310 6.9529-310
1 4 4 Matrix: z0
1 z0 1.0000E+00 1.0000E+00 5.4347-323 4.9407-324
2 z0 -2.1213E+00 -1.4142E+00 6.9423-310 6.9423-310
3 z0 2.1213E+00 1.4142E+00 6.9529-310 6.9529-310
4 z0 6.9423-310 6.9529-310 2.3477-313 3.1632-317
2 4 4 Matrix: z1
1 z1 1.0000E+00 1.0000E+00 6.3240-322 0.0000E+00
2 z1 1.4142E+00 2.1213E+00 0.0000E+00 6.9423-310
3 z1 1.4142E+00 2.1213E+00 6.9423-310 6.9529-310
4 z1 6.9423-310 6.9423-310 6.9423-310 6.9529-310
1 3 0 1.00000000E-03 2.12132000E+00 abort: FACE MATCH FAIL
2 4 3 2.12132000E-03 1.00000031E+00 abort: FACE MATCH FAIL
0 quit
The mesh within my .rea file is as follows:
20.0000 20.0000 -8.00000 -4.00000 XFAC,YFAC,XZERO,YZERO
**MESH DATA** 1st line is X of corner 1,2,3,4. 2nd line is Y.
2 2 2 NEL,NDIM,NELV
ELEMENT 1 [ 1A] GROUP 0
-2.121320 -1.414213 -0.8742276E-07 -0.1311341E-06
2.121320 1.414213 2.000000 3.000000
ELEMENT 2 [ 1B] GROUP 0
-0.1311341E-06 -0.8742276E-07 1.414213 2.121320
3.000000 2.000000 1.414213 2.121320
***** CURVED SIDE DATA *****
4 Curved sides follow IEDGE,IEL,CURVE(I),I=1,5, CCURVE
2 1 -2.00000 0.00000 0.00000 0.00000 0.00000 C
4 1 3.00000 0.00000 0.00000 0.00000 0.00000 C
2 2 -2.00000 0.00000 0.00000 0.00000 0.00000 C
4 2 3.00000 0.00000 0.00000 0.00000 0.00000 C
***** BOUNDARY CONDITIONS *****
***** FLUID BOUNDARY CONDITIONS *****
P 1 1 2.00000 3.00000 0.00000 0.00000 0.00000
W 1 2 0.00000 0.00000 0.00000 0.00000 0.00000
E 1 3 2.00000 1.00000 0.00000 0.00000 0.00000
W 1 4 0.00000 0.00000 0.00000 0.00000 0.00000
E 2 1 1.00000 3.00000 0.00000 0.00000 0.00000
W 2 2 0.00000 0.00000 0.00000 0.00000 0.00000
P 2 3 1.00000 1.00000 0.00000 0.00000 0.00000
W 2 4 0.00000 0.00000 0.00000 0.00000 0.00000
***** NO THERMAL BOUNDARY CONDITIONS *****
<mime-attachment.png>