Method test_elasticity_tensor

python/pymesh/tests/test_assembler.py:62–70

(self)

        self.assert_PSD(K)

    def test_elasticity_tensor(self):
        assembler = self.single_triangle_2D_setup()
        E = assembler.assemble("elasticity_tensor")
        self.assert_symmetry(E)
        self.assert_PD(E)

        # The default material should yield identity elasticity tensor.
        I = scipy.sparse.eye(assembler.mesh.dim * (assembler.mesh.dim + 1) / 2)
        self.assertAlmostEqual(0.0, scipy.sparse.linalg.norm(E-I))