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

Reply via email to