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.698E+02 m^2. Volume: 2.500E+02 m^3.
 G =    0.0000000000000000       ; curpol =    1.0000000000000000     
 Done initializing plasma surface. Took    4.69999984E-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   0.00000000     sec.
 Building inductance matrix.
   Number of OpenMP threads:           4
 Done. Took   3.21099997     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  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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.42522178585173E+25, method 2 =  2.42522178585173E+25, relative difference =  4.42739641489285E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.94985408229689E-18, method 2 = -6.25277607468888E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.950E-18,  chi2_M: 2.425E+25
   max(B_n): 2.200E-10,  max(M): 4.398E+14,  min(M): 3.139E+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 =  5.00575838200666E+20, method 2 =  5.00575838200666E+20, relative difference =  4.58224273917212E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.00493690263592E-16, method 2 =  2.50111042987555E-12, relative difference = -9.99919643969793E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.005E-16,  chi2_M: 5.006E+20
   max(B_n): 6.381E-10,  max(M): 1.061E+12,  min(M): 2.174E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.30026909342494E+20, method 2 =  4.30026909342494E+20, relative difference = -7.61998825843291E-17
   2 methods of computing chi2_B that should agree: method 1 =  8.92628832234054E-17, method 2 =  4.83169060316868E-12, relative difference = -9.99963051759104E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 8.926E-17,  chi2_M: 4.300E+20
   max(B_n): 5.977E-10,  max(M): 1.065E+12,  min(M): 3.318E+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 =  5.80956476396187E+19, method 2 =  5.80956476396187E+19, relative difference =  3.52522105047221E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.03673437988903E-19, method 2 =  4.54747350886464E-13, relative difference = -9.99996465413087E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 8.037E-19,  chi2_M: 5.810E+19
   max(B_n): 5.171E-11,  max(M): 3.541E+11,  min(M): 6.388E+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 =  3.45532524570283E+19, method 2 =  3.45532524570284E+19, relative difference = -6.51979145176467E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.11970211121135E-19, method 2 =  3.12638803734444E-12, relative difference = -9.99999864399304E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.120E-19,  chi2_M: 3.455E+19
   max(B_n): 2.571E-11,  max(M): 2.166E+11,  min(M): 8.510E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.456E-32 (  6 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.86869076928987E+18, method 2 =  2.86869076928987E+18, relative difference = -5.35435891676895E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.94770532977777E-19, method 2 = -1.81898940354586E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.948E-19,  chi2_M: 2.869E+18
   max(B_n): 3.371E-11,  max(M): 9.048E+10,  min(M): 4.588E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.637E-31 (  7 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.46948026677564E+16, method 2 =  3.46948026677564E+16, relative difference = -1.15291043396462E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.16921235549116E-19, method 2 =  1.70530256582424E-13, relative difference = -9.99997455924081E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.169E-19,  chi2_M: 3.469E+16
   max(B_n): 2.504E-11,  max(M): 1.190E+10,  min(M): 3.882E+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 =  2.34082530034712E+16, method 2 =  2.34082530034712E+16, relative difference =  8.54399514437673E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.44522378880373E-21, method 2 =  3.41060513164848E-13, relative difference = -9.99999944612624E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.445E-21,  chi2_M: 2.341E+16
   max(B_n): 5.221E-12,  max(M): 7.461E+09,  min(M): 1.729E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.64834329323041E+16, method 2 =  3.64834329323040E+16, relative difference =  9.86749247714680E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.07065086788694E-20, method 2 = -1.70530256582424E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.071E-20,  chi2_M: 3.648E+16
   max(B_n): 5.564E-12,  max(M): 1.045E+10,  min(M): 2.032E+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 =  9.56096694069942E+14, method 2 =  9.56096694069941E+14, relative difference =  2.61479828923780E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.20141472245878E-22, method 2 =  2.27373675443232E-13, relative difference = -9.99999998943224E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.201E-22,  chi2_M: 9.561E+14
   max(B_n): 6.007E-13,  max(M): 1.449E+09,  min(M): 7.727E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.438E-28 ( 11 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.94698917570292E+14, method 2 =  4.94698917570292E+14, relative difference = -1.89509207864153E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.32977534402904E-24, method 2 =  3.41060513164848E-13, relative difference = -9.99999999962882E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.330E-24,  chi2_M: 4.947E+14
   max(B_n): 1.398E-13,  max(M): 8.548E+08,  min(M): 4.145E+08
 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 =  4.70057658518382E+14, method 2 =  4.70057658518382E+14, relative difference =  2.65924823763108E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.65750921027811E-25, method 2 =  2.27373675443232E-13, relative difference = -9.99999999997662E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.658E-25,  chi2_M: 4.701E+14
   max(B_n): 3.608E-14,  max(M): 6.908E+08,  min(M): 5.721E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.68375931427876E+14, method 2 =  4.68375931427875E+14, relative difference =  9.34078740267998E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.43301625413502E-26, method 2 =  2.84217094304040E-13, relative difference = -9.99999999999688E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 4.433E-26,  chi2_M: 4.684E+14
   max(B_n): 2.165E-14,  max(M): 6.458E+08,  min(M): 6.152E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.68285246611209E+14, method 2 =  4.68285246611209E+14, relative difference = -5.33862644209163E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.26244859729136E-25, method 2 =  1.70530256582424E-13, relative difference = -9.99999999995001E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.262E-25,  chi2_M: 4.683E+14
   max(B_n): 4.796E-14,  max(M): 6.347E+08,  min(M): 6.248E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.68282418877647E+14, method 2 =  4.68282418877645E+14, relative difference =  1.33466466987585E-15
   2 methods of computing chi2_B that should agree: method 1 =  9.61826306325308E-24, method 2 =  3.41060513164848E-13, relative difference = -9.99999999943598E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.618E-24,  chi2_M: 4.683E+14
   max(B_n): 1.948E-13,  max(M): 6.327E+08,  min(M): 6.247E+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 =  4.68282329910255E+14, method 2 =  4.68282329910254E+14, relative difference =  6.00599215549946E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.24016303177279E-22, method 2 =  2.84217094304040E-13, relative difference = -9.99999998423625E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.240E-22,  chi2_M: 4.683E+14
   max(B_n): 9.021E-13,  max(M): 6.326E+08,  min(M): 6.250E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.68282311894665E+14, method 2 =  4.68282311894664E+14, relative difference =  5.33865989916431E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.23142212442444E-21, method 2 =  2.84217094304040E-13, relative difference = -9.99999963187140E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 5.231E-21,  chi2_M: 4.683E+14
   max(B_n): 4.319E-12,  max(M): 6.325E+08,  min(M): 6.249E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.859E-24 ( 18 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.68282311422089E+14, method 2 =  4.68282311422090E+14, relative difference = -1.06773198091038E-15
   2 methods of computing chi2_B that should agree: method 1 =  1.22177620840000E-19, method 2 =  5.11590769747272E-13, relative difference = -9.99999522362017E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.222E-19,  chi2_M: 4.683E+14
   max(B_n): 2.084E-11,  max(M): 6.324E+08,  min(M): 6.249E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.68282311255118E+14, method 2 =  4.68282311255119E+14, relative difference = -1.26793172778317E-15
   2 methods of computing chi2_B that should agree: method 1 =  2.85218623515467E-18, method 2 =  3.97903932025656E-13, relative difference = -9.99985664048223E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.852E-18,  chi2_M: 4.683E+14
   max(B_n): 1.007E-10,  max(M): 6.324E+08,  min(M): 6.249E+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 =  4.68282310744758E+14, method 2 =  4.68282310744758E+14, relative difference =  6.67332489034230E-17
   2 methods of computing chi2_B that should agree: method 1 =  6.64435791224944E-17, method 2 =  3.97903932025656E-13, relative difference = -9.99666087813589E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.644E-17,  chi2_M: 4.683E+14
   max(B_n): 4.860E-10,  max(M): 6.324E+08,  min(M): 6.249E+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 =  4.68282308310301E+14, method 2 =  4.68282308310301E+14, relative difference =  2.66932997001395E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.53246127254670E-15, method 2 =  1.70530256582424E-13, relative difference = -9.82187178121429E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.532E-15,  chi2_M: 4.683E+14
   max(B_n): 2.338E-09,  max(M): 6.324E+08,  min(M): 6.249E+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 =  4.68282297149884E+14, method 2 =  4.68282297149883E+14, relative difference =  8.67532260930144E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.37753363757610E-14, method 2 =  3.97903932025656E-13, relative difference = -8.43516523270451E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.378E-14,  chi2_M: 4.683E+14
   max(B_n): 1.105E-08,  max(M): 6.324E+08,  min(M): 6.249E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.336E-20 ( 23 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.68282253287846E+14, method 2 =  4.68282253287846E+14, relative difference = -6.67332570914045E-17
   2 methods of computing chi2_B that should agree: method 1 =  6.30954485086082E-13, method 2 =  9.09494701772928E-13, relative difference = -1.80817529758831E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.310E-13,  chi2_M: 4.683E+14
   max(B_n): 4.887E-08,  max(M): 6.324E+08,  min(M): 6.249E+08
 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 =  4.68282118232254E+14, method 2 =  4.68282118232255E+14, relative difference = -1.33466552675415E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.37968893748233E-12, method 2 =  9.60653778747655E-12, relative difference = -1.19480744268165E-02
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.380E-12,  chi2_M: 4.683E+14
   max(B_n): 1.901E-07,  max(M): 6.324E+08,  min(M): 6.249E+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 =  4.68281613385482E+14, method 2 =  4.68281613385482E+14, relative difference = -7.34066831099410E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.74226716893043E-10, method 2 =  1.74622982740402E-10, relative difference = -1.13592142339724E-03
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.742E-10,  chi2_M: 4.683E+14
   max(B_n): 7.430E-07,  max(M): 6.324E+08,  min(M): 6.249E+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 =  4.68279229590759E+14, method 2 =  4.68279229590758E+14, relative difference =  8.67537943878127E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.96632852005052E-09, method 2 =  3.96664745494490E-09, relative difference = -4.02036884251946E-05
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.966E-09,  chi2_M: 4.683E+14
   max(B_n): 3.301E-06,  max(M): 6.324E+08,  min(M): 6.249E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.274E-17 ( 27 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.68267714705017E+14, method 2 =  4.68267714705017E+14, relative difference =  2.00205986994122E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.25117879009840E-08, method 2 =  9.25120389183576E-08, relative difference = -1.35667593681226E-06
   Diagnostics:    4.00000019E-03  sec.
   chi2_B: 9.251E-08,  chi2_M: 4.683E+14
   max(B_n): 1.563E-05,  max(M): 6.324E+08,  min(M): 6.249E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.68212070349001E+14, method 2 =  4.68212070349000E+14, relative difference =  6.00689341029502E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.16042628110822E-06, method 2 =  2.16042661804750E-06, relative difference = -7.79798080910216E-08
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.160E-06,  chi2_M: 4.682E+14
   max(B_n): 7.517E-05,  max(M): 6.323E+08,  min(M): 6.249E+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 =  4.67943284822642E+14, method 2 =  4.67943284822643E+14, relative difference = -5.34252778292895E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.04324357663318E-05, method 2 =  5.04324361827457E-05, relative difference = -4.12843400078257E-09
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.043E-05,  chi2_M: 4.679E+14
   max(B_n): 3.635E-04,  max(M): 6.322E+08,  min(M): 6.247E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.66647519787840E+14, method 2 =  4.66647519787840E+14, relative difference =  4.01802199838641E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.17469919652481E-03, method 2 =  1.17469919672430E-03, relative difference = -8.49116304610792E-11
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 1.175E-03,  chi2_M: 4.666E+14
   max(B_n): 1.754E-03,  max(M): 6.313E+08,  min(M): 6.238E+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 =  4.60460402855430E+14, method 2 =  4.60460402855430E+14, relative difference =  2.03600568949323E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.70739091190901E-02, method 2 =  2.70739091192809E-02, relative difference = -3.52231933796718E-12
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.707E-02,  chi2_M: 4.605E+14
   max(B_n): 8.423E-03,  max(M): 6.271E+08,  min(M): 6.197E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.32225177481968E+14, method 2 =  4.32225177481968E+14, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  5.93593877920285E-01, method 2 =  5.93593877920739E-01, relative difference = -3.82484754275460E-13
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.936E-01,  chi2_M: 4.322E+14
   max(B_n): 3.943E-02,  max(M): 6.077E+08,  min(M): 6.003E+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.26528929088661E+14, method 2 =  3.26528929088661E+14, relative difference =  1.91407236640370E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.04741942614484E+01, method 2 =  1.04741942614485E+01, relative difference = -6.86854283966291E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.047E+01,  chi2_M: 3.265E+14
   max(B_n): 1.655E-01,  max(M): 5.284E+08,  min(M): 5.215E+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.22555692796166E+14, method 2 =  1.22555692796166E+14, relative difference = -5.09972230371622E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.18227641769586E+01, method 2 =  9.18227641769590E+01, relative difference = -1.77978555415599E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.182E+01,  chi2_M: 1.226E+14
   max(B_n): 4.894E-01,  max(M): 3.242E+08,  min(M): 3.191E+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.48574055354331E+13, method 2 =  1.48574055354331E+13, relative difference =  5.91561055464170E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.60002964213981E+02, method 2 =  2.60002964213982E+02, relative difference = -2.18626041563227E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.600E+02,  chi2_M: 1.486E+13
   max(B_n): 8.224E-01,  max(M): 1.130E+08,  min(M): 1.110E+08
 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 =  8.62593441331052E+11, method 2 =  8.62593441331051E+11, relative difference =  5.66061862523024E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.52583190451681E+02, method 2 =  3.52583190451681E+02, relative difference = -4.03049694371335E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.526E+02,  chi2_M: 8.626E+11
   max(B_n): 9.572E-01,  max(M): 2.724E+07,  min(M): 2.673E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.95790859588560E+10, method 2 =  3.95790859588560E+10, relative difference = -9.63816412938522E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.77869678665884E+02, method 2 =  3.77869678665884E+02, relative difference = -8.27372031643964E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.779E+02,  chi2_M: 3.958E+10
   max(B_n): 9.908E-01,  max(M): 5.835E+06,  min(M): 5.725E+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.71949550967202E+09, method 2 =  1.71949550967202E+09, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  3.83440811876165E+02, method 2 =  3.83440811876165E+02, relative difference = -5.18859651478834E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.834E+02,  chi2_M: 1.719E+09
   max(B_n): 9.981E-01,  max(M): 1.216E+06,  min(M): 1.193E+06
 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 =  7.38415893686688E+07, method 2 =  7.38415893686689E+07, relative difference = -1.00899515579565E-15
   2 methods of computing chi2_B that should agree: method 1 =  3.84608934056748E+02, method 2 =  3.84608934056748E+02, relative difference = -5.17283787233797E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.846E+02,  chi2_M: 7.384E+07
   max(B_n): 9.996E-01,  max(M): 2.521E+05,  min(M): 2.473E+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 =  3.16339611003375E+06, method 2 =  3.16339611003375E+06, relative difference = -5.88811860558014E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.84851301120021E+02, method 2 =  3.84851301120021E+02, relative difference =  2.21553436465117E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.849E+02,  chi2_M: 3.163E+06
   max(B_n): 9.999E-01,  max(M): 5.217E+04,  min(M): 5.118E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   3.3410000801086426      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
