No, I have the same variable for the entire mesh. I use the subdomains in order to impose specific BCs and calculate some properties for the surface only.
Viviana Palacio Betancur PhD Student | 2016 Cohort de Pablo Group Institute for Molecular Engineering University of Chicago On Fri, Jan 5, 2018 at 10:52 AM, John Peterson <jwpeter...@gmail.com> wrote: > > > On Fri, Jan 5, 2018 at 8:39 AM, Viviana Palacio Betancur < > vpalac...@uchicago.edu> 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 Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users