Segmentation error when using my own .inp files on the Monte Carlo 01_disclination simulation

[edpaja]

Hi all,

I am trying to simulate stable configurations of mobius strips using a very basic mobius strip mesh.

vertex file:
0 0.000 0.000 0.000 1
1 0.000 1.000 0.000 1
2 1.000 0.000 0.000 1
3 1.000 1.000 0.000 1
4 2.000 0.500 -0.500 1
5 2.000 0.500 0.500 1
6 0.500 2.000 -0.500 1
7 0.500 2.000 0.500 1
8 -0.500 0.500 -0.500 1
9 -0.500 0.500 0.500 1

face file:
0 0 2 3 1 1
1 0 3 1 1 1
2 2 4 3 1 1
3 2 5 4 1 1
4 4 5 6 1 1
5 5 7 6 1 1
6 6 7 8 1 1
7 7 9 8 1 1
8 8 1 0 1 1
9 8 9 1 1 1

I think I am getting the error because the mesh is not reading correctly but I do not understand the documentation well enough to understand why. Could it be because two adjacent faces have to have opposite normal vectors?

Please let me know if you have any ideas or advice.

Thank you

[fdmatoz]

Hi, perhaps you are having issues with the orientation of the surface. " two adjacent faces have to have opposite normal vectors?" Yes you can but you must be careful that the order of vertices are correct otherwise the "circulation of the half edges won't work". Is there anyway for you to attach the files so I can test it?

[edpaja]

mobius_inp.zip

[edpaja]

faces_twist1.txt]
vertices_twist1.txt]

I used the Monte Carlo disclination.py example script.