uflacs + discontinuous Galerkin doesn't work together
Issue #45
resolved
from dolfin import *
parameters['form_compiler']['representation'] = 'uflacs'
mesh = RectangleMesh(0,0,2,1,100,50,"crossed")
norm = FacetNormal(mesh)
M = VectorFunctionSpace(mesh, "Lagrange", 2)
P0 = FunctionSpace(mesh, "DG", 0)
fi = TestFunction(P0)
phi = Function(P0)
phi_k = Function(P0)
vm0 = Function(M)
dt = 1.e-5
bcs = []
vmn = abs(dot(vm0("+"),norm("+")))
phi_star = conditional(ge(dot(vm0,norm),0.),phi,0.1)
F = phi*fi*dx - phi_k*fi*dx - dt*phi*dot(vm0,nabla_grad(fi))*dx \
+ avg(dt)*(dot(vm0("+"),jump(fi,norm))*avg(phi) \
+ 0.5*vmn*dot(jump(phi,norm),jump(fi,norm)))*dS \
+ dt*phi_star*dot(vm0,norm)*fi*ds
J = derivative(F, phi)
problem = NonlinearVariationalProblem(F, phi, bcs, J)
solver = NonlinearVariationalSolver(problem)
fails with message
In instant.recompile: The module did not compile with command 'make VERBOSE=1', see '/home/klara/.instant/error/ffc_form_bc63ceaa445ac7280e961c8add8f8c283034afb3/compile.log'
Traceback (most recent call last):
File "uflacs+dg.py", line 35, in <module>
problem = NonlinearVariationalProblem(F, phi, bcs, J)
File "/usr/lib/python2.7/dist-packages/dolfin/fem/solving.py", line 117, in __init__
F = Form(F, form_compiler_parameters=form_compiler_parameters)
File "/usr/lib/python2.7/dist-packages/dolfin/fem/form.py", line 70, in __init__
mpi_comm=mesh.mpi_comm())
File "/usr/lib/python2.7/dist-packages/dolfin/compilemodules/jit.py", line 68, in mpi_jit
output = local_jit(*args, **kwargs)
File "/usr/lib/python2.7/dist-packages/dolfin/compilemodules/jit.py", line 128, in jit
return form_compiler.jit(form, parameters=p)
File "/usr/lib/python2.7/dist-packages/ffc/jitcompiler.py", line 74, in jit
return jit_form(ufl_object, parameters)
File "/usr/lib/python2.7/dist-packages/ffc/jitcompiler.py", line 148, in jit_form
cache_dir = cache_dir)
File "/usr/lib/python2.7/dist-packages/ffc/backends/ufc/build.py", line 73, in build_ufc_module
**kwargs)
File "/usr/lib/python2.7/dist-packages/instant/build.py", line 563, in build_module
recompile(modulename, module_path, new_compilation_checksum, build_system)
File "/usr/lib/python2.7/dist-packages/instant/build.py", line 165, in recompile
instant_error(msg % (cmd, compile_log_filename_dest))
File "/usr/lib/python2.7/dist-packages/instant/output.py", line 85, in instant_error
raise RuntimeError(text)
RuntimeError: In instant.recompile: The module did not compile with command 'make VERBOSE=1', see '/home/klara/.instant/error/ffc_form_bc63ceaa445ac7280e961c8add8f8c283034afb3/compile.log'
Comments (5)
-
-
I tested this and it works fine with latest dev, but fails with 1.5.
-
That was my hunch, 1.5 had some similar issues that are fixed. Please use the latest development version, or wait for the 1.6 release in late June.
-
- changed status to resolved
-
reporter I am using version 1.5 and so I will wait for the release of 1.6. Tanks very much for your help!
- Log in to comment
Which version are you using? 1.5 or recent dev?
Can you post the first error message in the compile.log file mentioned at that last line?