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: 2.961E+02 m^2. Volume: 3.698E+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.30000003E-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.60000008E-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.66900003     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  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  1.000E-03 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 2 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 3 dgemm took  0.000E+00 sec.
 Total time for matrix_B:   1.00000005E-03 sec.
 Regularization term:   0.00000000     sec.
 Optimal LWORK:        9408
 Solving system for lambda= 0.000E+00 (  1 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.73062226982330E+30, method 2 =  3.73062226982330E+30, relative difference = -3.77249365323695E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.46359508026209E-16, method 2 =  1.13686837721616E-13, relative difference = -9.93925286382549E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.464E-16,  chi2_M: 3.731E+30
   max(B_n): 3.481E-09,  max(M): 1.548E+17,  min(M): 5.720E+15
 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 =  3.07242395976696E+18, method 2 =  3.07242395976695E+18, relative difference =  2.49965502826717E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.12093046313393E-24, method 2 =  5.68434188608080E-14, relative difference = -9.99999999960561E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.121E-24,  chi2_M: 3.072E+18
   max(B_n): 1.594E-13,  max(M): 9.841E+10,  min(M): 4.745E+09
 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.75027098906632E+18, method 2 =  1.75027098906632E+18, relative difference = -5.11920728616985E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.85635646303186E-25, method 2 =  1.13686837721616E-13, relative difference = -9.99999999993216E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.856E-25,  chi2_M: 1.750E+18
   max(B_n): 1.109E-13,  max(M): 6.082E+10,  min(M): 7.156E+09
 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 =  2.01945260682685E+19, method 2 =  2.01945260682685E+19, relative difference = -2.02827240716285E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.89382114625454E-24, method 2 = -8.52651282912120E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.894E-24,  chi2_M: 2.019E+19
   max(B_n): 4.197E-13,  max(M): 1.872E+11,  min(M): 7.947E+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 =  1.26095618522329E+18, method 2 =  1.26095618522329E+18, relative difference =  6.09061606580724E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.53332806520096E-26, method 2 =  1.98951966012828E-13, relative difference = -9.99999999999444E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.533E-26,  chi2_M: 1.261E+18
   max(B_n): 4.685E-14,  max(M): 5.816E+10,  min(M): 1.490E+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 =  9.59963651692479E+16, method 2 =  9.59963651692478E+16, relative difference =  5.83355420814823E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.25565364858787E-26, method 2 =  1.98951966012828E-13, relative difference = -9.99999999999773E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.256E-26,  chi2_M: 9.600E+16
   max(B_n): 2.487E-14,  max(M): 1.624E+10,  min(M): 1.703E+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 =  4.90629849705632E+16, method 2 =  4.90629849705632E+16, relative difference =  8.15278565378751E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.72649315490821E-26, method 2 =  1.70530256582424E-13, relative difference = -9.99999999999680E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.726E-26,  chi2_M: 4.906E+16
   max(B_n): 2.948E-14,  max(M): 1.084E+10,  min(M): 1.832E+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 =  9.49153125547023E+16, method 2 =  9.49153125547023E+16, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  3.64409432389069E-26, method 2 =  1.42108547152020E-13, relative difference = -9.99999999999487E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.644E-26,  chi2_M: 9.492E+16
   max(B_n): 3.364E-14,  max(M): 1.714E+10,  min(M): 1.158E+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 =  2.15493586842861E+16, method 2 =  2.15493586842861E+16, relative difference = -1.85620373144414E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.30582017829666E-26, method 2 =  5.68434188608080E-14, relative difference = -9.99999999999189E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.306E-26,  chi2_M: 2.155E+16
   max(B_n): 2.731E-14,  max(M): 4.690E+09,  min(M): 3.808E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.14108922793194E+16, method 2 =  2.14108922793194E+16, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  2.53634947609148E-26, method 2 =  1.13686837721616E-13, relative difference = -9.99999999999554E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.536E-26,  chi2_M: 2.141E+16
   max(B_n): 2.826E-14,  max(M): 4.415E+09,  min(M): 4.098E+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 =  2.13986548111527E+16, method 2 =  2.13986548111527E+16, relative difference = -6.54246732963017E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.88547079945143E-26, method 2 =  5.68434188608080E-14, relative difference = -9.99999999996874E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 8.885E-26,  chi2_M: 2.140E+16
   max(B_n): 4.485E-14,  max(M): 4.315E+09,  min(M): 4.168E+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 =  2.13966853311134E+16, method 2 =  2.13966853311134E+16, relative difference = -2.80417265905914E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.61629475082403E-24, method 2 =  2.84217094304040E-14, relative difference = -9.99999999886263E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.616E-24,  chi2_M: 2.140E+16
   max(B_n): 1.384E-13,  max(M): 4.298E+09,  min(M): 4.204E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.13966447423927E+16, method 2 =  2.13966447423927E+16, relative difference =  1.86945198565403E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.73820836304146E-23, method 2 =  1.13686837721616E-13, relative difference = -9.99999999342367E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.738E-23,  chi2_M: 2.140E+16
   max(B_n): 5.895E-13,  max(M): 4.298E+09,  min(M): 4.202E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.13966424327170E+16, method 2 =  2.13966424327170E+16, relative difference =  7.47780874981341E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.72948988872386E-22, method 2 =  8.52651282912120E-14, relative difference = -9.99999979523892E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 8.729E-22,  chi2_M: 2.140E+16
   max(B_n): 2.812E-12,  max(M): 4.299E+09,  min(M): 4.202E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.13966423504804E+16, method 2 =  2.13966423504804E+16, relative difference =  4.67363048659617E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.03900516127587E-20, method 2 =  1.70530256582424E-13, relative difference = -9.99999760862974E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.039E-20,  chi2_M: 2.140E+16
   max(B_n): 1.356E-11,  max(M): 4.299E+09,  min(M): 4.202E+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 =  2.13966423449544E+16, method 2 =  2.13966423449544E+16, relative difference =  2.80417829268193E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.76246974450360E-19, method 2 =  1.13686837721616E-13, relative difference = -9.99991621809720E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 4.762E-19,  chi2_M: 2.140E+16
   max(B_n): 6.552E-11,  max(M): 4.299E+09,  min(M): 4.202E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.833E-24 ( 17 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.13966423351931E+16, method 2 =  2.13966423351931E+16, relative difference =  4.67363048993535E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.11231330152826E-17, method 2 =  1.13686837721616E-13, relative difference = -9.99804338918101E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.112E-17,  chi2_M: 2.140E+16
   max(B_n): 3.166E-10,  max(M): 4.299E+09,  min(M): 4.202E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.13966422886806E+16, method 2 =  2.13966422886807E+16, relative difference = -4.67363050009498E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.59727515194377E-16, method 2 =  1.70530256582424E-13, relative difference = -9.96958515845450E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.597E-16,  chi2_M: 2.140E+16
   max(B_n): 1.530E-09,  max(M): 4.299E+09,  min(M): 4.202E+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 =  2.13966420642072E+16, method 2 =  2.13966420642073E+16, relative difference = -1.86945221965052E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.05786686335467E-15, method 2 =  8.52651282912120E-14, relative difference = -8.67330964055623E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 6.058E-15,  chi2_M: 2.140E+16
   max(B_n): 7.380E-09,  max(M): 4.299E+09,  min(M): 4.202E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.13966409855691E+16, method 2 =  2.13966409855691E+16, relative difference =  1.86945231389254E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.40606211406898E-13, method 2 =  2.84217094304040E-13, relative difference = -3.38048503852232E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.406E-13,  chi2_M: 2.140E+16
   max(B_n): 3.539E-08,  max(M): 4.299E+09,  min(M): 4.202E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.13966358732738E+16, method 2 =  2.13966358732739E+16, relative difference = -5.60835828168155E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.21681149948803E-12, method 2 =  3.24007487506606E-12, relative difference = -3.60287826493390E-03
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.217E-12,  chi2_M: 2.140E+16
   max(B_n): 1.661E-07,  max(M): 4.299E+09,  min(M): 4.202E+09
 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 =  2.13966118827919E+16, method 2 =  2.13966118827919E+16, relative difference = -8.41254685489546E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.30134867354264E-11, method 2 =  7.31574800738599E-11, relative difference = -9.85102182579370E-04
   Diagnostics:    0.00000000      sec.
   chi2_B: 7.301E-11,  chi2_M: 2.140E+16
   max(B_n): 7.546E-07,  max(M): 4.299E+09,  min(M): 4.202E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.13964971895767E+16, method 2 =  2.13964971895767E+16, relative difference =  2.80419731642940E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.68916706766791E-09, method 2 =  1.68927272170549E-09, relative difference = -3.12730266542633E-05
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.689E-09,  chi2_M: 2.140E+16
   max(B_n): 3.475E-06,  max(M): 4.298E+09,  min(M): 4.202E+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 =  2.13959433286519E+16, method 2 =  2.13959433286518E+16, relative difference =  1.40213495330333E-15
   2 methods of computing chi2_B that should agree: method 1 =  3.94180531174493E-08, method 2 =  3.94180972307367E-08, relative difference = -5.59556589543324E-07
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.942E-08,  chi2_M: 2.140E+16
   max(B_n): 1.652E-05,  max(M): 4.298E+09,  min(M): 4.202E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.13932669016017E+16, method 2 =  2.13932669016016E+16, relative difference =  9.34873579242947E-17
   2 methods of computing chi2_B that should agree: method 1 =  9.20537522102283E-07, method 2 =  9.20537672755017E-07, relative difference = -8.18286697396910E-08
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.205E-07,  chi2_M: 2.139E+16
   max(B_n): 7.993E-05,  max(M): 4.298E+09,  min(M): 4.201E+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 =  2.13803390122701E+16, method 2 =  2.13803390122701E+16, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  2.14881214580771E-05, method 2 =  2.14881215754303E-05, relative difference = -2.73065282634218E-09
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.149E-05,  chi2_M: 2.138E+16
   max(B_n): 3.865E-04,  max(M): 4.297E+09,  min(M): 4.200E+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 =  2.13180247033716E+16, method 2 =  2.13180247033716E+16, relative difference = -6.56721257940271E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.00436151540485E-04, method 2 =  5.00436151668282E-04, relative difference = -1.27686131134177E-10
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.004E-04,  chi2_M: 2.132E+16
   max(B_n): 1.865E-03,  max(M): 4.291E+09,  min(M): 4.194E+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 =  2.10206784787952E+16, method 2 =  2.10206784787952E+16, relative difference =  5.70866445253188E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.15253934955142E-02, method 2 =  1.15253934955604E-02, relative difference = -2.00257768698382E-12
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.153E-02,  chi2_M: 2.102E+16
   max(B_n): 8.943E-03,  max(M): 4.262E+09,  min(M): 4.163E+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 =  1.96678355389334E+16, method 2 =  1.96678355389334E+16, relative difference =  5.08444357296182E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.51851162950640E-01, method 2 =  2.51851162950743E-01, relative difference = -2.04542824972752E-13
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.519E-01,  chi2_M: 1.967E+16
   max(B_n): 4.167E-02,  max(M): 4.127E+09,  min(M): 4.023E+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 =  1.46647959698443E+16, method 2 =  1.46647959698443E+16, relative difference = -1.36381031424690E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.38561449931726E+00, method 2 =  4.38561449931737E+00, relative difference = -1.26575537452959E-14
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.386E+00,  chi2_M: 1.466E+16
   max(B_n): 1.728E-01,  max(M): 3.570E+09,  min(M): 3.467E+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 =  5.32379255973918E+15, method 2 =  5.32379255973919E+15, relative difference = -6.57426066234907E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.71858062628770E+01, method 2 =  3.71858062628771E+01, relative difference = -1.43309263769300E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.719E+01,  chi2_M: 5.324E+15
   max(B_n): 5.013E-01,  max(M): 2.153E+09,  min(M): 2.087E+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 =  6.23973255630202E+14, method 2 =  6.23973255630202E+14, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  1.01797724606473E+02, method 2 =  1.01797724606473E+02, relative difference =  2.09398413915507E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.018E+02,  chi2_M: 6.240E+14
   max(B_n): 8.298E-01,  max(M): 7.365E+08,  min(M): 7.152E+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 =  3.57401710482472E+13, method 2 =  3.57401710482472E+13, relative difference = -1.09295783579961E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.36189931686153E+02, method 2 =  1.36189931686153E+02, relative difference = -7.30420977342543E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.362E+02,  chi2_M: 3.574E+13
   max(B_n): 9.594E-01,  max(M): 1.760E+08,  min(M): 1.714E+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 =  1.63440426286371E+12, method 2 =  1.63440426286371E+12, relative difference = -2.24063865850635E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.45467525315076E+02, method 2 =  1.45467525315076E+02, relative difference = -8.79218177114085E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.455E+02,  chi2_M: 1.634E+12
   max(B_n): 9.913E-01,  max(M): 3.761E+07,  min(M): 3.668E+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 =  7.09546742520179E+10, method 2 =  7.09546742520179E+10, relative difference =  1.07524903914741E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.47505356785868E+02, method 2 =  1.47505356785868E+02, relative difference = -8.67071509968859E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.475E+02,  chi2_M: 7.095E+10
   max(B_n): 9.982E-01,  max(M): 7.836E+06,  min(M): 7.644E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.04660066077979E+09, method 2 =  3.04660066077978E+09, relative difference =  6.26057972535301E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.47932357727613E+02, method 2 =  1.47932357727613E+02, relative difference = -9.60631935669437E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.479E+02,  chi2_M: 3.047E+09
   max(B_n): 9.996E-01,  max(M): 1.624E+06,  min(M): 1.584E+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 =  1.30513219324691E+08, method 2 =  1.30513219324691E+08, relative difference =  9.13388621992475E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.48020941647788E+02, method 2 =  1.48020941647788E+02, relative difference = -9.60057040375836E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.480E+02,  chi2_M: 1.305E+08
   max(B_n): 9.999E-01,  max(M): 3.360E+05,  min(M): 3.279E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.58839714207925E+06, method 2 =  5.58839714207925E+06, relative difference =  1.66652897232885E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.48039280820518E+02, method 2 =  1.48039280820518E+02, relative difference = -7.67950486462100E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.480E+02,  chi2_M: 5.588E+06
   max(B_n): 1.000E+00,  max(M): 6.954E+04,  min(M): 6.784E+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 =  2.39263991830210E+05, method 2 =  2.39263991830210E+05, relative difference =  4.86555962459855E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.48043075874211E+02, method 2 =  1.48043075874211E+02, relative difference = -6.71939450183654E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.480E+02,  chi2_M: 2.393E+05
   max(B_n): 1.000E+00,  max(M): 1.439E+04,  min(M): 1.404E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.02437411447240E+04, method 2 =  1.02437411447240E+04, relative difference =  3.55141618251984E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.48043861141446E+02, method 2 =  1.48043861141447E+02, relative difference = -4.79954204302481E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.480E+02,  chi2_M: 1.024E+04
   max(B_n): 1.000E+00,  max(M): 2.977E+03,  min(M): 2.905E+03
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   1.7519999742507935      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
