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.70000009E-02  sec.
 Initializing coil surface.
   Building a plain circular torus.
   Evaluating coil surface & derivatives:   8.00000038E-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    9.99999978E-03  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.57200003     sec.
 Computing g and RHS_B:
  js=   1 j_RZetaZ= 1 dgemm took  2.000E-03 sec, rest took  0.000E+00 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  0.000E+00 sec, rest took  0.000E+00 sec.
 Total time to compute g and RHS_B:   3.00000003E-03 sec.
 About to form matrix_B.
  js=   1, ks=   1, j_RZetaZ= 1, k_RZetaZ= 1 dgemm took  1.000E-03 sec.
  js=   1, ks=   1, j_RZetaZ= 2, k_RZetaZ= 1 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 2, k_RZetaZ= 2 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 1 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 2 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 3 dgemm took  0.000E+00 sec.
 Total time for matrix_B:   1.00000005E-03 sec.
 Regularization term:   1.00000005E-03 sec.
 Optimal LWORK:        9408
 Solving system for lambda= 0.000E+00 (  1 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  8.99848295275107E+31, method 2 =  8.99848295275108E+31, relative difference = -5.00484320636917E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.39617705644606E-14, method 2 =  4.12114786740858E-12, relative difference = -9.93247206586133E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.396E-14,  chi2_M: 8.998E+31
   max(B_n): 1.801E-08,  max(M): 6.493E+17,  min(M): 6.212E+16
 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.02596019151622E+22, method 2 =  1.02596019151622E+22, relative difference =  3.06613066083106E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.18513892477015E-21, method 2 = -6.56541487842333E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.185E-21,  chi2_M: 1.026E+22
   max(B_n): 4.901E-12,  max(M): 4.592E+12,  min(M): 1.531E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.19380383307070E+19, method 2 =  5.19380383307069E+19, relative difference =  5.52042412873504E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.31683384872256E-24, method 2 =  3.12638803734444E-13, relative difference = -9.99999999985179E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.317E-24,  chi2_M: 5.194E+19
   max(B_n): 1.898E-13,  max(M): 3.375E+11,  min(M): 5.203E+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 =  2.35597928187297E+18, method 2 =  2.35597928187297E+18, relative difference =  8.69277593295246E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.94014528415174E-25, method 2 =  3.41060513164848E-13, relative difference = -9.99999999994171E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 9.940E-25,  chi2_M: 2.356E+18
   max(B_n): 1.571E-13,  max(M): 8.973E+10,  min(M): 5.341E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.129E-32 (  5 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.34025906009917E+17, method 2 =  2.34025906009917E+17, relative difference = -2.73473997350054E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.51588573941500E-25, method 2 =  3.41060513164848E-13, relative difference = -9.99999999999111E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.516E-25,  chi2_M: 2.340E+17
   max(B_n): 6.573E-14,  max(M): 2.197E+10,  min(M): 2.155E+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 =  2.56133094654444E+16, method 2 =  2.56133094654444E+16, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  1.19270888013277E-26, method 2 =  3.69482222595252E-13, relative difference = -9.99999999999935E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.193E-26,  chi2_M: 2.561E+16
   max(B_n): 1.632E-14,  max(M): 7.770E+09,  min(M): 1.011E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.637E-31 (  7 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.05594192975674E+20, method 2 =  1.05594192975674E+20, relative difference = -1.55160047520552E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.26987578275696E-22, method 2 = -2.01993088921881E-10, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.270E-22,  chi2_M: 1.056E+20
   max(B_n): 1.387E-12,  max(M): 6.065E+11,  min(M): 3.506E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.40827314459905E+16, method 2 =  1.40827314459905E+16, relative difference =  1.42017903818611E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.06258069477740E-26, method 2 =  2.27373675443232E-13, relative difference = -9.99999999999907E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.063E-26,  chi2_M: 1.408E+16
   max(B_n): 1.477E-14,  max(M): 4.925E+09,  min(M): 1.060E+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.28099720177681E+15, method 2 =  2.28099720177681E+15, relative difference =  1.09601186623666E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.00581978110327E-26, method 2 =  3.41060513164848E-13, relative difference = -9.99999999999941E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.006E-26,  chi2_M: 2.281E+15
   max(B_n): 1.477E-14,  max(M): 2.176E+09,  min(M): 1.782E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.976E-29 ( 10 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  7.72018651904902E+15, method 2 =  7.72018651904902E+15, relative difference = -6.47652746169131E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.41574036520761E-26, method 2 =  2.27373675443232E-13, relative difference = -9.99999999999876E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.416E-26,  chi2_M: 7.720E+15
   max(B_n): 1.554E-14,  max(M): 3.748E+09,  min(M): 7.086E+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 =  8.78891415016851E+18, method 2 =  8.78891415016850E+18, relative difference =  2.91276027534177E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.22088925886796E-24, method 2 = -1.18683374239481E-09, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.221E-24,  chi2_M: 8.789E+18
   max(B_n): 1.809E-13,  max(M): 9.048E+10,  min(M): 8.123E+10
 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.69135463170676E+15, method 2 =  1.69135463170676E+15, relative difference =  5.91242061986072E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.83086588346105E-25, method 2 =  2.84217094304040E-13, relative difference = -9.99999999995897E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.831E-25,  chi2_M: 1.691E+15
   max(B_n): 5.451E-14,  max(M): 1.241E+09,  min(M): 1.151E+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.69049637310699E+15, method 2 =  1.69049637310699E+15, relative difference = -8.13370570841784E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.36066847192890E-23, method 2 =  3.41060513164848E-13, relative difference = -9.99999999920210E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.361E-23,  chi2_M: 1.690E+15
   max(B_n): 2.232E-13,  max(M): 1.213E+09,  min(M): 1.168E+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.69038331326726E+15, method 2 =  1.69038331326726E+15, relative difference = -2.21843174300616E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.18449554807101E-22, method 2 =  3.97903932025656E-13, relative difference = -9.99999998399365E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.184E-22,  chi2_M: 1.690E+15
   max(B_n): 1.041E-12,  max(M): 1.205E+09,  min(M): 1.176E+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.69037938454619E+15, method 2 =  1.69037938454619E+15, relative difference =  1.47895793267212E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.44237157648392E-21, method 2 =  3.12638803734444E-13, relative difference = -9.99999952389970E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.442E-21,  chi2_M: 1.690E+15
   max(B_n): 4.997E-12,  max(M): 1.205E+09,  min(M): 1.177E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.69037929931428E+15, method 2 =  1.69037929931427E+15, relative difference =  6.65531103259707E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.73850910931665E-19, method 2 =  2.84217094304040E-13, relative difference = -9.99998776633721E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.739E-19,  chi2_M: 1.690E+15
   max(B_n): 2.412E-11,  max(M): 1.205E+09,  min(M): 1.178E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.69037928587883E+15, method 2 =  1.69037928587883E+15, relative difference = -7.39479009499414E-17
   2 methods of computing chi2_B that should agree: method 1 =  4.06054415547791E-18, method 2 =  2.84217094304040E-13, relative difference = -9.99971426868928E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.061E-18,  chi2_M: 1.690E+15
   max(B_n): 1.165E-10,  max(M): 1.205E+09,  min(M): 1.178E+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.69037926878986E+15, method 2 =  1.69037926878986E+15, relative difference = -1.47895803395042E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.48196756223193E-17, method 2 =  2.84217094304040E-13, relative difference = -9.99332988376779E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.482E-17,  chi2_M: 1.690E+15
   max(B_n): 5.630E-10,  max(M): 1.205E+09,  min(M): 1.178E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.69037918680765E+15, method 2 =  1.69037918680765E+15, relative difference = -7.39479052839425E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.21204466517316E-15, method 2 =  3.41060513164848E-13, relative difference = -9.87112021542552E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.212E-15,  chi2_M: 1.690E+15
   max(B_n): 2.720E-09,  max(M): 1.205E+09,  min(M): 1.178E+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 =  1.69037879276633E+15, method 2 =  1.69037879276633E+15, relative difference =  2.21843767565438E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.13661738850562E-14, method 2 =  2.55795384873636E-13, relative difference = -6.65542953404467E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.137E-14,  chi2_M: 1.690E+15
   max(B_n): 1.311E-08,  max(M): 1.205E+09,  min(M): 1.178E+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.69037693607148E+15, method 2 =  1.69037693607148E+15, relative difference = -8.87376044946564E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.16664942692232E-12, method 2 =  1.53477230924182E-12, relative difference = -1.36270052687965E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.167E-12,  chi2_M: 1.690E+15
   max(B_n): 6.259E-08,  max(M): 1.205E+09,  min(M): 1.178E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.69036896042904E+15, method 2 =  1.69036896042904E+15, relative difference = -4.43690115919805E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.39317006126928E-11, method 2 =  2.42152964347042E-11, relative difference = -5.89020789255556E-03
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.393E-11,  chi2_M: 1.690E+15
   max(B_n): 2.859E-07,  max(M): 1.205E+09,  min(M): 1.178E+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.69034520117298E+15, method 2 =  1.69034520117298E+15, relative difference =  1.47898784122035E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.30138033737966E-10, method 2 =  3.30402372128447E-10, relative difference = -4.00185042630566E-04
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.301E-10,  chi2_M: 1.690E+15
   max(B_n): 1.103E-06,  max(M): 1.206E+09,  min(M): 1.179E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.69031706685238E+15, method 2 =  1.69031706685237E+15, relative difference =  6.65555606141349E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.83219209312091E-09, method 2 =  1.83254655894416E-09, relative difference = -9.67233565332655E-05
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 1.832E-09,  chi2_M: 1.690E+15
   max(B_n): 3.010E-06,  max(M): 1.207E+09,  min(M): 1.180E+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.69029794258635E+15, method 2 =  1.69029794258635E+15, relative difference =  5.91611676737820E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.32539267838045E-09, method 2 =  7.32569560568663E-09, relative difference = -2.06760958850857E-05
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.325E-09,  chi2_M: 1.690E+15
   max(B_n): 7.436E-06,  max(M): 1.208E+09,  min(M): 1.181E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.69023687173721E+15, method 2 =  1.69023687173721E+15, relative difference = -2.21862394715468E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.04140935821221E-07, method 2 =  1.04141292922577E-07, relative difference = -1.71450708251992E-06
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.041E-07,  chi2_M: 1.690E+15
   max(B_n): 2.556E-05,  max(M): 1.208E+09,  min(M): 1.181E+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.68994431676488E+15, method 2 =  1.68994431676488E+15, relative difference =  2.95867736611089E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.35364892507666E-06, method 2 =  2.35364925060821E-06, relative difference = -6.91546495681062E-08
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.354E-06,  chi2_M: 1.690E+15
   max(B_n): 1.121E-04,  max(M): 1.208E+09,  min(M): 1.181E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.68853161715629E+15, method 2 =  1.68853161715629E+15, relative difference = -7.40288181340167E-17
   2 methods of computing chi2_B that should agree: method 1 =  5.48476798978807E-05, method 2 =  5.48476801895958E-05, relative difference = -2.65932042189998E-09
   Diagnostics:    0.00000000      sec.
   chi2_B: 5.485E-05,  chi2_M: 1.689E+15
   max(B_n): 5.302E-04,  max(M): 1.207E+09,  min(M): 1.180E+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.68172905652988E+15, method 2 =  1.68172905652988E+15, relative difference =  2.97313052931196E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.27584495166797E-03, method 2 =  1.27584495203337E-03, relative difference = -1.43196408118029E-10
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.276E-03,  chi2_M: 1.682E+15
   max(B_n): 2.559E-03,  max(M): 1.205E+09,  min(M): 1.178E+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.64942529364986E+15, method 2 =  1.64942529364986E+15, relative difference = -1.51567943672549E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.92275121097180E-02, method 2 =  2.92275121100829E-02, relative difference = -6.24227171484457E-12
   Diagnostics:    4.00000019E-03  sec.
   chi2_B: 2.923E-02,  chi2_M: 1.649E+15
   max(B_n): 1.226E-02,  max(M): 1.193E+09,  min(M): 1.167E+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 =  1.50569494431359E+15, method 2 =  1.50569494431359E+15, relative difference =  6.64145153556237E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.23171948448379E-01, method 2 =  6.23171948448743E-01, relative difference = -2.92266183744615E-13
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.232E-01,  chi2_M: 1.506E+15
   max(B_n): 5.653E-02,  max(M): 1.141E+09,  min(M): 1.114E+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.01935835310099E+15, method 2 =  1.01935835310099E+15, relative difference =  6.13130797524431E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.85374428369538E+00, method 2 =  9.85374428369556E+00, relative difference = -9.10374958056756E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.854E+00,  chi2_M: 1.019E+15
   max(B_n): 2.242E-01,  max(M): 9.397E+08,  min(M): 9.155E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.624E-13 ( 33 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.95799898598100E+14, method 2 =  2.95799898598100E+14, relative difference = -1.05645742774439E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.67855690917328E+01, method 2 =  6.67855690917331E+01, relative difference = -2.02144148250030E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.679E+01,  chi2_M: 2.958E+14
   max(B_n): 5.821E-01,  max(M): 5.066E+08,  min(M): 4.928E+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 =  2.83615634990848E+13, method 2 =  2.83615634990848E+13, relative difference =  4.82056463510595E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.49565931930945E+02, method 2 =  1.49565931930945E+02, relative difference = -4.75069908358648E-16
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 1.496E+02,  chi2_M: 2.836E+13
   max(B_n): 8.705E-01,  max(M): 1.569E+08,  min(M): 1.526E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.793E-12 ( 35 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.50822361369626E+12, method 2 =  1.50822361369626E+12, relative difference =  4.04682407142657E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.85775504676349E+02, method 2 =  1.85775504676349E+02, relative difference = -3.05979084593732E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.858E+02,  chi2_M: 1.508E+12
   max(B_n): 9.701E-01,  max(M): 3.617E+07,  min(M): 3.519E+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 =  6.77438484452246E+10, method 2 =  6.77438484452246E+10, relative difference = -2.81553037890181E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.94900713398323E+02, method 2 =  1.94900713398323E+02, relative difference =  1.45826605428159E-16
   Diagnostics:    3.00000003E-03  sec.
   chi2_B: 1.949E+02,  chi2_M: 6.774E+10
   max(B_n): 9.937E-01,  max(M): 7.666E+06,  min(M): 7.459E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.92968473703522E+09, method 2 =  2.92968473703522E+09, relative difference =  3.25521140329709E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.96872683803809E+02, method 2 =  1.96872683803809E+02, relative difference = -7.94012652476460E-16
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 1.969E+02,  chi2_M: 2.930E+09
   max(B_n): 9.987E-01,  max(M): 1.594E+06,  min(M): 1.551E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.281E-10 ( 38 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.25691883251554E+08, method 2 =  1.25691883251554E+08, relative difference =  5.92765451848040E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.97284447596809E+02, method 2 =  1.97284447596809E+02, relative difference = -7.20323111543225E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.973E+02,  chi2_M: 1.257E+08
   max(B_n): 9.997E-01,  max(M): 3.302E+05,  min(M): 3.213E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.38361571542758E+06, method 2 =  5.38361571542758E+06, relative difference =  3.45984046352577E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97369808504495E+02, method 2 =  1.97369808504495E+02, relative difference = -7.20011577397785E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.974E+02,  chi2_M: 5.384E+06
   max(B_n): 9.999E-01,  max(M): 6.834E+04,  min(M): 6.649E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-08 ( 40 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.30511118405684E+05, method 2 =  2.30511118405684E+05, relative difference =  5.05031265442267E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97387477772499E+02, method 2 =  1.97387477772499E+02, relative difference =  0.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 2.305E+05
   max(B_n): 1.000E+00,  max(M): 1.414E+04,  min(M): 1.376E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   1.6649999618530273      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
