suzanne.soy/geosolver
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

bug fix in make_hcs_3d.py

Browse Source
This commit is contained in:
kwikrick 2014-02-09 07:50:44 +00:00
parent 4ca358c1ef
commit bac3e1e103
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
geosolver/intersections.py
View File

@ -343,13 +343,13 @@ def make_hcs_3d (a, b, c, righthanded=True):
v = vector.vector([0.0,1.0,0.0]) v = vector.vector([0.0,1.0,0.0])
elif tol_eq(nu, 0.0): elif tol_eq(nu, 0.0):
# determine u perpendicular from v # determine u perpendicular from v
u,dummy = perp_3d(v)[0] u,dummy = perp_3d(v)
elif tol_eq(nv, 0.0): elif tol_eq(nv, 0.0):
# determine v perpendicular from u # determine v perpendicular from u
dummy,v = perp_3d(u)[0] dummy,v = perp_3d(u)
# ensure that u and v are different # ensure that u and v are different
if tol_eq(vector.norm(u-v),0.0): if tol_eq(vector.norm(u-v),0.0):
dummy,v = perp_3d(u)[0] dummy,v = perp_3d(u)
# make the basis vectors orthogonal # make the basis vectors orthogonal
w = vector.cross(u,v) w = vector.cross(u,v)
v = vector.cross(w,u) v = vector.cross(w,u)