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

[John P. <jwp...@gm...>]

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