Function test_code_generation_volume_terms

tests/pytest/test_code_generation.py:7–38

(unit_mesh_3d)

@pytest.mark.slow
def test_code_generation_volume_terms(unit_mesh_3d):
    fes = L2(unit_mesh_3d, order=5)
    gfu = GridFunction(fes)

    cfvec = CF((1,)*10)
    cf10 = InnerProduct(cfvec, cfvec)

    cf_tensor3 = CF(tuple(range(1,25)), dims=(2,3,4))
    cf_tensor = InnerProduct(cf_tensor3, cf_tensor3)

    # piecewise polynomials - also test interpolation on L2 space and the resulting GridFunction
    functions = [x,y,x*y, specialcf.mesh_size, CoefficientFunction((x,y)).Norm()**2, Id(3)[:,2][2], cf10, cf_tensor]

    for cf in functions:
        gfu.Set(cf)
        # should all give the same results
        cfs = [ cf.Compile(), cf.Compile(True, wait=True), gfu, gfu.Compile(), gfu.Compile(True, wait=True) ]

        for f in cfs:
            assert Integrate( (cf-f)*(cf-f), unit_mesh_3d) == approx(0)

    cf8x8 = CoefficientFunction( tuple(range(8*8)), dims=(8,8) )
    functions = [sin(x)*y, exp(x)+y*y*y, (1+x)**(1+y), cf8x8[1::3, 2:6], cf8x8[:,3], cf8x8[:,1:8:2]]
    for cf in functions:
        cfs = [ cf.Compile(), cf.Compile(True, wait=True)]
        for f in cfs:
            assert Integrate( Norm(cf-f), unit_mesh_3d) == approx(0)

    cf = atan2(1+x,1+y)
    cfs = [ cf.Compile(), cf.Compile(True, wait=True, maxderiv=0)]
    for f in cfs:
        assert Integrate( (cf-f)*(cf-f), unit_mesh_3d) == approx(0)

@pytest.mark.slow
def test_code_generation_boundary_terms(unit_mesh_3d):