Change swapping of swappable BCs

The swapping should happen at the beginning of a controller
step before the Schwarz loop for a multidomain simulation.

N/A for a single domain simulation.

src/ics_bcs.jl

@@ -381,9 +381,6 @@ function apply_bc_detail(model::SolidMechanics, bc::SMContactSchwarzBC)
     else
         apply_sm_schwarz_contact_neumann(model, bc)
     end
-    if (bc.swap_bcs == true)
-        bc.is_dirichlet = !bc.is_dirichlet
-    end
 end
 
 function apply_bc_detail(model::SolidMechanics, bc::CouplingSchwarzBoundaryCondition)
@@ -392,9 +389,6 @@ function apply_bc_detail(model::SolidMechanics, bc::CouplingSchwarzBoundaryCondi
     else
         apply_sm_schwarz_coupling_neumann(model, bc)
     end
-    if (bc.swap_bcs == true)
-        bc.is_dirichlet = !bc.is_dirichlet
-    end
 end
 
 function apply_bc_detail(model::LinearOpInfRom, bc::CouplingSchwarzBoundaryCondition)

src/schwarz.jl

@@ -110,6 +110,7 @@ function schwarz(sim::MultiDomainSimulation)
     save_schwarz_solutions(sim)
     resize_histories(sim)
     set_subcycle_times(sim)
+    swap_swappable_bcs(sim)
     is_schwarz = true
     while true
         println("Schwarz iteration=", iteration_number)
@@ -189,6 +190,22 @@ function set_subcycle_times(sim::MultiDomainSimulation)
     end
 end
 
+function swap_swappable_bcs(sim::MultiDomainSimulation)
+    for subsim ∈ sim.subsims
+        swap_swappable_bcs(subsim)
+    end
+end
+
+function swap_swappable_bcs(sim::SingleDomainSimulation)
+    for bc ∈ sim.model.boundary_conditions
+        if typeof(bc) == ContactSchwarzBoundaryCondition || typeof(bc) == CouplingSchwarzBoundaryCondition
+            if (bc.swap_bcs == true)
+                bc.is_dirichlet = !bc.is_dirichlet
+            end
+        end
+    end
+end
+
 function subcycle(sim::MultiDomainSimulation, is_schwarz::Bool)
     subsim_index = 1
     for subsim ∈ sim.subsims