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 ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
