Skip to content
Snippets Groups Projects
Commit 20a217a8 authored by Wolfgang Merkt's avatar Wolfgang Merkt
Browse files

test_geometry: More output

parent 04bc85a1
No related branches found
No related tags found
No related merge requests found
......@@ -13,15 +13,21 @@ def isapprox(a,b,epsilon=1e-6):
return abs(a-b)<epsilon
# --- Quaternion ---------------------------------------------------------------
# Coefficient initialisation
verbose and print("[Quaternion] Coefficient initialisation")
q = Quaternion(1,2,3,4)
q.normalize()
assert(isapprox(np.linalg.norm(q.coeffs()),q.norm()))
assert(isapprox(np.linalg.norm(q.coeffs()),1))
# Coefficient-vector initialisation
verbose and print("[Quaternion] Coefficient-vector initialisation")
v = np.array([0.5,-0.5,0.5,0.5])
qv = Quaternion(v)
assert(isapprox(qv.coeffs(), v))
# Angle axis initialisation
verbose and print("[Quaternion] AngleAxis initialisation")
r = AngleAxis(q)
q2 = Quaternion(r)
assert(q==q)
......@@ -34,6 +40,8 @@ Rr = r.matrix()
assert(isapprox(Rq.dot(Rq.T),np.eye(3)))
assert(isapprox(Rr,Rq))
# Rotation matrix initialisation
verbose and print("[Quaternion] Rotation Matrix initialisation")
qR = Quaternion(Rr)
assert(q.isApprox(qR))
assert(isapprox(q.coeffs(),qR.coeffs()))
......@@ -43,7 +51,7 @@ try:
qR[5]
print("Error, this message should not appear.")
except RuntimeError as e:
if verbose: print("As expected, catched exception: ",e)
if verbose: print("As expected, caught exception: ", e)
# --- Angle Vector ------------------------------------------------
r = AngleAxis(.1,np.array([1,0,0],np.double))
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment