Adaptive mesh refinement, problem with boundary terms in linear form
Hello,
when running the attached file the following problem occurs. When I consider
L = inner(f, v)*dx
I can solve a linear elasticity equation with adaptive mesh refinement.
Whenever I include terms with ds in the linear form, like
L = inner(f, v)*dx + inner(g, v)*ds(2)
I receive the following error message:
In instant.recompile: The module did not compile with command 'make VERBOSE=1 ', see '/home/
Traceback (most recent call last):
File "Test_MinimalEx
solve(a == L, u, bc, tol=5.0E-7, M=MF) #Solve adaptively
File "/usr/lib/
_solve_
File "/usr/lib/
solver = AdaptiveLinearV
File "/usr/lib/
ec = generate_
File "/usr/lib/
forms = [Form(form, form_compiler_
File "/usr/lib/
common_cell)
File "/usr/lib/
return local_jit(*args, **kwargs)
File "/usr/lib/
return jit_compile(form, parameters=p, common_
File "/usr/lib/
return jit_form(
File "/usr/lib/
cache_dir = cache_dir)
File "/usr/lib/
**kwargs)
File "/usr/lib/
recompile(
File "/usr/lib/
instant_
File "/usr/lib/
raise RuntimeError(text)
RuntimeError: In instant.recompile: The module did not compile with command 'make VERBOSE=1 ', see '/home/
By choosing a fine mesh manually, I can solve the problem by FEniCS accurately. The PDE is well-posed and numerical locking effects have been taken care of.
It seams to me that there is a bug or what might be the issue?
I am working with FEniCS v1.2 in python.
Thanks in advance,
SJ
Question information
- Language:
- English Edit question
- Status:
- Expired
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
This question was originally filed as bug #1188623.