Re: [Libmesh-users] Solve on part of domain

2010-09-14 Thread Tim Kroeger
On Wed, 15 Sep 2010, Tim Kroeger wrote: > Just to keep this clear: NumericVector::solve_only_on() will only > accept a list of dofs, nothing else, because in particular a > NumericVector does not know anything about subdomain ids. On the > other hand, System::solve_only_on() will only accept (a r

Re: [Libmesh-users] Solve on part of domain

2010-09-14 Thread Tim Kroeger
On Tue, 14 Sep 2010, Roy Stogner wrote: > On Tue, 14 Sep 2010, Tim Kroeger wrote: > >> I suggest to make two methods, one of which takes a vector and the other >> takes a set. Then, the set-taking method will by default just convert the >> set to a vector and call the vector-taking method. Roy

Re: [Libmesh-users] Largest subdomain_id around a dof

2010-09-14 Thread Tim Kroeger
On Tue, 14 Sep 2010, Roy Stogner wrote: > On Tue, 14 Sep 2010, Tim Kroeger wrote: > >> For a dof of a vector the corresponds to a System, I need to know the >> maximal subdomain_id of all elements which contain this dof. > >> While it is principally possible to determine this on the fly (via >> El

Re: [Libmesh-users] Largest subdomain_id around a dof

2010-09-14 Thread Roy Stogner
On Tue, 14 Sep 2010, Tim Kroeger wrote: > For a dof of a vector the corresponds to a System, I need to know the > maximal subdomain_id of all elements which contain this dof. > While it is principally possible to determine this on the fly (via > Elem::find_point_neighbor() and then a loop over t

Re: [Libmesh-users] Solve on part of domain

2010-09-14 Thread Roy Stogner
On Tue, 14 Sep 2010, Tim Kroeger wrote: > I suggest to make two methods, one of which takes a vector and the other > takes a set. Then, the set-taking method will by default just convert the > set to a vector and call the vector-taking method. Roy, what do you think? I think this is one of t

[Libmesh-users] Largest subdomain_id around a dof

2010-09-14 Thread Tim Kroeger
Dear all, For a dof of a vector the corresponds to a System, I need to know the maximal subdomain_id of all elements which contain this dof. While it is principally possible to determine this on the fly (via Elem::find_point_neighbor() and then a loop over the resulting elements and a loop ov

Re: [Libmesh-users] Solve on part of domain

2010-09-14 Thread Tim Kroeger
Dear Jed, Thank you for your help so far. I still have some questions: On Mon, 13 Sep 2010, Jed Brown wrote: > On Mon, 13 Sep 2010 15:18:30 +0200 (CEST), Tim Kroeger > wrote: > >> + >> + /** >> + * Set of dofs to which any solve is limited (\p NULL means solve on >> + * all dofs). >> +