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: 8.883E+02 m^2. Volume: 3.329E+02 m^3.
 G =    0.0000000000000000       ; curpol =    1.0000000000000000     
 Done initializing plasma surface. Took    4.30000015E-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.09999997E-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   3.30200005     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  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:   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  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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.04203675010815E+22, method 2 =  3.04203675010815E+22, relative difference =  8.96207976416812E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.67840868298371E-18, method 2 = -3.97903932025656E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.678E-18,  chi2_M: 3.042E+22
   max(B_n): 7.713E-11,  max(M): 1.371E+13,  min(M): 1.352E+11
 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 =  2.09244022941906E+19, method 2 =  2.09244022941905E+19, relative difference =  5.87256917891110E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.88535055510507E-18, method 2 = -1.64845914696343E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.885E-18,  chi2_M: 2.092E+19
   max(B_n): 7.741E-11,  max(M): 2.103E+11,  min(M): 2.155E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.09243888288470E+19, method 2 =  2.09243888288470E+19, relative difference =  4.89381079837459E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.88545872874595E-18, method 2 = -1.36424205265939E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.885E-18,  chi2_M: 2.092E+19
   max(B_n): 7.742E-11,  max(M): 2.103E+11,  min(M): 2.155E+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 =  1.63818522118639E+19, method 2 =  1.63818522118638E+19, relative difference =  3.12540971178586E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.59928919316333E-18, method 2 = -1.70530256582424E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.599E-18,  chi2_M: 1.638E+19
   max(B_n): 6.949E-11,  max(M): 1.877E+11,  min(M): 1.470E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  7.84346043913656E+19, method 2 =  7.84346043913656E+19, relative difference =  1.04443696294104E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.62881812283987E-18, method 2 = -2.10320649784990E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.629E-18,  chi2_M: 7.843E+19
   max(B_n): 1.108E-10,  max(M): 4.076E+11,  min(M): 1.053E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.07835093289645E+17, method 2 =  3.07835093289646E+17, relative difference = -7.27662325975408E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.45981210983950E-20, method 2 = -7.38964445190504E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 7.460E-20,  chi2_M: 3.078E+17
   max(B_n): 1.389E-11,  max(M): 2.720E+10,  min(M): 1.803E+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.18414588809589E+17, method 2 =  3.18414588809590E+17, relative difference = -2.00995815673106E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.42462020473911E-19, method 2 =  3.01270119962282E-12, relative difference = -9.99999440726625E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 8.425E-19,  chi2_M: 3.184E+17
   max(B_n): 5.648E-11,  max(M): 2.869E+10,  min(M): 6.749E+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 =  1.70222696157850E+16, method 2 =  1.70222696157849E+16, relative difference =  9.98691736396638E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.11178449031264E-20, method 2 = -3.97903932025656E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.112E-20,  chi2_M: 1.702E+16
   max(B_n): 1.445E-11,  max(M): 6.904E+09,  min(M): 8.924E+07
 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 =  4.36985459633266E+15, method 2 =  4.36985459633265E+15, relative difference =  8.00941981670816E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.53397020931685E-20, method 2 = -1.13686837721616E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.534E-20,  chi2_M: 4.370E+15
   max(B_n): 1.259E-11,  max(M): 3.849E+09,  min(M): 2.372E+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 =  1.91957294845986E+15, method 2 =  1.91957294845985E+15, relative difference =  4.55830553718756E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.20894846555787E-20, method 2 = -6.82121026329696E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.209E-20,  chi2_M: 1.920E+15
   max(B_n): 1.128E-11,  max(M): 2.200E+09,  min(M): 1.097E+08
 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 =  3.53947463235284E+15, method 2 =  3.53947463235284E+15, relative difference =  2.82527805358294E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.73212878965153E-20, method 2 = -9.66338120633736E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.732E-20,  chi2_M: 3.539E+15
   max(B_n): 9.565E-12,  max(M): 2.592E+09,  min(M): 4.055E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.26163348756280E+14, method 2 =  3.26163348756280E+14, relative difference =  2.87432663288152E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.28385743049559E-20, method 2 = -6.82121026329696E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.284E-20,  chi2_M: 3.262E+14
   max(B_n): 7.507E-12,  max(M): 7.456E+08,  min(M): 2.015E+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 =  2.65344829484907E+14, method 2 =  2.65344829484906E+14, relative difference =  5.29970756441663E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.32896938636319E-21, method 2 = -4.54747350886464E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 7.329E-21,  chi2_M: 2.653E+14
   max(B_n): 4.252E-12,  max(M): 5.364E+08,  min(M): 4.079E+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 =  2.61944990229166E+14, method 2 =  2.61944990229165E+14, relative difference =  8.35099002308172E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.55721301829987E-22, method 2 = -6.25277607468888E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 7.557E-22,  chi2_M: 2.619E+14
   max(B_n): 1.375E-12,  max(M): 4.843E+08,  min(M): 4.557E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.61805190499930E+14, method 2 =  2.61805190499930E+14, relative difference = -6.56499589147931E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.21159368845457E-23, method 2 = -4.54747350886464E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 4.212E-23,  chi2_M: 2.618E+14
   max(B_n): 3.572E-13,  max(M): 4.744E+08,  min(M): 4.653E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.61798615113288E+14, method 2 =  2.61798615113288E+14, relative difference =  5.96832798112342E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.39261614498305E-23, method 2 = -5.68434188608080E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.393E-23,  chi2_M: 2.618E+14
   max(B_n): 2.923E-13,  max(M): 4.725E+08,  min(M): 4.661E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.61798441125247E+14, method 2 =  2.61798441125247E+14, relative difference = -5.37149875284106E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.17884927531869E-22, method 2 = -6.82121026329696E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.179E-22,  chi2_M: 2.618E+14
   max(B_n): 1.109E-12,  max(M): 4.725E+08,  min(M): 4.659E+08
 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 =  2.61798428480361E+14, method 2 =  2.61798428480361E+14, relative difference = -5.96833223587213E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.21114946626049E-20, method 2 = -6.82121026329696E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.211E-20,  chi2_M: 2.618E+14
   max(B_n): 5.311E-12,  max(M): 4.725E+08,  min(M): 4.658E+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 =  2.61798428201631E+14, method 2 =  2.61798428201631E+14, relative difference =  6.56516546644909E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.82978938277241E-19, method 2 = -6.25277607468888E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.830E-19,  chi2_M: 2.618E+14
   max(B_n): 2.566E-11,  max(M): 4.725E+08,  min(M): 4.658E+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 =  2.61798428137694E+14, method 2 =  2.61798428137693E+14, relative difference =  4.17783257057884E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.60986425442543E-18, method 2 = -5.68434188608080E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.610E-18,  chi2_M: 2.618E+14
   max(B_n): 1.240E-10,  max(M): 4.725E+08,  min(M): 4.658E+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 =  2.61798427891490E+14, method 2 =  2.61798427891490E+14, relative difference =  1.19366644985937E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.54389921643879E-16, method 2 = -5.68434188608080E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.544E-16,  chi2_M: 2.618E+14
   max(B_n): 5.993E-10,  max(M): 4.725E+08,  min(M): 4.658E+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 =  2.61798426707909E+14, method 2 =  2.61798426707909E+14, relative difference =  1.79049968288384E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.60612623532239E-15, method 2 = -5.11590769747272E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.606E-15,  chi2_M: 2.618E+14
   max(B_n): 2.897E-09,  max(M): 4.725E+08,  min(M): 4.658E+08
 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 =  2.61798420987969E+14, method 2 =  2.61798420987968E+14, relative difference =  7.75883212868327E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.42290782799063E-14, method 2 = -3.97903932025656E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.423E-14,  chi2_M: 2.618E+14
   max(B_n): 1.400E-08,  max(M): 4.725E+08,  min(M): 4.658E+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 =  2.61798393343930E+14, method 2 =  2.61798393343930E+14, relative difference = -9.54933285902827E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.96735445692464E-12, method 2 =  1.42108547152020E-12, relative difference =  1.61215484689191E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.967E-12,  chi2_M: 2.618E+14
   max(B_n): 6.765E-08,  max(M): 4.725E+08,  min(M): 4.658E+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 =  2.61798259742935E+14, method 2 =  2.61798259742935E+14, relative difference = -2.98416804132741E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.59516495614199E-11, method 2 =  4.53610482509248E-11, relative difference =  6.46789903972482E-03
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.595E-11,  chi2_M: 2.618E+14
   max(B_n): 3.269E-07,  max(M): 4.725E+08,  min(M): 4.658E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.61797614074636E+14, method 2 =  2.61797614074637E+14, relative difference = -7.75885604297717E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.07327332536698E-09, method 2 =  1.07269215732231E-09, relative difference =  2.70818915475025E-04
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.073E-09,  chi2_M: 2.618E+14
   max(B_n): 1.579E-06,  max(M): 4.725E+08,  min(M): 4.658E+08
 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 =  2.61794493912086E+14, method 2 =  2.61794493912086E+14, relative difference =  6.56526412880624E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.50660051419769E-08, method 2 =  2.50653329203487E-08, relative difference =  1.34092097712550E-05
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.507E-08,  chi2_M: 2.618E+14
   max(B_n): 7.620E-06,  max(M): 4.725E+08,  min(M): 4.658E+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 =  2.61779418050177E+14, method 2 =  2.61779418050176E+14, relative difference =  7.75939535326899E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.85325133369170E-07, method 2 =  5.85324528401543E-07, relative difference =  5.16779397778062E-07
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.853E-07,  chi2_M: 2.618E+14
   max(B_n): 3.666E-05,  max(M): 4.725E+08,  min(M): 4.658E+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 =  2.61706584794954E+14, method 2 =  2.61706584794954E+14, relative difference = -1.19408535419482E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.36662392404912E-05, method 2 =  1.36662386012176E-05, relative difference =  2.33887895316947E-08
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.367E-05,  chi2_M: 2.617E+14
   max(B_n): 1.762E-04,  max(M): 4.724E+08,  min(M): 4.657E+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 =  2.61355020515592E+14, method 2 =  2.61355020515592E+14, relative difference =  7.17414953920175E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.18771529006690E-04, method 2 =  3.18771528384332E-04, relative difference =  9.76182782534934E-10
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.188E-04,  chi2_M: 2.614E+14
   max(B_n): 8.493E-04,  max(M): 4.721E+08,  min(M): 4.654E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.59665880086091E+14, method 2 =  2.59665880086091E+14, relative difference = -1.80520444135590E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.39748645029543E-03, method 2 =  7.39748644957672E-03, relative difference =  4.85783847176479E-11
   Diagnostics:    0.00000000      sec.
   chi2_B: 7.397E-03,  chi2_M: 2.597E+14
   max(B_n): 4.094E-03,  max(M): 4.706E+08,  min(M): 4.639E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.51728720883608E+14, method 2 =  2.51728720883608E+14, relative difference =  4.96566301855546E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.67503135079714E-01, method 2 =  1.67503135079130E-01, relative difference =  1.74392888484254E-12
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.675E-01,  chi2_M: 2.517E+14
   max(B_n): 1.948E-02,  max(M): 4.633E+08,  min(M): 4.567E+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.18057502800435E+14, method 2 =  2.18057502800435E+14, relative difference =  7.16554110697117E-17
   2 methods of computing chi2_B that should agree: method 1 =  3.38907826175672E+00, method 2 =  3.38907826175614E+00, relative difference =  8.62212901303434E-14
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.389E+00,  chi2_M: 2.181E+14
   max(B_n): 8.761E-02,  max(M): 4.311E+08,  min(M): 4.250E+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.22383700524896E+14, method 2 =  1.22383700524896E+14, relative difference =  7.02197266722771E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.44276363153964E+01, method 2 =  4.44276363153960E+01, relative difference =  4.23821388200151E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.443E+01,  chi2_M: 1.224E+14
   max(B_n): 3.169E-01,  max(M): 3.230E+08,  min(M): 3.182E+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 =  2.50061337835792E+13, method 2 =  2.50061337835792E+13, relative difference =  4.68635019768444E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.12028671992821E+02, method 2 =  2.12028671992821E+02, relative difference =  6.70232689835608E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.120E+02,  chi2_M: 2.501E+13
   max(B_n): 6.916E-01,  max(M): 1.460E+08,  min(M): 1.437E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.87870885163700E+12, method 2 =  1.87870885163699E+12, relative difference =  1.29951282652057E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.72072250443695E+02, method 2 =  3.72072250443695E+02, relative difference = -3.81938042900421E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.721E+02,  chi2_M: 1.879E+12
   max(B_n): 9.155E-01,  max(M): 4.005E+07,  min(M): 3.936E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.24372715026728E+10, method 2 =  9.24372715026729E+10, relative difference = -4.12679561351473E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.27597957478802E+02, method 2 =  4.27597957478802E+02, relative difference =  1.32936600530011E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.276E+02,  chi2_M: 9.244E+10
   max(B_n): 9.813E-01,  max(M): 8.886E+06,  min(M): 8.729E+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 =  4.07819318439124E+09, method 2 =  4.07819318439124E+09, relative difference = -1.75385447663989E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.40633259318579E+02, method 2 =  4.40633259318579E+02, relative difference =  1.29003922556173E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.406E+02,  chi2_M: 4.078E+09
   max(B_n): 9.961E-01,  max(M): 1.867E+06,  min(M): 1.833E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.069E-09 ( 39 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.75699239218633E+08, method 2 =  1.75699239218633E+08, relative difference =  3.39242475041232E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.43404808003855E+02, method 2 =  4.43404808003855E+02, relative difference =  5.12790279534486E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.434E+02,  chi2_M: 1.757E+08
   max(B_n): 9.992E-01,  max(M): 3.874E+05,  min(M): 3.806E+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 =  7.53205271805080E+06, method 2 =  7.53205271805080E+06, relative difference = -3.70943729210844E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.43981544216617E+02, method 2 =  4.43981544216617E+02, relative difference =  0.00000000000000E+00
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 4.440E+02,  chi2_M: 7.532E+06
   max(B_n): 9.998E-01,  max(M): 8.022E+04,  min(M): 7.879E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   3.4230000972747803      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
