Constrained linear operator on distributed data structures

[Doug Shi-Dong]

Hello everyone,

I see that some functions within the constrained_linear_operator.h have a note saying

Currently, this function may not work correctly for distributed data structures.

I am currently using it in parallel and seems to be working fine for a very simple linear elasticity problem with hanging nodes and Dirichlet boundary condition. I was wondering what I should be careful about. In which scenarios would those functions fail?

I couldn't find any issues on GitHub regarding those, or any mention of a specific problem when git blaming.

Happy holidays!

Doug

[David Wells]

Hi Doug,

If I recall correctly, we added that comment since this is a part of
the library that is not very well tested. The code is probably fine.

Thanks,
David

> --
> The deal.II project is located at http://www.dealii.org/
> For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en
> ---
> You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an email to dealii+un...@googlegroups.com.
> To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/3d00528c-a107-4594-9068-61c7532dfb75%40googlegroups.com.

[Doug Shi-Dong]

Hello Dr. Wells,

Funnily enough, I've got most of it working and it seems to work great in parallel, with Dirichlet  boundary conditions (inhomogeneous constraints) and hanging nodes (the linear operator part of the constraint). In serial, it works fine with the Dirichlet BC.

The one scenario where it doesn't work well for me is actually in serial with hanging nodes. 

I'll report back here with the issue in the upcoming month or on Github with an issue if I can pinpoint where it comes from. 

Doug

On Tuesday, December 31, 2019 at 3:47:23 PM UTC-5, David Wells wrote:

Hi Doug,

If I recall correctly, we added that comment since this is a part of
the library that is not very well tested. The code is probably fine.

Thanks,
David

On Tue, Dec 24, 2019 at 2:33 PM Doug Shi-Dong <doug....@gmail.com> wrote:
>
> Hello everyone,
>
> I see that some functions within the constrained_linear_operator.h have a note saying
>
> Currently, this function may not work correctly for distributed data structures.
>
> I am currently using it in parallel and seems to be working fine for a very simple linear elasticity problem with hanging nodes and Dirichlet boundary condition. I was wondering what I should be careful about. In which scenarios would those functions fail?
>
> I couldn't find any issues on GitHub regarding those, or any mention of a specific problem when git blaming.
>
> Happy holidays!
>
> Doug
>
> --
> The deal.II project is located at http://www.dealii.org/
> For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en
> ---
> You received this message because you are subscribed to the Google Groups "deal.II User Group" group.

> To unsubscribe from this group and stop receiving emails from it, send an email to dea...@googlegroups.com.