[Libmesh-users] Solving a system with 2 subdomains

[Viviana P. B. <vpa...@uc...>]

Hello,

I am solving a system on a mesh with two subdomains: bulk with elements
TET10, and surface with elements TRI6. So far, I've sucessfully asembled
the system's matrix and rhs, and initialized the system.

My problem is when I execute system.solve(). Apparently, the solution is
only done for the bulk subdomain and thus the surface solution is almost 0.
Is the default setting to solve only for the first subdomain? Should I be
using the system.restrict_solve_to?

Thank you. Best,
Viviana.

Viviana Palacio Betancur
PhD Student | 2016 Cohort
de Pablo Group
Institute for Molecular Engineering
University of Chicago