[CIG-SHORT] Issue with Normal Orientation

[Brad Aagaard]

Romain,

Matt and I discussed this issue and we believe we have mapped out how to 
fix this problem. This problem is occurring because Matt changed the way 
we create the fault mesh when we switched to DMPlex.

We are planning to go ahead and release ASAP and do a bugfix release 
when Matt fixes this issue (hopefully in the next week or so).

Brad


On 02/18/2015 10:15 AM, Romain Jolivet wrote:
> Wonderful. Let me know when you will have time to update the git repository (no rush, I will stick to 4 processes until then).
> R
>
> Romain Jolivet
> Cambridge, UK
> jolivetinsar at gmail.com
> UK: +44 7596 703 148
> France: +33 6 52 91 76 39
> USA: +1 (626) 560 6356
>
> On 18 Feb 2015, at 18:05, Brad Aagaard <baagaard at usgs.gov> wrote:
>
>> On 02/18/2015 07:37 AM, Romain Jolivet wrote:
>>>
>>> On 18 Feb 2015, at 15:30, Matthew Knepley <knepley at mcs.anl.gov> wrote:
>>>
>>>> On Wed, Feb 18, 2015 at 9:21 AM, Romain Jolivet <jolivetinsar at gmail.com> wrote:
>>>> Hello,
>>>>
>>>> I am having some more problems with the orientation of the vectors normal to the fault.
>>>>
>>>> My problem is in 2D, with a fault slicing the entire medium (see tar ball attached).
>>>> The fault is FaultCohesiveKin class. I enforce the slip.
>>>>
>>>> When I run with 1, 2, 3, 5, 7 or 9 nodes, I have one orientation for the vectors normal to the fault.
>>>> When I run with 4, 6 or 8 nodes, I have the opposite orientation.
>>>> I haven’t tried with more nodes.
>>>>
>>>> That should not happen, right?
>>>>
>>>> This can happen. We have to choose some orientation for the fault, and different partitions
>>>> produce different starting points, and thus the fault can have the opposite normal.
>>>>
>>>>    Thanks,
>>>>
>>>>       Matt
>>>
>>> Hmmm, that is not very convenient…
>>> Isn’t it possible to add a constraint on the normal vector for that starting point?
>>>
>>> Something like:
>>> If it points down, take the opposite,
>>> If it is horizontal and points to the right, make it point to the left.
>>
>> Romain,
>>
>> I know I originally setup the fault slip orientation information so that it didn't matter which direction the normal pointed in a 2-D simulation. Positive tangential slip was left lateral and negative tangential slip was right lateral. I can see from your example that this is not happening now. I will fix this.
>>
>> Brad
>>
>>
>> _______________________________________________
>> CIG-SHORT mailing list
>> CIG-SHORT at geodynamics.org
>> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>
>
>
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>