Hi!

Yes, it now works fine for me too.

Thanks!

-Torquil
On 18 Mar 2016 22:40, "David Knezevic" <david.kneze...@akselos.com> 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 CondensedEigenSystem, not EigenSystem. You used
> CondensedEigenSystem in the first version of your example, I'm not sure why
> you switched that?
> - You need to call eigen_system.initialize_condensed_dofs(); after
> equation_systems.init(), as mentioned in my previous email.
>
> After making those changes, it worked fine for me. I have merged Roy's PR
> now. Let me know if you have any further issues with this.
>
> Thanks,
> David
>
>
>
> On Fri, Mar 18, 2016 at 10:55 AM, David Knezevic <
> david.kneze...@akselos.com> wrote:
>
>> On Fri, Mar 18, 2016 at 10:41 AM, Torquil Macdonald Sørensen <
>> torq...@gmail.com> 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 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/libMesh/libmesh/pull/872
>>> >
>>> > I don't know if there's any MOOSE test coverage for
>>> > CondensedEigenSystem and libMesh only has the one example IIRC.
>>> > ---
>>> > Roy
>>>
>>> Hi!
>>>
>>> After merging the change, and changing to use DirichletBoundary, the
>>> results seem to be independent of the number of MPI nodes, as one would
>>> expect.
>>>
>>> The output of my program is now:
>>>
>>> 1 1 1 1 1 1 1 1 1 1 1 1 1 32.9925 73.0811 73.0811 73.0811 120.25 120.25
>>> 120.25 175.419 175.731
>>>
>>> where the unit values seem to be spurious eigenvalues. Is there a
>>> convenient way in libmesh for projecting out those parts of the vector
>>> space?
>>>
>>> The program is available here:
>>>
>>>
>>> https://drive.google.com/file/d/0B5vWRGzG5kNlbHhWN1M3WHVDeU0/view?usp=sharing
>>>
>>> Best regards,
>>> Torquil Macdonald Sørensen
>>>
>>
>>
>> It sounds like the Dirichlet dofs aren't getting condensed out, since
>> they produce these spurious eigenvalues. You need to
>> call initialize_condensed_dofs() (no argument required) after
>> equation_systems.init().
>>
>> David
>>
>>
>>
>
------------------------------------------------------------------------------
Transform Data into Opportunity.
Accelerate data analysis in your applications with
Intel Data Analytics Acceleration Library.
Click to learn more.
http://pubads.g.doubleclick.net/gampad/clk?id=278785231&iu=/4140
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to