[CIG-SHORT] Bizarre normal_dir inversion

[Romain Jolivet]

Dear Pylith developers, 

I am running quite a simple problem and something very bizarre happens.

I have an elastic 2D medium with a fault that cuts in two.
I impose dirichlet conditions on one side.
I impose slip on the fault (a combination of creep and an earthquake).

When I compute on one node, everything goes fine. Same with two nodes.
However, when I get to 3 nodes (and more), some normal_dir of the fault are flipped to the other direction.
As a result, some nodes of the fault slip in one direction, and the other ones slip in the other direction.

I looked into FaultCohesiveLagrange, in _calcOrientation and, from what I printed to the screen, flipGlobal is always = 0, whatever I do (I mean run with 1, 2, 3 or 4 cores).
Please find attached the files I use to run the problem.
Let me know if you find anything that could help me.

Cheers, 
R

PS: I run on an iMac (Maverick). I compiled things manually because I never got the installer to work properly and it is possible that some links are not correctly set up.

-------------- next part --------------
A non-text attachment was scrubbed...
Name: NormalReversed.tgz
Type: application/octet-stream
Size: 81920 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20150202/47ade79d/attachment-0001.obj>
-------------- next part --------------


—————————————————————————————————————
—————————————————————————————————————	
!!!!!! ADDRESS CHANGED !!!!!

Romain Jolivet
Postdoctoral Fellow

University of Cambridge
Department of Earth Sciences
Bullard Labs
Madingley Rise 
Madingley Road
Cambridge CB3 0EZ
United Kingdom

email: rpj29 at cam.ac.uk
Phone: +44 1223 748 938
Mobile: +44 7596 703 148

France: +33 6 52 91 76 39
US: +1 (626) 560 6356
—————————————————————————————————————
—————————————————————————————————————