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)
