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

2018-01-05 Thread John Peterson
On Fri, Jan 5, 2018 at 10:02 AM, Viviana Palacio Betancur < [email protected]> wrote: > 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. > OK, the normal way this would be handled wo

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

2018-01-05 Thread Viviana Palacio Betancur
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, 201

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

2018-01-05 Thread John Peterson
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 avoi

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

2018-01-05 Thread Viviana Palacio Betancur
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. Best, Viviana. Viviana Palacio Betancur PhD Student | 2016 Cohort de Pablo Grou

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

2018-01-04 Thread John Peterson
On Thu, Jan 4, 2018 at 9:16 AM, Viviana Palacio Betancur < [email protected]> wrote: > Hello, > > I am solving a system on a mesh with two subdomains: bulk with elements > TET10, and surface with elements TRI6. So far, I've sucessfully asembled > the system's matrix and rhs, and initialized t

[Libmesh-users] Solving a system with 2 subdomains

2018-01-04 Thread Viviana Palacio Betancur
Hello, I am solving a system on a mesh with two subdomains: bulk with elements TET10, and surface with elements TRI6. So far, I've sucessfully asembled the system's matrix and rhs, and initialized the system. My problem is when I execute system.solve(). Apparently, the solution is only done for t