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
 No ports specified.
 We will use the following values of the regularization weight lambda:
 0.000E+00 1.000E-34 4.833E-34 2.336E-33 1.129E-32 5.456E-32 2.637E-31 1.274E-30 6.158E-30 2.976E-29 1.438E-28 6.952E-28 3.360E-27 1.624E-26 7.848E-26 3.793E-25 1.833E-24 8.859E-24 4.281E-23 2.069E-22 1.000E-21 4.833E-21 2.336E-20 1.129E-19 5.456E-19 2.637E-18 1.274E-17 6.158E-17 2.976E-16 1.438E-15 6.952E-15 3.360E-14 1.624E-13 7.848E-13 3.793E-12 1.833E-11 8.859E-11 4.281E-10 2.069E-09 1.000E-08
 Initializing plasma surface.
 Reading FOCUS format data from file:plasma
 Number of modes for Bnormal read from FOCUS file:           1
 Plasma surface area: 3.948E+02 m^2. Volume: 6.575E+01 m^3.
 G =    0.0000000000000000       ; curpol =    1.0000000000000000     
 Done initializing plasma surface. Took    1.79999992E-02  sec.
 Initializing coil surface.
   Building a plain circular torus.
   Evaluating coil surface & derivatives:   8.99999961E-03  sec.
 s_integration:  0.21132486540518713       0.78867513459481287     
 s_weights:  0.50000000000000011       0.50000000000000011     
 s_magnetization  0.50000000000000000     
 s_magnetization_weights   1.0000000000000000     
 interpolate_magnetization_to_integration:
   1.0000000000000000     
   1.0000000000000000     
 Maximum allowed uniform d before singularity:  2.9000E+01
 Coil surface area: 1.184E+03 m^2. Volume: 5.918E+02 m^3.
 Done initializing coil surface. Took    1.20000001E-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.58399999     sec.
 Computing g and RHS_B:
  js=   1 j_RZetaZ= 1 dgemm took  3.000E-03 sec, rest took  1.000E-03 sec.
  js=   1 j_RZetaZ= 2 dgemm took  1.000E-03 sec, rest took  0.000E+00 sec.
  js=   1 j_RZetaZ= 3 dgemm took  0.000E+00 sec, rest took  0.000E+00 sec.
 Total time to compute g and RHS_B:   8.00000038E-03 sec.
 About to form matrix_B.
  js=   1, ks=   1, j_RZetaZ= 1, k_RZetaZ= 1 dgemm took  1.000E-03 sec.
  js=   1, ks=   1, j_RZetaZ= 2, k_RZetaZ= 1 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 2, k_RZetaZ= 2 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 1 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 2 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 3 dgemm took  1.000E-03 sec.
 Total time for matrix_B:   2.00000009E-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 =  3.61486051948918E+32, method 2 =  3.61486051948918E+32, relative difference =  5.98011405829654E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.45760207680438E-14, method 2 =  7.72502062318381E-11, relative difference = -9.98071096520533E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.458E-14,  chi2_M: 3.615E+32
   max(B_n): 3.805E-08,  max(M): 1.319E+18,  min(M): 1.553E+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.27281399294976E+22, method 2 =  6.27281399294976E+22, relative difference =  2.00594374616279E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.60699683870860E-21, method 2 = -8.98410235095071E-11, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 5.607E-21,  chi2_M: 6.273E+22
   max(B_n): 1.079E-11,  max(M): 1.122E+13,  min(M): 8.911E+11
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.833E-34 (  3 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  7.83213730702725E+20, method 2 =  7.83213730702725E+20, relative difference = -3.34703018759382E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.13191417181755E-24, method 2 = -4.83169060316868E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 6.132E-24,  chi2_M: 7.832E+20
   max(B_n): 3.182E-13,  max(M): 1.214E+12,  min(M): 4.474E+09
 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 =  3.35974508739356E+19, method 2 =  3.35974508739356E+19, relative difference = -3.65742033409233E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.08914697468425E-24, method 2 = -6.82121026329696E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.089E-24,  chi2_M: 3.360E+19
   max(B_n): 1.588E-13,  max(M): 2.827E+11,  min(M): 1.768E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.129E-32 (  5 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  7.03836994290807E+18, method 2 =  7.03836994290806E+18, relative difference =  4.36464696359898E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.17059228916476E-25, method 2 = -9.37916411203332E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.171E-25,  chi2_M: 7.038E+18
   max(B_n): 7.694E-14,  max(M): 1.437E+11,  min(M): 9.377E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.456E-32 (  6 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.27432259143524E+18, method 2 =  1.27432259143524E+18, relative difference =  2.00891047306689E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.59416527355986E-26, method 2 = -6.25277607468888E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.594E-26,  chi2_M: 1.274E+18
   max(B_n): 3.264E-14,  max(M): 5.091E+10,  min(M): 8.053E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.637E-31 (  7 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.25810766825875E+21, method 2 =  1.25810766825875E+21, relative difference =  9.37636761750810E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.48030178009578E-21, method 2 =  5.57420776203799E-09, relative difference = -9.99999999999469E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.480E-21,  chi2_M: 1.258E+21
   max(B_n): 4.307E-12,  max(M): 2.099E+12,  min(M): 1.100E+11
 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.56631348173135E+17, method 2 =  1.56631348173136E+17, relative difference = -9.19356193249558E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.25035043921953E-26, method 2 = -7.95807864051312E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.250E-26,  chi2_M: 1.566E+17
   max(B_n): 2.476E-14,  max(M): 1.804E+10,  min(M): 3.324E+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 =  4.60282365262762E+16, method 2 =  4.60282365262763E+16, relative difference = -6.08322241153331E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.06770531587611E-26, method 2 = -3.97903932025656E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.068E-26,  chi2_M: 4.603E+16
   max(B_n): 2.409E-14,  max(M): 8.671E+09,  min(M): 2.718E+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 =  7.58485599395220E+15, method 2 =  7.58485599395221E+15, relative difference = -6.59208296635659E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.02370827466972E-26, method 2 = -5.68434188608080E-14, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.024E-26,  chi2_M: 7.585E+15
   max(B_n): 2.420E-14,  max(M): 3.400E+09,  min(M): 1.379E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.18602050911634E+17, method 2 =  4.18602050911635E+17, relative difference = -3.82224596490989E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.12856329938558E-26, method 2 = -5.45981038158061E-11, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 8.129E-26,  chi2_M: 4.186E+17
   max(B_n): 4.008E-14,  max(M): 2.179E+10,  min(M): 1.559E+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 =  6.76615158941064E+15, method 2 =  6.76615158941065E+15, relative difference = -5.17280754613549E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.35503244563758E-25, method 2 =  5.68434188608080E-14, relative difference = -9.99999999995232E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.355E-25,  chi2_M: 6.766E+15
   max(B_n): 4.741E-14,  max(M): 2.457E+09,  min(M): 2.332E+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 =  6.76404707637133E+15, method 2 =  6.76404707637132E+15, relative difference =  1.10880363713014E-15
   2 methods of computing chi2_B that should agree: method 1 =  2.67716377056740E-24, method 2 =  5.68434188608080E-14, relative difference = -9.99999999905806E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.677E-24,  chi2_M: 6.764E+15
   max(B_n): 1.444E-13,  max(M): 2.423E+09,  min(M): 2.366E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.624E-26 ( 14 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396529973921E+15, method 2 =  6.76396529973920E+15, relative difference =  5.17447953219829E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.18459477475784E-23, method 2 =  1.70530256582424E-13, relative difference = -9.99999999274663E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 6.185E-23,  chi2_M: 6.764E+15
   max(B_n): 6.136E-13,  max(M): 2.416E+09,  min(M): 2.364E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.848E-26 ( 15 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396293090802E+15, method 2 =  6.76396293090802E+15, relative difference =  1.47842324124883E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.44372915747353E-21, method 2 =  5.68434188608080E-14, relative difference = -9.99999949203297E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.444E-21,  chi2_M: 6.764E+15
   max(B_n): 2.909E-12,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.793E-25 ( 16 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396282229215E+15, method 2 =  6.76396282229216E+15, relative difference = -3.69605816247347E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.37151707219825E-20, method 2 =  5.68434188608080E-14, relative difference = -9.99998813753592E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.372E-20,  chi2_M: 6.764E+15
   max(B_n): 1.403E-11,  max(M): 2.415E+09,  min(M): 2.363E+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 =  6.76396281087117E+15, method 2 =  6.76396281087117E+15, relative difference =  5.91369306994285E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.87424768733811E-19, method 2 =  1.13686837721616E-13, relative difference = -9.99986147572918E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 7.874E-19,  chi2_M: 6.764E+15
   max(B_n): 6.778E-11,  max(M): 2.415E+09,  min(M): 2.363E+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 =  6.76396280734398E+15, method 2 =  6.76396280734398E+15, relative difference = -7.39211634128331E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.83917997893761E-17, method 2 =  5.68434188608080E-14, relative difference = -9.99353105376853E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.839E-17,  chi2_M: 6.764E+15
   max(B_n): 3.276E-10,  max(M): 2.415E+09,  min(M): 2.363E+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 =  6.76396279141703E+15, method 2 =  6.76396279141702E+15, relative difference =  1.10881745380341E-15
   2 methods of computing chi2_B that should agree: method 1 =  4.29574878232154E-16, method 2 =  2.84217094304040E-14, relative difference = -9.70221438072788E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 4.296E-16,  chi2_M: 6.764E+15
   max(B_n): 1.583E-09,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.069E-22 ( 20 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396271450416E+15, method 2 =  6.76396271450416E+15, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  1.00330253107407E-14, method 2 =  1.13686837721616E-13, relative difference = -8.37810597832351E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.003E-14,  chi2_M: 6.764E+15
   max(B_n): 7.649E-09,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-21 ( 21 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396234289406E+15, method 2 =  6.76396234289406E+15, relative difference = -5.17448179420596E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.34272684018534E-13, method 2 =  2.55795384873636E-13, relative difference = -4.39177784093444E-02
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.343E-13,  chi2_M: 6.764E+15
   max(B_n): 3.693E-08,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.833E-21 ( 22 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396054906167E+15, method 2 =  6.76396054906167E+15, relative difference =  1.47842376185758E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.46486248000183E-12, method 2 =  5.62749846721999E-12, relative difference = -1.46619811591057E-02
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.465E-12,  chi2_M: 6.764E+15
   max(B_n): 1.777E-07,  max(M): 2.415E+09,  min(M): 2.363E+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 =  6.76395191121606E+15, method 2 =  6.76395191121605E+15, relative difference =  7.39212824932854E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.27154579151708E-10, method 2 =  1.27187149701058E-10, relative difference = -1.28058221107278E-04
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.272E-10,  chi2_M: 6.764E+15
   max(B_n): 8.460E-07,  max(M): 2.415E+09,  min(M): 2.363E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76391035067710E+15, method 2 =  6.76391035067709E+15, relative difference =  4.43530420195426E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.95739042722131E-09, method 2 =  2.95744939649012E-09, relative difference = -9.96971525242530E-06
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.957E-09,  chi2_M: 6.764E+15
   max(B_n): 3.981E-06,  max(M): 2.415E+09,  min(M): 2.363E+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 =  6.76370974668116E+15, method 2 =  6.76370974668116E+15, relative difference =  1.47847858269005E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.89985272227303E-08, method 2 =  6.89985881763278E-08, relative difference = -4.41701968312857E-07
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.900E-08,  chi2_M: 6.764E+15
   max(B_n): 1.891E-05,  max(M): 2.415E+09,  min(M): 2.363E+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 =  6.76274043443493E+15, method 2 =  6.76274043443494E+15, relative difference = -4.43607148475553E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.61123917793163E-06, method 2 =  1.61123927000517E-06, relative difference = -2.85722755831801E-08
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.611E-06,  chi2_M: 6.763E+15
   max(B_n): 9.104E-05,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.274E-17 ( 27 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.75805875905638E+15, method 2 =  6.75805875905636E+15, relative difference =  1.25775763470675E-15
   2 methods of computing chi2_B that should agree: method 1 =  3.76078215441162E-05, method 2 =  3.76078216106634E-05, relative difference = -8.84752273230293E-10
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.761E-05,  chi2_M: 6.758E+15
   max(B_n): 4.404E-04,  max(M): 2.414E+09,  min(M): 2.362E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.158E-17 ( 28 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.73550097214234E+15, method 2 =  6.73550097214234E+15, relative difference =  2.96934111994325E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.75479806290646E-04, method 2 =  8.75479806353496E-04, relative difference = -3.58949786071805E-11
   Diagnostics:    0.00000000      sec.
   chi2_B: 8.755E-04,  chi2_M: 6.736E+15
   max(B_n): 2.125E-03,  max(M): 2.410E+09,  min(M): 2.358E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.976E-16 ( 29 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.62805758356087E+15, method 2 =  6.62805758356086E+15, relative difference =  2.26310646986585E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.01223758388828E-02, method 2 =  2.01223758389801E-02, relative difference = -2.41901558372314E-12
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.012E-02,  chi2_M: 6.628E+15
   max(B_n): 1.018E-02,  max(M): 2.391E+09,  min(M): 2.339E+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 =  6.14335280104900E+15, method 2 =  6.14335280104901E+15, relative difference = -1.13944294373013E-15
   2 methods of computing chi2_B that should agree: method 1 =  4.35616184540566E-01, method 2 =  4.35616184540550E-01, relative difference =  1.89872686713442E-14
   Diagnostics:    0.00000000      sec.
   chi2_B: 4.356E-01,  chi2_M: 6.143E+15
   max(B_n): 4.731E-02,  max(M): 2.303E+09,  min(M): 2.250E+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 =  4.41162927503261E+15, method 2 =  4.41162927503261E+15, relative difference = -3.40010437524561E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.30618648915259E+00, method 2 =  7.30618648915265E+00, relative difference = -4.19400127893645E-15
   Diagnostics:    0.00000000      sec.
   chi2_B: 7.306E+00,  chi2_M: 4.412E+15
   max(B_n): 1.930E-01,  max(M): 1.955E+09,  min(M): 1.904E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.360E-14 ( 32 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.46150008573861E+15, method 2 =  1.46150008573860E+15, relative difference =  5.13171369142253E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.65321565588320E+01, method 2 =  5.65321565588320E+01, relative difference =  3.77064725111265E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 5.653E+01,  chi2_M: 1.462E+15
   max(B_n): 5.353E-01,  max(M): 1.126E+09,  min(M): 1.095E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.624E-13 ( 33 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.56983574601765E+14, method 2 =  1.56983574601765E+14, relative difference = -6.96728942995864E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.41829307622223E+02, method 2 =  1.41829307622223E+02, relative difference = -4.00787536890586E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.418E+02,  chi2_M: 1.570E+14
   max(B_n): 8.480E-01,  max(M): 3.688E+08,  min(M): 3.592E+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 =  8.69531229175062E+12, method 2 =  8.69531229175061E+12, relative difference =  2.24618154525969E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.83490844582143E+02, method 2 =  1.83490844582143E+02, relative difference =  2.32341642127658E-16
   Diagnostics:    3.00000003E-03  sec.
   chi2_B: 1.835E+02,  chi2_M: 8.695E+12
   max(B_n): 9.643E-01,  max(M): 8.672E+07,  min(M): 8.460E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.793E-12 ( 35 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.94387149389084E+11, method 2 =  3.94387149389084E+11, relative difference =  7.73797477232018E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.94389000944839E+02, method 2 =  1.94389000944839E+02, relative difference =  2.92420963040693E-16
   Diagnostics:    3.00000003E-03  sec.
   chi2_B: 1.944E+02,  chi2_M: 3.944E+11
   max(B_n): 9.924E-01,  max(M): 1.846E+07,  min(M): 1.802E+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 =  1.70914238318475E+10, method 2 =  1.70914238318475E+10, relative difference = -5.57984124546258E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.96765024188282E+02, method 2 =  1.96765024188282E+02, relative difference =  7.22224631833137E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.968E+02,  chi2_M: 1.709E+10
   max(B_n): 9.984E-01,  max(M): 3.843E+06,  min(M): 3.752E+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 =  7.33589315091546E+08, method 2 =  7.33589315091547E+08, relative difference = -6.50005593584224E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97262094128866E+02, method 2 =  1.97262094128866E+02, relative difference =  5.04283316294051E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.973E+02,  chi2_M: 7.336E+08
   max(B_n): 9.997E-01,  max(M): 7.962E+05,  min(M): 7.774E+05
 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 =  3.14238110463200E+07, method 2 =  3.14238110463200E+07, relative difference = -4.14924721428518E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97365179948588E+02, method 2 =  1.97365179948588E+02, relative difference =  2.88011385167411E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 3.142E+07
   max(B_n): 9.999E-01,  max(M): 1.648E+05,  min(M): 1.609E+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 =  1.34550327844362E+06, method 2 =  1.34550327844362E+06, relative difference =  3.46087069996874E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97386519918373E+02, method 2 =  1.97386519918374E+02, relative difference = -2.15985185630894E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 1.346E+06
   max(B_n): 1.000E+00,  max(M): 3.410E+04,  min(M): 3.329E+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 =  5.76067566285368E+04, method 2 =  5.76067566285368E+04, relative difference =  6.31519464036196E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.97390935884973E+02, method 2 =  1.97390935884973E+02, relative difference =  4.31960707359426E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 5.761E+04
   max(B_n): 1.000E+00,  max(M): 7.055E+03,  min(M): 6.889E+03
 Volume of magnetization region:  1.18E+00 meters^3.
 Evaluating the outer magnetization surface in terms of the standard toroidal angle.
   Computing offset points:   1.00000005E-03 sec
 Done. Took    1.00000005E-03 sec
 REGCOIL_PM complete. Total time=   1.6779999732971191      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
