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-26 1.600E-26 2.560E-26 4.095E-26 6.551E-26 1.048E-25 1.677E-25 2.683E-25 4.292E-25 6.866E-25 1.099E-24 1.758E-24 2.812E-24 4.498E-24 7.197E-24 1.151E-23 1.842E-23 2.947E-23 4.715E-23 7.543E-23 1.207E-22 1.931E-22 3.089E-22 4.942E-22 7.906E-22 1.265E-21 2.024E-21 3.237E-21 5.179E-21 8.286E-21 1.326E-20 2.121E-20 3.393E-20 5.429E-20 8.685E-20 1.389E-19 2.223E-19 3.556E-19 5.690E-19 9.103E-19 1.456E-18 2.330E-18 3.728E-18 5.964E-18 9.541E-18 1.526E-17 2.442E-17 3.907E-17 6.251E-17 1.000E-16 1.600E-16 2.560E-16 4.095E-16 6.551E-16 1.048E-15 1.677E-15 2.683E-15 4.292E-15 6.866E-15 1.099E-14 1.758E-14 2.812E-14 4.498E-14 7.197E-14 1.151E-13 1.842E-13 2.947E-13 4.715E-13 7.543E-13 1.207E-12 1.931E-12 3.089E-12 4.942E-12 7.906E-12 1.265E-11 2.024E-11 3.237E-11 5.179E-11 8.286E-11 1.326E-10 2.121E-10 3.393E-10 5.429E-10 8.685E-10 1.389E-09 2.223E-09 3.556E-09 5.690E-09 9.103E-09 1.456E-08 2.330E-08 3.728E-08 5.964E-08 9.541E-08 1.526E-07 2.442E-07 3.907E-07 6.251E-07 1.000E-06
 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.70000009E-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.09999999E-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.53600001     sec.
 Computing g and RHS_B:
  js=   1 j_RZetaZ= 1 dgemm took  0.000E+00 sec, rest took  1.000E-03 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  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:   1.00000005E-03 sec.
 Regularization term:   0.00000000     sec.
 Optimal LWORK:        9408
 Solving system for lambda= 0.000E+00 (  1 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      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-26 (  2 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76397136445425E+15, method 2 =  6.76397136445425E+15, relative difference =  1.47842139790118E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.35025395388844E-23, method 2 =  5.68434188608080E-14, relative difference = -9.99999999173078E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.350E-23,  chi2_M: 6.764E+15
   max(B_n): 3.866E-13,  max(M): 2.417E+09,  min(M): 2.365E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.600E-26 (  3 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396721839485E+15, method 2 =  6.76396721839485E+15, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  6.01135815206808E-23, method 2 =  0.00000000000000E+00, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 6.011E-23,  chi2_M: 6.764E+15
   max(B_n): 6.054E-13,  max(M): 2.417E+09,  min(M): 2.365E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.560E-26 (  4 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396455099592E+15, method 2 =  6.76396455099592E+15, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  1.53718143824450E-22, method 2 =  1.13686837721616E-13, relative difference = -9.99999997295762E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.537E-22,  chi2_M: 6.764E+15
   max(B_n): 9.555E-13,  max(M): 2.416E+09,  min(M): 2.364E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.095E-26 (  5 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396370109354E+15, method 2 =  6.76396370109354E+15, relative difference =  3.69605768226672E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.93191327452098E-22, method 2 =  8.52651282912120E-14, relative difference = -9.99999990777207E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.932E-22,  chi2_M: 6.764E+15
   max(B_n): 1.522E-12,  max(M): 2.416E+09,  min(M): 2.364E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.551E-26 (  6 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396306893355E+15, method 2 =  6.76396306893355E+15, relative difference = -1.47842321108011E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.00606183186343E-21, method 2 =  0.00000000000000E+00, relative difference =  1.00000000000000E+00
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.006E-21,  chi2_M: 6.764E+15
   max(B_n): 2.430E-12,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.048E-25 (  7 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396289573471E+15, method 2 =  6.76396289573471E+15, relative difference =  7.39211624468394E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.57531599008739E-21, method 2 =  2.84217094304040E-14, relative difference = -9.99999818778264E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.575E-21,  chi2_M: 6.764E+15
   max(B_n): 3.883E-12,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.677E-25 (  8 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396287185669E+15, method 2 =  6.76396287185669E+15, relative difference = -5.17448138954562E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.59083243878404E-21, method 2 =  1.13686837721616E-13, relative difference = -9.99999884052856E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.591E-21,  chi2_M: 6.764E+15
   max(B_n): 6.208E-12,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.683E-25 (  9 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396283132615E+15, method 2 =  6.76396283132615E+15, relative difference = -7.39211631507397E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.68675002442568E-20, method 2 =  8.52651282912120E-14, relative difference = -9.99999604351808E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.687E-20,  chi2_M: 6.764E+15
   max(B_n): 9.927E-12,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.292E-25 ( 10 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396281817764E+15, method 2 =  6.76396281817763E+15, relative difference =  7.39211632944356E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.31738726549547E-20, method 2 =  5.68434188608080E-14, relative difference = -9.99998480955554E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.317E-20,  chi2_M: 6.764E+15
   max(B_n): 1.588E-11,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.866E-25 ( 11 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396281376781E+15, method 2 =  6.76396281376780E+15, relative difference =  2.95684653370517E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.10504935196088E-19, method 2 =  5.68434188608080E-14, relative difference = -9.99996111960801E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.105E-19,  chi2_M: 6.764E+15
   max(B_n): 2.540E-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.099E-24 ( 12 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396281172282E+15, method 2 =  6.76396281172282E+15, relative difference =  7.39211633649782E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.82832999776112E-19, method 2 =  1.13686837721616E-13, relative difference = -9.99995024361627E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.828E-19,  chi2_M: 6.764E+15
   max(B_n): 4.063E-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.758E-24 ( 13 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396281120755E+15, method 2 =  6.76396281120754E+15, relative difference =  5.91369306964876E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.23925163324673E-19, method 2 =  5.68434188608080E-14, relative difference = -9.99974529472067E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.239E-19,  chi2_M: 6.764E+15
   max(B_n): 6.499E-11,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.812E-24 ( 14 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396281028581E+15, method 2 =  6.76396281028579E+15, relative difference =  9.60975123948878E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.85289150131789E-18, method 2 =  8.52651282912120E-14, relative difference = -9.99956539061765E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.853E-18,  chi2_M: 6.764E+15
   max(B_n): 1.040E-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.498E-24 ( 15 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396280942818E+15, method 2 =  6.76396280942818E+15, relative difference =  2.95684653560222E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.74255608435179E-18, method 2 =  1.13686837721616E-13, relative difference = -9.99916571551336E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 4.743E-18,  chi2_M: 6.764E+15
   max(B_n): 1.663E-10,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.197E-24 ( 16 of 100)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396280813542E+15, method 2 =  6.76396280813542E+15, relative difference =  7.39211634041837E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.21387435319732E-17, method 2 =  8.52651282912120E-14, relative difference = -9.99715311150227E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.214E-17,  chi2_M: 6.764E+15
   max(B_n): 2.661E-10,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.151E-23 ( 17 of 100)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396280609853E+15, method 2 =  6.76396280609854E+15, relative difference = -8.87053961117331E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.10696521328114E-17, method 2 =  8.52651282912120E-14, relative difference = -9.99271487994026E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.107E-17,  chi2_M: 6.764E+15
   max(B_n): 4.257E-10,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.842E-23 ( 18 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396280285388E+15, method 2 =  6.76396280285388E+15, relative difference =  5.17448144233328E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.95239150565176E-17, method 2 =  5.68434188608080E-14, relative difference = -9.97205909913349E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.952E-17,  chi2_M: 6.764E+15
   max(B_n): 6.811E-10,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.947E-23 ( 19 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396279767178E+15, method 2 =  6.76396279767178E+15, relative difference = -2.95684654074150E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.03544321556312E-16, method 2 =  8.52651282912120E-14, relative difference = -9.95236984140879E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.035E-16,  chi2_M: 6.764E+15
   max(B_n): 1.090E-09,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.715E-23 ( 20 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396278938490E+15, method 2 =  6.76396278938490E+15, relative difference = -2.21763490827306E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.20977642879921E-16, method 2 =  2.84217094304040E-14, relative difference = -9.63999359039415E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 5.210E-16,  chi2_M: 6.764E+15
   max(B_n): 1.743E-09,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.543E-23 ( 21 of 100)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396277612845E+15, method 2 =  6.76396277612846E+15, relative difference = -1.47842327507955E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.33345201672967E-15, method 2 =  2.84217094304040E-14, relative difference = -9.10371716913798E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.333E-15,  chi2_M: 6.764E+15
   max(B_n): 2.789E-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.207E-22 ( 22 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396275492061E+15, method 2 =  6.76396275492061E+15, relative difference =  1.47842327971502E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.41297355553838E-15, method 2 =  5.68434188608080E-14, relative difference = -8.86718291000358E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.413E-15,  chi2_M: 6.764E+15
   max(B_n): 4.461E-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.931E-22 ( 23 of 100)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396272099244E+15, method 2 =  6.76396272099245E+15, relative difference = -1.47842328713082E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.73542200821412E-15, method 2 =  8.52651282912120E-14, relative difference = -8.14141045336679E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 8.735E-15,  chi2_M: 6.764E+15
   max(B_n): 7.137E-09,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.089E-22 ( 24 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396266671521E+15, method 2 =  6.76396266671522E+15, relative difference = -2.95684659798878E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.23577235693919E-14, method 2 =  8.52651282912120E-14, relative difference = -5.84517169302478E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.236E-14,  chi2_M: 6.764E+15
   max(B_n): 1.142E-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.942E-22 ( 25 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396257988686E+15, method 2 =  6.76396257988686E+15, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  5.72215802286322E-14, method 2 =  5.68434188608080E-14, relative difference =  3.31531469638344E-03
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.722E-14,  chi2_M: 6.764E+15
   max(B_n): 1.826E-08,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.906E-22 ( 26 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396244099280E+15, method 2 =  6.76396244099279E+15, relative difference =  6.65290506749103E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.46444863871542E-13, method 2 =  2.84217094304040E-13, relative difference = -3.19908057391797E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.464E-13,  chi2_M: 6.764E+15
   max(B_n): 2.920E-08,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.265E-21 ( 27 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396221882951E+15, method 2 =  6.76396221882951E+15, relative difference =  1.47842339689036E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.74766409550587E-13, method 2 =  4.54747350886464E-13, relative difference = -9.64190651807112E-02
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.748E-13,  chi2_M: 6.764E+15
   max(B_n): 4.670E-08,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.024E-21 ( 28 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396186351712E+15, method 2 =  6.76396186351713E+15, relative difference = -8.87054084731359E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.58970055986289E-13, method 2 =  9.66338120633736E-13, relative difference = -3.82695338695429E-03
   Diagnostics:    0.00000000      sec.
   chi2_B: 9.590E-13,  chi2_M: 6.764E+15
   max(B_n): 7.464E-08,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.237E-21 ( 29 of 100)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396129535162E+15, method 2 =  6.76396129535162E+15, relative difference = -2.95684719747651E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.45350855307058E-12, method 2 =  2.44426701101474E-12, relative difference =  1.88688557385103E-03
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.454E-12,  chi2_M: 6.764E+15
   max(B_n): 1.192E-07,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.179E-21 ( 30 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76396038703240E+15, method 2 =  6.76396038703239E+15, relative difference =  8.13133088500102E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.27601313651472E-12, method 2 =  6.39488462184090E-12, relative difference = -9.38145722530117E-03
   Diagnostics:    0.00000000      sec.
   chi2_B: 6.276E-12,  chi2_M: 6.764E+15
   max(B_n): 1.904E-07,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.286E-21 ( 31 of 100)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76395893533379E+15, method 2 =  6.76395893533378E+15, relative difference =  1.10881808593208E-15
   2 methods of computing chi2_B that should agree: method 1 =  1.60497914729570E-11, method 2 =  1.61151092470391E-11, relative difference = -2.03071586169923E-03
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.605E-11,  chi2_M: 6.764E+15
   max(B_n): 3.036E-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.326E-20 ( 32 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76395661590066E+15, method 2 =  6.76395661590066E+15, relative difference = -5.91369848617425E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.10327983708385E-11, method 2 =  4.11546352552250E-11, relative difference = -1.48242716691728E-03
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 4.103E-11,  chi2_M: 6.764E+15
   max(B_n): 4.837E-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.121E-20 ( 33 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76395291090772E+15, method 2 =  6.76395291090772E+15, relative difference =  2.21763814703834E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.04878147233760E-10, method 2 =  1.05046638054773E-10, relative difference = -8.02624715237405E-04
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.049E-10,  chi2_M: 6.764E+15
   max(B_n): 7.693E-07,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.393E-20 ( 34 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76394699287981E+15, method 2 =  6.76394699287980E+15, relative difference =  7.39213362444050E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.68033736403394E-10, method 2 =  2.68101985057001E-10, relative difference = -1.27297344451440E-04
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.680E-10,  chi2_M: 6.764E+15
   max(B_n): 1.222E-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.429E-20 ( 35 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76393753787262E+15, method 2 =  6.76393753787261E+15, relative difference =  6.65292956181753E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.85070156520806E-10, method 2 =  6.85133727529319E-10, relative difference = -4.63952914258713E-05
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.851E-10,  chi2_M: 6.764E+15
   max(B_n): 1.939E-06,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.685E-20 ( 36 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76392242590194E+15, method 2 =  6.76392242590194E+15, relative difference = -7.39216047311966E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.75147509851770E-09, method 2 =  1.75148784364865E-09, relative difference = -3.63838589268773E-06
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.751E-09,  chi2_M: 6.764E+15
   max(B_n): 3.076E-06,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.389E-19 ( 37 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76389826242038E+15, method 2 =  6.76389826242036E+15, relative difference =  8.13140556911910E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.47949979196580E-09, method 2 =  4.47960246674484E-09, relative difference = -1.14603869976257E-05
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.479E-09,  chi2_M: 6.764E+15
   max(B_n): 4.885E-06,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.223E-19 ( 38 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76385961459136E+15, method 2 =  6.76385961459136E+15, relative difference =  3.69611455951402E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.14602030123317E-08, method 2 =  1.14602585199464E-08, relative difference = -2.42174943256335E-06
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.146E-08,  chi2_M: 6.764E+15
   max(B_n): 7.767E-06,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.556E-19 ( 39 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76379779052627E+15, method 2 =  6.76379779052627E+15, relative difference = -2.21768900616896E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.93256916099034E-08, method 2 =  2.93258324290946E-08, relative difference = -2.40094683866040E-06
   Diagnostics:    4.00000019E-03  sec.
   chi2_B: 2.933E-08,  chi2_M: 6.764E+15
   max(B_n): 1.237E-05,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.690E-19 ( 40 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76369888595453E+15, method 2 =  6.76369888595452E+15, relative difference =  8.87088574043351E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.50509173026313E-08, method 2 =  7.50510196212417E-08, relative difference = -6.81660827035484E-07
   Diagnostics:    3.00000003E-03  sec.
   chi2_B: 7.505E-08,  chi2_M: 6.764E+15
   max(B_n): 1.972E-05,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 9.103E-19 ( 41 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76354065880558E+15, method 2 =  6.76354065880558E+15, relative difference =  7.39257772257258E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.92082242268452E-07, method 2 =  1.92082268313243E-07, relative difference = -6.77959310393704E-08
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.921E-07,  chi2_M: 6.764E+15
   max(B_n): 3.147E-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.456E-18 ( 42 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76328753003580E+15, method 2 =  6.76328753003580E+15, relative difference = -7.39285440371264E-17
   2 methods of computing chi2_B that should agree: method 1 =  4.91614756854870E-07, method 2 =  4.91614912334626E-07, relative difference = -1.58131676311790E-07
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.916E-07,  chi2_M: 6.763E+15
   max(B_n): 5.028E-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.330E-18 ( 43 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76288258966445E+15, method 2 =  6.76288258966445E+15, relative difference =  7.39329706483650E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.25822501868393E-06, method 2 =  1.25822504060125E-06, relative difference = -8.70961544383012E-09
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.258E-06,  chi2_M: 6.763E+15
   max(B_n): 8.044E-05,  max(M): 2.415E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.728E-18 ( 44 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76223481809616E+15, method 2 =  6.76223481809616E+15, relative difference = -2.21820158623582E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.22016704823693E-06, method 2 =  3.22016708764750E-06, relative difference = -6.11933521430880E-09
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.220E-06,  chi2_M: 6.762E+15
   max(B_n): 1.288E-04,  max(M): 2.414E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.964E-18 ( 45 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.76119866882628E+15, method 2 =  6.76119866882626E+15, relative difference =  1.18322214622763E-15
   2 methods of computing chi2_B that should agree: method 1 =  8.24089636388021E-06, method 2 =  8.24089639195336E-06, relative difference = -1.70328236243252E-09
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.241E-06,  chi2_M: 6.761E+15
   max(B_n): 2.061E-04,  max(M): 2.414E+09,  min(M): 2.363E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 9.541E-18 ( 46 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.75954147193985E+15, method 2 =  6.75954147193985E+15, relative difference =  2.95878057454397E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.10877838455813E-05, method 2 =  2.10877839208479E-05, relative difference = -1.78460159161806E-09
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.109E-05,  chi2_M: 6.760E+15
   max(B_n): 3.297E-04,  max(M): 2.414E+09,  min(M): 2.362E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.526E-17 ( 47 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.75689145873632E+15, method 2 =  6.75689145873632E+15, relative difference =  1.47997049546660E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.39540068502843E-05, method 2 =  5.39540069723898E-05, relative difference = -1.13157087802273E-09
   Diagnostics:    0.00000000      sec.
   chi2_B: 5.395E-05,  chi2_M: 6.757E+15
   max(B_n): 5.275E-04,  max(M): 2.414E+09,  min(M): 2.362E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.442E-17 ( 48 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.75265505394312E+15, method 2 =  6.75265505394313E+15, relative difference = -9.62584338763819E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.38011212985898E-04, method 2 =  1.38011213010714E-04, relative difference = -8.99032046016281E-11
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.380E-04,  chi2_M: 6.753E+15
   max(B_n): 8.437E-04,  max(M): 2.413E+09,  min(M): 2.361E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.907E-17 ( 49 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.74588569145824E+15, method 2 =  6.74588569145824E+15, relative difference =  6.67073265960907E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.52891948755383E-04, method 2 =  3.52891948836032E-04, relative difference = -1.14268345234738E-10
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.529E-04,  chi2_M: 6.746E+15
   max(B_n): 1.349E-03,  max(M): 2.412E+09,  min(M): 2.360E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.251E-17 ( 50 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.73507683211567E+15, method 2 =  6.73507683211566E+15, relative difference =  9.65096636909215E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.01795643349913E-04, method 2 =  9.01795643414971E-04, relative difference = -3.60714771460177E-11
   Diagnostics:    0.00000000      sec.
   chi2_B: 9.018E-04,  chi2_M: 6.735E+15
   max(B_n): 2.157E-03,  max(M): 2.410E+09,  min(M): 2.358E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-16 ( 51 of 100)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.71783816809640E+15, method 2 =  6.71783816809640E+15, relative difference =  1.48857411413852E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.30227747510036E-03, method 2 =  2.30227747528033E-03, relative difference = -3.90863816737313E-11
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.302E-03,  chi2_M: 6.718E+15
   max(B_n): 3.446E-03,  max(M): 2.407E+09,  min(M): 2.355E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.600E-16 ( 52 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.69039615385338E+15, method 2 =  6.69039615385338E+15, relative difference =  3.73669950554439E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.86870269126903E-03, method 2 =  5.86870269137307E-03, relative difference = -8.86397266516575E-12
   Diagnostics:    0.00000000      sec.
   chi2_B: 5.869E-03,  chi2_M: 6.690E+15
   max(B_n): 5.501E-03,  max(M): 2.402E+09,  min(M): 2.350E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.560E-16 ( 53 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.64684143325641E+15, method 2 =  6.64684143325641E+15, relative difference = -4.51342194653536E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.49233778703814E-02, method 2 =  1.49233778704740E-02, relative difference = -3.10499851514022E-12
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.492E-02,  chi2_M: 6.647E+15
   max(B_n): 8.771E-03,  max(M): 2.394E+09,  min(M): 2.342E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.095E-16 ( 54 of 100)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.57803989486457E+15, method 2 =  6.57803989486456E+15, relative difference =  3.80052422903627E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.78014957769990E-02, method 2 =  3.78014957770461E-02, relative difference = -6.23007249469114E-13
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.780E-02,  chi2_M: 6.578E+15
   max(B_n): 1.396E-02,  max(M): 2.382E+09,  min(M): 2.330E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.551E-16 ( 55 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.47016937336078E+15, method 2 =  6.47016937336079E+15, relative difference = -4.63666378248413E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.51672658181530E-02, method 2 =  9.51672658183043E-02, relative difference = -7.94528730239026E-13
   Diagnostics:    0.00000000      sec.
   chi2_B: 9.517E-02,  chi2_M: 6.470E+15
   max(B_n): 2.214E-02,  max(M): 2.363E+09,  min(M): 2.310E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.048E-15 ( 56 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.30303198952532E+15, method 2 =  6.30303198952531E+15, relative difference =  4.75961411109057E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.37289985627434E-01, method 2 =  2.37289985627541E-01, relative difference = -2.24580506594436E-13
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.373E-01,  chi2_M: 6.303E+15
   max(B_n): 3.493E-02,  max(M): 2.333E+09,  min(M): 2.280E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.677E-15 ( 57 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.04880219306743E+15, method 2 =  6.04880219306743E+15, relative difference =  4.95965962226095E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.82848415835721E-01, method 2 =  5.82848415835713E-01, relative difference =  7.14308597094622E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.828E-01,  chi2_M: 6.049E+15
   max(B_n): 5.471E-02,  max(M): 2.286E+09,  min(M): 2.233E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.683E-15 ( 58 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.67294210905606E+15, method 2 =  5.67294210905606E+15, relative difference = -4.40688438545688E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.39910272931027E+00, method 2 =  1.39910272931041E+00, relative difference = -5.31661765012438E-14
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.399E+00,  chi2_M: 5.673E+15
   max(B_n): 8.468E-02,  max(M): 2.215E+09,  min(M): 2.161E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.292E-15 ( 59 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.14060854013821E+15, method 2 =  5.14060854013820E+15, relative difference =  3.89058996495040E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.24496694174026E+00, method 2 =  3.24496694174030E+00, relative difference = -5.74790040964166E-15
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 3.245E+00,  chi2_M: 5.141E+15
   max(B_n): 1.288E-01,  max(M): 2.109E+09,  min(M): 2.056E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.866E-15 ( 60 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.43256577384699E+15, method 2 =  4.43256577384700E+15, relative difference = -4.51205938510916E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.16153291843059E+00, method 2 =  7.16153291843057E+00, relative difference =  1.86031069705974E-15
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 7.162E+00,  chi2_M: 4.433E+15
   max(B_n): 1.911E-01,  max(M): 1.959E+09,  min(M): 1.909E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.099E-14 ( 61 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.57013266657521E+15, method 2 =  3.57013266657521E+15, relative difference =  4.20152453729103E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.47635796066682E+01, method 2 =  1.47635796066682E+01, relative difference = -4.21120697252363E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.476E+01,  chi2_M: 3.570E+15
   max(B_n): 2.741E-01,  max(M): 1.759E+09,  min(M): 1.712E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.758E-14 ( 62 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.63475998406749E+15, method 2 =  2.63475998406749E+15, relative difference =  1.89770606439874E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.78873106262029E+01, method 2 =  2.78873106262028E+01, relative difference =  5.09581397277137E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.789E+01,  chi2_M: 2.635E+15
   max(B_n): 3.763E-01,  max(M): 1.512E+09,  min(M): 1.471E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.812E-14 ( 63 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.75443488365711E+15, method 2 =  1.75443488365710E+15, relative difference =  4.98736093400098E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.75294289173847E+01, method 2 =  4.75294289173847E+01, relative difference =  2.24242985433875E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 4.753E+01,  chi2_M: 1.754E+15
   max(B_n): 4.909E-01,  max(M): 1.234E+09,  min(M): 1.200E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.498E-14 ( 64 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.04721103880775E+15, method 2 =  1.04721103880775E+15, relative difference = -1.79047004903108E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.26139984064391E+01, method 2 =  7.26139984064390E+01, relative difference =  3.91408131409049E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 7.261E+01,  chi2_M: 1.047E+15
   max(B_n): 6.066E-01,  max(M): 9.532E+08,  min(M): 9.270E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.197E-14 ( 65 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.63011167531837E+14, method 2 =  5.63011167531836E+14, relative difference =  6.66061388522626E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.99227847652137E+01, method 2 =  9.99227847652136E+01, relative difference =  7.11091806968397E-17
   Diagnostics:    0.00000000      sec.
   chi2_B: 9.992E+01,  chi2_M: 5.630E+14
   max(B_n): 7.117E-01,  max(M): 6.988E+08,  min(M): 6.798E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.151E-13 ( 66 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.76564660630135E+14, method 2 =  2.76564660630135E+14, relative difference =  3.38980402580707E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.25633618891565E+02, method 2 =  1.25633618891565E+02, relative difference = -5.09010619789625E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.256E+02,  chi2_M: 2.766E+14
   max(B_n): 7.981E-01,  max(M): 4.896E+08,  min(M): 4.766E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.842E-13 ( 67 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.26508433143286E+14, method 2 =  1.26508433143286E+14, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  1.47092485591262E+02, method 2 =  1.47092485591262E+02, relative difference =  1.93223394901230E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.471E+02,  chi2_M: 1.265E+14
   max(B_n): 8.636E-01,  max(M): 3.310E+08,  min(M): 3.224E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.947E-13 ( 68 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.49135682302521E+13, method 2 =  5.49135682302520E+13, relative difference =  6.40210627956103E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.63422805716366E+02, method 2 =  1.63422805716366E+02, relative difference = -5.21745590631927E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.634E+02,  chi2_M: 5.491E+13
   max(B_n): 9.102E-01,  max(M): 2.180E+08,  min(M): 2.125E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.715E-13 ( 69 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.29810120682338E+13, method 2 =  2.29810120682338E+13, relative difference = -3.39954566700700E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.75050899391444E+02, method 2 =  1.75050899391444E+02, relative difference =  1.21771908324425E-15
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.751E+02,  chi2_M: 2.298E+13
   max(B_n): 9.419E-01,  max(M): 1.410E+08,  min(M): 1.375E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.543E-13 ( 70 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.38427173015714E+12, method 2 =  9.38427173015714E+12, relative difference = -1.04063749226457E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.82960865419200E+02, method 2 =  1.82960865419200E+02, relative difference =  7.76715538737867E-17
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.830E+02,  chi2_M: 9.384E+12
   max(B_n): 9.629E-01,  max(M): 9.009E+07,  min(M): 8.789E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.207E-12 ( 71 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.77097980582483E+12, method 2 =  3.77097980582483E+12, relative difference =  1.94225880994819E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.88180448077848E+02, method 2 =  1.88180448077847E+02, relative difference =  2.26551507242503E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.882E+02,  chi2_M: 3.771E+12
   max(B_n): 9.765E-01,  max(M): 5.710E+07,  min(M): 5.572E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.931E-12 ( 72 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.49973809284873E+12, method 2 =  1.49973809284873E+12, relative difference =  4.06972100935734E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.91557235141042E+02, method 2 =  1.91557235141042E+02, relative difference =  6.67673514616391E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.916E+02,  chi2_M: 1.500E+12
   max(B_n): 9.852E-01,  max(M): 3.601E+07,  min(M): 3.514E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.089E-12 ( 73 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.92537072299477E+11, method 2 =  5.92537072299478E+11, relative difference = -6.18038861397876E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.93714285643284E+02, method 2 =  1.93714285643284E+02, relative difference =  7.33598695006450E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.937E+02,  chi2_M: 5.925E+11
   max(B_n): 9.907E-01,  max(M): 2.263E+07,  min(M): 2.209E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.942E-12 ( 74 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.33134224611829E+11, method 2 =  2.33134224611829E+11, relative difference =  6.54506608281389E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.95081118766217E+02, method 2 =  1.95081118766217E+02, relative difference =  9.46996369848677E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.951E+02,  chi2_M: 2.331E+11
   max(B_n): 9.942E-01,  max(M): 1.419E+07,  min(M): 1.386E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.906E-12 ( 75 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.14865003545717E+10, method 2 =  9.14865003545719E+10, relative difference = -7.50542981916770E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.95942825454050E+02, method 2 =  1.95942825454050E+02, relative difference = -7.25255169831902E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.959E+02,  chi2_M: 9.149E+10
   max(B_n): 9.963E-01,  max(M): 8.892E+06,  min(M): 8.681E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.265E-11 ( 76 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.58420141366818E+10, method 2 =  3.58420141366818E+10, relative difference = -4.25723537865683E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.96484342767688E+02, method 2 =  1.96484342767688E+02, relative difference =  1.44651268544122E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.965E+02,  chi2_M: 3.584E+10
   max(B_n): 9.977E-01,  max(M): 5.565E+06,  min(M): 5.434E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.024E-11 ( 77 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.40274655247229E+10, method 2 =  1.40274655247229E+10, relative difference =  8.83820822188487E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.96823961936478E+02, method 2 =  1.96823961936478E+02, relative difference =  2.16602509807037E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.968E+02,  chi2_M: 1.403E+10
   max(B_n): 9.986E-01,  max(M): 3.482E+06,  min(M): 3.399E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.237E-11 ( 78 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.48637019359584E+09, method 2 =  5.48637019359583E+09, relative difference =  3.47652193619549E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97036690085498E+02, method 2 =  1.97036690085498E+02, relative difference =  0.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.970E+02,  chi2_M: 5.486E+09
   max(B_n): 9.991E-01,  max(M): 2.177E+06,  min(M): 2.126E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.179E-11 ( 79 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.14494037925582E+09, method 2 =  2.14494037925582E+09, relative difference = -2.77884855690349E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97169832131981E+02, method 2 =  1.97169832131981E+02, relative difference =  2.16222551313371E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.972E+02,  chi2_M: 2.145E+09
   max(B_n): 9.994E-01,  max(M): 1.361E+06,  min(M): 1.329E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.286E-11 ( 80 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  8.38369345237452E+08, method 2 =  8.38369345237453E+08, relative difference = -3.55479629079882E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97253121813385E+02, method 2 =  1.97253121813385E+02, relative difference =  2.88175002444756E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.973E+02,  chi2_M: 8.384E+08
   max(B_n): 9.996E-01,  max(M): 8.511E+05,  min(M): 8.310E+05
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.326E-10 ( 81 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.27632366322320E+08, method 2 =  3.27632366322320E+08, relative difference = -9.09626931008891E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.97305209269166E+02, method 2 =  1.97305209269165E+02, relative difference =  6.48222583228093E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.973E+02,  chi2_M: 3.276E+08
   max(B_n): 9.998E-01,  max(M): 5.321E+05,  min(M): 5.195E+05
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.121E-10 ( 82 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.28025127891491E+08, method 2 =  1.28025127891491E+08, relative difference =  4.65569890513284E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97337777283514E+02, method 2 =  1.97337777283514E+02, relative difference =  4.32077068389760E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.973E+02,  chi2_M: 1.280E+08
   max(B_n): 9.999E-01,  max(M): 3.326E+05,  min(M): 3.248E+05
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.393E-10 ( 83 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.00238096623839E+07, method 2 =  5.00238096623839E+07, relative difference =  5.95762749555350E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97358138185782E+02, method 2 =  1.97358138185781E+02, relative difference =  4.32032492173940E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 5.002E+07
   max(B_n): 9.999E-01,  max(M): 2.079E+05,  min(M): 2.030E+05
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.429E-10 ( 84 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.95452625165799E+07, method 2 =  1.95452625165799E+07, relative difference = -4.76495301009877E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97370866474024E+02, method 2 =  1.97370866474024E+02, relative difference =  4.32004630746422E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 1.955E+07
   max(B_n): 9.999E-01,  max(M): 1.300E+05,  min(M): 1.269E+05
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.685E-10 ( 85 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  7.63652452519920E+06, method 2 =  7.63652452519920E+06, relative difference =  6.09781695548987E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.97378822982095E+02, method 2 =  1.97378822982095E+02, relative difference =  7.19978693787788E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 7.637E+06
   max(B_n): 1.000E+00,  max(M): 8.123E+04,  min(M): 7.931E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.389E-09 ( 86 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.98361944819648E+06, method 2 =  2.98361944819648E+06, relative difference = -6.24290457134832E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97383796483146E+02, method 2 =  1.97383796483147E+02, relative difference = -2.15988165721831E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 2.984E+06
   max(B_n): 1.000E+00,  max(M): 5.078E+04,  min(M): 4.958E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.223E-09 ( 87 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.16570055866673E+06, method 2 =  1.16570055866673E+06, relative difference = -5.99203565416929E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97386905291282E+02, method 2 =  1.97386905291281E+02, relative difference =  2.15984763947201E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 1.166E+06
   max(B_n): 1.000E+00,  max(M): 3.174E+04,  min(M): 3.099E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.556E-09 ( 88 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.55436690609799E+05, method 2 =  4.55436690609799E+05, relative difference =  5.75128178100650E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97388848505253E+02, method 2 =  1.97388848505253E+02, relative difference =  7.19942125546359E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 4.554E+05
   max(B_n): 1.000E+00,  max(M): 1.984E+04,  min(M): 1.937E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.690E-09 ( 89 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.77937474534689E+05, method 2 =  1.77937474534689E+05, relative difference = -8.99591656735746E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97390063135811E+02, method 2 =  1.97390063135810E+02, relative difference =  8.63925234499236E-16
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 1.779E+05
   max(B_n): 1.000E+00,  max(M): 1.240E+04,  min(M): 1.211E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 9.103E-09 ( 90 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.95193699627261E+04, method 2 =  6.95193699627261E+04, relative difference =  2.09321736318512E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97390822352642E+02, method 2 =  1.97390822352642E+02, relative difference =  7.19934926346986E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 6.952E+04
   max(B_n): 1.000E+00,  max(M): 7.751E+03,  min(M): 7.568E+03
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.456E-08 ( 91 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.71608655931174E+04, method 2 =  2.71608655931174E+04, relative difference =  8.03651590841843E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97391296907289E+02, method 2 =  1.97391296907289E+02, relative difference =  7.91926515081575E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.974E+02,  chi2_M: 2.716E+04
   max(B_n): 1.000E+00,  max(M): 4.845E+03,  min(M): 4.730E+03
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.330E-08 ( 92 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.06116029983343E+04, method 2 =  1.06116029983343E+04, relative difference =  5.14245417162152E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97391593531005E+02, method 2 =  1.97391593531005E+02, relative difference =  3.59966056836404E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 1.061E+04
   max(B_n): 1.000E+00,  max(M): 3.028E+03,  min(M): 2.957E+03
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.728E-08 ( 93 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.14589344944927E+03, method 2 =  4.14589344944926E+03, relative difference =  4.38744850953053E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97391778937538E+02, method 2 =  1.97391778937538E+02, relative difference =  1.43986287490715E-16
   Diagnostics:    8.00000038E-03  sec.
   chi2_B: 1.974E+02,  chi2_M: 4.146E+03
   max(B_n): 1.000E+00,  max(M): 1.893E+03,  min(M): 1.848E+03
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.964E-08 ( 94 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.61977662028013E+03, method 2 =  1.61977662028013E+03, relative difference =  6.31680644530844E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97391894826987E+02, method 2 =  1.97391894826987E+02, relative difference =  2.87972405911757E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.974E+02,  chi2_M: 1.620E+03
   max(B_n): 1.000E+00,  max(M): 1.183E+03,  min(M): 1.155E+03
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 9.541E-08 ( 95 of 100)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.32837228936542E+02, method 2 =  6.32837228936540E+02, relative difference =  9.88054398316046E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97391967264341E+02, method 2 =  1.97391967264341E+02, relative difference =  7.91923825643240E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.974E+02,  chi2_M: 6.328E+02
   max(B_n): 1.000E+00,  max(M): 7.395E+02,  min(M): 7.220E+02
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.526E-07 ( 96 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.47245761055072E+02, method 2 =  2.47245761055072E+02, relative difference = -2.87383182113214E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97392012541705E+02, method 2 =  1.97392012541705E+02, relative difference = -3.59965292724282E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.974E+02,  chi2_M: 2.472E+02
   max(B_n): 1.000E+00,  max(M): 4.622E+02,  min(M): 4.513E+02
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.442E-07 ( 97 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.65974476637915E+01, method 2 =  9.65974476637916E+01, relative difference = -2.94228368531290E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97392040842565E+02, method 2 =  1.97392040842565E+02, relative difference = -1.43986096445867E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.974E+02,  chi2_M: 9.660E+01
   max(B_n): 1.000E+00,  max(M): 2.889E+02,  min(M): 2.821E+02
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.907E-07 ( 98 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.77400459918734E+01, method 2 =  3.77400459918734E+01, relative difference =  4.70682213737300E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97392058532167E+02, method 2 =  1.97392058532167E+02, relative difference =  7.19930417711622E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.974E+02,  chi2_M: 3.774E+01
   max(B_n): 1.000E+00,  max(M): 1.806E+02,  min(M): 1.763E+02
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.251E-07 ( 99 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.47448101171810E+01, method 2 =  1.47448101171810E+01, relative difference = -2.40946722987011E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97392069589146E+02, method 2 =  1.97392069589146E+02, relative difference =  7.19930377384493E-17
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.974E+02,  chi2_M: 1.474E+01
   max(B_n): 1.000E+00,  max(M): 1.129E+02,  min(M): 1.102E+02
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-06 (100 of 100)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.76070895282086E+00, method 2 =  5.76070895282087E+00, relative difference = -3.08357331354228E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.97392076500369E+02, method 2 =  1.97392076500369E+02, relative difference = -4.31958211306687E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.974E+02,  chi2_M: 5.761E+00
   max(B_n): 1.000E+00,  max(M): 7.055E+01,  min(M): 6.889E+01
 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    2.00000009E-03 sec
 REGCOIL_PM complete. Total time=   1.6900000572204590      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
