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: 1.184E+02 m^2. Volume: 5.918E+00 m^3.
 G =    0.0000000000000000       ; curpol =    1.0000000000000000     
 Done initializing plasma surface. Took    2.19999999E-02  sec.
 Initializing coil surface.
   Building a plain circular torus.
   Evaluating coil surface & derivatives:   1.20000001E-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.49999997E-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.79900002     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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.71219902190540E+33, method 2 =  2.71219902190540E+33, relative difference =  1.06271838395260E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.64028298673724E-14, method 2 =  2.47837306233123E-11, relative difference = -9.95458733748457E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.640E-14,  chi2_M: 2.712E+33
   max(B_n): 6.634E-08,  max(M): 3.585E+18,  min(M): 2.640E+17
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-34 (  2 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.91525815751022E+19, method 2 =  6.91525815751022E+19, relative difference = -5.92313389710780E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.56377674699508E-24, method 2 =  0.00000000000000E+00, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.564E-24,  chi2_M: 6.915E+19
   max(B_n): 2.846E-13,  max(M): 3.692E+11,  min(M): 1.069E+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 =  2.18970668593072E+20, method 2 =  2.18970668593072E+20, relative difference =  4.48936839950400E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.71767954909436E-23, method 2 =  3.58113538823090E-12, relative difference = -9.99999999973653E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.718E-23,  chi2_M: 2.190E+20
   max(B_n): 1.735E-12,  max(M): 7.426E+11,  min(M): 2.475E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.336E-33 (  4 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.48152369949408E+18, method 2 =  1.48152369949408E+18, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  8.28063327187523E-26, method 2 =  4.97379915032070E-14, relative difference = -9.99999999996670E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.281E-26,  chi2_M: 1.482E+18
   max(B_n): 8.515E-14,  max(M): 5.728E+10,  min(M): 1.204E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.69912146138542E+19, method 2 =  4.69912146138542E+19, relative difference =  8.71652293659248E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.11615959716679E-24, method 2 = -5.68434188608080E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.116E-24,  chi2_M: 4.699E+19
   max(B_n): 1.698E-13,  max(M): 2.513E+11,  min(M): 1.485E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.84316203053373E+18, method 2 =  1.84316203053373E+18, relative difference = -2.77783500049498E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.29423751468618E-24, method 2 = -7.81597009336110E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    3.00000003E-03  sec.
   chi2_B: 2.294E-24,  chi2_M: 1.843E+18
   max(B_n): 2.034E-13,  max(M): 7.926E+10,  min(M): 4.055E+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 =  8.74435889636028E+17, method 2 =  8.74435889636026E+17, relative difference =  8.78280511015702E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.29556580681661E-24, method 2 = -3.55271367880050E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 5.296E-24,  chi2_M: 8.744E+17
   max(B_n): 2.661E-13,  max(M): 3.458E+10,  min(M): 1.941E+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 =  8.37641726762776E+17, method 2 =  8.37641726762776E+17, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  2.80477299143543E-25, method 2 =  0.00000000000000E+00, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.805E-25,  chi2_M: 8.376E+17
   max(B_n): 1.058E-13,  max(M): 2.833E+10,  min(M): 2.478E+10
 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 =  8.37567776722500E+17, method 2 =  8.37567776722500E+17, relative difference =  1.52823453286227E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.84529526418303E-25, method 2 =  7.10542735760100E-15, relative difference = -9.99999999835469E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.845E-25,  chi2_M: 8.376E+17
   max(B_n): 1.527E-13,  max(M): 2.789E+10,  min(M): 2.502E+10
 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 =  8.36174329261148E+17, method 2 =  8.36174329261148E+17, relative difference =  7.65390633990774E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.34822669338195E-23, method 2 =  7.10542735760100E-15, relative difference = -9.99999996205079E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.348E-23,  chi2_M: 8.362E+17
   max(B_n): 6.087E-13,  max(M): 2.689E+10,  min(M): 2.619E+10
 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 =  8.36170461254582E+17, method 2 =  8.36170461254582E+17, relative difference = -3.06157669831939E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.14326521678119E-22, method 2 =  2.84217094304040E-14, relative difference = -9.99999977881238E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.143E-22,  chi2_M: 8.362E+17
   max(B_n): 2.863E-12,  max(M): 2.689E+10,  min(M): 2.624E+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 =  8.36170300081111E+17, method 2 =  8.36170300081110E+17, relative difference =  7.65394322111080E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.31194205842803E-21, method 2 =  2.13162820728030E-14, relative difference = -9.99999313957328E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.312E-21,  chi2_M: 8.362E+17
   max(B_n): 1.373E-11,  max(M): 2.689E+10,  min(M): 2.625E+10
 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 =  8.36170292007851E+17, method 2 =  8.36170292007851E+17, relative difference = -7.65394329500995E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.70619764658958E-19, method 2 =  1.42108547152020E-14, relative difference = -9.99975987691157E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.706E-19,  chi2_M: 8.362E+17
   max(B_n): 6.625E-11,  max(M): 2.690E+10,  min(M): 2.625E+10
 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 =  8.36170291311565E+17, method 2 =  8.36170291311565E+17, relative difference =  2.29618299041503E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.98248315177958E-18, method 2 =  0.00000000000000E+00, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.982E-18,  chi2_M: 8.362E+17
   max(B_n): 3.200E-10,  max(M): 2.690E+10,  min(M): 2.625E+10
 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 =  8.36170289422004E+17, method 2 =  8.36170289422004E+17, relative difference = -7.65394331867968E-17
   2 methods of computing chi2_B that should agree: method 1 =  9.28000182791337E-17, method 2 =  7.10542735760100E-15, relative difference = -9.74215869148538E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.280E-17,  chi2_M: 8.362E+17
   max(B_n): 1.544E-09,  max(M): 2.690E+10,  min(M): 2.625E+10
 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 =  8.36170280442371E+17, method 2 =  8.36170280442372E+17, relative difference = -3.82697170043769E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.14426581539413E-15, method 2 =  1.42108547152020E-14, relative difference = -7.37786608006615E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.144E-15,  chi2_M: 8.362E+17
   max(B_n): 7.413E-09,  max(M): 2.690E+10,  min(M): 2.625E+10
 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 =  8.36170238750745E+17, method 2 =  8.36170238750744E+17, relative difference =  6.88854940425237E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.79549794558945E-14, method 2 =  6.39488462184090E-14, relative difference = -1.42925111506596E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.795E-14,  chi2_M: 8.362E+17
   max(B_n): 3.481E-08,  max(M): 2.690E+10,  min(M): 2.625E+10
 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 =  8.36170056458342E+17, method 2 =  8.36170056458341E+17, relative difference =  1.53078909022590E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.01071465466235E-12, method 2 =  1.03028696685215E-12, relative difference = -9.58956229308226E-03
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.011E-12,  chi2_M: 8.362E+17
   max(B_n): 1.543E-07,  max(M): 2.689E+10,  min(M): 2.625E+10
 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 =  8.36169242602951E+17, method 2 =  8.36169242602952E+17, relative difference = -1.53079058016464E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.19575331569116E-11, method 2 =  2.19699813897023E-11, relative difference = -2.83381222888834E-04
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.196E-11,  chi2_M: 8.362E+17
   max(B_n): 6.663E-07,  max(M): 2.689E+10,  min(M): 2.625E+10
 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 =  8.36165357781141E+17, method 2 =  8.36165357781141E+17, relative difference =  3.82699423053301E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.06951372153332E-10, method 2 =  5.06979347392189E-10, relative difference = -2.75908780721578E-05
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.070E-10,  chi2_M: 8.362E+17
   max(B_n): 3.047E-06,  max(M): 2.689E+10,  min(M): 2.625E+10
 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 =  8.36146590717103E+17, method 2 =  8.36146590717103E+17, relative difference =  7.65416025258344E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.18319932632075E-08, method 2 =  1.18319931630140E-08, relative difference =  4.23400937772081E-09
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.183E-08,  chi2_M: 8.361E+17
   max(B_n): 1.468E-05,  max(M): 2.689E+10,  min(M): 2.625E+10
 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 =  8.36055901276028E+17, method 2 =  8.36055901276027E+17, relative difference =  4.59299431310659E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.76319295282497E-07, method 2 =  2.76319312320084E-07, relative difference = -3.08295269217179E-08
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.763E-07,  chi2_M: 8.361E+17
   max(B_n): 7.117E-05,  max(M): 2.689E+10,  min(M): 2.625E+10
 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 =  8.35617837665626E+17, method 2 =  8.35617837665627E+17, relative difference = -5.36130249746138E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.45026529054476E-06, method 2 =  6.45026531032045E-06, relative difference = -1.53293591571724E-09
   Diagnostics:    0.00000000      sec.
   chi2_B: 6.450E-06,  chi2_M: 8.356E+17
   max(B_n): 3.440E-04,  max(M): 2.689E+10,  min(M): 2.624E+10
 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 =  8.33506049076388E+17, method 2 =  8.33506049076386E+17, relative difference =  1.22854537304762E-15
   2 methods of computing chi2_B that should agree: method 1 =  1.50240798719668E-04, method 2 =  1.50240798717505E-04, relative difference =  7.19920054355308E-12
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.502E-04,  chi2_M: 8.335E+17
   max(B_n): 1.657E-03,  max(M): 2.686E+10,  min(M): 2.620E+10
 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 =  8.23419059018715E+17, method 2 =  8.23419059018714E+17, relative difference =  2.33174102417316E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.46412772242956E-03, method 2 =  3.46412772245941E-03, relative difference = -4.30879358514967E-12
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 3.464E-03,  chi2_M: 8.234E+17
   max(B_n): 7.882E-03,  max(M): 2.673E+10,  min(M): 2.601E+10
 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 =  7.77195911714799E+17, method 2 =  7.77195911714798E+17, relative difference =  5.76431235994971E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.63076453612680E-02, method 2 =  7.63076453612967E-02, relative difference = -1.88504402077883E-13
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.631E-02,  chi2_M: 7.772E+17
   max(B_n): 3.640E-02,  max(M): 2.603E+10,  min(M): 2.521E+10
 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 =  6.00578097231274E+17, method 2 =  6.00578097231273E+17, relative difference =  8.52511942011160E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.37699927902657E+00, method 2 =  1.37699927902659E+00, relative difference = -6.12759579156490E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.377E+00,  chi2_M: 6.006E+17
   max(B_n): 1.530E-01,  max(M): 2.291E+10,  min(M): 2.213E+10
 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 =  2.39220065615717E+17, method 2 =  2.39220065615718E+17, relative difference = -1.00326032175552E-15
   2 methods of computing chi2_B that should agree: method 1 =  1.28107239551293E+01, method 2 =  1.28107239551293E+01, relative difference = -4.85315971187681E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.281E+01,  chi2_M: 2.392E+17
   max(B_n): 4.652E-01,  max(M): 1.446E+10,  min(M): 1.397E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.976E-16 ( 29 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    2.00000009E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.08921040579141E+16, method 2 =  3.08921040579141E+16, relative difference =  3.88448775696966E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.86405787434593E+01, method 2 =  3.86405787434594E+01, relative difference = -5.51655352118948E-16
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 3.864E+01,  chi2_M: 3.089E+16
   max(B_n): 8.082E-01,  max(M): 5.196E+09,  min(M): 5.020E+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.84114947767846E+15, method 2 =  1.84114947767846E+15, relative difference = -6.78923691506108E-17
   2 methods of computing chi2_B that should agree: method 1 =  5.37904099036591E+01, method 2 =  5.37904099036591E+01, relative difference = -3.96284060875933E-16
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 5.379E+01,  chi2_M: 1.841E+15
   max(B_n): 9.536E-01,  max(M): 1.267E+09,  min(M): 1.227E+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 =  8.50383671637896E+13, method 2 =  8.50383671637895E+13, relative difference =  9.18702964386958E-17
   2 methods of computing chi2_B that should agree: method 1 =  5.80292386796052E+01, method 2 =  5.80292386796052E+01, relative difference =  4.28559745353749E-16
   Diagnostics:    4.99999989E-03  sec.
   chi2_B: 5.803E+01,  chi2_M: 8.504E+13
   max(B_n): 9.902E-01,  max(M): 2.718E+08,  min(M): 2.642E+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:    2.00000009E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.69976013757221E+12, method 2 =  3.69976013757221E+12, relative difference =  1.97964691700424E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.89687536514121E+01, method 2 =  5.89687536514122E+01, relative difference = -6.62721323530460E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.897E+01,  chi2_M: 3.700E+12
   max(B_n): 9.980E-01,  max(M): 5.663E+07,  min(M): 5.516E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.624E-13 ( 33 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.58929457011977E+11, method 2 =  1.58929457011977E+11, relative difference =  1.24812770106898E-15
   2 methods of computing chi2_B that should agree: method 1 =  5.91660019451950E+01, method 2 =  5.91660019451950E+01, relative difference =  1.20093079200834E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.917E+01,  chi2_M: 1.589E+11
   max(B_n): 9.996E-01,  max(M): 1.173E+07,  min(M): 1.144E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.848E-13 ( 34 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    3.00000003E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.80901216950044E+09, method 2 =  6.80901216950043E+09, relative difference =  2.80121196046040E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.92069390432899E+01, method 2 =  5.92069390432899E+01, relative difference =  0.00000000000000E+00
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 5.921E+01,  chi2_M: 6.809E+09
   max(B_n): 9.999E-01,  max(M): 2.429E+06,  min(M): 2.367E+06
 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 =  2.91558230006736E+08, method 2 =  2.91558230006736E+08, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  5.92154148057818E+01, method 2 =  5.92154148057818E+01, relative difference = -5.39967898799259E-16
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 5.922E+01,  chi2_M: 2.916E+08
   max(B_n): 1.000E+00,  max(M): 5.025E+05,  min(M): 4.899E+05
 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 =  1.24829471811737E+07, method 2 =  1.24829471811737E+07, relative difference = -7.46075875431135E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.92171687856229E+01, method 2 =  5.92171687856230E+01, relative difference = -9.59914498219118E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.922E+01,  chi2_M: 1.248E+07
   max(B_n): 1.000E+00,  max(M): 1.040E+05,  min(M): 1.014E+05
 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 =  5.34439741720851E+05, method 2 =  5.34439741720851E+05, relative difference =  2.17826843213581E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.92175317180145E+01, method 2 =  5.92175317180145E+01, relative difference = -1.19988576886927E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.922E+01,  chi2_M: 5.344E+05
   max(B_n): 1.000E+00,  max(M): 2.152E+04,  min(M): 2.097E+04
 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 =  2.28811708682099E+04, method 2 =  2.28811708682099E+04, relative difference =  7.94972169047993E-17
   2 methods of computing chi2_B that should agree: method 1 =  5.92176068141548E+01, method 2 =  5.92176068141549E+01, relative difference = -6.59936335986212E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.922E+01,  chi2_M: 2.288E+04
   max(B_n): 1.000E+00,  max(M): 4.452E+03,  min(M): 4.340E+03
 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 =  9.79619238969574E+02, method 2 =  9.79619238969571E+02, relative difference =  1.21854670528173E-15
   2 methods of computing chi2_B that should agree: method 1 =  5.92176223526011E+01, method 2 =  5.92176223526012E+01, relative difference = -2.39976786480651E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.922E+01,  chi2_M: 9.796E+02
   max(B_n): 1.000E+00,  max(M): 9.211E+02,  min(M): 8.979E+02
 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 =  4.19407581421601E+01, method 2 =  4.19407581421600E+01, relative difference =  5.08247418908130E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.92176255677210E+01, method 2 =  5.92176255677210E+01, relative difference =  1.79982580088641E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.922E+01,  chi2_M: 4.194E+01
   max(B_n): 1.000E+00,  max(M): 1.906E+02,  min(M): 1.858E+02
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   1.9090000391006470      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
