[CIG-SHORT] Principle stress

[Brad Aagaard]

Xiao,

The file Charles sent is included in the binary. It is
src/pylith/playpen/postproc/princaxes.py. This is a python script, and 
the parameters are in princaxes.cfg, which is included in the same 
directory. Note that this script requires TVTK that is part of the 
Enthought Mayavi distribution. You probably don't have this software 
installed.

A better approach is to have PyLith output to HDF5 files, which can be 
accessed from Matlab or Python (via h5py) with just a few lines of code 
(see src/pylith/examples/2d/greensfns/invert_slip.py for an example of 
using h5py). The Python script princaxes.py should still help you get 
started, but it just solves the eigenvalue problem for the principal 
stresses.

Regards,
Brad


On 12/04/2013 08:17 AM, Ma, Xiao wrote:
> Hi Charles,
> I didn't get your email with the attachment the other day, can you resent me the script for the VTK files for the principle stress purpose ?
> Best ,
> Xiao
>
>
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>