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: 7.106E+02 m^2. Volume: 2.130E+02 m^3.
 G =    0.0000000000000000       ; curpol =    1.0000000000000000     
 Done initializing plasma surface. Took    4.10000011E-02  sec.
 Initializing coil surface.
   Building a plain circular torus.
   Evaluating coil surface & derivatives:   1.70000009E-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.19999999E-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.80299997     sec.
 Computing g and RHS_B:
  js=   1 j_RZetaZ= 1 dgemm took  2.000E-03 sec, rest took  1.000E-03 sec.
  js=   1 j_RZetaZ= 2 dgemm took  1.000E-03 sec, rest took  0.000E+00 sec.
  js=   1 j_RZetaZ= 3 dgemm took  1.000E-03 sec, rest took  0.000E+00 sec.
 Total time to compute g and RHS_B:   4.99999989E-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  1.000E-03 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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.97558979815603E+26, method 2 =  6.97558979815605E+26, relative difference = -8.86627953363156E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.41298008219861E-17, method 2 = -4.32009983342141E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 2.413E-17,  chi2_M: 6.976E+26
   max(B_n): 3.672E-10,  max(M): 2.033E+15,  min(M): 1.576E+12
 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.35892442932071E+21, method 2 =  1.35892442932071E+21, relative difference = -4.82263756438316E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.99171773896644E-18, method 2 = -3.12638803734444E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.992E-18,  chi2_M: 1.359E+21
   max(B_n): 1.395E-10,  max(M): 1.638E+12,  min(M): 1.626E+11
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.833E-34 (  3 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.65324060564061E+19, method 2 =  3.65324060564062E+19, relative difference = -1.06513652399242E-15
   2 methods of computing chi2_B that should agree: method 1 =  7.59224398509062E-19, method 2 = -1.76214598468505E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    4.99999989E-03  sec.
   chi2_B: 7.592E-19,  chi2_M: 3.653E+19
   max(B_n): 4.812E-11,  max(M): 3.153E+11,  min(M): 1.788E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.08036999229622E+22, method 2 =  1.08036999229622E+22, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  1.71662728195955E-17, method 2 = -1.26192389870994E-11, relative difference =  1.00000000000000E+00
   Diagnostics:    3.00000003E-03  sec.
   chi2_B: 1.717E-17,  chi2_M: 1.080E+22
   max(B_n): 2.335E-10,  max(M): 4.365E+12,  min(M): 1.467E+12
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.129E-32 (  5 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.30348607632749E+21, method 2 =  3.30348607632750E+21, relative difference = -1.11095246512435E-15
   2 methods of computing chi2_B that should agree: method 1 =  4.73780506587477E-16, method 2 =  5.45242073712870E-10, relative difference = -9.99998262129363E-01
   Diagnostics:    4.00000019E-03  sec.
   chi2_B: 4.738E-16,  chi2_M: 3.303E+21
   max(B_n): 1.222E-09,  max(M): 2.396E+12,  min(M): 1.931E+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:    2.00000009E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.67607273988729E+18, method 2 =  2.67607273988729E+18, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  3.12563472413802E-20, method 2 = -9.66338120633736E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.126E-20,  chi2_M: 2.676E+18
   max(B_n): 9.898E-12,  max(M): 7.945E+10,  min(M): 5.154E+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 =  7.17492315659724E+16, method 2 =  7.17492315659724E+16, relative difference = -1.11499451985686E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.66999296902465E-23, method 2 = -8.52651282912120E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    7.00000022E-03  sec.
   chi2_B: 3.670E-23,  chi2_M: 7.175E+16
   max(B_n): 3.575E-13,  max(M): 1.458E+10,  min(M): 8.253E+08
 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.88281677472549E+16, method 2 =  1.88281677472548E+16, relative difference =  1.38090972786191E-15
   2 methods of computing chi2_B that should agree: method 1 =  1.78110106480159E-23, method 2 = -6.82121026329696E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    4.99999989E-03  sec.
   chi2_B: 1.781E-23,  chi2_M: 1.883E+16
   max(B_n): 2.554E-13,  max(M): 7.262E+09,  min(M): 7.745E+08
 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.78087552276553E+19, method 2 =  2.78087552276553E+19, relative difference = -2.94583483976054E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.95159596013288E-21, method 2 =  2.59262833424145E-10, relative difference = -9.99999999938660E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.952E-21,  chi2_M: 2.781E+19
   max(B_n): 6.302E-12,  max(M): 2.487E+11,  min(M): 5.991E+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 =  7.98219563261637E+15, method 2 =  7.98219563261637E+15, relative difference =  1.87918220629771E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.72281873293432E-26, method 2 = -1.81898940354586E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 9.723E-26,  chi2_M: 7.982E+15
   max(B_n): 2.681E-14,  max(M): 4.233E+09,  min(M): 6.995E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.438E-28 ( 11 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.48261593003492E+15, method 2 =  3.48261593003492E+15, relative difference = -3.58925596480421E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.69147940617505E-26, method 2 = -1.19371179607697E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.691E-26,  chi2_M: 3.483E+15
   max(B_n): 2.010E-14,  max(M): 2.918E+09,  min(M): 4.913E+07
 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 =  6.75474820362477E+14, method 2 =  6.75474820362476E+14, relative difference =  9.25275052687544E-17
   2 methods of computing chi2_B that should agree: method 1 =  8.26024787644500E-26, method 2 = -6.25277607468888E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.260E-26,  chi2_M: 6.755E+14
   max(B_n): 2.642E-14,  max(M): 9.425E+08,  min(M): 5.384E+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:    2.00000009E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.46563784015592E+14, method 2 =  6.46563784015592E+14, relative difference = -1.93329727229179E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.02675892511418E-24, method 2 = -6.25277607468888E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.027E-24,  chi2_M: 6.466E+14
   max(B_n): 5.662E-14,  max(M): 7.806E+08,  min(M): 6.943E+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 =  6.45262480189674E+14, method 2 =  6.45262480189672E+14, relative difference =  1.25917750519319E-15
   2 methods of computing chi2_B that should agree: method 1 =  2.28505560053441E-23, method 2 = -6.25277607468888E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.285E-23,  chi2_M: 6.453E+14
   max(B_n): 2.180E-13,  max(M): 7.463E+08,  min(M): 7.258E+08
 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 =  6.45220986275077E+14, method 2 =  6.45220986275078E+14, relative difference = -7.74928296871662E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.30892621786284E-22, method 2 = -6.82121026329696E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.309E-22,  chi2_M: 6.452E+14
   max(B_n): 1.006E-12,  max(M): 7.420E+08,  min(M): 7.310E+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 =  6.45217623636654E+14, method 2 =  6.45217623636653E+14, relative difference =  9.68665419393381E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.23846191501432E-20, method 2 = -6.25277607468888E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 1.238E-20,  chi2_M: 6.452E+14
   max(B_n): 4.841E-12,  max(M): 7.429E+08,  min(M): 7.325E+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 =  6.45217453724617E+14, method 2 =  6.45217453724616E+14, relative difference =  2.90599702344733E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.88896320751782E-19, method 2 = -6.82121026329696E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.889E-19,  chi2_M: 6.452E+14
   max(B_n): 2.337E-11,  max(M): 7.431E+08,  min(M): 7.329E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.859E-24 ( 18 of  40)
   Additions:    4.99999989E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.45217449345017E+14, method 2 =  6.45217449345016E+14, relative difference =  5.81199408634525E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.71204757268092E-18, method 2 = -5.11590769747272E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.712E-18,  chi2_M: 6.452E+14
   max(B_n): 1.126E-10,  max(M): 7.431E+08,  min(M): 7.329E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.45217443513625E+14, method 2 =  6.45217443513625E+14, relative difference =  1.93733137962443E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.52858956865634E-16, method 2 = -5.11590769747272E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.529E-16,  chi2_M: 6.452E+14
   max(B_n): 5.384E-10,  max(M): 7.431E+08,  min(M): 7.330E+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 =  6.45217418833994E+14, method 2 =  6.45217418833996E+14, relative difference = -1.35613201760928E-15
   2 methods of computing chi2_B that should agree: method 1 =  3.17388118170861E-15, method 2 = -6.25277607468888E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.174E-15,  chi2_M: 6.452E+14
   max(B_n): 2.473E-09,  max(M): 7.431E+08,  min(M): 7.330E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.45217342363434E+14, method 2 =  6.45217342363434E+14, relative difference =  2.90599752500741E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.56836068313845E-14, method 2 = -6.82121026329696E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.568E-14,  chi2_M: 6.452E+14
   max(B_n): 9.734E-09,  max(M): 7.431E+08,  min(M): 7.330E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.45217244204905E+14, method 2 =  6.45217244204904E+14, relative difference =  1.93733197806944E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.72828588208602E-13, method 2 = -3.41060513164848E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 2.728E-13,  chi2_M: 6.452E+14
   max(B_n): 2.738E-08,  max(M): 7.431E+08,  min(M): 7.330E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.336E-20 ( 23 of  40)
   Additions:    4.00000019E-03  sec.
   DSYSV:    8.00000038E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.45217175503805E+14, method 2 =  6.45217175503806E+14, relative difference = -6.78066264523083E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.11139275385732E-12, method 2 =  3.97903932025656E-13, relative difference =  4.72729337117873E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.111E-12,  chi2_M: 6.452E+14
   max(B_n): 6.828E-08,  max(M): 7.431E+08,  min(M): 7.330E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.129E-19 ( 24 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.45216963552170E+14, method 2 =  6.45216963552169E+14, relative difference =  9.68666410379437E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.54886065116648E-11, method 2 =  1.48361323226709E-11, relative difference =  2.15162344038110E-02
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 1.549E-11,  chi2_M: 6.452E+14
   max(B_n): 2.335E-07,  max(M): 7.430E+08,  min(M): 7.330E+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 =  6.45215949511092E+14, method 2 =  6.45215949511093E+14, relative difference = -1.93733586553026E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.49317304570186E-10, method 2 =  3.48734374711057E-10, relative difference =  8.35081236004375E-04
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 3.493E-10,  chi2_M: 6.452E+14
   max(B_n): 1.021E-06,  max(M): 7.430E+08,  min(M): 7.330E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.45211049280427E+14, method 2 =  6.45211049280426E+14, relative difference =  7.74940231661603E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.14608666568002E-09, method 2 =  8.14549139249721E-09, relative difference =  3.65387061147950E-05
   Diagnostics:    7.00000022E-03  sec.
   chi2_B: 8.146E-09,  chi2_M: 6.452E+14
   max(B_n): 4.827E-06,  max(M): 7.430E+08,  min(M): 7.330E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.45187367617676E+14, method 2 =  6.45187367617675E+14, relative difference =  5.81226506936536E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.90249878890077E-07, method 2 =  1.90249238585238E-07, relative difference =  1.68280242839969E-06
   Diagnostics:    2.89999992E-02  sec.
   chi2_B: 1.902E-07,  chi2_M: 6.452E+14
   max(B_n): 2.322E-05,  max(M): 7.430E+08,  min(M): 7.330E+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:    4.99999989E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.45072934170359E+14, method 2 =  6.45072934170360E+14, relative difference = -8.71994421411343E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.44290603748107E-06, method 2 =  4.44290537870984E-06, relative difference =  7.41374315729513E-08
   Diagnostics:    3.00000003E-03  sec.
   chi2_B: 4.443E-06,  chi2_M: 6.451E+14
   max(B_n): 1.123E-04,  max(M): 7.430E+08,  min(M): 7.330E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.44520314201080E+14, method 2 =  6.44520314201080E+14, relative difference = -9.69713422880586E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.03684993982779E-04, method 2 =  1.03684993291608E-04, relative difference =  3.33302894873600E-09
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 1.037E-04,  chi2_M: 6.445E+14
   max(B_n): 5.425E-04,  max(M): 7.426E+08,  min(M): 7.326E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.41859523841757E+14, method 2 =  6.41859523841757E+14, relative difference = -3.89493324806745E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.41179444492500E-03, method 2 =  2.41179444441286E-03, relative difference =  1.06173927296870E-10
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.412E-03,  chi2_M: 6.419E+14
   max(B_n): 2.617E-03,  max(M): 7.411E+08,  min(M): 7.311E+08
 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 =  6.29229462926928E+14, method 2 =  6.29229462926928E+14, relative difference =  9.93278345697205E-17
   2 methods of computing chi2_B that should agree: method 1 =  5.52242956631225E-02, method 2 =  5.52242956623559E-02, relative difference =  6.94086714283724E-12
   Diagnostics:    9.99999978E-03  sec.
   chi2_B: 5.522E-02,  chi2_M: 6.292E+14
   max(B_n): 1.252E-02,  max(M): 7.338E+08,  min(M): 7.239E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.73144536837741E+14, method 2 =  5.73144536837742E+14, relative difference = -2.18095073695848E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.17491215149630E+00, method 2 =  1.17491215149568E+00, relative difference =  2.66095472190298E-13
   Diagnostics:    8.00000038E-03  sec.
   chi2_B: 1.175E+00,  chi2_M: 5.731E+14
   max(B_n): 5.773E-02,  max(M): 7.005E+08,  min(M): 6.907E+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:    3.00000003E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.84815086255338E+14, method 2 =  3.84815086255338E+14, relative difference =  6.49662679373533E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.84251633592601E+01, method 2 =  1.84251633592594E+01, relative difference =  1.65823971499984E-14
   Diagnostics:    8.99999961E-03  sec.
   chi2_B: 1.843E+01,  chi2_M: 3.848E+14
   max(B_n): 2.284E-01,  max(M): 5.744E+08,  min(M): 5.656E+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.09711422548928E+14, method 2 =  1.09711422548928E+14, relative difference = -2.84838162462649E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.22695248809668E+02, method 2 =  1.22695248809667E+02, relative difference =  2.20062506257018E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.227E+02,  chi2_M: 1.097E+14
   max(B_n): 5.885E-01,  max(M): 3.071E+08,  min(M): 3.016E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.03716630866194E+13, method 2 =  1.03716630866194E+13, relative difference = -2.82470369075102E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.70921227613896E+02, method 2 =  2.70921227613896E+02, relative difference =  8.39261203139365E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.709E+02,  chi2_M: 1.037E+13
   max(B_n): 8.736E-01,  max(M): 9.449E+07,  min(M): 9.267E+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 =  5.48878866875414E+11, method 2 =  5.48878866875414E+11, relative difference =  2.22399366903860E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.34881811181137E+02, method 2 =  3.34881811181136E+02, relative difference =  1.01845039586330E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.349E+02,  chi2_M: 5.489E+11
   max(B_n): 9.709E-01,  max(M): 2.174E+07,  min(M): 2.131E+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 =  2.46254163078581E+10, method 2 =  2.46254163078581E+10, relative difference = -1.54908945210713E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.50929229870112E+02, method 2 =  3.50929229870111E+02, relative difference =  1.21484790997277E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.509E+02,  chi2_M: 2.463E+10
   max(B_n): 9.938E-01,  max(M): 4.606E+06,  min(M): 4.514E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.06470370442094E+09, method 2 =  1.06470370442094E+09, relative difference = -5.03841398081985E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.54393502028999E+02, method 2 =  3.54393502028998E+02, relative difference =  2.40594502447269E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.544E+02,  chi2_M: 1.065E+09
   max(B_n): 9.987E-01,  max(M): 9.577E+05,  min(M): 9.386E+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 =  4.56765273379626E+07, method 2 =  4.56765273379626E+07, relative difference =  4.07790446819471E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.55116711201543E+02, method 2 =  3.55116711201543E+02, relative difference =  3.20139363019428E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.551E+02,  chi2_M: 4.568E+07
   max(B_n): 9.997E-01,  max(M): 1.984E+05,  min(M): 1.944E+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 =  1.95638954278332E+06, method 2 =  1.95638954278332E+06, relative difference = -2.38020740309852E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.55266629571005E+02, method 2 =  3.55266629571005E+02, relative difference =  8.80011736852304E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.553E+02,  chi2_M: 1.956E+06
   max(B_n): 9.999E-01,  max(M): 4.105E+04,  min(M): 4.023E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   4.0809998512268066      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
