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: 2.369E+02 m^2. Volume: 2.367E+01 m^3.
 G =    0.0000000000000000       ; curpol =    1.0000000000000000     
 Done initializing plasma surface. Took    2.19999999E-02  sec.
 Initializing coil surface.
   Building a plain circular torus.
   Evaluating coil surface & derivatives:   1.20000001E-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    1.49999997E-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   1.65699995     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  1.000E-03 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  1.000E-03 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:   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 =  2.09824155962433E+31, method 2 =  2.09824155962434E+31, relative difference = -3.21955277745295E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.58538728843869E-15, method 2 = -1.42108547152020E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.585E-15,  chi2_M: 2.098E+31
   max(B_n): 1.097E-08,  max(M): 2.687E+17,  min(M): 3.658E+16
 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 =  2.09885447895416E+20, method 2 =  2.09885447895416E+20, relative difference = -4.68369774968791E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.50634647795702E-24, method 2 =  1.70530256582424E-13, relative difference = -9.99999999947149E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.506E-24,  chi2_M: 2.099E+20
   max(B_n): 3.431E-13,  max(M): 6.054E+11,  min(M): 2.039E+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 =  8.24261785083042E+19, method 2 =  8.24261785083041E+19, relative difference =  3.97543603173338E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.59719864181520E-24, method 2 =  2.27373675443232E-13, relative difference = -9.99999999985951E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.597E-24,  chi2_M: 8.243E+19
   max(B_n): 2.319E-13,  max(M): 3.604E+11,  min(M): 1.386E+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 =  5.93694522405396E+18, method 2 =  5.93694522405396E+18, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  2.68120372733284E-25, method 2 =  5.68434188608080E-14, relative difference = -9.99999999990566E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.681E-25,  chi2_M: 5.937E+18
   max(B_n): 1.077E-13,  max(M): 1.312E+11,  min(M): 3.451E+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 =  4.28983839303804E+17, method 2 =  4.28983839303803E+17, relative difference =  4.47569307765990E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.18588587840974E-26, method 2 =  5.68434188608080E-14, relative difference = -9.99999999998879E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.186E-26,  chi2_M: 4.290E+17
   max(B_n): 3.642E-14,  max(M): 3.223E+10,  min(M): 3.037E+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 =  1.12664969403594E+17, method 2 =  1.12664969403594E+17, relative difference =  2.13020960526124E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.91983613096193E-26, method 2 =  5.68434188608080E-14, relative difference = -9.99999999998973E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.920E-26,  chi2_M: 1.127E+17
   max(B_n): 3.331E-14,  max(M): 1.726E+10,  min(M): 9.469E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  7.84541736363501E+16, method 2 =  7.84541736363502E+16, relative difference = -7.13792490627338E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.88222785715482E-26, method 2 =  1.13686837721616E-13, relative difference = -9.99999999999493E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.882E-26,  chi2_M: 7.845E+16
   max(B_n): 3.275E-14,  max(M): 1.274E+10,  min(M): 8.153E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.274E-30 (  8 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.37707534123856E+16, method 2 =  5.37707534123855E+16, relative difference =  7.43898819739922E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.83461498344594E-26, method 2 =  9.94759830064140E-14, relative difference = -9.99999999999430E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.835E-26,  chi2_M: 5.377E+16
   max(B_n): 3.253E-14,  max(M): 8.428E+09,  min(M): 5.135E+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 =  5.23960526947306E+16, method 2 =  5.23960526947305E+16, relative difference =  6.10733029574540E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.92062561657436E-26, method 2 =  1.13686837721616E-13, relative difference = -9.99999999999486E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.921E-26,  chi2_M: 5.240E+16
   max(B_n): 3.464E-14,  max(M): 7.193E+09,  min(M): 6.171E+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 =  5.22606328689088E+16, method 2 =  5.22606328689088E+16, relative difference = -6.12315585237347E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.85919177211981E-26, method 2 =  9.94759830064140E-14, relative difference = -9.99999999999023E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.859E-26,  chi2_M: 5.226E+16
   max(B_n): 4.197E-14,  max(M): 6.810E+09,  min(M): 6.531E+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 =  5.22507635619861E+16, method 2 =  5.22507635619860E+16, relative difference =  9.95200767512449E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.68973292413186E-25, method 2 =  1.13686837721616E-13, relative difference = -9.99999999989991E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.690E-25,  chi2_M: 5.225E+16
   max(B_n): 1.030E-13,  max(M): 6.717E+09,  min(M): 6.571E+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 =  5.22501797381774E+16, method 2 =  5.22501797381774E+16, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  1.29009248276115E-23, method 2 =  1.13686837721616E-13, relative difference = -9.99999999773045E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.290E-23,  chi2_M: 5.225E+16
   max(B_n): 4.098E-13,  max(M): 6.714E+09,  min(M): 6.565E+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 =  5.22501414193538E+16, method 2 =  5.22501414193538E+16, relative difference =  3.82774083604525E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.01710928948163E-22, method 2 =  5.68434188608080E-14, relative difference = -9.99999989384490E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.017E-22,  chi2_M: 5.225E+16
   max(B_n): 1.927E-12,  max(M): 6.719E+09,  min(M): 6.564E+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 =  5.22501398539948E+16, method 2 =  5.22501398539947E+16, relative difference =  6.12438552115253E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.05102077776680E-21, method 2 =  5.68434188608080E-14, relative difference = -9.99999751914292E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.051E-21,  chi2_M: 5.225E+16
   max(B_n): 9.274E-12,  max(M): 6.720E+09,  min(M): 6.564E+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 =  5.22501397861343E+16, method 2 =  5.22501397861343E+16, relative difference = -6.12438552910664E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.64706903242156E-19, method 2 =  8.52651282912120E-14, relative difference = -9.99996136601484E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.647E-19,  chi2_M: 5.225E+16
   max(B_n): 4.478E-11,  max(M): 6.720E+09,  min(M): 6.564E+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 =  5.22501397665858E+16, method 2 =  5.22501397665858E+16, relative difference =  4.59328914854848E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.84660442692319E-18, method 2 =  9.94759830064140E-14, relative difference = -9.99922665641026E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.847E-18,  chi2_M: 5.225E+16
   max(B_n): 2.163E-10,  max(M): 6.720E+09,  min(M): 6.564E+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 =  5.22501396886478E+16, method 2 =  5.22501396886478E+16, relative difference = -1.53109638513333E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.97780062956896E-17, method 2 =  1.27897692436818E-13, relative difference = -9.98597081323900E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.978E-17,  chi2_M: 5.225E+16
   max(B_n): 1.045E-09,  max(M): 6.720E+09,  min(M): 6.564E+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 =  5.22501393143615E+16, method 2 =  5.22501393143615E+16, relative difference = -7.65548198050556E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.08940743503313E-15, method 2 =  1.42108547152020E-13, relative difference = -9.71020290252845E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.089E-15,  chi2_M: 5.225E+16
   max(B_n): 5.035E-09,  max(M): 6.720E+09,  min(M): 6.564E+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 =  5.22501375313787E+16, method 2 =  5.22501375313787E+16, relative difference =  6.12438579339288E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.80403028924822E-14, method 2 =  1.84741111297626E-13, relative difference = -5.87249668882511E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.804E-14,  chi2_M: 5.225E+16
   max(B_n): 2.402E-08,  max(M): 6.720E+09,  min(M): 6.564E+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 =  5.22501293063385E+16, method 2 =  5.22501293063386E+16, relative difference = -4.59329006810487E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.06821370054831E-12, method 2 =  1.15107923193136E-12, relative difference = -3.73387082751902E-02
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.068E-12,  chi2_M: 5.225E+16
   max(B_n): 1.107E-07,  max(M): 6.720E+09,  min(M): 6.564E+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 =  5.22500919667657E+16, method 2 =  5.22500919667657E+16, relative difference = -7.65548891769272E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.34949567736265E-11, method 2 =  2.35900188272353E-11, relative difference = -2.01894664690161E-03
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.349E-11,  chi2_M: 5.225E+16
   max(B_n): 4.885E-07,  max(M): 6.720E+09,  min(M): 6.564E+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 =  5.22499149674002E+16, method 2 =  5.22499149674002E+16, relative difference =  7.65551485106853E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.39485648845408E-10, method 2 =  5.39571942681505E-10, relative difference = -7.99714832406381E-05
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.395E-10,  chi2_M: 5.225E+16
   max(B_n): 2.213E-06,  max(M): 6.720E+09,  min(M): 6.564E+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 =  5.22490605392088E+16, method 2 =  5.22490605392089E+16, relative difference = -4.59338402496058E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.25825788666245E-08, method 2 =  1.25826886687719E-08, relative difference = -4.36324180717479E-06
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.258E-08,  chi2_M: 5.225E+16
   max(B_n): 1.047E-05,  max(M): 6.719E+09,  min(M): 6.564E+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 =  5.22449315274410E+16, method 2 =  5.22449315274410E+16, relative difference =  7.65624508072912E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.93849231016606E-07, method 2 =  2.93849325316842E-07, relative difference = -1.60456810042864E-07
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.938E-07,  chi2_M: 5.224E+16
   max(B_n): 5.073E-05,  max(M): 6.719E+09,  min(M): 6.564E+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 =  5.22249832494153E+16, method 2 =  5.22249832494153E+16, relative difference =  3.06366780886983E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.86083911006389E-06, method 2 =  6.86083917855740E-06, relative difference = -4.99162755718031E-09
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.861E-06,  chi2_M: 5.222E+16
   max(B_n): 2.454E-04,  max(M): 6.718E+09,  min(M): 6.563E+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 =  5.21287365221488E+16, method 2 =  5.21287365221489E+16, relative difference = -9.20797303030841E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.59952797340056E-04, method 2 =  1.59952797460505E-04, relative difference = -3.76515918227646E-10
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.600E-04,  chi2_M: 5.213E+16
   max(B_n): 1.185E-03,  max(M): 6.712E+09,  min(M): 6.556E+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 =  5.16673340673043E+16, method 2 =  5.16673340673043E+16, relative difference =  6.19346838339197E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.70278802793268E-03, method 2 =  3.70278802803625E-03, relative difference = -1.39862181090426E-11
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.703E-03,  chi2_M: 5.167E+16
   max(B_n): 5.694E-03,  max(M): 6.685E+09,  min(M): 6.525E+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 =  4.95217961885036E+16, method 2 =  4.95217961885037E+16, relative difference = -1.05004268831573E-15
   2 methods of computing chi2_B that should agree: method 1 =  8.28818069325647E-02, method 2 =  8.28818069326616E-02, relative difference = -5.84200603889248E-13
   Diagnostics:    0.00000000      sec.
   chi2_B: 8.288E-02,  chi2_M: 4.952E+16
   max(B_n): 2.683E-02,  max(M): 6.552E+09,  min(M): 6.382E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.08239348844515E+16, method 2 =  4.08239348844515E+16, relative difference =  9.79817357469741E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.59554901003805E+00, method 2 =  1.59554901003817E+00, relative difference = -3.69483119833390E-14
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.596E+00,  chi2_M: 4.082E+16
   max(B_n): 1.168E-01,  max(M): 5.960E+09,  min(M): 5.783E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.95550347905508E+16, method 2 =  1.95550347905508E+16, relative difference =  5.11377254354571E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.78502974180933E+01, method 2 =  1.78502974180934E+01, relative difference = -1.19416957451905E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.785E+01,  chi2_M: 1.956E+16
   max(B_n): 3.885E-01,  max(M): 4.130E+09,  min(M): 3.998E+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 =  3.15899267250103E+15, method 2 =  3.15899267250102E+15, relative difference =  3.16556606384365E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.73521749399644E+01, method 2 =  6.73521749399645E+01, relative difference = -4.21986512770199E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.735E+01,  chi2_M: 3.159E+15
   max(B_n): 7.546E-01,  max(M): 1.659E+09,  min(M): 1.607E+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 =  2.08713852535346E+14, method 2 =  2.08713852535346E+14, relative difference =  7.48632628366337E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.03937280944699E+02, method 2 =  1.03937280944699E+02, relative difference =  3.41813221060761E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.039E+02,  chi2_M: 2.087E+14
   max(B_n): 9.373E-01,  max(M): 4.259E+08,  min(M): 4.138E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.90382154435242E+12, method 2 =  9.90382154435243E+12, relative difference = -4.93023069744667E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.15196478939460E+02, method 2 =  1.15196478939460E+02, relative difference =  2.46723768747659E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.152E+02,  chi2_M: 9.904E+12
   max(B_n): 9.864E-01,  max(M): 9.265E+07,  min(M): 9.025E+07
 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 =  4.33431250552175E+11, method 2 =  4.33431250552174E+11, relative difference =  1.26736686741021E-15
   2 methods of computing chi2_B that should agree: method 1 =  1.17753991056211E+02, method 2 =  1.17753991056210E+02, relative difference =  3.01706434485485E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.178E+02,  chi2_M: 4.334E+11
   max(B_n): 9.972E-01,  max(M): 1.937E+07,  min(M): 1.889E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.86417230023428E+10, method 2 =  1.86417230023428E+10, relative difference =  7.16212789344072E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.18293805984161E+02, method 2 =  1.18293805984161E+02, relative difference = -7.80857077683179E-16
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 1.183E+02,  chi2_M: 1.864E+10
   max(B_n): 9.994E-01,  max(M): 4.017E+06,  min(M): 3.918E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  7.98871051863758E+08, method 2 =  7.98871051863759E+08, relative difference = -5.96888768332096E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.18405964697106E+02, method 2 =  1.18405964697106E+02, relative difference =  4.20063225956917E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.184E+02,  chi2_M: 7.989E+08
   max(B_n): 9.999E-01,  max(M): 8.316E+05,  min(M): 8.110E+05
 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 =  3.42090362904601E+07, method 2 =  3.42090362904601E+07, relative difference = -6.53387064195221E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.18429191805669E+02, method 2 =  1.18429191805669E+02, relative difference =  2.39989051660855E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.184E+02,  chi2_M: 3.421E+07
   max(B_n): 1.000E+00,  max(M): 1.721E+05,  min(M): 1.678E+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 =  1.46466201619432E+06, method 2 =  1.46466201619432E+06, relative difference = -3.97413603069397E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.18433998668763E+02, method 2 =  1.18433998668763E+02, relative difference =  2.99974139076125E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.184E+02,  chi2_M: 1.465E+06
   max(B_n): 1.000E+00,  max(M): 3.561E+04,  min(M): 3.473E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.27075762063740E+04, method 2 =  6.27075762063740E+04, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  1.18434993311662E+02, method 2 =  1.18434993311662E+02, relative difference = -4.19960267759050E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.184E+02,  chi2_M: 6.271E+04
   max(B_n): 1.000E+00,  max(M): 7.368E+03,  min(M): 7.186E+03
 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 =  2.68472444313827E+03, method 2 =  2.68472444313828E+03, relative difference = -3.38766499518210E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.18435199118577E+02, method 2 =  1.18435199118577E+02, relative difference = -2.39976878849574E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.184E+02,  chi2_M: 2.685E+03
   max(B_n): 1.000E+00,  max(M): 1.524E+03,  min(M): 1.487E+03
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   1.7480000257492065      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
