- Timestamp:
- 06/07/17 10:50:54 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified issm/branches/trunk-larour-NatGeoScience2016/src/py3/geometry/slope.py ¶
r19895 r21759 1 import numpy as np y1 import numpy as np 2 2 from GetNodalFunctionsCoeff import GetNodalFunctionsCoeff 3 3 … … 33 33 alpha,beta=GetNodalFunctionsCoeff(index,x,y)[0:2] 34 34 35 summation=np y.array([[1],[1],[1]])36 sx=np y.dot(surf[index-1]*alpha,summation).reshape(-1,)37 sy=np y.dot(surf[index-1]*beta,summation).reshape(-1,)35 summation=np.array([[1],[1],[1]]) 36 sx=np.dot(surf[index-1]*alpha,summation).reshape(-1,) 37 sy=np.dot(surf[index-1]*beta,summation).reshape(-1,) 38 38 39 s=np y.sqrt(sx**2+sy**2)39 s=np.sqrt(sx**2+sy**2) 40 40 41 41 if md.mesh.dimension()==3: 42 42 sx=project3d(md,'vector',sx,'type','element') 43 43 sy=project3d(md,'vector',sy,'type','element') 44 s=np y.sqrt(sx**2+sy**2)44 s=np.sqrt(sx**2+sy**2) 45 45 46 46 return (sx,sy,s)
Note:
See TracChangeset
for help on using the changeset viewer.