nate-sime/dolfin_dg

Utility functions for the automatic generation of nonlinear DG FEM formulations using UFL

dolfin_dg's Issues

Problem with CompressibleNavierStokesOperator.

Upon loading and using the CompressibleNavierStokesOperator I get the following error message. I trying to figure out whether this is a dune-fem bug, or a ufl problem. Would it be possible to double check with FEinCS (or FEniCSx) whether this still works?

Traceback (most recent call last):
  File "/home/robertk/work/Dune/stable/spacetimelobattocode/navierstokes/LobattoNavierStokes.py", line 118, in <module>
    F = SpatialOperator(L_h, V_h, time, dt, cfl=1)
  File "/home/robertk/work/Dune/stable/spacetimelobattocode/navierstokes/../common/spatialoperator.py", line 6, in __init__
    self._op = molGalerkin( L_h )
  File "/home/robertk/work/Dune/stable/dune-fem/build-cmake/python/dune/fem/operator/__init__.py", line 180, in molGalerkin
    return _galerkin(integrands, domainSpace=domainSpace, rangeSpace=rangeSpace,
  File "/home/robertk/work/Dune/stable/dune-fem/build-cmake/python/dune/fem/operator/__init__.py", line 105, in _galerkin
    integrands = makeIntegrands(domainSpace.gridView,integrands)
  File "/home/robertk/work/Dune/stable/dune-fem/build-cmake/python/dune/fem/model/_models.py", line 73, in integrands
    Integrands = integrands.load(view, form, *args,
  File "/home/robertk/work/Dune/stable/dune-fem/build-cmake/python/dune/models/integrands/load.py", line 293, in load
    module = builder.load(source.name(), source, "Integrands")
  File "/home/robertk/work/Dune/stable/dune-common/build-cmake/python/dune/generator/cmakebuilder.py", line 382, in load
    self._buildModule( moduleName, source, pythonName, extraCMake )
  File "/home/robertk/work/Dune/stable/dune-common/build-cmake/python/dune/generator/cmakebuilder.py", line 654, in _buildModule
    compilationMessage = self._configureWithMake( moduleName, source, pythonName )
  File "/home/robertk/work/Dune/stable/dune-common/build-cmake/python/dune/generator/cmakebuilder.py", line 630, in _configureWithMake
    code = str(source)
  File "/home/robertk/work/Dune/stable/dune-fem/build-cmake/python/dune/models/integrands/load.py", line 135, in __str__
    integrands = _compileUFL(self.integrands, self.form, *self.args, tempVars=self.tempVars)
  File "/home/robertk/work/Dune/stable/dune-fem/build-cmake/python/dune/models/integrands/ufl.py", line 236, in _compileUFL
    integrals = splitForm(form, [phi], boundaryId)
  File "/home/robertk/work/Dune/stable/dune-fem/build-cmake/python/dune/ufl/linear.py", line 186, in splitForm
    form = applyRestrictions(form)
  File "/home/robertk/work/Dune/stable/dune-fem/build-cmake/python/dune/ufl/applyrestrictions.py", line 42, in applyRestrictions
    return map_integrand_dags(RestrictionPropagator(), form, only_integral_type=integral_types)
  File "/home/robertk/work/Dune/stable/ufl/ufl/algorithms/map_integrands.py", line 46, in map_integrand_dags
    return map_integrands(lambda expr: map_expr_dag(function, expr, compress),
  File "/home/robertk/work/Dune/stable/ufl/ufl/algorithms/map_integrands.py", line 27, in map_integrands
    mapped_integrals = [map_integrands(function, itg, only_integral_type)
  File "/home/robertk/work/Dune/stable/ufl/ufl/algorithms/map_integrands.py", line 27, in <listcomp>
    mapped_integrals = [map_integrands(function, itg, only_integral_type)
  File "/home/robertk/work/Dune/stable/ufl/ufl/algorithms/map_integrands.py", line 35, in map_integrands
    return itg.reconstruct(function(itg.integrand()))
  File "/home/robertk/work/Dune/stable/ufl/ufl/algorithms/map_integrands.py", line 46, in <lambda>
    return map_integrands(lambda expr: map_expr_dag(function, expr, compress),
  File "/home/robertk/work/Dune/stable/ufl/ufl/corealg/map_dag.py", line 36, in map_expr_dag
    result, = map_expr_dags(function, [expression], compress=compress,
  File "/home/robertk/work/Dune/stable/ufl/ufl/corealg/map_dag.py", line 97, in map_expr_dags
    r = handlers[v._ufl_typecode_](v)
  File "/home/robertk/work/Dune/stable/dune-fem/build-cmake/python/dune/ufl/applyrestrictions.py", line 17, in restricted
    return map_expr_dag(RestrictionPropagator(expr.side()), expr.ufl_operands[0])
  File "/home/robertk/work/Dune/stable/ufl/ufl/corealg/map_dag.py", line 36, in map_expr_dag
    result, = map_expr_dags(function, [expression], compress=compress,
  File "/home/robertk/work/Dune/stable/ufl/ufl/corealg/map_dag.py", line 99, in map_expr_dags
    r = handlers[v._ufl_typecode_](v, *[vcache[u] for u in v.ufl_operands])
TypeError: RestrictionPropagator.variable() takes 3 positional arguments but 4 were given