On Fri, Jan 5, 2018 at 8:39 AM, Viviana Palacio Betancur < [email protected]> wrote:
> Hi John, > > Thank you for the clarification, it's really helpful. > > My assembly function is separated, one for the matrix and another one for > the rhs. I loop over the elements of each subdomain to avoid the if > statement. > OK... and are you using subdomain restricted variables? In other words, are you using a variable that only has DOFs on the boundary subdomain (system.add_variable(name, order, family, std::set<subdomain_id>))? -- John ------------------------------------------------------------------------------ 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 [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
