This is REGCOIL_PM,
a regularized least-squares method for designing permanent magnets for stellarators.
 Successfully read parameters from regcoil_nml namelist in regcoil_in.benchmark.
 Resolution parameters:
   ntheta_plasma  =   96
   ntheta_coil    =   96
   nzeta_plasma   =    4
   nzeta_coil     =    4
   mpol_magnetization =   24
   ntor_magnetization =    0
   ns_magnetization =    1
   ns_integration   =    2
 Not assuming stellarator symmetry
 We will use the following values of the regularization weight lambda:
 0.000E+00 1.000E-34 4.833E-34 2.336E-33 1.129E-32 5.456E-32 2.637E-31 1.274E-30 6.158E-30 2.976E-29 1.438E-28 6.952E-28 3.360E-27 1.624E-26 7.848E-26 3.793E-25 1.833E-24 8.859E-24 4.281E-23 2.069E-22 1.000E-21 4.833E-21 2.336E-20 1.129E-19 5.456E-19 2.637E-18 1.274E-17 6.158E-17 2.976E-16 1.438E-15 6.952E-15 3.360E-14 1.624E-13 7.848E-13 3.793E-12 1.833E-11 8.859E-11 4.281E-10 2.069E-09 1.000E-08
 Initializing plasma surface.
 Reading FOCUS format data from file:plasma
 Number of modes for Bnormal read from FOCUS file:           1
 Plasma surface area: 3.948E+02 m^2. Volume: 6.575E+01 m^3.
 G =    0.0000000000000000       ; curpol =    1.0000000000000000     
 Done initializing plasma surface. Took    1.89999994E-02  sec.
 Initializing coil surface.
   Building a plain circular torus.
   Evaluating coil surface & derivatives:   9.99999978E-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.30000003E-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.62500000     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:   2.00000009E-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 =  5.33619979756588E+32, method 2 =  5.33619979756587E+32, relative difference =  5.40141649649604E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.20323180620744E-14, method 2 =  3.01270119962282E-12, relative difference = -9.46985859949944E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.203E-14,  chi2_M: 5.336E+32
   max(B_n): 4.437E-08,  max(M): 1.602E+18,  min(M): 1.649E+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 =  1.36408852117093E+23, method 2 =  1.36408852117094E+23, relative difference = -2.45984270662998E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.33396320210660E-20, method 2 = -2.09467998502078E-11, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.334E-20,  chi2_M: 1.364E+23
   max(B_n): 1.657E-11,  max(M): 1.658E+13,  min(M): 1.018E+12
 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 =  6.24986498240133E+20, method 2 =  6.24986498240134E+20, relative difference = -3.14579595805058E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.63080781678665E-23, method 2 = -1.10844666778576E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.631E-23,  chi2_M: 6.250E+20
   max(B_n): 5.008E-13,  max(M): 1.155E+12,  min(M): 2.332E+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 =  1.79989896879977E+19, method 2 =  1.79989896879977E+19, relative difference =  2.84460410764843E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.39650430841733E-25, method 2 = -3.97903932025656E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.397E-25,  chi2_M: 1.800E+19
   max(B_n): 1.690E-13,  max(M): 1.763E+11,  min(M): 4.538E+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.07434953383250E+18, method 2 =  1.07434953383250E+18, relative difference =  8.93563937776767E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.74392272224285E-26, method 2 = -3.12638803734444E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 9.744E-26,  chi2_M: 1.074E+18
   max(B_n): 5.618E-14,  max(M): 5.954E+10,  min(M): 8.996E+08
 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 =  6.76219334872461E+17, method 2 =  6.76219334872461E+17, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  6.21712862885310E-26, method 2 = -1.42108547152020E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 6.217E-26,  chi2_M: 6.762E+17
   max(B_n): 4.585E-14,  max(M): 4.137E+10,  min(M): 2.205E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.637E-31 (  7 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  7.92553905130628E+19, method 2 =  7.92553905130629E+19, relative difference = -8.26896436642986E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.55272281945558E-23, method 2 = -3.61012553184992E-10, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 9.553E-23,  chi2_M: 7.926E+19
   max(B_n): 1.172E-12,  max(M): 5.344E+11,  min(M): 2.019E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.274E-30 (  8 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.08499833365801E+17, method 2 =  2.08499833365800E+17, relative difference =  4.60432022655739E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.91494197462830E-26, method 2 = -1.13686837721616E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 5.915E-26,  chi2_M: 2.085E+17
   max(B_n): 4.419E-14,  max(M): 1.698E+10,  min(M): 8.666E+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.97895770488872E+17, method 2 =  1.97895770488872E+17, relative difference = -8.08506415295001E-17
   2 methods of computing chi2_B that should agree: method 1 =  6.55614389546098E-26, method 2 = -1.42108547152020E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 6.556E-26,  chi2_M: 1.979E+17
   max(B_n): 4.707E-14,  max(M): 1.399E+10,  min(M): 1.167E+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.97288063690049E+17, method 2 =  1.97288063690049E+17, relative difference =  4.05498429573949E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.34235393395281E-25, method 2 = -8.52651282912120E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.342E-25,  chi2_M: 1.973E+17
   max(B_n): 6.950E-14,  max(M): 1.367E+10,  min(M): 1.232E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.07307611238950E+17, method 2 =  3.07307611238950E+17, relative difference =  1.04130190173253E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.15452208515085E-24, method 2 = -1.57172053150134E-11, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.155E-24,  chi2_M: 3.073E+17
   max(B_n): 1.799E-13,  max(M): 2.298E+10,  min(M): 2.932E+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 =  1.97093187110604E+17, method 2 =  1.97093187110604E+17, relative difference =  1.62359747026884E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.52048634168657E-23, method 2 = -8.52651282912120E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.520E-23,  chi2_M: 1.971E+17
   max(B_n): 7.243E-13,  max(M): 1.304E+10,  min(M): 1.277E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.97093023753222E+17, method 2 =  1.97093023753221E+17, relative difference =  5.68259585586521E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.22227805442444E-21, method 2 = -5.68434188608080E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.222E-21,  chi2_M: 1.971E+17
   max(B_n): 3.410E-12,  max(M): 1.302E+10,  min(M): 1.276E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.97093017868346E+17, method 2 =  1.97093017868346E+17, relative difference =  4.05899716109874E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.18980368141718E-20, method 2 = -5.68434188608080E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.190E-20,  chi2_M: 1.971E+17
   max(B_n): 1.640E-11,  max(M): 1.302E+10,  min(M): 1.275E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.97093017642724E+17, method 2 =  1.97093017642724E+17, relative difference =  3.24719773259621E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.21218726087358E-18, method 2 = -2.84217094304040E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.212E-18,  chi2_M: 1.971E+17
   max(B_n): 7.921E-11,  max(M): 1.302E+10,  min(M): 1.275E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.97093017512867E+17, method 2 =  1.97093017512867E+17, relative difference = -2.43539830105176E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.83132473438426E-17, method 2 = -5.68434188608080E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.831E-17,  chi2_M: 1.971E+17
   max(B_n): 3.827E-10,  max(M): 1.302E+10,  min(M): 1.275E+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.97093016940191E+17, method 2 =  1.97093016940191E+17, relative difference =  1.21769915406404E-15
   2 methods of computing chi2_B that should agree: method 1 =  6.61317746628636E-16, method 2 = -1.13686837721616E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.613E-16,  chi2_M: 1.971E+17
   max(B_n): 1.850E-09,  max(M): 1.302E+10,  min(M): 1.275E+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.97093014174419E+17, method 2 =  1.97093014174419E+17, relative difference =  6.49439557947627E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.54465379040531E-14, method 2 = -5.68434188608080E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.545E-14,  chi2_M: 1.971E+17
   max(B_n): 8.939E-09,  max(M): 1.302E+10,  min(M): 1.275E+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.97093000807737E+17, method 2 =  1.97093000807737E+17, relative difference =  4.87079701494054E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.60788053769984E-13, method 2 =  3.12638803734444E-13, relative difference =  7.14988561845757E-02
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.608E-13,  chi2_M: 1.971E+17
   max(B_n): 4.320E-08,  max(M): 1.302E+10,  min(M): 1.275E+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.97092936207519E+17, method 2 =  1.97092936207519E+17, relative difference =  4.05899884284885E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.42700128531438E-12, method 2 =  8.38440428196918E-12, relative difference =  2.53381569879513E-03
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.427E-12,  chi2_M: 1.971E+17
   max(B_n): 2.088E-07,  max(M): 1.302E+10,  min(M): 1.275E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.97092623999618E+17, method 2 =  1.97092623999618E+17, relative difference = -2.43540316354472E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.96830965548824E-10, method 2 =  1.96735072677257E-10, relative difference =  2.43651286578648E-04
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.968E-10,  chi2_M: 1.971E+17
   max(B_n): 1.009E-06,  max(M): 1.302E+10,  min(M): 1.275E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.97091115131081E+17, method 2 =  1.97091115131081E+17, relative difference = -8.11807269412362E-17
   2 methods of computing chi2_B that should agree: method 1 =  4.59738788297902E-09, method 2 =  4.59729676549614E-09, relative difference =  9.90979966862198E-06
   Diagnostics:    0.00000000      sec.
   chi2_B: 4.597E-09,  chi2_M: 1.971E+17
   max(B_n): 4.877E-06,  max(M): 1.302E+10,  min(M): 1.275E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.336E-20 ( 23 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.97083823119201E+17, method 2 =  1.97083823119200E+17, relative difference =  4.87102383547416E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.07378198310670E-07, method 2 =  1.07378184566187E-07, relative difference =  6.40003486420858E-08
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.074E-07,  chi2_M: 1.971E+17
   max(B_n): 2.357E-05,  max(M): 1.302E+10,  min(M): 1.275E+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.97048587044444E+17, method 2 =  1.97048587044444E+17, relative difference =  4.87189486815998E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.50760671048470E-06, method 2 =  2.50760660946980E-06, relative difference =  2.01416947423700E-08
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.508E-06,  chi2_M: 1.970E+17
   max(B_n): 1.139E-04,  max(M): 1.302E+10,  min(M): 1.275E+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.96878426825203E+17, method 2 =  1.96878426825203E+17, relative difference = -2.43805280111348E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.85200723943804E-05, method 2 =  5.85200724003698E-05, relative difference = -5.11739272183616E-11
   Diagnostics:    0.00000000      sec.
   chi2_B: 5.852E-05,  chi2_M: 1.969E+17
   max(B_n): 5.502E-04,  max(M): 1.301E+10,  min(M): 1.275E+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.96059156331242E+17, method 2 =  1.96059156331242E+17, relative difference =  8.16080222898031E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.36117163641820E-03, method 2 =  1.36117163629024E-03, relative difference =  4.70052695313961E-11
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.361E-03,  chi2_M: 1.961E+17
   max(B_n): 2.653E-03,  max(M): 1.299E+10,  min(M): 1.272E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.274E-17 ( 27 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.92170924036864E+17, method 2 =  1.92170924036864E+17, relative difference =  1.66518427074126E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.11620286674755E-02, method 2 =  3.11620286673531E-02, relative difference =  1.96412913472775E-12
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.116E-02,  chi2_M: 1.922E+17
   max(B_n): 1.269E-02,  max(M): 1.286E+10,  min(M): 1.259E+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.74917097809018E+17, method 2 =  1.74917097809018E+17, relative difference =  3.65887616486056E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.62466345622882E-01, method 2 =  6.62466345622875E-01, relative difference =  5.61424313411315E-15
   Diagnostics:    0.00000000      sec.
   chi2_B: 6.625E-01,  chi2_M: 1.749E+17
   max(B_n): 5.835E-02,  max(M): 1.228E+10,  min(M): 1.200E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.976E-16 ( 29 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.17127689781956E+17, method 2 =  1.17127689781956E+17, relative difference = -6.83015264357451E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.03603875852612E+01, method 2 =  1.03603875852611E+01, relative difference =  6.08670933203502E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.036E+01,  chi2_M: 1.171E+17
   max(B_n): 2.297E-01,  max(M): 1.007E+10,  min(M): 9.801E+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 =  3.32062172254134E+16, method 2 =  3.32062172254135E+16, relative difference = -7.22756218724976E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.86028105977784E+01, method 2 =  6.86028105977784E+01, relative difference = -1.03573414787048E-16
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 6.860E+01,  chi2_M: 3.321E+16
   max(B_n): 5.898E-01,  max(M): 5.364E+09,  min(M): 5.216E+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 =  3.12527953601168E+15, method 2 =  3.12527953601169E+15, relative difference = -2.39978533554509E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.50808272429696E+02, method 2 =  1.50808272429696E+02, relative difference =  9.42312678624895E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.508E+02,  chi2_M: 3.125E+15
   max(B_n): 8.747E-01,  max(M): 1.643E+09,  min(M): 1.603E+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 =  1.65144775160013E+14, method 2 =  1.65144775160013E+14, relative difference = -1.89227906058312E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.86129530958017E+02, method 2 =  1.86129530958017E+02, relative difference =  3.05397099365331E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.861E+02,  chi2_M: 1.651E+14
   max(B_n): 9.713E-01,  max(M): 3.771E+08,  min(M): 3.690E+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 =  7.40661990310833E+12, method 2 =  7.40661990310832E+12, relative difference =  5.27399819499401E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.94979328976520E+02, method 2 =  1.94979328976520E+02, relative difference = -1.45767808206103E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.950E+02,  chi2_M: 7.407E+12
   max(B_n): 9.939E-01,  max(M): 7.982E+07,  min(M): 7.819E+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 =  3.20208658675102E+11, method 2 =  3.20208658675101E+11, relative difference =  4.76526435157466E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.96889193372882E+02, method 2 =  1.96889193372882E+02, relative difference =  0.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.969E+02,  chi2_M: 3.202E+11
   max(B_n): 9.987E-01,  max(M): 1.659E+07,  min(M): 1.626E+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.37369619843175E+10, method 2 =  1.37369619843175E+10, relative difference =  1.04135941865639E-15
   2 methods of computing chi2_B that should agree: method 1 =  1.97287874178609E+02, method 2 =  1.97287874178609E+02, relative difference = -1.44062120131485E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.973E+02,  chi2_M: 1.374E+10
   max(B_n): 9.997E-01,  max(M): 3.437E+06,  min(M): 3.368E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.833E-11 ( 36 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.88371415745284E+08, method 2 =  5.88371415745284E+08, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  1.97370517963258E+02, method 2 =  1.97370517963258E+02, relative difference =  1.08001348391715E-15
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.974E+02,  chi2_M: 5.884E+08
   max(B_n): 9.999E-01,  max(M): 7.113E+05,  min(M): 6.970E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.51923197553519E+07, method 2 =  2.51923197553519E+07, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  1.97387624588661E+02, method 2 =  1.97387624588661E+02, relative difference =  2.87978635840341E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.974E+02,  chi2_M: 2.519E+07
   max(B_n): 1.000E+00,  max(M): 1.472E+05,  min(M): 1.442E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.07858629394060E+06, method 2 =  1.07858629394060E+06, relative difference =  4.31732991531398E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97391164463417E+02, method 2 =  1.97391164463417E+02, relative difference =  2.87973471433382E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.974E+02,  chi2_M: 1.079E+06
   max(B_n): 1.000E+00,  max(M): 3.045E+04,  min(M): 2.984E+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 =  4.61780358355548E+04, method 2 =  4.61780358355548E+04, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  1.97391896924276E+02, method 2 =  1.97391896924276E+02, relative difference =  7.19931007130127E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 4.618E+04
   max(B_n): 1.000E+00,  max(M): 6.301E+03,  min(M): 6.175E+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 =  1.97703672756714E+03, method 2 =  1.97703672756715E+03, relative difference = -5.75036549075717E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.97392048481051E+02, method 2 =  1.97392048481050E+02, relative difference =  1.00790263611823E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 1.977E+03
   max(B_n): 1.000E+00,  max(M): 1.304E+03,  min(M): 1.278E+03
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   1.7059999704360962      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
