def test_fenics_vector_copy(): mesh = UnitSquare(3, 3) fs = FunctionSpace(mesh, "CG", 1) vec1 = FEniCSVector(Function(fs)) vec1.coeffs = np.array(range(fs.dim())) assert_equal(vec1.coeffs[0], 0) vec2 = vec1.copy() vec2.coeffs[0] = 5 assert_equal(vec1.coeffs[0], 0)
def test_fenics_vector_inner(): mesh = UnitSquare(3, 3) fs = FunctionSpace(mesh, "CG", 1) vec = FEniCSVector(Function(fs)) vec.coeffs = np.array(range(fs.dim())) assert_equal(vec.__inner__(vec), 1240)