Re: [Libmesh-users] Problem with parallel generalized eigenproblem

2016-03-19 Thread David Knezevic
On Thu, Mar 17, 2016 at 3:31 PM, Roy Stogner wrote: > > On Thu, 17 Mar 2016, David Knezevic wrote: > > Do you mean just condensing out all constrained dofs? I use that already >> with DirichletBoundary and that works fine. And I think it would work >> automatically with AMR too. >> > > Would it m

Re: [Libmesh-users] Problem with parallel generalized eigenproblem

2016-03-19 Thread Roy Stogner
On Thu, 17 Mar 2016, Torquil Macdonald Sørensen wrote: I have attached the program that produces these results I'm afraid our mailing lists aren't great about attachments. Is the program short enough to post as message text or as a gist? --- Roy---

Re: [Libmesh-users] Problem with parallel generalized eigenproblem

2016-03-19 Thread Torquil Macdonald Sørensen
On 17/03/16 19:06, Roy Stogner wrote: > > On Thu, 17 Mar 2016, Torquil Macdonald Sørensen wrote: > >> I have attached the program that produces these results > > I'm afraid our mailing lists aren't great about attachments. Is the > program short enough to post as message text or as a gist? > --- >

Re: [Libmesh-users] Problem with parallel generalized eigenproblem

2016-03-19 Thread David Knezevic
On Thu, Mar 17, 2016 at 2:06 PM, Roy Stogner wrote: > > On Thu, 17 Mar 2016, Torquil Macdonald Sørensen wrote: > > I have attached the program that produces these results >> > > I'm afraid our mailing lists aren't great about attachments. Is the > program short enough to post as message text or

Re: [Libmesh-users] Problem with parallel generalized eigenproblem

2016-03-19 Thread Roy Stogner
On Thu, 17 Mar 2016, David Knezevic wrote: > Torquil: Once this change is in, it'd be good if you can modify your > example to use DirichletBoundary and re-run it to check if it fixes > your problem. If you could try that out *before* the change is in, that would be helpful. https://github.com/

Re: [Libmesh-users] Problem with parallel generalized eigenproblem

2016-03-19 Thread David Knezevic
On Thu, Mar 17, 2016 at 4:15 PM, Roy Stogner wrote: > > On Thu, 17 Mar 2016, David Knezevic wrote: > > This would add all constrained dofs to the list of condensed dofs. >> That should do what we want, I agree. >> > > I'll put in a PR now. > > Looks like a few other bits to fix while I'm at it; e

Re: [Libmesh-users] Problem with parallel generalized eigenproblem

2016-03-19 Thread David Knezevic
On Thu, Mar 17, 2016 at 2:57 PM, Roy Stogner wrote: > > On Thu, 17 Mar 2016, David Knezevic wrote: > > I'd be interested to see the code, if you can send it through. In >> particular, I'd be interested to see how you're imposing the >> boundary conditions. Are you using CondensedEigenSystem? If s

Re: [Libmesh-users] Problem with parallel generalized eigenproblem

2016-03-19 Thread Torquil Macdonald Sørensen
On 17/03/16 21:48, Roy Stogner wrote: > > On Thu, 17 Mar 2016, David Knezevic wrote: > >> Torquil: Once this change is in, it'd be good if you can modify your >> example to use DirichletBoundary and re-run it to check if it fixes >> your problem. > > If you could try that out *before* the change is

Re: [Libmesh-users] Problem with parallel generalized eigenproblem

2016-03-19 Thread David Knezevic
On Fri, Mar 18, 2016 at 10:41 AM, Torquil Macdonald Sørensen < [email protected]> wrote: > On 17/03/16 21:48, Roy Stogner wrote: > > > > On Thu, 17 Mar 2016, David Knezevic wrote: > > > >> Torquil: Once this change is in, it'd be good if you can modify your > >> example to use DirichletBoundary an

Re: [Libmesh-users] Problem with parallel generalized eigenproblem

2016-03-19 Thread Roy Stogner
On Thu, 17 Mar 2016, David Knezevic wrote: > Do you mean just condensing out all constrained dofs? I use that already with > DirichletBoundary and that works fine. And I think it would work > automatically with AMR too. Would it make sense to just add that to initialize_condensed_dofs()? Just

Re: [Libmesh-users] Problem with parallel generalized eigenproblem

2016-03-19 Thread David Knezevic
Torquil, I ran your example with Roy's PR, and it worked fine for me for 1 up to 8 MPI processes. Note that there are some issues with your example: - You need to use CondensedEigenSystem, not EigenSystem. You used CondensedEigenSystem in the first version of your example, I'm not sure why you sw

Re: [Libmesh-users] Problem with parallel generalized eigenproblem

2016-03-18 Thread Roy Stogner
On Thu, 17 Mar 2016, David Knezevic wrote: > I'd be interested to see the code, if you can send it through. In > particular, I'd be interested to see how you're imposing the > boundary conditions. Are you using CondensedEigenSystem? If so, > perhaps there is an issue with how the "indices to cond

Re: [Libmesh-users] Problem with parallel generalized eigenproblem

2016-03-18 Thread Roy Stogner
On Thu, 17 Mar 2016, David Knezevic wrote: > This would add all constrained dofs to the list of condensed dofs. > That should do what we want, I agree. I'll put in a PR now. Looks like a few other bits to fix while I'm at it; e.g. unsigned int instead of dof_id_type. > Would you want to get ri

Re: [Libmesh-users] Problem with parallel generalized eigenproblem

2016-03-18 Thread Torquil Macdonald Sørensen
Hi! Yes, it now works fine for me too. Thanks! -Torquil On 18 Mar 2016 22:40, "David Knezevic" wrote: > Torquil, > > I ran your example with Roy's PR, and it worked fine for me for 1 up to 8 > MPI processes. Note that there are some issues with your example: > > - You need to use CondensedEige