Merge branch 'master' of https://gitlab.onelab.info/doc/models

DiffractionGratings/grating3D.pro

@@ -500,6 +500,15 @@ PostOperation {
 
 DefineConstant[
   R_ = {"helmholtz_vector", Name "GetDP/1ResolutionChoices", Visible 1},
-  C_ = {"-solve -pos -petsc_prealloc 500 -ksp_type preonly -pc_type lu -pc_factor_mat_solver_type mumps", Name "GetDP/9ComputeCommand", Visible 1},
+  C_ = {"-solve -pos -petsc_prealloc 500 -ksp_type preonly -pc_type lu -pc_factor_mat_solver_type mumps -ksp_error_if_not_converged", Name "GetDP/9ComputeCommand", Visible 1},
   P_ = {"postop_helmholtz_vector", Name "GetDP/2PostOperationChoices", Visible 1}
 ];
+
+// Notes on possible mumps craches (NaN's when large nDOF)
+// -ksp_view
+// -ksp_error_if_not_converged
+// -ksp_monitor_true_residual
+// related to? https://lists.mcs.anl.gov/mailman/htdig/petsc-users/2018-October/036377.html
+
+// N.B. pastix works with petsc 3.12
+// -pc_factor_mat_solver_type pastix
\ No newline at end of file

DiffractionGratings/grating3D_data_hole.geo

+// note : invcreasing paramaille to 7 leads to:
+// Ndofs 1116912
+// solve wall time mumps  : 634.318s - R00=0.2432916553175037
+// solve wall time pastix : 966.046s - R00=0.2432916553175064
+// beware : very near a Rayleigh anomaly
+
 nm  = 1000;
 pp1 = "1Incident Plane Wave";
 pp2 = "2Layers Thicknesses";