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: 3.948E+02 m^2. Volume: 6.575E+01 m^3.
 G =    0.0000000000000000       ; curpol =    1.0000000000000000     
 Done initializing plasma surface. Took    1.89999994E-02  sec.
 Initializing coil surface.
   Building a plain circular torus.
   Evaluating coil surface & derivatives:   8.00000038E-03  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    1.09999999E-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   1.66299999     sec.
 Computing g and RHS_B:
  js=   1 j_RZetaZ= 1 dgemm took  0.000E+00 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  1.000E-03 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:   1.00000005E-03 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 =  6.16155118200228E+32, method 2 =  6.16155118200228E+32, relative difference = -1.16947164617258E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.99873293832602E-14, method 2 =  8.64019966684282E-12, relative difference = -9.77120093785995E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.999E-14,  chi2_M: 6.162E+32
   max(B_n): 4.835E-08,  max(M): 1.736E+18,  min(M): 1.533E+17
 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 =  4.52748265173057E+21, method 2 =  4.52748265173057E+21, relative difference =  2.89503041938547E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.54203348686308E-22, method 2 = -3.15480974677484E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.542E-22,  chi2_M: 4.527E+21
   max(B_n): 3.355E-12,  max(M): 2.906E+12,  min(M): 1.318E+11
 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 =  7.70224285821993E+19, method 2 =  7.70224285821992E+19, relative difference =  3.19075890651412E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.29031681523394E-24, method 2 = -1.13686837721616E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.290E-24,  chi2_M: 7.702E+19
   max(B_n): 2.490E-13,  max(M): 4.173E+11,  min(M): 1.236E+11
 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 =  3.20355424213311E+18, method 2 =  3.20355424213311E+18, relative difference = -2.39733727588961E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.86692500033171E-25, method 2 =  5.68434188608080E-14, relative difference = -9.99999999989913E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.867E-25,  chi2_M: 3.204E+18
   max(B_n): 9.703E-14,  max(M): 1.064E+11,  min(M): 4.017E+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 =  3.18853116332257E+17, method 2 =  3.18853116332257E+17, relative difference =  1.00359690280257E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.87554156728190E-26, method 2 = -1.13686837721616E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.876E-26,  chi2_M: 3.189E+17
   max(B_n): 4.352E-14,  max(M): 3.123E+10,  min(M): 1.776E+09
 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 =  3.32069767592616E+17, method 2 =  3.32069767592616E+17, relative difference = -2.89095874930033E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.94229365537407E-26, method 2 =  2.84217094304040E-14, relative difference = -9.99999999997226E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.942E-26,  chi2_M: 3.321E+17
   max(B_n): 3.364E-14,  max(M): 2.701E+10,  min(M): 4.264E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.02946459793393E+20, method 2 =  1.02946459793393E+20, relative difference =  4.77452066818470E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.19307544367104E-22, method 2 =  1.53676182890194E-10, relative difference = -9.99999999998447E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.193E-22,  chi2_M: 1.029E+20
   max(B_n): 1.193E-12,  max(M): 6.013E+11,  min(M): 4.597E+10
 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 =  4.80723753929309E+16, method 2 =  4.80723753929310E+16, relative difference = -2.49623612353564E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.53178918685960E-26, method 2 = -1.13686837721616E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.532E-26,  chi2_M: 4.807E+16
   max(B_n): 3.086E-14,  max(M): 9.713E+09,  min(M): 1.546E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.63345849690714E+16, method 2 =  3.63345849690715E+16, relative difference = -6.60527704401445E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.51500779048234E-26, method 2 =  5.68434188608080E-14, relative difference = -9.99999999998763E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.515E-26,  chi2_M: 3.633E+16
   max(B_n): 3.197E-14,  max(M): 7.100E+09,  min(M): 3.608E+09
 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 =  3.46247400113778E+16, method 2 =  3.46247400113779E+16, relative difference = -1.27076766455261E-15
   2 methods of computing chi2_B that should agree: method 1 =  4.09825873795119E-26, method 2 =  5.68434188608080E-14, relative difference = -9.99999999998558E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 4.098E-26,  chi2_M: 3.462E+16
   max(B_n): 3.664E-14,  max(M): 6.122E+09,  min(M): 4.643E+09
 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 =  4.45100963557100E+16, method 2 =  4.45100963557100E+16, relative difference = -2.69601752916910E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.65792578526259E-25, method 2 = -1.42108547152020E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.658E-25,  chi2_M: 4.451E+16
   max(B_n): 5.529E-14,  max(M): 8.524E+09,  min(M): 2.310E+09
 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 =  3.42411449000184E+16, method 2 =  3.42411449000184E+16, relative difference =  5.84092618935451E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.94081028611423E-24, method 2 =  0.00000000000000E+00, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.941E-24,  chi2_M: 3.424E+16
   max(B_n): 1.489E-13,  max(M): 5.448E+09,  min(M): 5.321E+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 =  3.42409395307435E+16, method 2 =  3.42409395307435E+16, relative difference = -8.17734571063974E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.72216451772815E-23, method 2 = -2.84217094304040E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.722E-23,  chi2_M: 3.424E+16
   max(B_n): 6.054E-13,  max(M): 5.433E+09,  min(M): 5.318E+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 =  3.42409298770436E+16, method 2 =  3.42409298770437E+16, relative difference = -8.76144430298111E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.56706831352811E-21, method 2 =  0.00000000000000E+00, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.567E-21,  chi2_M: 3.424E+16
   max(B_n): 2.861E-12,  max(M): 5.431E+09,  min(M): 5.317E+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 =  3.42409293438609E+16, method 2 =  3.42409293438608E+16, relative difference =  1.22660222151740E-15
   2 methods of computing chi2_B that should agree: method 1 =  3.65885838118403E-20, method 2 =  0.00000000000000E+00, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.659E-20,  chi2_M: 3.424E+16
   max(B_n): 1.376E-11,  max(M): 5.430E+09,  min(M): 5.317E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.42409293125279E+16, method 2 =  3.42409293125279E+16, relative difference =  1.16819259299031E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.54549269972038E-19, method 2 =  5.68434188608080E-14, relative difference = -9.99969933672516E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 8.545E-19,  chi2_M: 3.424E+16
   max(B_n): 6.646E-11,  max(M): 5.430E+09,  min(M): 5.317E+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 =  3.42409292945381E+16, method 2 =  3.42409292945381E+16, relative difference =  3.50457778081221E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.99596275430712E-17, method 2 = -2.84217094304040E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.996E-17,  chi2_M: 3.424E+16
   max(B_n): 3.211E-10,  max(M): 5.430E+09,  min(M): 5.317E+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 =  3.42409292110220E+16, method 2 =  3.42409292110220E+16, relative difference = -4.08867408758681E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.66200053112103E-16, method 2 =  5.68434188608080E-14, relative difference = -9.83730477991405E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.662E-16,  chi2_M: 3.424E+16
   max(B_n): 1.552E-09,  max(M): 5.430E+09,  min(M): 5.317E+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 =  3.42409288075943E+16, method 2 =  3.42409288075944E+16, relative difference = -2.33638522043411E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.08891291999054E-14, method 2 = -2.84217094304040E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.089E-14,  chi2_M: 3.424E+16
   max(B_n): 7.501E-09,  max(M): 5.430E+09,  min(M): 5.317E+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 =  3.42409268578617E+16, method 2 =  3.42409268578617E+16, relative difference = -5.25686704531107E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.54339691510085E-13, method 2 =  2.84217094304040E-13, relative difference = -5.54767920132882E-02
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.543E-13,  chi2_M: 3.424E+16
   max(B_n): 3.625E-08,  max(M): 5.430E+09,  min(M): 5.317E+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 =  3.42409174349428E+16, method 2 =  3.42409174349429E+16, relative difference = -5.25686849197300E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.94066489318860E-12, method 2 =  6.02540239924565E-12, relative difference = -7.08148333000111E-03
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.941E-12,  chi2_M: 3.424E+16
   max(B_n): 1.752E-07,  max(M): 5.430E+09,  min(M): 5.317E+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 =  3.42408718946890E+16, method 2 =  3.42408718946890E+16, relative difference = -2.33638910381860E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.38757197412731E-10, method 2 =  1.38726363729802E-10, relative difference =  1.11118953505495E-04
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.388E-10,  chi2_M: 3.424E+16
   max(B_n): 8.467E-07,  max(M): 5.430E+09,  min(M): 5.317E+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 =  3.42406518031142E+16, method 2 =  3.42406518031141E+16, relative difference =  7.00921236488180E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.24096058966653E-09, method 2 =  3.24092752634897E-09, relative difference =  5.10087754813320E-06
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.241E-09,  chi2_M: 3.424E+16
   max(B_n): 4.092E-06,  max(M): 5.430E+09,  min(M): 5.317E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.42395881459023E+16, method 2 =  3.42395881459023E+16, relative difference =  1.16823835116098E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.56974532553293E-08, method 2 =  7.56974145588174E-08, relative difference =  2.55599892525711E-07
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.570E-08,  chi2_M: 3.424E+16
   max(B_n): 1.978E-05,  max(M): 5.430E+09,  min(M): 5.317E+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 =  3.42344482640886E+16, method 2 =  3.42344482640886E+16, relative difference =  5.84206873898408E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.76781596693643E-06, method 2 =  1.76781588834274E-06, relative difference =  2.22290372576063E-08
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.768E-06,  chi2_M: 3.423E+16
   max(B_n): 9.560E-05,  max(M): 5.429E+09,  min(M): 5.316E+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 =  3.42096238949592E+16, method 2 =  3.42096238949591E+16, relative difference =  5.84630806272823E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.12612871415931E-05, method 2 =  4.12612871798501E-05, relative difference = -4.63594876631251E-10
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.126E-05,  chi2_M: 3.421E+16
   max(B_n): 4.618E-04,  max(M): 5.428E+09,  min(M): 5.314E+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 =  3.40900295479273E+16, method 2 =  3.40900295479273E+16, relative difference =  4.10677262110241E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.60376405884141E-04, method 2 =  9.60376405942043E-04, relative difference = -3.01454331266888E-11
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.604E-04,  chi2_M: 3.409E+16
   max(B_n): 2.228E-03,  max(M): 5.418E+09,  min(M): 5.305E+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 =  3.35207919592637E+16, method 2 =  3.35207919592637E+16, relative difference = -5.96644614611287E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.20568300733589E-02, method 2 =  2.20568300733817E-02, relative difference = -5.16244848350136E-13
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.206E-02,  chi2_M: 3.352E+16
   max(B_n): 1.067E-02,  max(M): 5.374E+09,  min(M): 5.259E+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 =  3.09609902582580E+16, method 2 =  3.09609902582581E+16, relative difference = -8.39766421652653E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.75820394975781E-01, method 2 =  4.75820394975756E-01, relative difference =  2.62494192323746E-14
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.758E-01,  chi2_M: 3.096E+16
   max(B_n): 4.947E-02,  max(M): 5.169E+09,  min(M): 5.050E+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 =  2.19305163146447E+16, method 2 =  2.19305163146446E+16, relative difference =  1.82394246565408E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.87161340771317E+00, method 2 =  7.87161340771320E+00, relative difference = -1.86174589197839E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.872E+00,  chi2_M: 2.193E+16
   max(B_n): 2.003E-01,  max(M): 4.358E+09,  min(M): 4.244E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  7.03657642840334E+15, method 2 =  7.03657642840334E+15, relative difference =  3.55286413135325E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.89902632028115E+01, method 2 =  5.89902632028114E+01, relative difference =  3.01127125555119E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.899E+01,  chi2_M: 7.037E+15
   max(B_n): 5.468E-01,  max(M): 2.470E+09,  min(M): 2.402E+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 =  7.34425121052715E+14, method 2 =  7.34425121052715E+14, relative difference = -8.51005748692438E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.43807203604324E+02, method 2 =  1.43807203604324E+02, relative difference =  8.89369163930900E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.438E+02,  chi2_M: 7.344E+14
   max(B_n): 8.541E-01,  max(M): 7.972E+08,  min(M): 7.768E+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.02510787343694E+13, method 2 =  4.02510787343694E+13, relative difference =  6.79329619473076E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.84088261637265E+02, method 2 =  1.84088261637265E+02, relative difference =  2.31587629577441E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.841E+02,  chi2_M: 4.025E+13
   max(B_n): 9.660E-01,  max(M): 1.864E+08,  min(M): 1.821E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.82097889724299E+12, method 2 =  1.82097889724299E+12, relative difference = -3.35177724148307E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.94523525938807E+02, method 2 =  1.94523525938807E+02, relative difference =  5.11382787898400E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.945E+02,  chi2_M: 1.821E+12
   max(B_n): 9.928E-01,  max(M): 3.962E+07,  min(M): 3.875E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.793E-12 ( 35 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  7.88719864945620E+10, method 2 =  7.88719864945621E+10, relative difference = -2.90194080446142E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.96793359239910E+02, method 2 =  1.96793359239910E+02, relative difference = -7.22120643201053E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.968E+02,  chi2_M: 7.887E+10
   max(B_n): 9.985E-01,  max(M): 8.245E+06,  min(M): 8.066E+06
 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 =  3.38491693031732E+09, method 2 =  3.38491693031733E+09, relative difference = -6.33920198364485E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97267978756064E+02, method 2 =  1.97267978756064E+02, relative difference =  8.64459896926794E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.973E+02,  chi2_M: 3.385E+09
   max(B_n): 9.997E-01,  max(M): 1.708E+06,  min(M): 1.671E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.44991856851784E+08, method 2 =  1.44991856851784E+08, relative difference =  2.05544801168802E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97366398491613E+02, method 2 =  1.97366398491613E+02, relative difference =  5.04016812216600E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 1.450E+08
   max(B_n): 9.999E-01,  max(M): 3.535E+05,  min(M): 3.459E+05
 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 =  6.20822421725131E+06, method 2 =  6.20822421725131E+06, relative difference = -7.50071632422308E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.97386772091695E+02, method 2 =  1.97386772091695E+02, relative difference = -5.03964789292987E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 6.208E+06
   max(B_n): 1.000E+00,  max(M): 7.315E+04,  min(M): 7.157E+04
 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 =  2.65800392622253E+05, method 2 =  2.65800392622253E+05, relative difference =  2.18990123901699E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97390988064826E+02, method 2 =  1.97390988064826E+02, relative difference =  7.19934321952682E-17
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.974E+02,  chi2_M: 2.658E+05
   max(B_n): 1.000E+00,  max(M): 1.513E+04,  min(M): 1.481E+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 =  1.13798486241897E+04, method 2 =  1.13798486241897E+04, relative difference = -1.59843022839452E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97391860424749E+02, method 2 =  1.97391860424749E+02, relative difference =  1.43986228050367E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 1.138E+04
   max(B_n): 1.000E+00,  max(M): 3.132E+03,  min(M): 3.064E+03
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   1.7430000305175781      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
