Definition at line 75 of file general_optics.py. 00075 : 00076 "eulermat(theta, eta, phi) returns r(eta)*r(theta)*r(phi) phi=orientation around normal, theta=yaw, eta=normal roll" 00077 pm=Numeric.identity(3).astype(numeric_float) 00078 em=Numeric.identity(3).astype(numeric_float) 00079 tm=Numeric.identity(3).astype(numeric_float) 00080 sp, cp=sincosdeg(phi) 00081 st, ct=sincosdeg(theta) 00082 se, ce=sincosdeg(eta) 00083 pm[0,0]=pm[1,1]=cp; pm[0,1]=-sp; pm[1,0]=sp 00084 em[0,0]=em[1,1]=ce; em[0,1]=-se; em[1,0]=se 00085 tm[0,0]=tm[2,2]=ct; tm[2,0]=-st; tm[0,2]=st 00086 return Numeric.dot(em, Numeric.dot(tm, pm)) 00087 00088 class euler:
|