GalerkinScheme fails with non virtualized integrands
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)