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.553E+02 m^2. Volume: 5.326E+01 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:   1.09999999E-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.64199996     sec.
 Computing g and RHS_B:
  js=   1 j_RZetaZ= 1 dgemm took  0.000E+00 sec, rest took  0.000E+00 sec.
  js=   1 j_RZetaZ= 2 dgemm took  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  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:   0.00000000     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 =  1.57594801183654E+33, method 2 =  1.57594801183654E+33, relative difference = -2.74339991535464E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.24511900801168E-13, method 2 = -5.40012479177676E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.245E-13,  chi2_M: 1.576E+33
   max(B_n): 8.082E-08,  max(M): 2.991E+18,  min(M): 1.414E+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 =  6.82911857320491E+19, method 2 =  6.82911857320490E+19, relative difference =  5.39806119996826E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.44142326021924E-23, method 2 = -6.25277607468888E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.441E-23,  chi2_M: 6.829E+19
   max(B_n): 5.441E-13,  max(M): 4.587E+11,  min(M): 3.292E+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 =  1.05746722942156E+19, method 2 =  1.05746722942156E+19, relative difference = -6.77846064688068E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.52532668482415E-24, method 2 = -8.52651282912120E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.525E-24,  chi2_M: 1.057E+19
   max(B_n): 1.875E-13,  max(M): 1.556E+11,  min(M): 1.860E+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 =  6.17812447023507E+19, method 2 =  6.17812447023507E+19, relative difference =  4.64089063568334E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.57395984972637E-24, method 2 = -1.42108547152020E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.574E-24,  chi2_M: 6.178E+19
   max(B_n): 2.682E-13,  max(M): 3.787E+11,  min(M): 3.815E+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 =  4.64829020345351E+17, method 2 =  4.64829020345351E+17, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  4.26310215605973E-26, method 2 = -2.27373675443232E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.263E-26,  chi2_M: 4.648E+17
   max(B_n): 3.364E-14,  max(M): 3.072E+10,  min(M): 7.230E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.456E-32 (  6 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.52657591784862E+18, method 2 =  1.52657591784863E+18, relative difference = -5.86934452144847E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.43985827296492E-25, method 2 = -5.68434188608080E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.440E-25,  chi2_M: 1.527E+18
   max(B_n): 7.483E-14,  max(M): 7.020E+10,  min(M): 6.662E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.51311380240948E+18, method 2 =  1.51311380240948E+18, relative difference =  4.22968846745608E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.02130605187392E-25, method 2 = -7.95807864051312E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.021E-25,  chi2_M: 1.513E+18
   max(B_n): 4.918E-14,  max(M): 6.033E+10,  min(M): 3.325E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.08333794075808E+16, method 2 =  6.08333794075808E+16, relative difference = -3.28766874284608E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.73248270004911E-26, method 2 =  2.84217094304040E-14, relative difference = -9.99999999998781E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.732E-26,  chi2_M: 6.083E+16
   max(B_n): 2.023E-14,  max(M): 1.036E+10,  min(M): 2.566E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.158E-30 (  9 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  8.45270647208569E+16, method 2 =  8.45270647208571E+16, relative difference = -7.57153938934875E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.96312422579638E-25, method 2 = -9.37916411203332E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.963E-25,  chi2_M: 8.453E+16
   max(B_n): 4.996E-14,  max(M): 1.415E+10,  min(M): 2.582E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.976E-29 ( 10 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.05792770261272E+16, method 2 =  1.05792770261272E+16, relative difference =  3.78097670580074E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.76928632809007E-26, method 2 = -2.84217094304040E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.769E-26,  chi2_M: 1.058E+16
   max(B_n): 2.010E-14,  max(M): 3.643E+09,  min(M): 2.276E+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 =  1.03340705664987E+16, method 2 =  1.03340705664987E+16, relative difference = -6.77371027704497E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.83337841309620E-26, method 2 =  0.00000000000000E+00, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.833E-26,  chi2_M: 1.033E+16
   max(B_n): 2.298E-14,  max(M): 3.120E+09,  min(M): 2.786E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.03149509121196E+16, method 2 =  1.03149509121196E+16, relative difference =  3.87786624878668E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.08566974616047E-25, method 2 =  5.68434188608080E-14, relative difference = -9.99999999989143E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.086E-25,  chi2_M: 1.031E+16
   max(B_n): 5.740E-14,  max(M): 2.996E+09,  min(M): 2.898E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.360E-27 ( 13 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.03142251517194E+16, method 2 =  1.03142251517194E+16, relative difference = -9.69534778706380E-17
   2 methods of computing chi2_B that should agree: method 1 =  6.91406654520130E-24, method 2 =  0.00000000000000E+00, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 6.914E-24,  chi2_M: 1.031E+16
   max(B_n): 2.294E-13,  max(M): 2.980E+09,  min(M): 2.914E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.624E-26 ( 14 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.03141961617056E+16, method 2 =  1.03141961617056E+16, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  1.61681326198528E-22, method 2 =  0.00000000000000E+00, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.617E-22,  chi2_M: 1.031E+16
   max(B_n): 1.065E-12,  max(M): 2.983E+09,  min(M): 2.917E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.848E-26 ( 15 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.03141951226497E+16, method 2 =  1.03141951226497E+16, relative difference =  3.87815040576079E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.77883765559152E-21, method 2 =  5.68434188608080E-14, relative difference = -9.99999867043979E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.779E-21,  chi2_M: 1.031E+16
   max(B_n): 5.100E-12,  max(M): 2.983E+09,  min(M): 2.918E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.793E-25 ( 16 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.03141950534086E+16, method 2 =  1.03141950534086E+16, relative difference = -9.69537607948890E-17
   2 methods of computing chi2_B that should agree: method 1 =  8.82718874233727E-20, method 2 =  0.00000000000000E+00, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 8.827E-20,  chi2_M: 1.031E+16
   max(B_n): 2.460E-11,  max(M): 2.983E+09,  min(M): 2.918E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.833E-24 ( 17 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.03141950486662E+16, method 2 =  1.03141950486662E+16, relative difference =  2.90861282518402E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.06181523180556E-18, method 2 =  0.00000000000000E+00, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.062E-18,  chi2_M: 1.031E+16
   max(B_n): 1.189E-10,  max(M): 2.983E+09,  min(M): 2.918E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.859E-24 ( 18 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.03141950399001E+16, method 2 =  1.03141950399001E+16, relative difference =  9.69537609218688E-17
   2 methods of computing chi2_B that should agree: method 1 =  4.81576732310062E-17, method 2 = -2.84217094304040E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.816E-17,  chi2_M: 1.031E+16
   max(B_n): 5.743E-10,  max(M): 2.983E+09,  min(M): 2.918E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.281E-23 ( 19 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.03141949982260E+16, method 2 =  1.03141949982260E+16, relative difference =  5.81722567881642E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.12472532732956E-15, method 2 =  2.84217094304040E-14, relative difference = -9.23867272884071E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.125E-15,  chi2_M: 1.031E+16
   max(B_n): 2.775E-09,  max(M): 2.983E+09,  min(M): 2.918E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.069E-22 ( 20 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.03141947969184E+16, method 2 =  1.03141947969184E+16, relative difference =  9.69537632059043E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.62588355531420E-14, method 2 =  5.68434188608080E-14, relative difference = -3.68035542758174E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.626E-14,  chi2_M: 1.031E+16
   max(B_n): 1.340E-08,  max(M): 2.983E+09,  min(M): 2.918E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-21 ( 21 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.03141938257683E+16, method 2 =  1.03141938257682E+16, relative difference =  1.16344526801698E-15
   2 methods of computing chi2_B that should agree: method 1 =  6.12121993440930E-13, method 2 =  6.53699316899292E-13, relative difference = -3.28461237922969E-02
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.121E-13,  chi2_M: 1.031E+16
   max(B_n): 6.449E-08,  max(M): 2.983E+09,  min(M): 2.918E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.833E-21 ( 22 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.03141891620646E+16, method 2 =  1.03141891620646E+16, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  1.42005999470332E-11, method 2 =  1.42392764246324E-11, relative difference = -1.35993831666999E-03
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.420E-11,  chi2_M: 1.031E+16
   max(B_n): 3.066E-07,  max(M): 2.983E+09,  min(M): 2.918E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.336E-20 ( 23 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.03141668617897E+16, method 2 =  1.03141668617897E+16, relative difference = -9.69540257977254E-17
   2 methods of computing chi2_B that should agree: method 1 =  3.28223258290914E-10, method 2 =  3.28213900502305E-10, relative difference =  1.42554218379577E-05
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.282E-10,  chi2_M: 1.031E+16
   max(B_n): 1.428E-06,  max(M): 2.983E+09,  min(M): 2.918E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.129E-19 ( 24 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.03140595800538E+16, method 2 =  1.03140595800538E+16, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  7.63530220876363E-09, method 2 =  7.63535012993088E-09, relative difference = -3.13812181627550E-06
   Diagnostics:    0.00000000      sec.
   chi2_B: 7.635E-09,  chi2_M: 1.031E+16
   max(B_n): 6.689E-06,  max(M): 2.983E+09,  min(M): 2.918E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.456E-19 ( 25 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.03135413193558E+16, method 2 =  1.03135413193558E+16, relative difference =  1.26047878196817E-15
   2 methods of computing chi2_B that should agree: method 1 =  1.78255936077682E-07, method 2 =  1.78255987748344E-07, relative difference = -1.44933894705646E-07
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.783E-07,  chi2_M: 1.031E+16
   max(B_n): 3.197E-05,  max(M): 2.983E+09,  min(M): 2.918E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.637E-18 ( 26 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.03110371735369E+16, method 2 =  1.03110371735369E+16, relative difference =  8.72851086513230E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.16245605903474E-06, method 2 =  4.16245609358157E-06, relative difference = -4.14981305536944E-09
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 4.162E-06,  chi2_M: 1.031E+16
   max(B_n): 1.547E-04,  max(M): 2.983E+09,  min(M): 2.917E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.274E-17 ( 27 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.02989476721902E+16, method 2 =  1.02989476721902E+16, relative difference = -9.70972988531882E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.71091624838568E-05, method 2 =  9.71091625103782E-05, relative difference = -1.36554437246677E-10
   Diagnostics:    3.00000003E-03  sec.
   chi2_B: 9.711E-05,  chi2_M: 1.030E+16
   max(B_n): 7.474E-04,  max(M): 2.981E+09,  min(M): 2.916E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.158E-17 ( 28 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.02408188950344E+16, method 2 =  1.02408188950344E+16, relative difference =  7.81187528262906E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.25538370972748E-03, method 2 =  2.25538370986555E-03, relative difference = -3.06104490143361E-11
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 2.255E-03,  chi2_M: 1.024E+16
   max(B_n): 3.602E-03,  max(M): 2.973E+09,  min(M): 2.907E+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 =  9.96671172644732E+15, method 2 =  9.96671172644733E+15, relative difference = -1.00333994545707E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.12681454707987E-02, method 2 =  5.12681454708854E-02, relative difference = -8.45907073507956E-13
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.127E-02,  chi2_M: 9.967E+15
   max(B_n): 1.716E-02,  max(M): 2.934E+09,  min(M): 2.867E+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 =  8.78536709792714E+15, method 2 =  8.78536709792714E+15, relative difference = -1.70738454441354E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.05546970281006E+00, method 2 =  1.05546970281009E+00, relative difference = -1.39899479712225E-14
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.055E+00,  chi2_M: 8.785E+15
   max(B_n): 7.760E-02,  max(M): 2.758E+09,  min(M): 2.689E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.23518099285832E+15, method 2 =  5.23518099285833E+15, relative difference = -3.82030726870444E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.46894807925295E+01, method 2 =  1.46894807925296E+01, relative difference = -9.67417086819641E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.469E+01,  chi2_M: 5.235E+15
   max(B_n): 2.881E-01,  max(M): 2.132E+09,  min(M): 2.072E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.18468406335948E+15, method 2 =  1.18468406335948E+15, relative difference =  1.05513363322817E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.76403146990880E+01, method 2 =  7.76403146990879E+01, relative difference =  1.09820688674016E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.764E+01,  chi2_M: 1.185E+15
   max(B_n): 6.613E-01,  max(M): 1.015E+09,  min(M): 9.856E+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 =  9.48636304382479E+13, method 2 =  9.48636304382479E+13, relative difference =  1.64710120494189E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.45211327387948E+02, method 2 =  1.45211327387948E+02, relative difference =  9.78632657026554E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.452E+02,  chi2_M: 9.486E+13
   max(B_n): 9.045E-01,  max(M): 2.868E+08,  min(M): 2.792E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.75767974867645E+12, method 2 =  4.75767974867644E+12, relative difference =  5.13150606801390E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.70103649372415E+02, method 2 =  1.70103649372414E+02, relative difference =  5.01253962544553E-16
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 1.701E+02,  chi2_M: 4.758E+12
   max(B_n): 9.787E-01,  max(M): 6.416E+07,  min(M): 6.258E+07
 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 =  2.10813330516475E+11, method 2 =  2.10813330516476E+11, relative difference = -6.51425131541990E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.76050148366596E+02, method 2 =  1.76050148366596E+02, relative difference = -1.61440985390256E-16
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 1.761E+02,  chi2_M: 2.108E+11
   max(B_n): 9.955E-01,  max(M): 1.350E+07,  min(M): 1.318E+07
 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 =  9.09069935934521E+09, method 2 =  9.09069935934521E+09, relative difference =  1.04906595049354E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.77319461087791E+02, method 2 =  1.77319461087791E+02, relative difference =  3.20570672345236E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.773E+02,  chi2_M: 9.091E+09
   max(B_n): 9.991E-01,  max(M): 2.804E+06,  min(M): 2.736E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.859E-11 ( 37 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.89783342184137E+08, method 2 =  3.89783342184136E+08, relative difference =  6.11669492507281E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.77583813281368E+02, method 2 =  1.77583813281368E+02, relative difference =  8.00233672912859E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.776E+02,  chi2_M: 3.898E+08
   max(B_n): 9.998E-01,  max(M): 5.805E+05,  min(M): 5.666E+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.66930680082618E+07, method 2 =  1.66930680082618E+07, relative difference = -5.57909770783024E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.77638585217467E+02, method 2 =  1.77638585217467E+02, relative difference =  2.39996080206420E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.776E+02,  chi2_M: 1.669E+07
   max(B_n): 1.000E+00,  max(M): 1.201E+05,  min(M): 1.173E+05
 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 =  7.14731341659699E+05, method 2 =  7.14731341659699E+05, relative difference =  1.62879833360327E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.77649921451622E+02, method 2 =  1.77649921451623E+02, relative difference = -4.79961531052134E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.776E+02,  chi2_M: 7.147E+05
   max(B_n): 1.000E+00,  max(M): 2.486E+04,  min(M): 2.426E+04
 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 =  3.06004285852313E+04, method 2 =  3.06004285852313E+04, relative difference =  1.78329796768640E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.77652267210509E+02, method 2 =  1.77652267210510E+02, relative difference = -3.19970129024310E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.777E+02,  chi2_M: 3.060E+04
   max(B_n): 1.000E+00,  max(M): 5.144E+03,  min(M): 5.021E+03
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   1.7389999628067017      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
