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
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---
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?
> ---
>
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
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/
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
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
Hi!
I'm implementing a generalized eigenvalue problem in libmesh with SLEPc,
and noticed a problem when running with certain numbers of MPI nodes. I
would get some erroneous eigenvalues in those cases. This seems to be
independent of PETSc/SLEPc settings, because I have reproduced the
problem with
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
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
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
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
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
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
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
15 matches
Mail list logo
