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: 8.290E+02 m^2. Volume: 2.900E+02 m^3.
 G =    0.0000000000000000       ; curpol =    1.0000000000000000     
 Done initializing plasma surface. Took    4.39999998E-02  sec.
 Initializing coil surface.
   Building a plain circular torus.
   Evaluating coil surface & derivatives:   1.99999996E-02  sec.
 s_integration:  0.21132486540518713       0.78867513459481287     
 s_weights:  0.50000000000000011       0.50000000000000011     
 s_magnetization  0.50000000000000000     
 s_magnetization_weights   1.0000000000000000     
 interpolate_magnetization_to_integration:
   1.0000000000000000     
   1.0000000000000000     
 Maximum allowed uniform d before singularity:  2.9000E+01
 Coil surface area: 1.184E+03 m^2. Volume: 5.918E+02 m^3.
 Done initializing coil surface. Took    2.40000002E-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   3.16700006     sec.
 Computing g and RHS_B:
  js=   1 j_RZetaZ= 1 dgemm took  0.000E+00 sec, rest took  1.000E-03 sec.
  js=   1 j_RZetaZ= 2 dgemm took  0.000E+00 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:   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 =  2.59625074728754E+22, method 2 =  2.59625074728754E+22, relative difference = -2.42328519561259E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.15290048228395E-19, method 2 = -1.70530256582424E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.153E-19,  chi2_M: 2.596E+22
   max(B_n): 2.084E-11,  max(M): 1.200E+13,  min(M): 2.885E+11
 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.12837259457997E+19, method 2 =  1.12837259457997E+19, relative difference = -4.53750828812524E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.29424862861441E-20, method 2 = -3.41060513164848E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 9.294E-20,  chi2_M: 1.128E+19
   max(B_n): 1.699E-11,  max(M): 1.073E+11,  min(M): 8.442E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.833E-34 (  3 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.50294510549881E+18, method 2 =  1.50294510549882E+18, relative difference = -2.55498353595924E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.99332970547495E-20, method 2 = -1.70530256582424E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 9.993E-20,  chi2_M: 1.503E+18
   max(B_n): 1.698E-11,  max(M): 5.316E+10,  min(M): 1.244E+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 =  3.05696160670587E+17, method 2 =  3.05696160670587E+17, relative difference =  1.04679103361336E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.74306247568849E-20, method 2 = -6.25277607468888E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.743E-20,  chi2_M: 3.057E+17
   max(B_n): 1.771E-11,  max(M): 2.723E+10,  min(M): 1.272E+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 =  6.98876823694211E+19, method 2 =  6.98876823694210E+19, relative difference =  7.03299899690281E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.20069901439162E-16, method 2 =  2.39879227592610E-11, relative difference = -9.99989989187339E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.201E-16,  chi2_M: 6.989E+19
   max(B_n): 6.894E-10,  max(M): 3.772E+11,  min(M): 4.681E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.456E-32 (  6 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.93237663235800E+16, method 2 =  6.93237663235800E+16, relative difference =  3.46201617032405E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.70225112119673E-19, method 2 = -2.27373675443232E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.702E-19,  chi2_M: 6.932E+16
   max(B_n): 2.773E-11,  max(M): 1.562E+10,  min(M): 4.711E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.637E-31 (  7 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.69945158665884E+18, method 2 =  1.69945158665884E+18, relative difference = -3.76592075363709E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.75774994486790E-18, method 2 =  2.33058017329313E-12, relative difference = -9.99995917119448E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.758E-18,  chi2_M: 1.699E+18
   max(B_n): 1.129E-10,  max(M): 5.844E+10,  min(M): 2.094E+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.41004258609200E+21, method 2 =  1.41004258609200E+21, relative difference = -1.85912115411028E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.71408538701409E-15, method 2 = -3.78452114091488E-09, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.714E-15,  chi2_M: 1.410E+21
   max(B_n): 2.126E-09,  max(M): 1.902E+12,  min(M): 2.220E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.13521998907136E+15, method 2 =  2.13521998907136E+15, relative difference = -7.61045704104117E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.94731495891077E-20, method 2 = -4.54747350886464E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 6.947E-20,  chi2_M: 2.135E+15
   max(B_n): 1.359E-11,  max(M): 2.891E+09,  min(M): 2.418E+08
 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 =  5.36599620589579E+16, method 2 =  5.36599620589578E+16, relative difference =  5.96347794000312E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.73612750841250E-20, method 2 = -8.81072992342524E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 7.736E-20,  chi2_M: 5.366E+16
   max(B_n): 1.431E-11,  max(M): 1.036E+10,  min(M): 1.219E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.75032995422313E+14, method 2 =  4.75032995422313E+14, relative difference =  5.26279231988391E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.50155627023425E-21, method 2 = -3.41060513164848E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.502E-21,  chi2_M: 4.750E+14
   max(B_n): 2.576E-12,  max(M): 9.760E+08,  min(M): 1.416E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.952E-28 ( 12 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.54463151847910E+14, method 2 =  3.54463151847910E+14, relative difference =  7.05291928643878E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.64851234443590E-22, method 2 = -1.13686837721616E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.649E-22,  chi2_M: 3.545E+14
   max(B_n): 8.449E-13,  max(M): 6.487E+08,  min(M): 4.401E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.360E-27 ( 13 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.47748658980779E+14, method 2 =  3.47748658980780E+14, relative difference = -8.08773787436935E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.22804378001109E-23, method 2 = -1.13686837721616E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.228E-23,  chi2_M: 3.477E+14
   max(B_n): 1.914E-13,  max(M): 5.572E+08,  min(M): 5.236E+08
 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.47603234493579E+14, method 2 =  3.47603234493578E+14, relative difference =  4.49506749347830E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.21248715361847E-25, method 2 = -2.27373675443232E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.212E-25,  chi2_M: 3.476E+14
   max(B_n): 5.951E-14,  max(M): 5.452E+08,  min(M): 5.362E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.848E-26 ( 15 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.47596559205226E+14, method 2 =  3.47596559205226E+14, relative difference = -4.49515381732383E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.86608015112045E-24, method 2 = -2.27373675443232E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.866E-24,  chi2_M: 3.476E+14
   max(B_n): 8.626E-14,  max(M): 5.445E+08,  min(M): 5.383E+08
 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 =  3.47596287478129E+14, method 2 =  3.47596287478128E+14, relative difference =  6.29322026386039E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.32410966692305E-23, method 2 = -3.41060513164848E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.324E-23,  chi2_M: 3.476E+14
   max(B_n): 3.642E-13,  max(M): 5.442E+08,  min(M): 5.384E+08
 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.47596276242888E+14, method 2 =  3.47596276242887E+14, relative difference =  5.39418897194922E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.01133779897356E-21, method 2 = -1.70530256582424E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.011E-21,  chi2_M: 3.476E+14
   max(B_n): 1.741E-12,  max(M): 5.442E+08,  min(M): 5.384E+08
 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.47596275424619E+14, method 2 =  3.47596275424619E+14, relative difference = -4.49515748720630E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.36282962867911E-20, method 2 = -2.27373675443232E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.363E-20,  chi2_M: 3.476E+14
   max(B_n): 8.403E-12,  max(M): 5.442E+08,  min(M): 5.384E+08
 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.47596275379095E+14, method 2 =  3.47596275379095E+14, relative difference = -3.59612599023602E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.51945148926014E-19, method 2 = -3.41060513164848E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.519E-19,  chi2_M: 3.476E+14
   max(B_n): 4.060E-11,  max(M): 5.442E+08,  min(M): 5.384E+08
 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.47596275278424E+14, method 2 =  3.47596275278424E+14, relative difference = -2.69709449345815E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.28920144793743E-17, method 2 = -1.70530256582424E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.289E-17,  chi2_M: 3.476E+14
   max(B_n): 1.962E-10,  max(M): 5.442E+08,  min(M): 5.384E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-21 ( 21 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.47596274800740E+14, method 2 =  3.47596274800740E+14, relative difference =  8.99031499054876E-17
   2 methods of computing chi2_B that should agree: method 1 =  3.01115899419170E-16, method 2 = -2.84217094304040E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.011E-16,  chi2_M: 3.476E+14
   max(B_n): 9.482E-10,  max(M): 5.442E+08,  min(M): 5.384E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.833E-21 ( 22 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.47596272492635E+14, method 2 =  3.47596272492634E+14, relative difference =  1.79806301004923E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.03250463938949E-15, method 2 = -1.13686837721616E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 7.033E-15,  chi2_M: 3.476E+14
   max(B_n): 4.582E-09,  max(M): 5.442E+08,  min(M): 5.384E+08
 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.47596261342735E+14, method 2 =  3.47596261342734E+14, relative difference =  7.19225227090394E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.64178439706619E-13, method 2 = -5.68434188608080E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.642E-13,  chi2_M: 3.476E+14
   max(B_n): 2.212E-08,  max(M): 5.442E+08,  min(M): 5.384E+08
 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.47596207565944E+14, method 2 =  3.47596207565944E+14, relative difference =  3.59612669181052E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.82616876626156E-12, method 2 =  3.52429196937010E-12, relative difference =  4.10690986250839E-02
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.826E-12,  chi2_M: 3.476E+14
   max(B_n): 1.065E-07,  max(M): 5.442E+08,  min(M): 5.384E+08
 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.47595949623057E+14, method 2 =  3.47595949623057E+14, relative difference = -8.09129106092851E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.86469790325193E-11, method 2 =  8.83915163285565E-11, relative difference =  1.44297828244518E-03
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.865E-11,  chi2_M: 3.476E+14
   max(B_n): 5.063E-07,  max(M): 5.442E+08,  min(M): 5.384E+08
 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.47594720690518E+14, method 2 =  3.47594720690518E+14, relative difference =  1.79807103732300E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.04125844036461E-09, method 2 =  2.04101979761617E-09, relative difference =  5.84582271293045E-05
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.041E-09,  chi2_M: 3.476E+14
   max(B_n): 2.350E-06,  max(M): 5.441E+08,  min(M): 5.384E+08
 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.47588824888948E+14, method 2 =  3.47588824888949E+14, relative difference = -7.19240614481414E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.73666215631911E-08, method 2 =  4.73663135380775E-08, relative difference =  3.25151029312282E-06
   Diagnostics:    4.00000019E-03  sec.
   chi2_B: 4.737E-08,  chi2_M: 3.476E+14
   max(B_n): 1.092E-05,  max(M): 5.441E+08,  min(M): 5.384E+08
 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 =  3.47560354267286E+14, method 2 =  3.47560354267286E+14, relative difference = -4.49562207201107E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.10542057343960E-06, method 2 =  1.10542032416561E-06, relative difference =  1.12750760263725E-07
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.105E-06,  chi2_M: 3.476E+14
   max(B_n): 5.200E-05,  max(M): 5.441E+08,  min(M): 5.384E+08
 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 =  3.47422809526091E+14, method 2 =  3.47422809526091E+14, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  2.58082204612508E-05, method 2 =  2.58082203004051E-05, relative difference =  3.11617189442762E-09
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.581E-05,  chi2_M: 3.474E+14
   max(B_n): 2.503E-04,  max(M): 5.440E+08,  min(M): 5.383E+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.46759214994381E+14, method 2 =  3.46759214994381E+14, relative difference = -2.70360515153200E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.01655564610094E-04, method 2 =  6.01655564480552E-04, relative difference =  1.07654522508479E-10
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.017E-04,  chi2_M: 3.468E+14
   max(B_n): 1.210E-03,  max(M): 5.435E+08,  min(M): 5.378E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.43578751879304E+14, method 2 =  3.43578751879304E+14, relative difference =  7.27635217930539E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.39241019721298E-02, method 2 =  1.39241019719520E-02, relative difference =  6.38425640254512E-12
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.392E-02,  chi2_M: 3.436E+14
   max(B_n): 5.820E-03,  max(M): 5.410E+08,  min(M): 5.353E+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 =  3.28807706384706E+14, method 2 =  3.28807706384707E+14, relative difference = -1.90080703056499E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.11245961850935E-01, method 2 =  3.11245961850659E-01, relative difference =  4.42935687464167E-13
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.112E-01,  chi2_M: 3.288E+14
   max(B_n): 2.752E-02,  max(M): 5.293E+08,  min(M): 5.236E+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 =  2.69274096225684E+14, method 2 =  2.69274096225685E+14, relative difference = -8.12369266357729E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.95355656664013E+00, method 2 =  5.95355656664003E+00, relative difference =  8.20514805506845E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.954E+00,  chi2_M: 2.693E+14
   max(B_n): 1.203E-01,  max(M): 4.792E+08,  min(M): 4.737E+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.26437345626268E+14, method 2 =  1.26437345626268E+14, relative difference =  2.47157988371338E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.52943872563592E+01, method 2 =  6.52943872563591E+01, relative difference =  3.26464233274920E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.529E+01,  chi2_M: 1.264E+14
   max(B_n): 3.979E-01,  max(M): 3.287E+08,  min(M): 3.242E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.793E-12 ( 35 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.98908342782013E+13, method 2 =  1.98908342782014E+13, relative difference = -1.47288316770639E-15
   2 methods of computing chi2_B that should agree: method 1 =  2.39922851616398E+02, method 2 =  2.39922851616398E+02, relative difference =  2.96154672626247E-16
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 2.399E+02,  chi2_M: 1.989E+13
   max(B_n): 7.615E-01,  max(M): 1.306E+08,  min(M): 1.284E+08
 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.29706530305611E+12, method 2 =  1.29706530305611E+12, relative difference =  3.76450783819076E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.65426640450497E+02, method 2 =  3.65426640450498E+02, relative difference = -1.24442856800438E-15
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 3.654E+02,  chi2_M: 1.297E+12
   max(B_n): 9.391E-01,  max(M): 3.336E+07,  min(M): 3.278E+07
 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.13302908929718E+10, method 2 =  6.13302908929717E+10, relative difference =  4.97593891707715E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.03584151031908E+02, method 2 =  4.03584151031909E+02, relative difference = -5.63386036993448E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.036E+02,  chi2_M: 6.133E+10
   max(B_n): 9.868E-01,  max(M): 7.256E+06,  min(M): 7.126E+06
 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.68196379551655E+09, method 2 =  2.68196379551655E+09, relative difference =  8.88970162461283E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.12223710736995E+02, method 2 =  4.12223710736995E+02, relative difference = -6.20525647145119E-16
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 4.122E+02,  chi2_M: 2.682E+09
   max(B_n): 9.972E-01,  max(M): 1.517E+06,  min(M): 1.490E+06
 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 =  1.15331375481743E+08, method 2 =  1.15331375481743E+08, relative difference = -5.81413558038499E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.14045974223905E+02, method 2 =  4.14045974223905E+02, relative difference = -6.17794643102384E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.140E+02,  chi2_M: 1.153E+08
   max(B_n): 9.994E-01,  max(M): 3.147E+05,  min(M): 3.090E+05
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-08 ( 40 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.94223411154483E+06, method 2 =  4.94223411154483E+06, relative difference = -9.42208071891974E-17
   2 methods of computing chi2_B that should agree: method 1 =  4.14424534269596E+02, method 2 =  4.14424534269596E+02, relative difference =  0.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.144E+02,  chi2_M: 4.942E+06
   max(B_n): 9.999E-01,  max(M): 6.514E+04,  min(M): 6.397E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   3.2980000972747803      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
