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.922E+02 m^2. Volume: 1.479E+02 m^3.
 G =    0.0000000000000000       ; curpol =    1.0000000000000000     
 Done initializing plasma surface. Took    3.59999985E-02  sec.
 Initializing coil surface.
   Building a plain circular torus.
   Evaluating coil surface & derivatives:   1.70000009E-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.09999997E-02  sec.
 Using Bn coefficients in FOCUS file.
 Initializing basis functions
 Done. Took   0.00000000     sec.
 Building inductance matrix.
   Number of OpenMP threads:           4
 Done. Took   3.27999997     sec.
 Computing g and RHS_B:
  js=   1 j_RZetaZ= 1 dgemm took  1.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:   2.00000009E-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  0.000E+00 sec.
 Total time for matrix_B:   0.00000000     sec.
 Regularization term:   1.00000005E-03 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 =  3.74925598115303E+27, method 2 =  3.74925598115302E+27, relative difference =  3.66576606566439E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.90821959546307E-18, method 2 = -4.54747350886464E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.908E-18,  chi2_M: 3.749E+27
   max(B_n): 2.891E-10,  max(M): 4.548E+15,  min(M): 1.999E+12
 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 =  1.88861819643753E+20, method 2 =  1.88861819643753E+20, relative difference = -2.60253766974789E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.30572438963340E-22, method 2 = -2.84217094304040E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.306E-22,  chi2_M: 1.889E+20
   max(B_n): 1.521E-12,  max(M): 6.283E+11,  min(M): 6.486E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.833E-34 (  3 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.74659303462003E+19, method 2 =  2.74659303462003E+19, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  3.68974423359611E-23, method 2 = -5.68434188608080E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.690E-23,  chi2_M: 2.747E+19
   max(B_n): 6.384E-13,  max(M): 2.775E+11,  min(M): 1.852E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.336E-33 (  4 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.70634983606830E+21, method 2 =  1.70634983606831E+21, relative difference = -1.53628520048398E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.14275954177007E-20, method 2 =  4.17799128626939E-11, relative difference = -9.99999998974263E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.143E-20,  chi2_M: 1.706E+21
   max(B_n): 1.211E-11,  max(M): 2.107E+12,  min(M): 8.590E+10
 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 =  1.52853607862507E+20, method 2 =  1.52853607862507E+20, relative difference = -2.14375051123917E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.39291958656108E-21, method 2 =  2.16004991671070E-11, relative difference = -9.99999999871029E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.393E-21,  chi2_M: 1.529E+20
   max(B_n): 3.432E-12,  max(M): 7.258E+11,  min(M): 3.364E+10
 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 =  7.64713879740219E+17, method 2 =  7.64713879740220E+17, relative difference = -2.51074297311335E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.78023540106734E-23, method 2 = -1.08002495835535E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.780E-23,  chi2_M: 7.647E+17
   max(B_n): 3.042E-13,  max(M): 4.598E+10,  min(M): 6.387E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.637E-31 (  7 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.77619555576432E+18, method 2 =  3.77619555576432E+18, relative difference =  6.77930992236932E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.56914688689670E-21, method 2 = -8.07176547823474E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.569E-21,  chi2_M: 3.776E+18
   max(B_n): 4.041E-12,  max(M): 1.218E+11,  min(M): 8.032E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.274E-30 (  8 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.84780315561190E+16, method 2 =  1.84780315561190E+16, relative difference = -5.41183186619709E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.97398349173530E-26, method 2 = -6.82121026329696E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.974E-26,  chi2_M: 1.848E+16
   max(B_n): 3.675E-14,  max(M): 7.115E+09,  min(M): 1.101E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.158E-30 (  9 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  8.19570157281763E+17, method 2 =  8.19570157281764E+17, relative difference = -3.90448574971704E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.07049224832039E-21, method 2 = -1.10844666778576E-11, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.070E-21,  chi2_M: 8.196E+17
   max(B_n): 1.420E-12,  max(M): 3.037E+10,  min(M): 2.129E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.976E-29 ( 10 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.06041221318162E+15, method 2 =  9.06041221318162E+15, relative difference =  2.20740508593007E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.51742862002722E-24, method 2 = -9.66338120633736E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.517E-24,  chi2_M: 9.060E+15
   max(B_n): 8.160E-14,  max(M): 4.648E+09,  min(M): 5.941E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.438E-28 ( 11 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.55458312073861E+15, method 2 =  1.55458312073862E+15, relative difference = -5.62851859335942E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.11845947538643E-23, method 2 = -8.52651282912120E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.118E-23,  chi2_M: 1.555E+15
   max(B_n): 2.018E-13,  max(M): 1.582E+09,  min(M): 5.298E+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 =  1.35089272706078E+15, method 2 =  1.35089272706078E+15, relative difference =  9.25314034904682E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.74619971946425E-22, method 2 = -7.95807864051312E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.746E-22,  chi2_M: 1.351E+15
   max(B_n): 6.951E-13,  max(M): 1.179E+09,  min(M): 9.403E+08
 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 =  1.33857307742584E+15, method 2 =  1.33857307742584E+15, relative difference = -1.86766045288140E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.50374860236731E-21, method 2 = -7.38964445190504E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.504E-21,  chi2_M: 1.339E+15
   max(B_n): 3.074E-12,  max(M): 1.093E+09,  min(M): 1.035E+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 =  1.33806533431606E+15, method 2 =  1.33806533431606E+15, relative difference =  9.34184578243276E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.23307274431535E-19, method 2 = -7.95807864051312E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.233E-19,  chi2_M: 1.338E+15
   max(B_n): 1.451E-11,  max(M): 1.075E+09,  min(M): 1.054E+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 =  1.33805493294207E+15, method 2 =  1.33805493294207E+15, relative difference = -6.53934288090908E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.74268803317283E-18, method 2 = -6.82121026329696E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.743E-18,  chi2_M: 1.338E+15
   max(B_n): 6.842E-11,  max(M): 1.073E+09,  min(M): 1.054E+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 =  1.33805392685480E+15, method 2 =  1.33805392685480E+15, relative difference = -1.86838508510374E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.28219932288367E-17, method 2 = -5.68434188608080E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.282E-17,  chi2_M: 1.338E+15
   max(B_n): 3.004E-10,  max(M): 1.072E+09,  min(M): 1.054E+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 =  1.33805333877469E+15, method 2 =  1.33805333877469E+15, relative difference =  9.34192953133448E-17
   2 methods of computing chi2_B that should agree: method 1 =  5.93772327172848E-16, method 2 = -6.82121026329696E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.938E-16,  chi2_M: 1.338E+15
   max(B_n): 1.010E-09,  max(M): 1.072E+09,  min(M): 1.054E+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 =  1.33805290619622E+15, method 2 =  1.33805290619622E+15, relative difference =  4.67096627574116E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.25575030571275E-15, method 2 = -8.52651282912120E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.256E-15,  chi2_M: 1.338E+15
   max(B_n): 1.992E-09,  max(M): 1.072E+09,  min(M): 1.054E+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 =  1.33805282748668E+15, method 2 =  1.33805282748668E+15, relative difference = -1.86838662020233E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.50305861004553E-15, method 2 = -6.82121026329696E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.503E-15,  chi2_M: 1.338E+15
   max(B_n): 2.622E-09,  max(M): 1.072E+09,  min(M): 1.054E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.33805282021025E+15, method 2 =  1.33805282021025E+15, relative difference =  3.73677326072550E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.13264854346229E-15, method 2 = -7.38964445190504E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.133E-15,  chi2_M: 1.338E+15
   max(B_n): 3.495E-09,  max(M): 1.072E+09,  min(M): 1.054E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.33805281036630E+15, method 2 =  1.33805281036630E+15, relative difference = -6.53935325437915E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.00075877425544E-14, method 2 = -7.95807864051312E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.001E-14,  chi2_M: 1.338E+15
   max(B_n): 7.117E-09,  max(M): 1.072E+09,  min(M): 1.054E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.33805276399946E+15, method 2 =  1.33805276399946E+15, relative difference =  4.67096677213134E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.45220920887854E-13, method 2 = -5.68434188608080E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.452E-13,  chi2_M: 1.338E+15
   max(B_n): 2.449E-08,  max(M): 1.072E+09,  min(M): 1.054E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.33805253996478E+15, method 2 =  1.33805253996478E+15, relative difference = -9.34193510841440E-17
   2 methods of computing chi2_B that should agree: method 1 =  3.30300311275342E-12, method 2 =  2.50111042987555E-12, relative difference =  1.38159372139823E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.303E-12,  chi2_M: 1.338E+15
   max(B_n): 1.084E-07,  max(M): 1.072E+09,  min(M): 1.054E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.129E-19 ( 24 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.33805145722382E+15, method 2 =  1.33805145722382E+15, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  7.70598498016539E-11, method 2 =  7.63975549489260E-11, relative difference =  4.31582206022797E-03
   Diagnostics:    0.00000000      sec.
   chi2_B: 7.706E-11,  chi2_M: 1.338E+15
   max(B_n): 5.140E-07,  max(M): 1.072E+09,  min(M): 1.054E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.33804622443096E+15, method 2 =  1.33804622443096E+15, relative difference = -9.34197920203839E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.79980735809707E-09, method 2 =  1.79920789378230E-09, relative difference =  1.66563427166466E-04
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.800E-09,  chi2_M: 1.338E+15
   max(B_n): 2.475E-06,  max(M): 1.072E+09,  min(M): 1.054E+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 =  1.33802093514079E+15, method 2 =  1.33802093514079E+15, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  4.20376032030784E-08, method 2 =  4.20367314291070E-08, relative difference =  1.03690855864143E-05
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.204E-08,  chi2_M: 1.338E+15
   max(B_n): 1.197E-05,  max(M): 1.072E+09,  min(M): 1.054E+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 =  1.33789872387152E+15, method 2 =  1.33789872387152E+15, relative difference = -1.86860182717393E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.81791543697399E-07, method 2 =  9.81790890364209E-07, relative difference =  3.32725114331491E-07
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.818E-07,  chi2_M: 1.338E+15
   max(B_n): 5.788E-05,  max(M): 1.072E+09,  min(M): 1.054E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.33730832129119E+15, method 2 =  1.33730832129119E+15, relative difference = -7.47770715308558E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.29217952349456E-05, method 2 =  2.29217945957316E-05, relative difference =  1.39433660868940E-08
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.292E-05,  chi2_M: 1.337E+15
   max(B_n): 2.797E-04,  max(M): 1.072E+09,  min(M): 1.054E+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 =  1.33446044835599E+15, method 2 =  1.33446044835599E+15, relative difference = -9.36708166615172E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.34249050124361E-04, method 2 =  5.34249049394475E-04, relative difference =  6.83094955234308E-10
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.342E-04,  chi2_M: 1.334E+15
   max(B_n): 1.350E-03,  max(M): 1.070E+09,  min(M): 1.053E+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.32082428247724E+15, method 2 =  1.32082428247724E+15, relative difference =  2.83913617409181E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.23510466105960E-02, method 2 =  1.23510466098082E-02, relative difference =  3.18921197560271E-11
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.235E-02,  chi2_M: 1.321E+15
   max(B_n): 6.491E-03,  max(M): 1.065E+09,  min(M): 1.047E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.952E-15 ( 31 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.25777849050046E+15, method 2 =  1.25777849050046E+15, relative difference =  3.97526276507601E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.74714891711350E-01, method 2 =  2.74714891710744E-01, relative difference =  1.10399976085381E-12
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.747E-01,  chi2_M: 1.258E+15
   max(B_n): 3.060E-02,  max(M): 1.039E+09,  min(M): 1.022E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.00853452588270E+15, method 2 =  1.00853452588270E+15, relative difference = -2.47884423967729E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.14498314517500E+00, method 2 =  5.14498314517436E+00, relative difference =  6.18015385630285E-14
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.145E+00,  chi2_M: 1.009E+15
   max(B_n): 1.323E-01,  max(M): 9.314E+08,  min(M): 9.143E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.624E-13 ( 33 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.45114579637836E+14, method 2 =  4.45114579637836E+14, relative difference =  2.10619926393512E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.30366838766528E+01, method 2 =  5.30366838766521E+01, relative difference =  6.36366708504986E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.304E+01,  chi2_M: 4.451E+14
   max(B_n): 4.239E-01,  max(M): 6.195E+08,  min(M): 6.067E+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 =  6.47349704978317E+13, method 2 =  6.47349704978317E+13, relative difference =  6.03421917081253E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.80160186662199E+02, method 2 =  1.80160186662199E+02, relative difference =  2.44524888840879E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.802E+02,  chi2_M: 6.473E+13
   max(B_n): 7.803E-01,  max(M): 2.364E+08,  min(M): 2.313E+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 =  4.06655650870169E+12, method 2 =  4.06655650870169E+12, relative difference = -2.40144824720949E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.64341868360977E+02, method 2 =  2.64341868360977E+02, relative difference =  7.52631307504967E-16
   Diagnostics:    3.00000003E-03  sec.
   chi2_B: 2.643E+02,  chi2_M: 4.067E+12
   max(B_n): 9.449E-01,  max(M): 5.924E+07,  min(M): 5.796E+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.90387293276133E+11, method 2 =  1.90387293276134E+11, relative difference = -4.80876285384308E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.89066476592063E+02, method 2 =  2.89066476592062E+02, relative difference =  9.83223989356380E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.891E+02,  chi2_M: 1.904E+11
   max(B_n): 9.881E-01,  max(M): 1.282E+07,  min(M): 1.254E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.859E-11 ( 37 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  8.30755574177205E+09, method 2 =  8.30755574177206E+09, relative difference = -5.16582079882947E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.94614416108903E+02, method 2 =  2.94614416108902E+02, relative difference =  1.64000481272535E-15
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.946E+02,  chi2_M: 8.308E+09
   max(B_n): 9.975E-01,  max(M): 2.678E+06,  min(M): 2.620E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.281E-10 ( 38 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.57084015069971E+08, method 2 =  3.57084015069971E+08, relative difference = -8.34602534136274E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.95782293941148E+02, method 2 =  2.95782293941147E+02, relative difference =  9.60899621532422E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.958E+02,  chi2_M: 3.571E+08
   max(B_n): 9.995E-01,  max(M): 5.551E+05,  min(M): 5.431E+05
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.069E-09 ( 39 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.53004886277069E+07, method 2 =  1.53004886277069E+07, relative difference =  3.04344062884695E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.96024810999483E+02, method 2 =  2.96024810999482E+02, relative difference =  1.15213489036040E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.960E+02,  chi2_M: 1.530E+07
   max(B_n): 9.999E-01,  max(M): 1.149E+05,  min(M): 1.124E+05
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-08 ( 40 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.55175902595066E+05, method 2 =  6.55175902595065E+05, relative difference =  7.10742390651604E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.96075028369645E+02, method 2 =  2.96075028369644E+02, relative difference =  1.15193947643244E-15
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.961E+02,  chi2_M: 6.552E+05
   max(B_n): 1.000E+00,  max(M): 2.378E+04,  min(M): 2.326E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   3.3889999389648438      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
