Skip to content

GalerkinScheme fails with non virtualized integrands

Andreas Dedner requested to merge bugfix/nonVirtIntgrands into master

In a recent change the GalerkinScheme only works with the VirtualizedIntegrands (see #136 ). I believe it is this commit git diff ed2c4105df~ ed2c4105df. Here is a minimal test case

from ufl import TestFunction, TrialFunction, SpatialCoordinate, dx, grad, dot
from dune.grid import structuredGrid as leafGridView
from dune.fem.space import lagrange
from dune.fem.scheme import galerkin as solutionScheme

gridView = leafGridView([0, 0], [1, 1], [4, 4])
space    = lagrange(gridView)

x     = SpatialCoordinate(space)
u     = TrialFunction(space)
v     = TestFunction(space)
b     = dot(x,x)*v*dx
a     = dot(grad(u), grad(v)) * dx

# this works
scheme = solutionScheme(a == b, solver='cg', virtualize=True)
# this fails
scheme = solutionScheme(a == b, solver='cg', virtualize=False)

Merge request reports