unittest/python/explog.py

Go to the documentation of this file.

import unittest
import math
 
import numpy as np
import pinocchio as pin
 
from pinocchio.utils import rand, zero, eye
from pinocchio.explog import exp, log
 
from test_case import PinocchioTestCase as TestCase
 
 
class TestExpLog(TestCase):
    def test_exp3(self):
        v = zero(3)
        m = pin.exp3(v)
        self.assertApprox(m, eye(3))
 
    def test_Jexp3(self):
        v = zero(3)
        m = pin.Jexp3(v)
        self.assertApprox(m, eye(3))
 
    def test_log3(self):
        m = eye(3)
        v = pin.log3(m)
        self.assertApprox(v, zero(3))
 
    def test_log3_quat(self):
        """Test log3 over the quaternions."""
        from eigenpy import Quaternion
 
        quat_v = np.array([0.0, 0.0, 0.0, 1.0])
        quat = Quaternion(quat_v)
        v0 = pin.log3(quat)
        self.assertApprox(v0, zero(3))
 
        v1 = pin.log3(quat_v)
        self.assertApprox(v1, zero(3))
 
    def test_exp3_quat(self):
        self.assertApprox(pin.exp3_quat(zero(3)), np.array([0.0, 0.0, 0.0, 1.0]))
 
    def test_Jlog3(self):
        m = eye(3)
        J = pin.Jlog3(m)
        self.assertApprox(J, eye(3))
 
    def test_exp6(self):
        v = pin.Motion.Zero()
        M = pin.exp6(v)
        self.assertTrue(M.isIdentity())
 
        M2 = pin.exp6(np.array(v))
        self.assertTrue(M2.isIdentity())
 
    def test_Jexp6(self):
        v = pin.Motion.Zero()
        J = pin.Jexp6(v)
        self.assertApprox(J, eye(6))
 
        J2 = pin.Jexp6(np.array(v))
        self.assertApprox(J, J2)
 
    def test_log6(self):
        m = pin.SE3.Identity()
        v = pin.log6(m)
        self.assertApprox(v.vector, zero(6))
 
    def test_log6_homogeneous(self):
        m = eye(4)
        v = pin.log6(m)
        self.assertApprox(v.vector, zero(6))
 
    def test_log6_quat(self):
        q0 = np.array([0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0])
        v = pin.log6_quat(q0)
        self.assertApprox(v.vector, zero(6))
 
    def test_Jlog6(self):
        m = pin.SE3.Identity()
        J = pin.Jlog6(m)
        self.assertApprox(J, eye(6))
 
    def test_skew(self):
        u = np.random.rand((3))
        v = np.random.rand((3))
 
        u_skew = pin.skew(u)
        u_unskew = pin.unSkew(u_skew)
 
        self.assertApprox(u, u_unskew)
 
        v_skew = pin.skew(v)
        u_v_square = pin.skewSquare(u, v)
 
        self.assertApprox(u_v_square, u_skew.dot(v_skew))
 
    def test_explog(self):
        self.assertApprox(exp(42), math.exp(42))
        self.assertApprox(log(42), math.log(42))
        self.assertApprox(exp(log(42)), 42)
        self.assertApprox(log(exp(42)), 42)
 
        m = rand(3)
        self.assertLess(np.linalg.norm(m), np.pi)  # necessary for next test
        self.assertApprox(log(exp(m)), m)
 
        m = np.random.rand(3)
        self.assertLess(np.linalg.norm(m), np.pi)  # necessary for next test
        self.assertApprox(log(exp(m)), m)
 
        m = pin.SE3.Random()
        self.assertApprox(exp(log(m)), m)
 
        m = rand(6)
        self.assertLess(
            np.linalg.norm(m), np.pi
        )  # necessary for next test (actually, only angular part)
        self.assertApprox(log(exp(m)), m)
 
        m = np.random.rand(6)
        self.assertLess(
            np.linalg.norm(m), np.pi
        )  # necessary for next test (actually, only angular part)
        self.assertApprox(log(exp(m)), m)
 
        m = eye(4)
        self.assertApprox(exp(log(m)).homogeneous, m)
 
        with self.assertRaises(ValueError):
            exp(eye(4))
        with self.assertRaises(ValueError):
            exp(list(range(3)))
        with self.assertRaises(ValueError):
            log(list(range(3)))
        with self.assertRaises(ValueError):
            log(zero(5))
        with self.assertRaises(ValueError):
            log(zero((3, 1)))
 
 
if __name__ == "__main__":
    unittest.main()