Example #1
0
    def testNonPlainDoubleCouples(self):
        """Convert random MTs to plain double couples and compare angles."""
        for i in range(100):
            ms = [random.random() * 1.0e20 - 0.5e20 for j in range(6)]
            m = num.matrix([[ms[0], ms[3], ms[4]], [ms[3], ms[1], ms[5]], [ms[4], ms[5], ms[2]]], dtype=num.float)

            m1 = MomentTensor(m=m)
            m_plain = m1.m_plain_double_couple()
            m2 = MomentTensor(m=m_plain)

            self.assertAnglesSame(m1, m2)
Example #2
0
    def testNonPlainDoubleCouples(self):
        for i in range(100):
            ms = [random.random() * 1.0e20 - 0.5e20 for j in range(6)]
            m = num.matrix([[ms[0], ms[3], ms[4]], [ms[3], ms[1], ms[5]],
                            [ms[4], ms[5], ms[2]]],
                           dtype=num.float)

            m1 = MomentTensor(m=m)
            m_plain = m1.m_plain_double_couple()
            m2 = MomentTensor(m=m_plain)

            self.assertAnglesSame(m1, m2)