This is REGCOIL_PM,
a regularized least-squares method for designing permanent magnets for stellarators.
 Successfully read parameters from regcoil_nml namelist in regcoil_in.benchmark.
 Resolution parameters:
   ntheta_plasma  =   96
   ntheta_coil    =   96
   nzeta_plasma   =    4
   nzeta_coil     =    4
   mpol_magnetization =   24
   ntor_magnetization =    0
   ns_magnetization =    1
   ns_integration   =    2
 Not assuming stellarator symmetry
 We will use the following values of the regularization weight lambda:
 0.000E+00 1.000E-34 4.833E-34 2.336E-33 1.129E-32 5.456E-32 2.637E-31 1.274E-30 6.158E-30 2.976E-29 1.438E-28 6.952E-28 3.360E-27 1.624E-26 7.848E-26 3.793E-25 1.833E-24 8.859E-24 4.281E-23 2.069E-22 1.000E-21 4.833E-21 2.336E-20 1.129E-19 5.456E-19 2.637E-18 1.274E-17 6.158E-17 2.976E-16 1.438E-15 6.952E-15 3.360E-14 1.624E-13 7.848E-13 3.793E-12 1.833E-11 8.859E-11 4.281E-10 2.069E-09 1.000E-08
 Initializing plasma surface.
 Reading FOCUS format data from file:plasma
 Number of modes for Bnormal read from FOCUS file:           1
 Plasma surface area: 5.330E+02 m^2. Volume: 1.198E+02 m^3.
 G =    0.0000000000000000       ; curpol =    1.0000000000000000     
 Done initializing plasma surface. Took    4.10000011E-02  sec.
 Initializing coil surface.
   Building a plain circular torus.
   Evaluating coil surface & derivatives:   1.60000008E-02  sec.
 s_integration:  0.21132486540518713       0.78867513459481287     
 s_weights:  0.50000000000000011       0.50000000000000011     
 s_magnetization  0.50000000000000000     
 s_magnetization_weights   1.0000000000000000     
 interpolate_magnetization_to_integration:
   1.0000000000000000     
   1.0000000000000000     
 Maximum allowed uniform d before singularity:  2.9000E+01
 Coil surface area: 1.184E+03 m^2. Volume: 5.918E+02 m^3.
 Done initializing coil surface. Took    2.19999999E-02  sec.
 Using Bn coefficients in FOCUS file.
 Initializing basis functions
 Done. Took   1.00000005E-03 sec.
 Building inductance matrix.
   Number of OpenMP threads:           4
 Done. Took   4.12099981     sec.
 Computing g and RHS_B:
  js=   1 j_RZetaZ= 1 dgemm took  2.000E-03 sec, rest took  0.000E+00 sec.
  js=   1 j_RZetaZ= 2 dgemm took  1.000E-03 sec, rest took  0.000E+00 sec.
  js=   1 j_RZetaZ= 3 dgemm took  0.000E+00 sec, rest took  0.000E+00 sec.
 Total time to compute g and RHS_B:   3.00000003E-03 sec.
 About to form matrix_B.
  js=   1, ks=   1, j_RZetaZ= 1, k_RZetaZ= 1 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 2, k_RZetaZ= 1 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 2, k_RZetaZ= 2 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 1 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 2 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 3 dgemm took  1.000E-03 sec.
 Total time for matrix_B:   1.00000005E-03 sec.
 Regularization term:   0.00000000     sec.
 Optimal LWORK:        9408
 Solving system for lambda= 0.000E+00 (  1 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.75047909407466E+29, method 2 =  5.75047909407465E+29, relative difference =  1.16251717248505E-15
   2 methods of computing chi2_B that should agree: method 1 =  7.64753342814351E-16, method 2 =  5.11590769747272E-13, relative difference = -9.97014755152040E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.648E-16,  chi2_M: 5.750E+29
   max(B_n): 2.873E-09,  max(M): 6.828E+16,  min(M): 3.780E+14
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-34 (  2 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  8.38418855757079E+21, method 2 =  8.38418855757080E+21, relative difference = -7.50395337223053E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.35509305719453E-21, method 2 =  4.32009983342141E-12, relative difference = -9.99999995669039E-01
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 9.355E-21,  chi2_M: 8.384E+21
   max(B_n): 1.068E-11,  max(M): 3.454E+12,  min(M): 1.656E+12
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.833E-34 (  3 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.25063093776508E+18, method 2 =  5.25063093776509E+18, relative difference = -3.90048362620536E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.10206534095946E-23, method 2 =  3.41060513164848E-13, relative difference = -9.99999999876734E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.102E-23,  chi2_M: 5.251E+18
   max(B_n): 3.231E-13,  max(M): 1.299E+11,  min(M): 6.927E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.336E-33 (  4 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.50464324206685E+18, method 2 =  1.50464324206685E+18, relative difference =  4.25349997997443E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.43730840404958E-23, method 2 =  3.41060513164848E-13, relative difference = -9.99999999739794E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.437E-23,  chi2_M: 1.505E+18
   max(B_n): 3.652E-13,  max(M): 6.190E+10,  min(M): 7.284E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.129E-32 (  5 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.01691304302317E+17, method 2 =  2.01691304302317E+17, relative difference =  2.37987453975964E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.07807159490961E-23, method 2 =  3.41060513164848E-13, relative difference = -9.99999999936781E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.078E-23,  chi2_M: 2.017E+17
   max(B_n): 1.700E-13,  max(M): 2.481E+10,  min(M): 9.830E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.456E-32 (  6 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.47434454376916E+17, method 2 =  4.47434454376917E+17, relative difference = -2.14556565907931E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.71120122540136E-23, method 2 =  3.41060513164848E-13, relative difference = -9.99999999782373E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.711E-23,  chi2_M: 4.474E+17
   max(B_n): 2.932E-13,  max(M): 3.073E+10,  min(M): 1.360E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.637E-31 (  7 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.54362934074320E+16, method 2 =  2.54362934074320E+16, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  2.61799154001442E-22, method 2 =  3.41060513164848E-13, relative difference = -9.99999998464794E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.618E-22,  chi2_M: 2.544E+16
   max(B_n): 7.142E-13,  max(M): 9.271E+09,  min(M): 2.123E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.274E-30 (  8 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.12303257738585E+17, method 2 =  9.12303257738585E+17, relative difference =  2.10456334964680E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.59515200229312E-19, method 2 =  3.63797880709171E-12, relative difference = -9.99999912305595E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.595E-19,  chi2_M: 9.123E+17
   max(B_n): 1.734E-11,  max(M): 5.068E+10,  min(M): 3.991E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.158E-30 (  9 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.35309426041171E+16, method 2 =  1.35309426041171E+16, relative difference =  5.17332768662404E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.14408326840180E-23, method 2 =  1.13686837721616E-13, relative difference = -9.99999999446887E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.144E-23,  chi2_M: 1.353E+16
   max(B_n): 2.558E-13,  max(M): 5.558E+09,  min(M): 4.176E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.976E-29 ( 10 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.47264115180108E+15, method 2 =  4.47264115180109E+15, relative difference = -7.26640007479979E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.50483844422971E-23, method 2 =  3.41060513164848E-13, relative difference = -9.99999999618552E-01
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 6.505E-23,  chi2_M: 4.473E+15
   max(B_n): 3.620E-13,  max(M): 2.967E+09,  min(M): 5.707E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.438E-28 ( 11 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.45523229028406E+15, method 2 =  2.45523229028406E+15, relative difference = -8.14586875512545E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.68859344507967E-22, method 2 =  3.41060513164848E-13, relative difference = -9.99999996077767E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.689E-22,  chi2_M: 2.455E+15
   max(B_n): 1.133E-12,  max(M): 2.158E+09,  min(M): 4.465E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.952E-28 ( 12 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.05411275126382E+15, method 2 =  2.05411275126382E+15, relative difference = -9.12802862864455E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.56332274408239E-20, method 2 =  3.97903932025656E-13, relative difference = -9.99999921422104E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.563E-20,  chi2_M: 2.054E+15
   max(B_n): 5.434E-12,  max(M): 1.458E+09,  min(M): 1.182E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.360E-27 ( 13 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.03982339826640E+15, method 2 =  2.03982339826641E+15, relative difference = -1.22559629530904E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.41293289174773E-19, method 2 =  5.11590769747272E-13, relative difference = -9.99998665757548E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.413E-19,  chi2_M: 2.040E+15
   max(B_n): 2.534E-11,  max(M): 1.347E+09,  min(M): 1.288E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.624E-26 ( 14 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.03944202078833E+15, method 2 =  2.03944202078833E+15, relative difference =  7.35495289746058E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.24112707984930E-18, method 2 =  5.11590769747272E-13, relative difference = -9.99975601393491E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.241E-18,  chi2_M: 2.039E+15
   max(B_n): 1.083E-10,  max(M): 1.329E+09,  min(M): 1.301E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.848E-26 ( 15 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.03942258342757E+15, method 2 =  2.03942258342757E+15, relative difference = -3.67751149808054E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.06616254529333E-17, method 2 =  5.11590769747272E-13, relative difference = -9.99762879094112E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.066E-17,  chi2_M: 2.039E+15
   max(B_n): 3.380E-10,  max(M): 1.326E+09,  min(M): 1.301E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.793E-25 ( 16 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.03942071866310E+15, method 2 =  2.03942071866310E+15, relative difference =  4.90335314753265E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.92300259091207E-16, method 2 =  3.41060513164848E-13, relative difference = -9.98872974806205E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.923E-16,  chi2_M: 2.039E+15
   max(B_n): 6.036E-10,  max(M): 1.325E+09,  min(M): 1.300E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.833E-24 ( 17 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.03942055811918E+15, method 2 =  2.03942055811918E+15, relative difference = -6.12919191690796E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.73559692822187E-16, method 2 =  3.41060513164848E-13, relative difference = -9.98397114647642E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.736E-16,  chi2_M: 2.039E+15
   max(B_n): 7.304E-10,  max(M): 1.324E+09,  min(M): 1.300E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.859E-24 ( 18 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.03942054900699E+15, method 2 =  2.03942054900698E+15, relative difference =  1.59358990551628E-15
   2 methods of computing chi2_B that should agree: method 1 =  2.97908820160234E-16, method 2 =  5.11590769747272E-13, relative difference = -9.98836040597757E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.979E-16,  chi2_M: 2.039E+15
   max(B_n): 8.141E-10,  max(M): 1.324E+09,  min(M): 1.300E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.281E-23 ( 19 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.03942054752073E+15, method 2 =  2.03942054752073E+15, relative difference =  1.83875758462804E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.30457690034909E-16, method 2 =  5.11590769747272E-13, relative difference = -9.98708951016979E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.305E-16,  chi2_M: 2.039E+15
   max(B_n): 1.081E-09,  max(M): 1.324E+09,  min(M): 1.300E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.069E-22 ( 20 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.03942054238050E+15, method 2 =  2.03942054238050E+15, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  9.71207261311291E-16, method 2 =  3.97903932025656E-13, relative difference = -9.95130269271495E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 9.712E-16,  chi2_M: 2.039E+15
   max(B_n): 2.342E-09,  max(M): 1.324E+09,  min(M): 1.300E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-21 ( 21 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.03942051762108E+15, method 2 =  2.03942051762108E+15, relative difference = -3.06459601930966E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.59122001339279E-14, method 2 =  5.11590769747272E-13, relative difference = -9.39669723044359E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.591E-14,  chi2_M: 2.039E+15
   max(B_n): 8.431E-09,  max(M): 1.324E+09,  min(M): 1.300E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.833E-21 ( 22 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.03942039796433E+15, method 2 =  2.03942039796433E+15, relative difference =  1.22583847964618E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.64886888303226E-13, method 2 =  7.95807864051312E-13, relative difference = -3.71261242349840E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.649E-13,  chi2_M: 2.039E+15
   max(B_n): 3.785E-08,  max(M): 1.324E+09,  min(M): 1.300E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.336E-20 ( 23 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.03941981967193E+15, method 2 =  2.03941981967194E+15, relative difference = -6.74211354982902E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.51595986959403E-12, method 2 =  9.03810359886847E-12, relative difference = -2.97448924126613E-02
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.516E-12,  chi2_M: 2.039E+15
   max(B_n): 1.801E-07,  max(M): 1.324E+09,  min(M): 1.300E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.129E-19 ( 24 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.03941702482858E+15, method 2 =  2.03941702482858E+15, relative difference =  2.45168101429390E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.98902049170743E-10, method 2 =  1.99293026525993E-10, relative difference = -9.81873908325359E-04
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 1.989E-10,  chi2_M: 2.039E+15
   max(B_n): 8.684E-07,  max(M): 1.324E+09,  min(M): 1.300E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.456E-19 ( 25 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.03940351762664E+15, method 2 =  2.03940351762664E+15, relative difference = -1.10326376342551E-15
   2 methods of computing chi2_B that should agree: method 1 =  4.64576030193283E-09, method 2 =  4.64626737084473E-09, relative difference = -5.45703187456479E-05
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.646E-09,  chi2_M: 2.039E+15
   max(B_n): 4.200E-06,  max(M): 1.324E+09,  min(M): 1.300E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.637E-18 ( 26 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.03933824015359E+15, method 2 =  2.03933824015359E+15, relative difference =  8.58121505076176E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.08508540123519E-07, method 2 =  1.08508857010747E-07, relative difference = -1.46019274196658E-06
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.085E-07,  chi2_M: 2.039E+15
   max(B_n): 2.030E-05,  max(M): 1.324E+09,  min(M): 1.300E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.274E-17 ( 27 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.03902280285860E+15, method 2 =  2.03902280285860E+15, relative difference =  1.22607750952816E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.53406517646055E-06, method 2 =  2.53406574302062E-06, relative difference = -1.11788759003360E-07
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.534E-06,  chi2_M: 2.039E+15
   max(B_n): 9.810E-05,  max(M): 1.324E+09,  min(M): 1.300E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.158E-17 ( 28 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.03749934765884E+15, method 2 =  2.03749934765884E+15, relative difference =  3.06748564468571E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.91444743851366E-05, method 2 =  5.91444747897185E-05, relative difference = -3.42028517967266E-09
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.914E-05,  chi2_M: 2.037E+15
   max(B_n): 4.739E-04,  max(M): 1.324E+09,  min(M): 1.300E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.976E-16 ( 29 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.03016060953432E+15, method 2 =  2.03016060953432E+15, relative difference =  5.54143349406259E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.37647363090603E-03, method 2 =  1.37647363123961E-03, relative difference = -1.21168864970787E-10
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.376E-03,  chi2_M: 2.030E+15
   max(B_n): 2.286E-03,  max(M): 1.321E+09,  min(M): 1.297E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.438E-15 ( 30 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.99524585721489E+15, method 2 =  1.99524585721490E+15, relative difference = -8.77084893408962E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.15975781860544E-02, method 2 =  3.15975781864495E-02, relative difference = -6.25174182461022E-12
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.160E-02,  chi2_M: 1.995E+15
   max(B_n): 1.095E-02,  max(M): 1.310E+09,  min(M): 1.286E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.952E-15 ( 31 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.83856591165580E+15, method 2 =  1.83856591165580E+15, relative difference =  2.03963315985924E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.80068816295623E-01, method 2 =  6.80068816296057E-01, relative difference = -3.19238446555982E-13
   Diagnostics:    0.00000000      sec.
   chi2_B: 6.801E-01,  chi2_M: 1.839E+15
   max(B_n): 5.078E-02,  max(M): 1.258E+09,  min(M): 1.234E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.360E-14 ( 32 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.29044372512231E+15, method 2 =  1.29044372512231E+15, relative difference =  1.93731811107303E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.11486337855462E+01, method 2 =  1.11486337855465E+01, relative difference = -1.41010623643171E-14
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.115E+01,  chi2_M: 1.290E+15
   max(B_n): 2.052E-01,  max(M): 1.055E+09,  min(M): 1.033E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.624E-13 ( 33 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.05470773082294E+14, method 2 =  4.05470773082293E+14, relative difference =  1.54141812799205E-15
   2 methods of computing chi2_B that should agree: method 1 =  8.18184797337736E+01, method 2 =  8.18184797337737E+01, relative difference = -6.07906571535523E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 8.182E+01,  chi2_M: 4.055E+14
   max(B_n): 5.545E-01,  max(M): 5.920E+08,  min(M): 5.784E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.848E-13 ( 34 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.15510255382837E+13, method 2 =  4.15510255382836E+13, relative difference =  9.40109166836549E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.95834719175744E+02, method 2 =  1.95834719175744E+02, relative difference =  1.45131106220741E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.958E+02,  chi2_M: 4.155E+13
   max(B_n): 8.573E-01,  max(M): 1.895E+08,  min(M): 1.852E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.793E-12 ( 35 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.26211690498781E+12, method 2 =  2.26211690498781E+12, relative difference =  6.47554397728128E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.49024711460192E+02, method 2 =  2.49024711460192E+02, relative difference = -4.56528338312278E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.490E+02,  chi2_M: 2.262E+12
   max(B_n): 9.667E-01,  max(M): 4.420E+07,  min(M): 4.321E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.833E-11 ( 36 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.02175205367870E+11, method 2 =  1.02175205367870E+11, relative difference =  7.46697254366289E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.62720263995861E+02, method 2 =  2.62720263995860E+02, relative difference =  2.16364805653909E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.627E+02,  chi2_M: 1.022E+11
   max(B_n): 9.929E-01,  max(M): 9.394E+06,  min(M): 9.184E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.859E-11 ( 37 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.42399110374823E+09, method 2 =  4.42399110374823E+09, relative difference =  4.31137538047149E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.65694927386866E+02, method 2 =  2.65694927386866E+02, relative difference =  0.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.657E+02,  chi2_M: 4.424E+09
   max(B_n): 9.985E-01,  max(M): 1.955E+06,  min(M): 1.911E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.281E-10 ( 38 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.89849038315187E+08, method 2 =  1.89849038315187E+08, relative difference =  4.70937161212545E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.66316732654131E+02, method 2 =  2.66316732654131E+02, relative difference =  2.13442911732592E-16
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 2.663E+02,  chi2_M: 1.898E+08
   max(B_n): 9.997E-01,  max(M): 4.049E+05,  min(M): 3.959E+05
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.069E-09 ( 39 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  8.13200352047100E+06, method 2 =  8.13200352047099E+06, relative difference =  6.87153594268147E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.66445665290484E+02, method 2 =  2.66445665290484E+02, relative difference = -2.13339627044922E-16
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 2.664E+02,  chi2_M: 8.132E+06
   max(B_n): 9.999E-01,  max(M): 8.380E+04,  min(M): 8.193E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-08 ( 40 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.48192976713133E+05, method 2 =  3.48192976713132E+05, relative difference =  1.67170692134389E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.66472354924380E+02, method 2 =  2.66472354924380E+02, relative difference =  2.13318259137760E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.665E+02,  chi2_M: 3.482E+05
   max(B_n): 1.000E+00,  max(M): 1.734E+04,  min(M): 1.695E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   4.2480001449584961      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
