Skip to content

ConstantFace fix broke FEM Nonhomogeneous Dirichlet BCs #551

@s-mayani

Description

@s-mayani

The Constant Face ix in commit a791690 (a791690) broke the FEM test for non-homogeneous BCs.

Before, running the TestNonhomDirichlet_2d test on 2 ranks would give:

{0}>       Size                  Spacing           Relative Error               Norm Error                  Residue     Iterations
{0}>          8       0.1428571428571428     0.001783276015931643    0.0002991477913988007    5.728386296687791e-13             17
{0}>         16      0.06666666666666667    0.0003953108198504651    6.977841205835258e-05    1.482933104244011e-12             38
{0}>         32      0.03225806451612903    9.290220252853989e-05    1.687185626847187e-05    1.635197199953021e-12             80
{0}>         64      0.01587301587301587    2.251359104334525e-05     4.14962982155807e-06    2.287139014011642e-12            161
{0}>        128     0.007874015748031496    5.541266888295135e-06    1.029061875381498e-06    3.959959305127847e-12            322
{0}>        256      0.00392156862745098    1.374543362823382e-06    2.562334274879125e-07    6.147130983156806e-12            637

After that commit, running the same test on 2 ranks gives:

{0}>       Size                  Spacing           Relative Error               Norm Error                  Residue     Iterations
{0}>          8       0.1428571428571428      0.07910225330188655      0.04827256219720535    3.080196504432488e-13             18
{0}>         16      0.06666666666666667      0.04473324400570541      0.02846813739570635    7.782857367937769e-13             40
{0}>         32      0.03225806451612903      0.02476940904838091      0.01604914107553235    1.413308950149204e-12             81
{0}>         64      0.01587301587301587      0.01352901535756036     0.008823145596464899    3.393894984935643e-12            161
{0}>        128     0.007874015748031496     0.007310660872388467     0.004773738533200396    4.106713412817932e-12            322
{0}>        256      0.00392156862745098     0.003915541587443943     0.002554252712022781    6.153892356801769e-12            637

Single rank is working fine, and 1D also. The problem arises on multi-rank for 2D and 3D (already on CPU).

Metadata

Metadata

Assignees

Labels

No labels
No labels

Type

No type
No fields configured for issues without a type.

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions