Re: [Libmesh-users] Solving a system with 2 subdomains

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

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 <jwp...@gm...> wrote:

>
>
> On Fri, Jan 5, 2018 at 8:39 AM, Viviana Palacio Betancur <
> vpa...@uc...> 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
>