| 72f91173 | 15-Jan-2023 |
Jed Brown <jed@jedbrown.org> |
DMPlex tests/ex5: update tri and tet numerical values for prior code
This commit is in preparation for a change to the quadrature that will
not affect these numerical values. I bisected on ex5_tet_0
DMPlex tests/ex5: update tri and tet numerical values for prior code
This commit is in preparation for a change to the quadrature that will
not affect these numerical values. I bisected on ex5_tet_0 to find two
commits that introduced numerical changes and summarize their impact
here. I'm expecting the other changes to be similar.
The first numerical change was caused by
commit e1a13dae02585224ddd4368ee030f3c5e77a306b
Author: Matthew Knepley <knepley@gmail.com>
Date: Thu Apr 28 10:27:14 2022 -0400
Plex: Fix fault face division
- Should not split a face that only touches the fault on unsplit points
src/dm/impls/plex/plexsubmesh.c | 20 +++++++++++++++++---
1 file changed, 17 insertions(+), 3 deletions(-)
which has the numerical change:
diff --git 1/base.out 2/here.out
index df7196b66a9..9534c0139ad 100644
--- 1/base.out
+++ 2/here.out
@@ -150,15 +150,15 @@ Vec Object: Local Solution 1 MPI process
0.
1.5
1.
-0.
-1.
-0.
0.5
0.
0.
-0.5
0.
0.
+0.
+1.
+0.
Discrete System with 2 fields
cell total dim 27 total comp 6
cohesive cell
@@ -222,23 +222,23 @@ Vec Object: Local Residual 1 MPI process
0.
0.
0.
--3.46945e-18
-0.0833333
-0.0833333
0.0208333
-0.0625
+0.0208333
0.0416667
-0.0208333
-0.0208333
0.0416667
+0.0833333
0.
+0.104167
+0.0416667
0.
+-1.34017e-17
0.
0.
+-9.26063e-18
0.
0.
-0.
-0.
+7.55411e-18
0.
Mat Object: Jacobian 1 MPI process
type: seqaij
The second commit
commit d6177c40aeff79f5090ba5db2353273bea8a20ec
Author: Toby Isaac <toby.isaac@gmail.com>
Date: Wed Jun 29 10:29:10 2022 -0400
Use DMLabelStratumHasPoint() in CheckPoint_Private()
plex test ex61 shows there are corner cases where checking the labelId is necessary
src/dm/impls/plex/plex.c | 12 ++++--------
src/dm/impls/plex/tests/ex61.c | 1 -
2 files changed, 4 insertions(+), 9 deletions(-)
which had the numerical change
diff --git 1/base2.out 2/here.out
index 9534c0139ad..5af4b0229e2 100644
--- 1/base2.out
+++ 2/here.out
@@ -142,14 +142,14 @@ Vec Object: Local Solution 1 MPI process
1.
1.5
0.
-0.5
-1.
-1.
-1.
-1.
0.
-1.5
-1.
+0.
+0.
+0.
+0.
+0.
+0.
+0.
0.5
0.
0.
@@ -231,15 +231,15 @@ Vec Object: Local Residual 1 MPI process
0.
0.104167
0.0416667
-0.
--1.34017e-17
-0.
-0.
--9.26063e-18
-0.
-0.
-7.55411e-18
-0.
+0.0416667
+0.145833
+0.166667
+0.0833333
+0.166667
+0.166667
+0.0416667
+0.1875
+0.166667
Mat Object: Jacobian 1 MPI process
type: seqaij
row 0: (0, 0.) (1, 0.) (2, 0.) (3, 0.) (4, 0.) (5, 0.) (6, 0.) (7, 0.) (8, 0.) (9, 0.) (10, 0.) (11, 0.)
show more ...
|