[CIG-SHORT] Bizarre normal_dir inversion

[Brad Aagaard]

Romain,

I can reproduce your problem on my machine, but I can't reproduce it 
using other similar examples. I am wondering if there are issues 
creating your mesh. Can you send me the journal files you used to create 
your mesh?

Thanks,
Brad


On 02/02/2015 01:15 PM, Romain Jolivet wrote:
> Please find the tarball attached.
> Thanks for the help.
> Romain
>
> ----------------------------------------------------------------
> !!!!!! ADDRESS CHANGED!!!!!!
> ----------------------------------------------------------------
> Romain Jolivet
> Postdoctoral Fellow
>
> University of Cambridge
> Department of Earth Sciences
> Bullard Laboratories
> Madingley Road,
> Cambridge, CB3 0EZ,
> United Kingdom
>
> email: rpj29 at cam.ac.uk <mailto:rpj29 at cam.ac.uk>
> Phone: +44 7 596 703 148
> Phone (Fr): +33 7 52 91 76 39
> ----------------------------------------------------------------
> ----------------------------------------------------------------
>
> On Feb 2, 2015, at 20: 00, cig-short-request at geodynamics.org
> <mailto:cig-short-request at geodynamics.org> wrote:
>
>> Send CIG-SHORT mailing list submissions to
>> cig-short at geodynamics.org <mailto:cig-short at geodynamics.org>
>>
>> To subscribe or unsubscribe via the World Wide Web, visit
>> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>> or, via email, send a message with subject or body 'help' to
>> cig-short-request at geodynamics.org
>>
>> You can reach the person managing the list at
>> cig-short-owner at geodynamics.org
>>
>> When replying, please edit your Subject line so it is more specific
>> than "Re: Contents of CIG-SHORT digest..."
>>
>>
>> Today's Topics:
>>
>>   1. Re: Bizarre normal_dir inversion (Brad Aagaard)
>>
>>
>> ----------------------------------------------------------------------
>>
>> Message: 1
>> Date: Mon, 2 Feb 2015 11:52:11 -0800
>> From: Brad Aagaard <baagaard at usgs.gov>
>> To: <cig-short at geodynamics.org>
>> Subject: Re: [CIG-SHORT] Bizarre normal_dir inversion
>> Message-ID: <54CFD56B.9070003 at usgs.gov>
>> Content-Type: text/plain; charset="utf-8"; format=flowed
>>
>> Romain,
>>
>> Matt fixed a related problem in Dec/Jan in our current development
>> version. I would like to see if it resolves your issue. Can you send me
>> a tarball of all of your input files?
>>
>> We are getting very close to a bugfix release, so I would really like to
>> see if we have something else that needs fixing or it is already fixed.
>>
>> Thanks,
>> Brad
>>
>>
>> On 2/2/15 10:51 AM, Romain Jolivet wrote:
>>> 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.
>>>
>>>
>>>
>>>
>>>
>>> —————————————————————————————————————
>>> —————————————————————————————————————
>>> !!!!!! 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
>>> —————————————————————————————————————
>>> —————————————————————————————————————
>>>
>>>
>>>
>>> _______________________________________________
>>> CIG-SHORT mailing list
>>> CIG-SHORT at geodynamics.org
>>> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>>>
>>
>>
>>
>> ------------------------------
>>
>> Subject: Digest Footer
>>
>> _______________________________________________
>> CIG-SHORT mailing list
>> CIG-SHORT at geodynamics.org
>> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>>
>> ------------------------------
>>
>> End of CIG-SHORT Digest, Vol 101, Issue 2
>> *****************************************
>
>
>
>
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>