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: 1.777E+02 m^2. Volume: 1.331E+01 m^3.
 G =    0.0000000000000000       ; curpol =    1.0000000000000000     
 Done initializing plasma surface. Took    1.60000008E-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.40000004E-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.66199994     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  1.000E-03 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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  7.89828247427443E+32, method 2 =  7.89828247427442E+32, relative difference =  4.56159894715259E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.21596318339783E-14, method 2 =  2.24531504500192E-12, relative difference = -9.71758516429559E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.216E-14,  chi2_M: 7.898E+32
   max(B_n): 4.439E-08,  max(M): 1.725E+18,  min(M): 1.514E+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 =  8.79091775980405E+23, method 2 =  8.79091775980405E+23, relative difference = -1.52677719968787E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.65136081744723E-19, method 2 =  4.11731093663548E-10, relative difference = -9.99999999197845E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.651E-19,  chi2_M: 8.791E+23
   max(B_n): 8.983E-11,  max(M): 3.866E+13,  min(M): 5.882E+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 =  1.49127613948924E+23, method 2 =  1.49127613948924E+23, relative difference = -1.12502410222604E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.26749432433713E-21, method 2 = -1.13473674900888E-10, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.267E-21,  chi2_M: 1.491E+23
   max(B_n): 1.417E-11,  max(M): 1.811E+13,  min(M): 3.536E+12
 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.00276072832982E+20, method 2 =  1.00276072832982E+20, relative difference = -2.45083391338363E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.30604808498312E-25, method 2 =  9.52127265918534E-13, relative difference = -9.99999999998045E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.306E-25,  chi2_M: 1.003E+20
   max(B_n): 2.308E-13,  max(M): 5.488E+11,  min(M): 5.179E+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 =  2.52785090020351E+18, method 2 =  2.52785090020351E+18, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  7.82535749830903E-26, method 2 =  3.69482222595252E-13, relative difference = -9.99999999999576E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.825E-26,  chi2_M: 2.528E+18
   max(B_n): 5.262E-14,  max(M): 6.947E+10,  min(M): 6.730E+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 =  5.35034594179762E+18, method 2 =  5.35034594179762E+18, relative difference = -2.87084240291934E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.52960188757300E-26, method 2 =  4.26325641456060E-13, relative difference = -9.99999999999600E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.530E-26,  chi2_M: 5.350E+18
   max(B_n): 5.440E-14,  max(M): 8.533E+10,  min(M): 4.438E+10
 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.48122329354701E+17, method 2 =  7.48122329354701E+17, relative difference = -1.71095013445739E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.68412885319043E-26, method 2 = -2.13162820728030E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.684E-26,  chi2_M: 7.481E+17
   max(B_n): 4.685E-14,  max(M): 3.451E+10,  min(M): 9.573E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.78985790611867E+17, method 2 =  1.78985790611867E+17, relative difference = -2.68177713079406E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.21325941302639E-26, method 2 =  0.00000000000000E+00, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.213E-26,  chi2_M: 1.790E+17
   max(B_n): 4.285E-14,  max(M): 1.556E+10,  min(M): 8.012E+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.65664986444893E+17, method 2 =  1.65664986444893E+17, relative difference = -4.82902281989511E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.65630650500226E-26, method 2 =  5.68434188608080E-14, relative difference = -9.99999999998362E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 4.656E-26,  chi2_M: 1.657E+17
   max(B_n): 5.418E-14,  max(M): 1.258E+10,  min(M): 1.101E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.65192928067090E+17, method 2 =  1.65192928067090E+17, relative difference =  2.90569339508928E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.85662105714750E-25, method 2 =  5.68434188608080E-14, relative difference = -9.99999999986431E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.857E-25,  chi2_M: 1.652E+17
   max(B_n): 1.132E-13,  max(M): 1.213E+10,  min(M): 1.154E+10
 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 =  1.65159264279590E+17, method 2 =  1.65159264279589E+17, relative difference =  9.68761883857417E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.14014959913608E-24, method 2 =  5.68434188608080E-14, relative difference = -9.99999999713594E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 8.140E-24,  chi2_M: 1.652E+17
   max(B_n): 3.969E-13,  max(M): 1.199E+10,  min(M): 1.168E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.952E-28 ( 12 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.65157758261824E+17, method 2 =  1.65157758261824E+17, relative difference = -9.68770717669541E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.88436173956053E-22, method 2 =  4.26325641456060E-14, relative difference = -9.99999991159989E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.884E-22,  chi2_M: 1.652E+17
   max(B_n): 1.798E-12,  max(M): 1.196E+10,  min(M): 1.167E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.360E-27 ( 13 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.65157707073840E+17, method 2 =  1.65157707073840E+17, relative difference =  3.87508407169799E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.39424682710568E-21, method 2 =  8.52651282912120E-14, relative difference = -9.99999896927462E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 4.394E-21,  chi2_M: 1.652E+17
   max(B_n): 8.633E-12,  max(M): 1.195E+10,  min(M): 1.167E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.624E-26 ( 14 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.65157704583924E+17, method 2 =  1.65157704583924E+17, relative difference =  1.93754206505936E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.02599450413678E-19, method 2 =  2.84217094304040E-14, relative difference = -9.99992780231584E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.026E-19,  chi2_M: 1.652E+17
   max(B_n): 4.166E-11,  max(M): 1.195E+10,  min(M): 1.167E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.848E-26 ( 15 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.65157704417389E+17, method 2 =  1.65157704417389E+17, relative difference =  1.93754206701306E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.39557548208622E-18, method 2 =  5.68434188608080E-14, relative difference = -9.99915716732856E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.396E-18,  chi2_M: 1.652E+17
   max(B_n): 2.013E-10,  max(M): 1.195E+10,  min(M): 1.167E+10
 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 =  1.65157704176593E+17, method 2 =  1.65157704176592E+17, relative difference =  1.93754206983796E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.58697621883262E-17, method 2 =  7.10542735760100E-14, relative difference = -9.98428641553354E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.587E-17,  chi2_M: 1.652E+17
   max(B_n): 9.718E-10,  max(M): 1.195E+10,  min(M): 1.167E+10
 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.65157703054209E+17, method 2 =  1.65157703054208E+17, relative difference =  9.68771041502583E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.29581144944397E-15, method 2 =  5.68434188608080E-14, relative difference = -9.55423852619684E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.296E-15,  chi2_M: 1.652E+17
   max(B_n): 4.675E-09,  max(M): 1.195E+10,  min(M): 1.167E+10
 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.65157697777573E+17, method 2 =  1.65157697777572E+17, relative difference =  7.75016857963139E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.93713959405486E-14, method 2 =  9.94759830064140E-14, relative difference = -5.44090129258450E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.937E-14,  chi2_M: 1.652E+17
   max(B_n): 2.215E-08,  max(M): 1.195E+10,  min(M): 1.167E+10
 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.65157674298449E+17, method 2 =  1.65157674298449E+17, relative difference =  5.81262726105736E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.29136088137996E-13, method 2 =  6.67910171614494E-13, relative difference = -2.98941407717388E-02
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.291E-13,  chi2_M: 1.652E+17
   max(B_n): 1.000E-07,  max(M): 1.195E+10,  min(M): 1.167E+10
 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.65157570982617E+17, method 2 =  1.65157570982617E+17, relative difference =  7.75017452959949E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.34508647065449E-11, method 2 =  1.35287336888723E-11, relative difference = -2.88621725149862E-03
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.345E-11,  chi2_M: 1.652E+17
   max(B_n): 4.309E-07,  max(M): 1.195E+10,  min(M): 1.167E+10
 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.65157082926144E+17, method 2 =  1.65157082926144E+17, relative difference = -3.87509871608836E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.07850308922832E-10, method 2 =  3.07878167404851E-10, relative difference = -4.52447516894804E-05
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.079E-10,  chi2_M: 1.652E+17
   max(B_n): 1.939E-06,  max(M): 1.195E+10,  min(M): 1.167E+10
 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 =  1.65154726741329E+17, method 2 =  1.65154726741329E+17, relative difference = -2.90636550022448E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.17951667865199E-09, method 2 =  7.17955117579550E-09, relative difference = -2.40246399362986E-06
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.180E-09,  chi2_M: 1.652E+17
   max(B_n): 9.183E-06,  max(M): 1.195E+10,  min(M): 1.167E+10
 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 =  1.65143340378858E+17, method 2 =  1.65143340378858E+17, relative difference =  3.87542118581207E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.67669250058159E-07, method 2 =  1.67669284678595E-07, relative difference = -1.03240255586275E-07
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.677E-07,  chi2_M: 1.651E+17
   max(B_n): 4.459E-05,  max(M): 1.195E+10,  min(M): 1.167E+10
 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 =  1.65088327861023E+17, method 2 =  1.65088327861023E+17, relative difference =  3.87671259556746E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.91493984936432E-06, method 2 =  3.91493993845415E-06, relative difference = -1.13781872855238E-08
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.915E-06,  chi2_M: 1.651E+17
   max(B_n): 2.157E-04,  max(M): 1.195E+10,  min(M): 1.167E+10
 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.64822851671227E+17, method 2 =  1.64822851671227E+17, relative difference =  1.94147836149751E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.12918440313390E-05, method 2 =  9.12918440576505E-05, relative difference = -1.44106370111019E-10
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.129E-05,  chi2_M: 1.648E+17
   max(B_n): 1.041E-03,  max(M): 1.194E+10,  min(M): 1.166E+10
 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.63548972004771E+17, method 2 =  1.63548972004771E+17, relative difference =  1.95660049756023E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.11554062431264E-03, method 2 =  2.11554062428831E-03, relative difference =  5.74926175922461E-12
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.116E-03,  chi2_M: 1.635E+17
   max(B_n): 5.002E-03,  max(M): 1.190E+10,  min(M): 1.160E+10
 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.57597783837718E+17, method 2 =  1.57597783837718E+17, relative difference =  1.72591259455834E-15
   2 methods of computing chi2_B that should agree: method 1 =  4.75977153791374E-02, method 2 =  4.75977153792115E-02, relative difference = -7.78694975010023E-13
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.760E-02,  chi2_M: 1.576E+17
   max(B_n): 2.355E-02,  max(M): 1.170E+10,  min(M): 1.138E+10
 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.32960063920032E+17, method 2 =  1.32960063920032E+17, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  9.37655700343551E-01, method 2 =  9.37655700343569E-01, relative difference = -9.70913822482027E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.377E-01,  chi2_M: 1.330E+17
   max(B_n): 1.034E-01,  max(M): 1.077E+10,  min(M): 1.043E+10
 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 =  6.84424372762292E+16, method 2 =  6.84424372762291E+16, relative difference =  5.25989450882738E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.12729905199775E+01, method 2 =  1.12729905199776E+01, relative difference = -1.65455180509987E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.127E+01,  chi2_M: 6.844E+16
   max(B_n): 3.566E-01,  max(M): 7.732E+09,  min(M): 7.474E+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.22319549516479E+16, method 2 =  1.22319549516479E+16, relative difference =  8.17530806770410E-17
   2 methods of computing chi2_B that should agree: method 1 =  4.70573059985824E+01, method 2 =  4.70573059985825E+01, relative difference = -2.26492800856950E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.706E+01,  chi2_M: 1.223E+16
   max(B_n): 7.283E-01,  max(M): 3.268E+09,  min(M): 3.160E+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 =  8.51657292671382E+14, method 2 =  8.51657292671383E+14, relative difference = -5.87090610627728E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.65270107045380E+01, method 2 =  7.65270107045380E+01, relative difference =  2.78545860821648E-16
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 7.653E+01,  chi2_M: 8.517E+14
   max(B_n): 9.288E-01,  max(M): 8.613E+08,  min(M): 8.349E+08
 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 =  4.09998095308107E+13, method 2 =  4.09998095308106E+13, relative difference =  6.66923829961981E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.60493423569926E+01, method 2 =  8.60493423569927E+01, relative difference = -6.60590973780856E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.605E+01,  chi2_M: 4.100E+13
   max(B_n): 9.845E-01,  max(M): 1.886E+08,  min(M): 1.835E+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 =  1.80005902756078E+12, method 2 =  1.80005902756078E+12, relative difference =  1.01721924640506E-15
   2 methods of computing chi2_B that should agree: method 1 =  8.82408921724305E+01, method 2 =  8.82408921724305E+01, relative difference = -1.61046135927918E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.824E+01,  chi2_M: 1.800E+12
   max(B_n): 9.968E-01,  max(M): 3.949E+07,  min(M): 3.848E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  7.74718703556809E+10, method 2 =  7.74718703556809E+10, relative difference =  2.95438634547845E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.87048042734600E+01, method 2 =  8.87048042734601E+01, relative difference = -4.00509725250975E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.870E+01,  chi2_M: 7.747E+10
   max(B_n): 9.993E-01,  max(M): 8.191E+06,  min(M): 7.985E+06
 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 =  3.32043649363265E+09, method 2 =  3.32043649363266E+09, relative difference = -3.59016923767040E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.88012510490764E+01, method 2 =  8.88012510490764E+01, relative difference =  2.40044839695135E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.880E+01,  chi2_M: 3.320E+09
   max(B_n): 9.999E-01,  max(M): 1.696E+06,  min(M): 1.653E+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 =  1.42190919099764E+08, method 2 =  1.42190919099764E+08, relative difference =  5.23984277202423E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.88212268690447E+01, method 2 =  8.88212268690448E+01, relative difference = -3.99984756350924E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.882E+01,  chi2_M: 1.422E+08
   max(B_n): 1.000E+00,  max(M): 3.509E+05,  min(M): 3.421E+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 =  6.08794938525818E+06, method 2 =  6.08794938525818E+06, relative difference =  3.05956083298147E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.88253609840915E+01, method 2 =  8.88253609840916E+01, relative difference = -2.39979684142468E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.883E+01,  chi2_M: 6.088E+06
   max(B_n): 1.000E+00,  max(M): 7.261E+04,  min(M): 7.079E+04
 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 =  2.60647860968575E+05, method 2 =  2.60647860968576E+05, relative difference = -4.46638316517664E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.88262164256376E+01, method 2 =  8.88262164256376E+01, relative difference = -2.39977373016313E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 8.883E+01,  chi2_M: 2.606E+05
   max(B_n): 1.000E+00,  max(M): 1.502E+04,  min(M): 1.465E+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 =  1.11592235671334E+04, method 2 =  1.11592235671334E+04, relative difference = -4.07508056587203E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.88263934298492E+01, method 2 =  8.88263934298493E+01, relative difference = -6.39938386170099E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 8.883E+01,  chi2_M: 1.116E+04
   max(B_n): 1.000E+00,  max(M): 3.109E+03,  min(M): 3.031E+03
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-08 ( 40 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.77763651121338E+02, method 2 =  4.77763651121339E+02, relative difference = -7.73357750695467E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.88264300545309E+01, method 2 =  8.88264300545309E+01, relative difference = -7.99922652890468E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.883E+01,  chi2_M: 4.778E+02
   max(B_n): 1.000E+00,  max(M): 6.432E+02,  min(M): 6.271E+02
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   1.7439999580383301      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
