def test_gradient(): msh = circular_mesh(200, 1.0) bf = basis_from_nodes([0.0, 1.0]) gradient = bf.get_gradient_basis() e = msh.elements[57] chain_rule = gradient.chain_rule(e.mapping.get_jacobian(0.0)) value = np.array(gradient.evaluate(0, 0.5)) np.testing.assert_almost_equal(chain_rule * value, -31.83229765 * np.ones(2))
def test_circular_mesh(): a = circular_mesh(4, 1.0) np.testing.assert_almost_equal(a.vertices[1].loc[0], 0.0) np.testing.assert_almost_equal(a.vertices[1].loc[1], 1.0)