Changeset 21275
- Timestamp:
- 10/15/16 14:29:50 (8 years ago)
- Location:
- issm/trunk-jpl/src/m/mech
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/mech/mechanicalproperties.m
r20214 r21275 93 93 %eigen values and vectors 94 94 [directions,value]=eig(stress); 95 %sort by a bsolute value of eigenvaluein descending order96 [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'); 97 97 %re-order eigenvalues and associated vectors 98 98 value=value(idx,idx); … … 102 102 [directions,value]=eig(strain); 103 103 %same for strainrate 104 [val,idx]=sort( abs(diag(value)),'descend');104 [val,idx]=sort(diag(value),'descend'); 105 105 value=value(idx,idx); 106 106 directions=directions(:,idx); -
issm/trunk-jpl/src/m/mech/mechanicalproperties.py
r21254 r21275 109 109 #eigenvalues and vectors for stress 110 110 value,directions=np.linalg.eig(stress); 111 idx= abs(value).argsort()[::-1] # sort in descendingorder111 idx=value.argsort()[::-1] # sort in descending algebraic (not absolute) order 112 112 value=value[idx] 113 113 directions=directions[:,idx] … … 117 117 #eigenvalues and vectors for strain 118 118 value,directions=np.linalg.eig(strain); 119 idx= abs(value).argsort()[::-1] # sort in descending order119 idx=value.argsort()[::-1] # sort in descending order 120 120 value=value[idx] 121 121 directions=directions[:,idx]
Note:
See TracChangeset
for help on using the changeset viewer.