Example #1
0
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)