Changeset 21275


Ignore:
Timestamp:
10/15/16 14:29:50 (8 years ago)
Author:
cborstad
Message:

BUG: principal stress convention is to sort by algebraic (not absolute) value

Location:
issm/trunk-jpl/src/m/mech
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/m/mech/mechanicalproperties.m

    r20214 r21275  
    9393        %eigen values and vectors
    9494        [directions,value]=eig(stress);
    95         %sort by absolute value of eigenvalue in descending order
    96         [val,idx]=sort(abs(diag(value)),'descend');
     95        %sort by algebraic value of eigenvalue (not absolute value) in descending order
     96        [val,idx]=sort(diag(value),'descend');
    9797        %re-order eigenvalues and associated vectors
    9898        value=value(idx,idx);
     
    102102        [directions,value]=eig(strain);
    103103        %same for strainrate
    104         [val,idx]=sort(abs(diag(value)),'descend');
     104        [val,idx]=sort(diag(value),'descend');
    105105        value=value(idx,idx);
    106106        directions=directions(:,idx);
  • issm/trunk-jpl/src/m/mech/mechanicalproperties.py

    r21254 r21275  
    109109                #eigenvalues and vectors for stress
    110110                value,directions=np.linalg.eig(stress);
    111                 idx=abs(value).argsort()[::-1] # sort in descending order
     111                idx=value.argsort()[::-1] # sort in descending algebraic (not absolute) order
    112112                value=value[idx]
    113113                directions=directions[:,idx]
     
    117117                #eigenvalues and vectors for strain
    118118                value,directions=np.linalg.eig(strain);
    119                 idx=abs(value).argsort()[::-1] # sort in descending order
     119                idx=value.argsort()[::-1] # sort in descending order
    120120                value=value[idx]
    121121                directions=directions[:,idx]
Note: See TracChangeset for help on using the changeset viewer.