phynexis-deps/mfem
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fixed bug for H1->RT

Browse Source
This commit is contained in:
Andrew Ho
2026-08-18 10:35:32 -07:00
parent c11172b842
commit fbd217e8a4
2 changed files with 3 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files
fem/integ/bilininteg_mixedvecgrad_pa.cpp
+2 -2
View File
@@ -1227,7 +1227,7 @@ static void PAHdivH1Apply3D(const int D1D,
auto Bct = Reshape(bct.Read(), D1D, Q1D);
auto op = Reshape(pa_data.Read(), Q1D, Q1D, Q1D, op_entries, NE);
auto X = Reshape(x.Read(), D1D, D1D, D1D, NE);
auto Y = Reshape(y.ReadWrite(), 3*(D1D-1)*D1D*D1D, NE);
auto Y = Reshape(y.ReadWrite(), 3*D1D*(D1D-1)*(D1D-1), NE);
mfem::forall(NE, [=] MFEM_HOST_DEVICE (int e)
{
@@ -1428,7 +1428,7 @@ static void PAHdivH1ApplyTranspose3D(const int D1D,
auto Bt = Reshape(bct.Read(), D1D, Q1D);
auto Gt = Reshape(gct.Read(), D1D, Q1D);
auto op = Reshape(pa_data.Read(), Q1D, Q1D, Q1D, op_entries, NE);
auto X = Reshape(x.Read(), 3*(D1D-1)*D1D*D1D, NE);
auto X = Reshape(x.Read(), 3*(D1D-1)*(D1D-1)*D1D, NE);
auto Y = Reshape(y.ReadWrite(), D1D, D1D, D1D, NE);
mfem::forall(NE, [=] MFEM_HOST_DEVICE (int e)
tests/unit/fem/test_pa_coeff.cpp
+1 -1
View File
@@ -762,7 +762,7 @@ TEST_CASE("Hcurl/Hdiv MixedVectorGradientPA",
auto vFEType = GENERATE(0, 1);
CAPTURE(dimension, coeffType, order, vFEType);
const int ne = 1;
const int ne = 3;
Mesh mesh = MakeCartesianNonaligned(dimension, ne);
H1_FECollection scalar_fec(order, dimension);