Add some debug logs for periodic examples

Author: pjaap

examples/Example212_PeriodicElasticity2D.jl

@@ -125,6 +125,9 @@ function main(;
         threads = 1,
         kwargs...
     )
+    ## print options for better logs
+    @info "selected options" use_LM_restrictions threads
+
     xgrid = create_grid(; h, width, height)
 
     ## create finite element space and solution vector
@@ -166,6 +169,10 @@ function main(;
     sol, SC = solve(PD, FES, return_config = true)
     residual(SC) < 1.0e-10 || error("Residual is not zero!")
 
+    if use_LM_restrictions
+        @info "Lagrange residuals" SC.statistics[:restriction_residuals]
+    end
+
     plt = GridVisualizer(; Plotter, size = (1300, 800))
 
     magnification = 1

examples/Example312_PeriodicElasticity3D.jl

@@ -139,6 +139,8 @@ function main(;
         threads = 1,
         kwargs...
     )
+    ## print options for better logs
+    @info "selected options" periodic_coupling threads
 
     xgrid = create_grid(; h, width, height, depth)
 
@@ -187,6 +189,10 @@ function main(;
     sol, SC = solve(PD, FES; return_config = true, kwargs...)
     residual(SC) < 1.0e-10 || error("Residual is not zero!")
 
+    if use_LM_restrictions
+        @info "Lagrange residuals" SC.statistics[:restriction_residuals]
+    end
+
     displace_mesh!(xgrid, sol[u])
     plt = plot([grid(u)], sol; Plotter, do_vector_plots = false, width = 1200, height = 800, title = "displaced mesh", scene3d = :LScene)
 

examples/Example313_PeriodicPoisson.jl

@@ -88,6 +88,8 @@ function main(;
         periodic = true,
         kwargs...
     )
+    ## print options for better logs
+    @info "selected options" periodic
 
     xgrid = create_grid(; h)
 
@@ -112,6 +114,8 @@ function main(;
     sol, SC = solve(PD, FES; return_config = true, kwargs...)
     residual(SC) < 1.0e-10 || error("Residual is not zero!")
 
+    @info "Lagrange residuals" SC.statistics[:restriction_residuals]
+
     function exact_error!(out, u, qpinfo)
         # exact solution here is u(x,y,z) = 2z - 1
         val = qpinfo.x[3] * 2 - 1