[CIG-SHORT] Issue with Normal Orientation

[Romain Jolivet]

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.

R


>  
> In terms of codes versions:
> 
> Pylith version: 	
> >> git log
> commit 30f0bbe7ea750d78005f76cff5f3ca9868a89fa7
> Merge: f047ab5 1327120
> Author: Brad Aagaard <baagaard at usgs.gov>
> Date:   Fri Feb 13 17:09:54 2015 -0800
> 
>     Merge branch 'baagaard/add-release-2.1.0'
>     
>     * baagaard/add-release-2.1.0:
>       Reset permissions for full-scale test scripts.
>       Reset permissions for pytest drivers.
>       Update script making petsc tarball.
>       Updated release notes.
>       Updated copyright.
>       Updated release information.
>       Add updates from README to release announcement.
>       Added a description of new 2D gravity example.
>       Added description of DP plasticity changes to README.
>       Updates to info on new features.
>       Started work on release notes for v2.1.0.
>       Updated version number in configure.ac.
>       Updated version version number in manual to v2.1.0.
> 
> PETSc version:
> >> git log
> commit 04fcc16b6fc1fef768ea84f574fa77a92547ade2
> Merge: a1925ac ae5cfb4
> Author: Matthew G. Knepley <knepley at gmail.com>
> Date:   Tue Feb 10 07:14:54 2015 -0600
> 
>     Merge branch 'knepley/fix-veclock-fftw' into knepley/pylith
>     
>     * knepley/fix-veclock-fftw:
>       DM: Fix more Vec lock issues
>       FFTW: Updated for vector locking
> 
> Cheers, 
> R
> 
> 
> 
> 
> —————————————————————————————————————
> —————————————————————————————————————	
> 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
> 
> 
> 
> -- 
> What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
> -- Norbert Wiener
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20150218/d333ba41/attachment-0001.html>