Some corrections related to bubbles

ElmerCSC · Jan 25, 2025 · 1903674 · 1903674
1 parent 8cb7258
commit 1903674
Show file tree

Hide file tree

Showing 2 changed files with 25 additions and 15 deletions.
diff --git a/fem/src/DefUtils.F90 b/fem/src/DefUtils.F90
@@ -2171,10 +2171,12 @@ FUNCTION GetElementNOFDOFs( UElement,USolver ) RESULT(n)
              IF (p > 1) BDOFs = GetBubbleDOFs(Element, p)
              BDOFs = MAX(nb, BDOFs)
            ELSE
-             Bubbles = ListGetLogical(Solver % Values, 'Bubbles', Found )
-             ! The following is not a right way to obtain the bubble count
-             ! in order to support solverwise definitions
-             IF (Bubbles) BDOFs = SIZE(Element % BubbleIndexes)
+             IF (ASSOCIATED(Solver % Values)) THEN
+               Bubbles = ListGetLogical(Solver % Values, 'Bubbles', Found )
+               ! The following is not a right way to obtain the bubble count
+               ! in order to support solverwise definitions
+               IF (Bubbles) BDOFs = SIZE(Element % BubbleIndexes)
+             END IF
            END IF
            n = n + BDOFs
          END IF
@@ -2247,10 +2249,12 @@ FUNCTION GetElementNOFBDOFs( Element, USolver, Update ) RESULT(n)
       ELSE
         ! The element command hasn't been given, so the only way to activate the bubbles
         ! should be the "Bubbles" command. The following is not a reliable way to obtain
-        ! the bubble count when solverwise definitions are used.   
-        Bubbles = ListGetLogical(Solver % Values, 'Bubbles', Found)
-        IF (Bubbles .AND. ASSOCIATED(Element % BubbleIndexes)) THEN
-          n = SIZE(Element % BubbleIndexes)
+        ! the bubble count when solverwise definitions are used.
+        IF (ASSOCIATED(Solver % Values)) THEN
+          Bubbles = ListGetLogical(Solver % Values, 'Bubbles', Found)
+          IF (Bubbles .AND. ASSOCIATED(CurrElement % BubbleIndexes)) THEN
+            n = SIZE(CurrElement % BubbleIndexes)
+          END IF
         END IF
       END IF
 
@@ -2272,9 +2276,13 @@ FUNCTION GetElementNOFBDOFs( Element, USolver, Update ) RESULT(n)
           n = MAX(k,n)
           CurrElement % BDOFs = n
         ELSE
-          Bubbles = ListGetLogical(Solver % Values, 'Bubbles', Found)
-          IF (Bubbles .AND. ASSOCIATED(Element % BubbleIndexes)) THEN
-            CurrElement % BDOFs = SIZE(Element % BubbleIndexes)
+          IF (ASSOCIATED(Solver % Values)) THEN
+            Bubbles = ListGetLogical(Solver % Values, 'Bubbles', Found)
+            IF (Bubbles .AND. ASSOCIATED(CurrElement % BubbleIndexes)) THEN
+              CurrElement % BDOFs = SIZE(CurrElement % BubbleIndexes)
+            ELSE
+              CurrElement % BDOFs = 0
+            END IF
           ELSE
             CurrElement % BDOFs = 0
           END IF

diff --git a/fem/src/ElementUtils.F90 b/fem/src/ElementUtils.F90
@@ -4300,10 +4300,12 @@ FUNCTION mGetElementDOFs( Indexes, UElement, USolver, NotDG, UMesh )  RESULT(nd)
            IF (p > 1) BDOFs = GetBubbleDOFs(Element, p)
            BDOFs = MAX(nb, BDOFs)
          ELSE
-           Bubbles = ListGetLogical(Solver % Values, 'Bubbles', Found )
-           ! The following is not a right way to obtain the bubble count
-           ! in order to support solverwise definitions
-           IF (Bubbles) BDOFs = SIZE(Element % BubbleIndexes)
+           IF (ASSOCIATED(Solver % Values)) THEN
+             Bubbles = ListGetLogical(Solver % Values, 'Bubbles', Found )
+             ! The following is not a right way to obtain the bubble count
+             ! in order to support solverwise definitions
+             IF (Bubbles) BDOFs = SIZE(Element % BubbleIndexes)
+           END IF
          END IF
          DO i=1,BDOFs
            nd = nd + 1