Further steps. Add test case.

Author: raback

fem/src/MortarUtils.F90

@@ -3379,7 +3379,6 @@ FUNCTION LevelProjector( BMesh1, BMesh2, Repeating, AntiRepeating, &
       ArcCoeff = 1.0_dp
     END IF
 
-
 
     ! We have a weak projector if it is requested 
     NitscheProjector = ListGetLogical( BC,'Nitsche Projector', Found )     
@@ -3454,7 +3453,6 @@ FUNCTION LevelProjector( BMesh1, BMesh2, Repeating, AntiRepeating, &
     Projector % FORMAT = MATRIX_LIST
 
     IF( NitscheProjector ) THEN
-
       Projector % ProjectorType = PROJECTOR_TYPE_NITSCHE
     ELSE      
       Projector % ProjectorType = PROJECTOR_TYPE_GALERKIN
@@ -3769,9 +3767,11 @@ FUNCTION LevelProjector( BMesh1, BMesh2, Repeating, AntiRepeating, &
     END IF
 
     CALL Info(Caller,'Max number of rows in projector: '//I2S(ProjectorRows),Level=10)
-    ALLOCATE( Projector % InvPerm(ProjectorRows) )
-    Projector % InvPerm = 0
-
+    IF(.NOT. NitscheProjector ) THEN     
+      ALLOCATE( Projector % InvPerm(ProjectorRows) )
+      Projector % InvPerm = 0
+    END IF
+      
     ! If after strong projectors there are still something undone they must 
     ! be dealt with the weak projectors. 
     SomethingUndone = .FALSE.
@@ -3860,11 +3860,13 @@ FUNCTION LevelProjector( BMesh1, BMesh2, Repeating, AntiRepeating, &
     IF( CreateDual ) DEALLOCATE( DualNodePerm )
     IF( DoEdges ) DEALLOCATE( EdgePerm )
 
-    m = COUNT( Projector % InvPerm  == 0 ) 
-    IF( m > 0 ) THEN
-      CALL Warn(Caller,'Projector % InvPerm not set in for dofs: '//I2S(m))
+    IF(ASSOCIATED(Projector % InvPerm) ) THEN     
+      m = COUNT( Projector % InvPerm  == 0 ) 
+      IF( m > 0 ) THEN
+        CALL Warn(Caller,'Projector % InvPerm not set in for dofs: '//I2S(m))
+      END IF
     END IF
-
+      
     CALL Info(Caller,'Projector created',Level=10)
 
   CONTAINS
@@ -5480,8 +5482,7 @@ SUBROUTINE AddProjectorWeakGeneric()
 
       IF( SaveErr ) CLOSE(11)
 
-      
-      
+           
       DEALLOCATE( Nodes % x, Nodes % y, Nodes % z, &
           NodesM % x, NodesM % y, NodesM % z, &
           NodesT % x, NodesT % y, NodesT % z, &

fem/src/SolverUtils.F90

@@ -26034,9 +26034,6 @@ SUBROUTINE GenerateAddMatrix( Model, Solver )
        END IF
 
        CALL Info(Caller,'Adding mortar projector of type Nitsche for BC: '//I2S(bc_ind),Level=8)
-       IF( .NOT. ASSOCIATED( Atmp % InvPerm ) ) THEN
-         CALL Fatal(Caller,'InvPerm is required for geometric projector!')
-       END IF
 
        ! Enable that the user can for vector valued cases either set some 
        ! or skip some field components. 

fem/tests/ContactBlunt3DNormalProj/case.sif

@@ -20,7 +20,7 @@ Header
 End
 
 Simulation
-  Max Output Level = 5
+  Max Output Level = 25
   Coordinate System = Cartesian 3D
   Coordinate Mapping(3) = 1 2 3
   Simulation Type = Steady State

fem/tests/RotatingBCMagnetoDynamics2/case.sif

@@ -21,7 +21,7 @@ Header
 End
 
 Simulation
-  Max Output Level = 5
+  Max Output Level = 25
   Coordinate System = Cartesian
   Coordinate Mapping(3) = 1 2 3
   Simulation Type = Transient

fem/tests/RotatingBCPoissonNitsche/CMakeLists.txt

@@ -0,0 +1,8 @@
+INCLUDE(test_macros)
+INCLUDE_DIRECTORIES(${CMAKE_BINARY_DIR}/fem/src)
+
+CONFIGURE_FILE( case.sif case.sif COPYONLY)
+
+file(COPY ELMERSOLVER_STARTINFO mortar.msh DESTINATION "${CMAKE_CURRENT_BINARY_DIR}/")
+
+ADD_ELMER_TEST(RotatingBCPoisson2DNitsche LABELS mortar)

fem/tests/RotatingBCPoissonNitsche/ELMERSOLVER_STARTINFO

@@ -0,0 +1 @@
+case.sif

fem/tests/RotatingBCPoisson2D/nitsche.sif …/tests/RotatingBCPoissonNitsche/case.siffem/tests/RotatingBCPoisson2D/nitsche.sif renamed to fem/tests/RotatingBCPoissonNitsche/case.sif fem/tests/RotatingBCPoisson2D/nitsche.sif renamed to fem/tests/RotatingBCPoissonNitsche/case.sif

@@ -1,5 +1,8 @@
 ! Test case for 2D nonconforming mortar elements
-! Some timing settings inhereted from benchmark computations
+! This test uses the Nitshce method to ensure continuity over the
+! interface.
+!
+! P.R. 26.2.2026
 
 Header
   CHECK KEYWORDS Warn
@@ -18,7 +21,7 @@ Simulation
   Steady State Max Iterations = 1
   Output Intervals = 1
 
-  Timestep Intervals = 20
+  Timestep Intervals = 3 !90
   Output Intervals = 1
 
 !  signs = Variable "time"
@@ -27,12 +30,9 @@ Simulation
 !      2 1 
 !    End 
 
-  Post File = case-a.vtu
-
+  Post File = case.vtu
   Ascii Output = True
 
-!  Output File = case.result
-
   Simulation Timing = Logical True
 End 
 
@@ -161,8 +161,6 @@ Body Force 3
     Real MATC "omega*tx+phi0"
 End 
 
-
-
 Boundary Condition 1
   Target Boundaries(4) = 10 11 12 13
   Name = "Ambient"
@@ -187,5 +185,5 @@ Boundary Condition 3
 End
 
 
-Solver 2 :: Reference Norm = 6.93769175
+Solver 2 :: Reference Norm = 7.3154459807490548