def Delta_Theta(v1,v2,n): if np.abs(Parallel(v1,n)) + np.abs(Parallel(v2,n)) == 0: vproj1 = vctr.Subtract(v1,vctr.ScalarMul((vctr.Dot(v1,n)),n)) vproj2 = vctr.Subtract(v2,vctr.ScalarMul((vctr.Dot(v2,n)),n)) deltheta = vctr.Angle(vproj1,vproj2) else: deltheta = 0 return deltheta
def Delta_Theta(v1, v2, n): vproj1 = vctr.Subtract(v1, vctr.ScalarMul((vctr.Dot(v1, n)), n)) vproj2 = vctr.Subtract(v2, vctr.ScalarMul((vctr.Dot(v2, n)), n)) deltheta = vctr.Angle(vproj1, vproj2) sign = np.sign(vctr.Dot(n, vctr.Cross(va1, va2))) nopt = vctr.ScalarMul(sign, n) return sign * deltheta