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: 4.737E+02 m^2. Volume: 9.468E+01 m^3.
 G =    0.0000000000000000       ; curpol =    1.0000000000000000     
 Done initializing plasma surface. Took    1.89999994E-02  sec.
 Initializing coil surface.
   Building a plain circular torus.
   Evaluating coil surface & derivatives:   1.09999999E-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.30000003E-02  sec.
 Using Bn coefficients in FOCUS file.
 Initializing basis functions
 Done. Took   1.00000005E-03 sec.
 Building inductance matrix.
   Number of OpenMP threads:           4
 Done. Took   1.56900001     sec.
 Computing g and RHS_B:
  js=   1 j_RZetaZ= 1 dgemm took  0.000E+00 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  1.000E-03 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  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  1.000E-03 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:   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 =  9.07592221419639E+29, method 2 =  9.07592221419639E+29, relative difference = -7.75334368419272E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.68761569875255E-15, method 2 = -4.60431692772545E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.688E-15,  chi2_M: 9.076E+29
   max(B_n): 4.067E-09,  max(M): 6.999E+16,  min(M): 9.542E+14
 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 =  9.61443282979016E+20, method 2 =  9.61443282979015E+20, relative difference =  3.40820936399586E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.02553841265196E-21, method 2 =  1.03455022326671E-11, relative difference = -9.99999999608421E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.026E-21,  chi2_M: 9.614E+20
   max(B_n): 5.411E-12,  max(M): 1.702E+12,  min(M): 2.005E+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 =  2.72360110890550E+20, method 2 =  2.72360110890550E+20, relative difference =  6.01556518185734E-17
   2 methods of computing chi2_B that should agree: method 1 =  6.79450911493069E-22, method 2 =  6.73594513500575E-12, relative difference = -9.99999999798261E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 6.795E-22,  chi2_M: 2.724E+20
   max(B_n): 3.078E-12,  max(M): 8.160E+11,  min(M): 1.059E+11
 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 =  8.40246201823638E+19, method 2 =  8.40246201823635E+19, relative difference =  1.26743804100352E-15
   2 methods of computing chi2_B that should agree: method 1 =  1.70798026004009E-22, method 2 =  1.25055521493778E-12, relative difference = -9.99999999726844E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.708E-22,  chi2_M: 8.402E+19
   max(B_n): 1.688E-12,  max(M): 4.299E+11,  min(M): 3.156E+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 =  2.83941505133133E+19, method 2 =  2.83941505133133E+19, relative difference = -7.21275320083884E-17
   2 methods of computing chi2_B that should agree: method 1 =  4.51743513243416E-24, method 2 =  5.00222085975111E-12, relative difference = -9.99999999998194E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.517E-24,  chi2_M: 2.839E+19
   max(B_n): 2.608E-13,  max(M): 2.489E+11,  min(M): 3.748E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.53958520995691E+21, method 2 =  3.53958520995691E+21, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  2.68311825419232E-21, method 2 = -2.34280150834820E-10, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.683E-21,  chi2_M: 3.540E+21
   max(B_n): 6.073E-12,  max(M): 2.421E+12,  min(M): 7.494E+11
 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 =  2.92638997311408E+18, method 2 =  2.92638997311408E+18, relative difference =  1.74959593459501E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.20387470102896E-26, method 2 =  3.86535248253495E-12, relative difference = -9.99999999999952E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.204E-26,  chi2_M: 2.926E+18
   max(B_n): 4.807E-14,  max(M): 8.117E+10,  min(M): 9.298E+09
 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.58203553125218E+22, method 2 =  1.58203553125217E+22, relative difference =  6.62801801404585E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.26978495948379E-22, method 2 = -2.66578581431531E-08, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.270E-22,  chi2_M: 1.582E+22
   max(B_n): 2.666E-12,  max(M): 4.880E+12,  min(M): 1.597E+12
 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.37617370033672E+19, method 2 =  2.37617370033672E+19, relative difference =  3.44755940983571E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.46231689475471E-24, method 2 = -8.06608113634866E-11, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.462E-24,  chi2_M: 2.376E+19
   max(B_n): 2.498E-13,  max(M): 2.264E+11,  min(M): 1.910E+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 =  2.39474098049503E+16, method 2 =  2.39474098049503E+16, relative difference =  8.35163391903274E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.42533603883824E-26, method 2 = -1.25055521493778E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.425E-26,  chi2_M: 2.395E+16
   max(B_n): 2.154E-14,  max(M): 7.038E+09,  min(M): 1.495E+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 =  4.01135944058716E+15, method 2 =  4.01135944058716E+15, relative difference =  6.23230113638000E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.49348466493383E-26, method 2 = -8.52651282912120E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.493E-26,  chi2_M: 4.011E+15
   max(B_n): 2.431E-14,  max(M): 2.532E+09,  min(M): 8.029E+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 =  3.29179709186741E+15, method 2 =  3.29179709186740E+15, relative difference =  6.07570863022246E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.68439574703636E-26, method 2 =  8.52651282912120E-14, relative difference = -9.99999999998901E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 4.684E-26,  chi2_M: 3.292E+15
   max(B_n): 3.197E-14,  max(M): 1.859E+09,  min(M): 1.476E+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 =  3.25724547078708E+15, method 2 =  3.25724547078708E+15, relative difference =  6.90767711607659E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.39438161219181E-25, method 2 =  1.42108547152020E-13, relative difference = -9.99999999992408E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.394E-25,  chi2_M: 3.257E+15
   max(B_n): 6.861E-14,  max(M): 1.708E+09,  min(M): 1.620E+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 =  3.25611685738753E+15, method 2 =  3.25611685738753E+15, relative difference = -3.83892856045378E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.18981274657678E-23, method 2 =  1.42108547152020E-13, relative difference = -9.99999999832549E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.190E-23,  chi2_M: 3.256E+15
   max(B_n): 2.482E-13,  max(M): 1.681E+09,  min(M): 1.641E+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 =  3.25605932101636E+15, method 2 =  3.25605932101636E+15, relative difference =  7.67799279289434E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.76307095742340E-22, method 2 =  1.70530256582424E-13, relative difference = -9.99999996759436E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.763E-22,  chi2_M: 3.256E+15
   max(B_n): 1.128E-12,  max(M): 1.675E+09,  min(M): 1.640E+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 =  3.25605715654650E+15, method 2 =  3.25605715654651E+15, relative difference = -7.67799789685385E-17
   2 methods of computing chi2_B that should agree: method 1 =  6.45138941489678E-21, method 2 =  1.42108547152020E-13, relative difference = -9.99999909204770E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 6.451E-21,  chi2_M: 3.256E+15
   max(B_n): 5.399E-12,  max(M): 1.674E+09,  min(M): 1.640E+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 =  3.25605704834166E+15, method 2 =  3.25605704834166E+15, relative difference =  4.60679889120482E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.50663205187367E-19, method 2 =  1.70530256582424E-13, relative difference = -9.99998233004804E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.507E-19,  chi2_M: 3.256E+15
   max(B_n): 2.604E-11,  max(M): 1.674E+09,  min(M): 1.639E+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 =  3.25605704353191E+15, method 2 =  3.25605704353191E+15, relative difference =  1.53559963266995E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.51895330057027E-18, method 2 =  1.42108547152020E-13, relative difference = -9.99950476361410E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.519E-18,  chi2_M: 3.256E+15
   max(B_n): 1.258E-10,  max(M): 1.674E+09,  min(M): 1.639E+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 =  3.25605704028324E+15, method 2 =  3.25605704028324E+15, relative difference =  3.07119926840413E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.21923719558422E-17, method 2 =  1.42108547152020E-13, relative difference = -9.98843913854995E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.219E-17,  chi2_M: 3.256E+15
   max(B_n): 6.080E-10,  max(M): 1.674E+09,  min(M): 1.639E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.25605702555996E+15, method 2 =  3.25605702555995E+15, relative difference =  6.91019838515592E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.91978095317643E-15, method 2 =  1.70530256582424E-13, relative difference = -9.77735221393847E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.920E-15,  chi2_M: 3.256E+15
   max(B_n): 2.938E-09,  max(M): 1.674E+09,  min(M): 1.639E+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 =  3.25605695443478E+15, method 2 =  3.25605695443478E+15, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  4.48404954047804E-14, method 2 =  2.84217094304040E-13, relative difference = -7.27461108285274E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 4.484E-14,  chi2_M: 3.256E+15
   max(B_n): 1.420E-08,  max(M): 1.674E+09,  min(M): 1.639E+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 =  3.25605661070116E+15, method 2 =  3.25605661070116E+15, relative difference = -7.67799918399345E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.04732357045927E-12, method 2 =  1.16529008664656E-12, relative difference = -5.33154605678399E-02
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.047E-12,  chi2_M: 3.256E+15
   max(B_n): 6.861E-08,  max(M): 1.674E+09,  min(M): 1.639E+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 =  3.25605494962813E+15, method 2 =  3.25605494962813E+15, relative difference =  5.37460217064170E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.44598531051444E-11, method 2 =  2.45847786572995E-11, relative difference = -2.54718095876665E-03
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.446E-11,  chi2_M: 3.256E+15
   max(B_n): 3.313E-07,  max(M): 1.674E+09,  min(M): 1.639E+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 =  3.25604692484870E+15, method 2 =  3.25604692484869E+15, relative difference =  5.37461541676437E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.71078929069325E-10, method 2 =  5.71191094422829E-10, relative difference = -9.81951300454616E-05
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.711E-10,  chi2_M: 3.256E+15
   max(B_n): 1.595E-06,  max(M): 1.674E+09,  min(M): 1.639E+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 =  3.25600817185951E+15, method 2 =  3.25600817185951E+15, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  1.33283633433063E-08, method 2 =  1.33285311676445E-08, relative difference = -6.29571978628361E-06
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.333E-08,  chi2_M: 3.256E+15
   max(B_n): 7.636E-06,  max(M): 1.674E+09,  min(M): 1.639E+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 =  3.25582097348054E+15, method 2 =  3.25582097348054E+15, relative difference = -6.91069938527579E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.11172900140380E-07, method 2 =  3.11173010913990E-07, relative difference = -1.77993633478913E-07
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.112E-07,  chi2_M: 3.256E+15
   max(B_n): 3.655E-05,  max(M): 1.674E+09,  min(M): 1.639E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.25491653007679E+15, method 2 =  3.25491653007680E+15, relative difference = -3.84034425598775E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.26575035863698E-06, method 2 =  7.26575055409739E-06, relative difference = -1.34508071058805E-08
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.266E-06,  chi2_M: 3.255E+15
   max(B_n): 1.760E-04,  max(M): 1.674E+09,  min(M): 1.639E+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 =  3.25055073289639E+15, method 2 =  3.25055073289638E+15, relative difference =  5.38370308234098E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.69479515039218E-04, method 2 =  1.69479515193416E-04, relative difference = -4.54916888675246E-10
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 1.695E-04,  chi2_M: 3.251E+15
   max(B_n): 8.509E-04,  max(M): 1.672E+09,  min(M): 1.638E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.22957447144137E+15, method 2 =  3.22957447144136E+15, relative difference =  1.00632452626167E-15
   2 methods of computing chi2_B that should agree: method 1 =  3.93301751489513E-03, method 2 =  3.93301751512354E-03, relative difference = -2.90373024179261E-11
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.933E-03,  chi2_M: 3.230E+15
   max(B_n): 4.099E-03,  max(M): 1.667E+09,  min(M): 1.633E+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 =  3.13100477499544E+15, method 2 =  3.13100477499544E+15, relative difference = -4.79079435451255E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.90599271980210E-02, method 2 =  8.90599271982353E-02, relative difference = -1.20336070012080E-12
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.906E-02,  chi2_M: 3.131E+15
   max(B_n): 1.950E-02,  max(M): 1.642E+09,  min(M): 1.607E+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 =  2.71277770873263E+15, method 2 =  2.71277770873263E+15, relative difference = -1.01372109891185E-15
   2 methods of computing chi2_B that should agree: method 1 =  1.80228606653281E+00, method 2 =  1.80228606653301E+00, relative difference = -5.75967681990381E-14
   Diagnostics:    3.00000003E-03  sec.
   chi2_B: 1.802E+00,  chi2_M: 2.713E+15
   max(B_n): 8.760E-02,  max(M): 1.529E+09,  min(M): 1.495E+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.52364646619917E+15, method 2 =  1.52364646619916E+15, relative difference =  6.56320230568030E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.36428043749313E+01, method 2 =  2.36428043749314E+01, relative difference = -2.17885947562246E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.364E+01,  chi2_M: 1.524E+15
   max(B_n): 3.164E-01,  max(M): 1.147E+09,  min(M): 1.119E+09
 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.11697174707424E+14, method 2 =  3.11697174707424E+14, relative difference =  9.02318092116159E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.12969108545900E+02, method 2 =  1.12969108545900E+02, relative difference = -7.54764991852340E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.130E+02,  chi2_M: 3.117E+14
   max(B_n): 6.907E-01,  max(M): 5.191E+08,  min(M): 5.061E+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.34348933694660E+13, method 2 =  2.34348933694661E+13, relative difference = -5.83398216687149E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.98384049388001E+02, method 2 =  1.98384049388000E+02, relative difference =  5.73064407508220E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.984E+02,  chi2_M: 2.343E+13
   max(B_n): 9.153E-01,  max(M): 1.422E+08,  min(M): 1.389E+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.15330627372223E+12, method 2 =  1.15330627372223E+12, relative difference =  4.23375178931526E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.28038333382817E+02, method 2 =  2.28038333382817E+02, relative difference = -3.11589163637363E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.280E+02,  chi2_M: 1.153E+12
   max(B_n): 9.812E-01,  max(M): 3.154E+07,  min(M): 3.081E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.833E-11 ( 36 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.08846274540253E+10, method 2 =  5.08846274540252E+10, relative difference =  4.49805466580839E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.35001392841955E+02, method 2 =  2.35001392841955E+02, relative difference = -6.04713637793594E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.350E+02,  chi2_M: 5.088E+10
   max(B_n): 9.961E-01,  max(M): 6.625E+06,  min(M): 6.473E+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 =  2.19226561571323E+09, method 2 =  2.19226561571324E+09, relative difference = -4.35017686529732E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.36481935072821E+02, method 2 =  2.36481935072821E+02, relative difference =  6.00927707684140E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.365E+02,  chi2_M: 2.192E+09
   max(B_n): 9.992E-01,  max(M): 1.375E+06,  min(M): 1.344E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.281E-10 ( 38 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.39804617367892E+07, method 2 =  9.39804617367892E+07, relative difference =  3.17111895780663E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.36790026469713E+02, method 2 =  2.36790026469712E+02, relative difference =  4.20102081533987E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.368E+02,  chi2_M: 9.398E+07
   max(B_n): 9.998E-01,  max(M): 2.847E+05,  min(M): 2.782E+05
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.069E-09 ( 39 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.02470050476737E+06, method 2 =  4.02470050476737E+06, relative difference =  5.20653845025904E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.36853850011179E+02, method 2 =  2.36853850011179E+02, relative difference =  1.19996822635826E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.369E+02,  chi2_M: 4.025E+06
   max(B_n): 1.000E+00,  max(M): 5.892E+04,  min(M): 5.757E+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.72320416656598E+05, method 2 =  1.72320416656598E+05, relative difference = -2.53340529997083E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.36867059205275E+02, method 2 =  2.36867059205275E+02, relative difference = -4.19965458009110E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.369E+02,  chi2_M: 1.723E+05
   max(B_n): 1.000E+00,  max(M): 1.219E+04,  min(M): 1.191E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   1.6579999923706055      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
