Add the possibility to use a multiplier to allow Robin type of lumped…

… conditions when having an integral constraint.
ElmerCSC · Aug 29, 2022 · 952c012 · 952c012
1 parent e328d86
commit 952c012
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/fem/src/SolverUtils.F90 b/fem/src/SolverUtils.F90
@@ -5290,6 +5290,12 @@ SUBROUTINE SetDirichletBoundaries( Model, A, b, Name, DOF, NDOFs, Perm, &
             t = Offset + Perm(ind)
             b(t) = b(t) + SingleVal
           END IF
+
+          SingleVal = ListGetCReal( ValueList,TRIM(DirName)//' Multiplier',GotIt)
+          IF( GotIt ) THEN                        
+            t = Offset + Perm(ind)
+            CALL AddToMatrixElement(A,t,t,SingleVal) 
+          END IF
         END IF
 
         n = COUNT( LumpedNodeSet )