Dear All,
From my FE results I have imported the individual components of a stress
tensor and compiled them into an M3 matrix using the M3fromelements macro
(the result of print (rd) on the data is an 1x9 matrix NOT a 3x3 matrix as
expected).
Now, I need to plot and work with the Principle stresses (both magnitudes
and orientations) which "fall out" of the stress tensor as the eigenvalues
and eigenvectors.
I hunted for an eigenvector/value module but didn't find one. What I did
notice though was that Autoglyph creates a plot similar to what I need and
so presumably must export the data as well. Thus, I connected the output of
the M3 macro to an autoglyph and tried to sent the output to image.
Unfortunately this doesn't work since autoglyph returns: invalid data: data
must be scalar, 2-vector, 3-vector, 2x2 or 3x3 matrix. I have tried making
the 3x3 matrix without the matrix macro, but always it ends up in the 1x9
format and with the autoglyph error. Is there a way to reconfigure the
entire data set from the 1x9 matrix to 3x3? or do I have to work on each
record individually then recombine into the total dataset?
I looked into writing a eigenvactor/value module, but saw in Numerical
recipies that it is highly non trivial and since I doubt either my maths or
programming is up to it. However in the course of that hunt, I noticed in
the opendx source code the eigen.c and eigen.h routines, which are
presumably used by, amongst others, autoglyph. So, I looked in the
documentation for an DX_eigen call or similar but didn't find one in the
programmers manual. I know the manual is not up to date so am wondering if
the call exists but isn't listed anywhere I found or if it just isn't part
of the exported command set.
I'd greatly appreciate pointers on both of the above queries, namely the
description of my matrix and the possible implementation of
eigenvector/eigenvector calls in the programming and/or user interface.
Many thanks for your time
Peter
Dr Peter Connolly
____________________________________
Geophysikalisches Institut
der Universitat Karlsruhe
Hertzstrasse 16
76185 Karlsruhe
Germany
tel +49 (0)721 6084593
fax +49 (0)721 71173
Email addresses: (both checked)
[EMAIL PROTECTED]
OR
[EMAIL PROTECTED]
