Issue with subdomains using dolfinx-master

22 views
Skip to first unread message

Ganesh Diwan

unread,
Jul 20, 2019, 5:35:34 PM7/20/19
to fenic...@googlegroups.com
Dear dolfinx developers;
Is there a way to solve for the weak form defined over subdomains with any of the branches other than the dolfinx master?

I was trying to set up my weak form as:
a = inner(u, v) * dx(1) + inner(u, v) * dx(2)
L = inner(g, v) * dx(1) + inner(g, v) * dx(2)
solve(a==L,u,[])

but end up in error in solve:

    solve(a==L,u,[])
  File "/usr/local/lib/python3.6/dist-packages/dolfin/fem/solving.py", line 113, in solve
    return _solve_varproblem(*args, **kwargs)
  File "/usr/local/lib/python3.6/dist-packages/dolfin/fem/solving.py", line 126, in _solve_varproblem
    a = fem.Form(eq.lhs, form_compiler_parameters=form_compiler_parameters)
  File "/usr/local/lib/python3.6/dist-packages/dolfin/fem/form.py", line 64, in __init__
    self._cpp_object = cpp.fem.Form(ufc_form, function_spaces)
RuntimeError: Cell integral subdomain not supported. Under development.

Perhaps there is another branch where this could work?

Thanks in advance.

Cheers,
Ganesh
Reply all
Reply to author
Forward
0 new messages