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+01 m^2. Volume: 1.479E+00 m^3.
 G =    0.0000000000000000       ; curpol =    1.0000000000000000     
 Done initializing plasma surface. Took    1.70000009E-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   0.00000000     sec.
 Building inductance matrix.
   Number of OpenMP threads:           4
 Done. Took   1.72300005     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  0.000E+00 sec, rest took  1.000E-03 sec.
  js=   1 j_RZetaZ= 3 dgemm took  0.000E+00 sec, rest took  0.000E+00 sec.
 Total time to compute g and RHS_B:   3.00000003E-03 sec.
 About to form matrix_B.
  js=   1, ks=   1, j_RZetaZ= 1, k_RZetaZ= 1 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 2, k_RZetaZ= 1 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 2, k_RZetaZ= 2 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 1 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 2 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 3 dgemm took  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 =  3.26704745947810E+31, method 2 =  3.26704745947809E+31, relative difference =  7.58170741556802E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.72351089531005E-16, method 2 = -3.55271367880050E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 2.724E-16,  chi2_M: 3.267E+31
   max(B_n): 7.353E-09,  max(M): 3.378E+17,  min(M): 2.881E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.48772252534041E+19, method 2 =  1.48772252534041E+19, relative difference =  2.75320157504692E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.65416665911293E-25, method 2 = -4.61852778244065E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.654E-25,  chi2_M: 1.488E+19
   max(B_n): 1.522E-13,  max(M): 1.573E+11,  min(M): 5.975E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.76438741400487E+19, method 2 =  4.76438741400487E+19, relative difference = -6.87769426635600E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.63777222713831E-25, method 2 = -3.90798504668055E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.638E-25,  chi2_M: 4.764E+19
   max(B_n): 1.392E-13,  max(M): 3.026E+11,  min(M): 1.921E+10
 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.34848654237883E+19, method 2 =  1.34848654237883E+19, relative difference = -7.59369832637400E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.04516865152199E-25, method 2 = -3.19744231092045E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.045E-25,  chi2_M: 1.348E+19
   max(B_n): 1.650E-13,  max(M): 1.213E+11,  min(M): 9.074E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.33857180291696E+19, method 2 =  1.33857180291696E+19, relative difference =  9.94492784846586E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.66076911911718E-25, method 2 = -3.90798504668055E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.661E-25,  chi2_M: 1.339E+19
   max(B_n): 1.749E-13,  max(M): 1.110E+11,  min(M): 1.037E+11
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.456E-32 (  6 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.33839668349541E+19, method 2 =  1.33839668349541E+19, relative difference =  4.59056726287911E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.05135515119129E-25, method 2 = -3.55271367880050E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.051E-25,  chi2_M: 1.338E+19
   max(B_n): 1.842E-13,  max(M): 1.096E+11,  min(M): 1.028E+11
 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 =  1.33792912719367E+19, method 2 =  1.33792912719366E+19, relative difference =  5.35753341063367E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.79838527180653E-25, method 2 = -2.84217094304040E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.798E-25,  chi2_M: 1.338E+19
   max(B_n): 2.889E-13,  max(M): 1.077E+11,  min(M): 1.047E+11
 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.33792179059252E+19, method 2 =  1.33792179059252E+19, relative difference = -2.29609833818426E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.23537688864242E-23, method 2 = -3.19744231092045E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.235E-23,  chi2_M: 1.338E+19
   max(B_n): 8.312E-13,  max(M): 1.075E+11,  min(M): 1.048E+11
 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 =  1.33792173665746E+19, method 2 =  1.33792173665746E+19, relative difference = -5.35756300507372E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.62482794184742E-22, method 2 = -3.55271367880050E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.625E-22,  chi2_M: 1.338E+19
   max(B_n): 3.673E-12,  max(M): 1.076E+11,  min(M): 1.049E+11
 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.33791991828158E+19, method 2 =  1.33791991828158E+19, relative difference = -4.59220310277713E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.06395817507761E-21, method 2 = -2.48689957516035E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 6.064E-21,  chi2_M: 1.338E+19
   max(B_n): 1.749E-11,  max(M): 1.076E+11,  min(M): 1.050E+11
 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.33791985810748E+19, method 2 =  1.33791985810748E+19, relative difference =  5.35757052753468E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.41364154598075E-19, method 2 = -3.55271367880050E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.414E-19,  chi2_M: 1.338E+19
   max(B_n): 8.428E-11,  max(M): 1.076E+11,  min(M): 1.050E+11
 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.33791985445414E+19, method 2 =  1.33791985445414E+19, relative difference =  4.59220332185494E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.30039331604339E-18, method 2 = -4.61852778244065E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.300E-18,  chi2_M: 1.338E+19
   max(B_n): 4.070E-10,  max(M): 1.076E+11,  min(M): 1.050E+11
 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.33791985068119E+19, method 2 =  1.33791985068119E+19, relative difference = -1.53073444493501E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.70328565265614E-17, method 2 = -3.90798504668055E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.703E-17,  chi2_M: 1.338E+19
   max(B_n): 1.966E-09,  max(M): 1.076E+11,  min(M): 1.050E+11
 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.33791983314732E+19, method 2 =  1.33791983314732E+19, relative difference = -1.07151412549704E-15
   2 methods of computing chi2_B that should agree: method 1 =  1.79369596747229E-15, method 2 = -2.48689957516035E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.794E-15,  chi2_M: 1.338E+19
   max(B_n): 9.472E-09,  max(M): 1.076E+11,  min(M): 1.050E+11
 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.33791974946219E+19, method 2 =  1.33791974946219E+19, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  4.13416772332822E-14, method 2 =  1.06581410364015E-14, relative difference =  5.90070065971162E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.134E-14,  chi2_M: 1.338E+19
   max(B_n): 4.521E-08,  max(M): 1.076E+11,  min(M): 1.050E+11
 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.33791936000684E+19, method 2 =  1.33791936000684E+19, relative difference =  4.59220501896958E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.28022243016250E-13, method 2 =  8.88178419700125E-13, relative difference =  2.19380072555050E-02
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.280E-13,  chi2_M: 1.338E+19
   max(B_n): 2.091E-07,  max(M): 1.076E+11,  min(M): 1.050E+11
 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.33791756299463E+19, method 2 =  1.33791756299463E+19, relative difference =  7.65368531158231E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.07406742236108E-11, method 2 =  2.07087680337281E-11, relative difference =  7.69761621510751E-04
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.074E-11,  chi2_M: 1.338E+19
   max(B_n): 9.340E-07,  max(M): 1.076E+11,  min(M): 1.050E+11
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.859E-24 ( 18 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.33790898870216E+19, method 2 =  1.33790898870216E+19, relative difference =  2.29612030858690E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.78986181105885E-10, method 2 =  4.78955541893811E-10, relative difference =  3.19844217438897E-05
   Diagnostics:    0.00000000      sec.
   chi2_B: 4.790E-10,  chi2_M: 1.338E+19
   max(B_n): 4.284E-06,  max(M): 1.076E+11,  min(M): 1.050E+11
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.281E-23 ( 19 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.33786758213384E+19, method 2 =  1.33786758213384E+19, relative difference = -1.53079424851115E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.11765642626312E-08, method 2 =  1.11765281474163E-08, relative difference =  1.61566973543876E-06
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.118E-08,  chi2_M: 1.338E+19
   max(B_n): 2.098E-05,  max(M): 1.076E+11,  min(M): 1.050E+11
 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.33766754771354E+19, method 2 =  1.33766754771354E+19, relative difference = -1.53102316304274E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.60929306967892E-07, method 2 =  2.60929276407751E-07, relative difference =  5.85601958075382E-08
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.609E-07,  chi2_M: 1.338E+19
   max(B_n): 1.018E-04,  max(M): 1.076E+11,  min(M): 1.050E+11
 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.33670245646418E+19, method 2 =  1.33670245646418E+19, relative difference = -6.89457848710622E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.08310138350099E-06, method 2 =  6.08310135774559E-06, relative difference =  2.11696246584516E-09
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.083E-06,  chi2_M: 1.337E+19
   max(B_n): 4.894E-04,  max(M): 1.076E+11,  min(M): 1.049E+11
 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.33206720830580E+19, method 2 =  1.33206720830580E+19, relative difference =  2.30618994360438E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.41136086291696E-04, method 2 =  1.41136086256211E-04, relative difference =  1.25715226124083E-10
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.411E-04,  chi2_M: 1.332E+19
   max(B_n): 2.314E-03,  max(M): 1.076E+11,  min(M): 1.046E+11
 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.31007388795812E+19, method 2 =  1.31007388795812E+19, relative difference =  1.56327060544044E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.23197448398343E-03, method 2 =  3.23197448394197E-03, relative difference =  6.41448858084598E-12
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.232E-03,  chi2_M: 1.310E+19
   max(B_n): 1.073E-02,  max(M): 1.069E+11,  min(M): 1.035E+11
 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.21116335785569E+19, method 2 =  1.21116335785569E+19, relative difference =  1.26820217110879E-15
   2 methods of computing chi2_B that should agree: method 1 =  6.97404921030637E-02, method 2 =  6.97404921030227E-02, relative difference =  2.94109917515978E-13
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.974E-02,  chi2_M: 1.211E+19
   max(B_n): 4.887E-02,  max(M): 1.029E+11,  min(M): 9.938E+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 =  8.61036861921595E+18, method 2 =  8.61036861921595E+18, relative difference =  1.18926383443644E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.15797770459021E+00, method 2 =  1.15797770459018E+00, relative difference =  1.39020240141514E-14
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.158E+00,  chi2_M: 8.610E+18
   max(B_n): 1.980E-01,  max(M): 8.680E+10,  min(M): 8.376E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.78601595451322E+18, method 2 =  2.78601595451322E+18, relative difference = -7.35099882210773E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.75148028146029E+00, method 2 =  8.75148028146029E+00, relative difference =  3.04466807146341E-16
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 8.751E+00,  chi2_M: 2.786E+18
   max(B_n): 5.437E-01,  max(M): 4.938E+10,  min(M): 4.764E+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 =  2.92957302032446E+17, method 2 =  2.92957302032446E+17, relative difference = -4.36923739780426E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.14942915767175E+01, method 2 =  2.14942915767175E+01, relative difference =  6.61145526219396E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.149E+01,  chi2_M: 2.930E+17
   max(B_n): 8.523E-01,  max(M): 1.601E+10,  min(M): 1.545E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.60996040894853E+16, method 2 =  1.60996040894853E+16, relative difference = -2.48453314613644E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.75901963294620E+01, method 2 =  2.75901963294620E+01, relative difference = -6.43836244653097E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.759E+01,  chi2_M: 1.610E+16
   max(B_n): 9.656E-01,  max(M): 3.751E+09,  min(M): 3.624E+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 =  7.28832702153690E+14, method 2 =  7.28832702153690E+14, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  2.91733460865589E+01, method 2 =  2.91733460865590E+01, relative difference = -6.08897187909025E-16
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 2.917E+01,  chi2_M: 7.288E+14
   max(B_n): 9.929E-01,  max(M): 7.969E+08,  min(M): 7.723E+08
 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 =  3.15727148699718E+13, method 2 =  3.15727148699717E+13, relative difference =  5.56750506644528E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.95179115558927E+01, method 2 =  2.95179115558927E+01, relative difference = -6.01789471466058E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.952E+01,  chi2_M: 3.157E+13
   max(B_n): 9.986E-01,  max(M): 1.656E+08,  min(M): 1.610E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.952E-15 ( 31 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.35504426770273E+12, method 2 =  1.35504426770273E+12, relative difference =  7.20686787344309E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.95899698493985E+01, method 2 =  2.95899698493985E+01, relative difference =  1.80097193249052E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.959E+01,  chi2_M: 1.355E+12
   max(B_n): 9.997E-01,  max(M): 3.427E+07,  min(M): 3.339E+07
 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 =  5.80433906216196E+10, method 2 =  5.80433906216195E+10, relative difference =  5.91493278096631E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.96049127692191E+01, method 2 =  2.96049127692191E+01, relative difference = -1.20004193442324E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.960E+01,  chi2_M: 5.804E+10
   max(B_n): 9.999E-01,  max(M): 7.091E+06,  min(M): 6.911E+06
 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 =  2.48529152446889E+09, method 2 =  2.48529152446889E+09, relative difference =  3.83727344263990E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.96080060860243E+01, method 2 =  2.96080060860243E+01, relative difference =  2.99979139804138E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.961E+01,  chi2_M: 2.485E+09
   max(B_n): 1.000E+00,  max(M): 1.467E+06,  min(M): 1.430E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.848E-13 ( 34 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.06405904110784E+08, method 2 =  1.06405904110784E+08, relative difference =  2.10061104950530E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.96086461968430E+01, method 2 =  2.96086461968430E+01, relative difference = -4.79956247263924E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.961E+01,  chi2_M: 1.064E+08
   max(B_n): 1.000E+00,  max(M): 3.036E+05,  min(M): 2.959E+05
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.793E-12 ( 35 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.55561093663468E+06, method 2 =  4.55561093663468E+06, relative difference =  1.02217088725257E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.96087786472161E+01, method 2 =  2.96087786472161E+01, relative difference =  3.59965575189425E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.961E+01,  chi2_M: 4.556E+06
   max(B_n): 1.000E+00,  max(M): 6.282E+04,  min(M): 6.123E+04
 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.95041028503541E+05, method 2 =  1.95041028503540E+05, relative difference =  5.96876066128936E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.96088060531387E+01, method 2 =  2.96088060531388E+01, relative difference = -4.19959449006003E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.961E+01,  chi2_M: 1.950E+05
   max(B_n): 1.000E+00,  max(M): 1.300E+04,  min(M): 1.267E+04
 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.35035635209070E+03, method 2 =  8.35035635209070E+03, relative difference =  3.26750618808699E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.96088117238074E+01, method 2 =  2.96088117238074E+01, relative difference = -3.59965173064736E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.961E+01,  chi2_M: 8.350E+03
   max(B_n): 1.000E+00,  max(M): 2.690E+03,  min(M): 2.621E+03
 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.57506526155090E+02, method 2 =  3.57506526155091E+02, relative difference = -3.17999335408768E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.96088128971472E+01, method 2 =  2.96088128971472E+01, relative difference = -4.79953545066686E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.961E+01,  chi2_M: 3.575E+02
   max(B_n): 1.000E+00,  max(M): 5.565E+02,  min(M): 5.424E+02
 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.53060428580016E+01, method 2 =  1.53060428580016E+01, relative difference = -5.80279584958702E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.96088131399274E+01, method 2 =  2.96088131399274E+01, relative difference =  1.79982577924223E-16
   Diagnostics:    3.00000003E-03  sec.
   chi2_B: 2.961E+01,  chi2_M: 1.531E+01
   max(B_n): 1.000E+00,  max(M): 1.151E+02,  min(M): 1.122E+02
 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 =  6.55302572959352E-01, method 2 =  6.55302572959352E-01, relative difference = -5.08264305881501E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.96088131901620E+01, method 2 =  2.96088131901620E+01, relative difference =  2.99970962698105E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.961E+01,  chi2_M: 6.553E-01
   max(B_n): 1.000E+00,  max(M): 2.383E+01,  min(M): 2.322E+01
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   1.8150000572204590      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
