Exemplo n.º 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)
Exemplo n.º 2
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)