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