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
