Definition at line 48 of file general_optics.py. 00048 : 00049 "planesolve(x0, u, x1, v) returns a, x such that x_=a*u_+x0_ solves (x_-x1_).v_=0" 00050 if vec_mag(x1-x0) < 1e-10: 00051 return 0, x1 00052 a=Numeric.dot((x1-x0), v) / Numeric.dot(u,v) 00053 #print "\nplanesolve: ", x0, u, x1, v, a, a*u+x0, "\n" 00054 return a, a*u+x0 00055 def cross(a,b):
|