This is REGCOIL_PM,
a regularized least-squares method for designing permanent magnets for stellarators.
 Successfully read parameters from regcoil_nml namelist in regcoil_in.benchmark.
 Resolution parameters:
   ntheta_plasma  =   96
   ntheta_coil    =   96
   nzeta_plasma   =    4
   nzeta_coil     =    4
   mpol_magnetization =   24
   ntor_magnetization =    0
   ns_magnetization =    1
   ns_integration   =    2
 Not assuming stellarator symmetry
 We will use the following values of the regularization weight lambda:
 0.000E+00 1.000E-34 4.833E-34 2.336E-33 1.129E-32 5.456E-32 2.637E-31 1.274E-30 6.158E-30 2.976E-29 1.438E-28 6.952E-28 3.360E-27 1.624E-26 7.848E-26 3.793E-25 1.833E-24 8.859E-24 4.281E-23 2.069E-22 1.000E-21 4.833E-21 2.336E-20 1.129E-19 5.456E-19 2.637E-18 1.274E-17 6.158E-17 2.976E-16 1.438E-15 6.952E-15 3.360E-14 1.624E-13 7.848E-13 3.793E-12 1.833E-11 8.859E-11 4.281E-10 2.069E-09 1.000E-08
 Initializing plasma surface.
 Reading FOCUS format data from file:plasma
 Number of modes for Bnormal read from FOCUS file:           1
 Plasma surface area: 6.514E+02 m^2. Volume: 1.790E+02 m^3.
 G =    0.0000000000000000       ; curpol =    1.0000000000000000     
 Done initializing plasma surface. Took    3.59999985E-02  sec.
 Initializing coil surface.
   Building a plain circular torus.
   Evaluating coil surface & derivatives:   1.70000009E-02  sec.
 s_integration:  0.21132486540518713       0.78867513459481287     
 s_weights:  0.50000000000000011       0.50000000000000011     
 s_magnetization  0.50000000000000000     
 s_magnetization_weights   1.0000000000000000     
 interpolate_magnetization_to_integration:
   1.0000000000000000     
   1.0000000000000000     
 Maximum allowed uniform d before singularity:  2.9000E+01
 Coil surface area: 1.184E+03 m^2. Volume: 5.918E+02 m^3.
 Done initializing coil surface. Took    1.99999996E-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   4.06599998     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:   3.00000003E-03 sec.
 About to form matrix_B.
  js=   1, ks=   1, j_RZetaZ= 1, k_RZetaZ= 1 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 2, k_RZetaZ= 1 dgemm took  1.000E-03 sec.
  js=   1, ks=   1, j_RZetaZ= 2, k_RZetaZ= 2 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 1 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 2 dgemm took  0.000E+00 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 3 dgemm took  0.000E+00 sec.
 Total time for matrix_B:   1.00000005E-03 sec.
 Regularization term:   1.00000005E-03 sec.
 Optimal LWORK:        9408
 Solving system for lambda= 0.000E+00 (  1 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.88972019627995E+26, method 2 =  4.88972019627997E+26, relative difference = -1.82700269485287E-15
   2 methods of computing chi2_B that should agree: method 1 =  1.28864875128574E-17, method 2 =  1.59161572810262E-12, relative difference = -9.99983807167835E-01
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 1.289E-17,  chi2_M: 4.890E+26
   max(B_n): 2.545E-10,  max(M): 1.716E+15,  min(M): 1.625E+12
 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 =  1.64171467842482E+20, method 2 =  1.64171467842482E+20, relative difference =  2.99394289677424E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.06810407781517E-19, method 2 =  4.54747350886464E-13, relative difference = -9.99996891419784E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.068E-19,  chi2_M: 1.642E+20
   max(B_n): 4.963E-11,  max(M): 7.040E+11,  min(M): 3.850E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.833E-34 (  3 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.78842148491129E+20, method 2 =  1.78842148491129E+20, relative difference =  4.58057570271604E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.10067779097618E-20, method 2 =  6.82121026329696E-13, relative difference = -9.99999967727786E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.101E-20,  chi2_M: 1.788E+20
   max(B_n): 6.488E-12,  max(M): 7.091E+11,  min(M): 7.845E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.336E-33 (  4 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.23820297827646E+19, method 2 =  6.23820297827646E+19, relative difference =  3.93959608008620E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.10277030069517E-21, method 2 =  6.82121026329696E-13, relative difference = -9.99999993834612E-01
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 2.103E-21,  chi2_M: 6.238E+19
   max(B_n): 2.917E-12,  max(M): 3.878E+11,  min(M): 2.958E+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.22765535269027E+17, method 2 =  7.22765535269026E+17, relative difference =  2.65646313542792E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.53012876923973E-23, method 2 =  9.09494701772928E-13, relative difference = -9.99999999856401E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.530E-23,  chi2_M: 7.228E+17
   max(B_n): 5.049E-13,  max(M): 4.428E+10,  min(M): 1.738E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.33813678877017E+17, method 2 =  5.33813678877018E+17, relative difference = -1.19892019505077E-15
   2 methods of computing chi2_B that should agree: method 1 =  6.05649909745192E-24, method 2 =  1.08002495835535E-12, relative difference = -9.99999999988785E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.056E-24,  chi2_M: 5.338E+17
   max(B_n): 1.902E-13,  max(M): 3.667E+10,  min(M): 8.076E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.637E-31 (  7 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.72126306405516E+16, method 2 =  2.72126306405516E+16, relative difference = -7.34952833637352E-17
   2 methods of computing chi2_B that should agree: method 1 =  7.83016043882300E-25, method 2 =  9.09494701772928E-13, relative difference = -9.99999999998278E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.830E-25,  chi2_M: 2.721E+16
   max(B_n): 7.783E-14,  max(M): 7.851E+09,  min(M): 3.299E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.274E-30 (  8 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.17525791270104E+15, method 2 =  3.17525791270104E+15, relative difference =  3.14935047008306E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.73753349691946E-25, method 2 =  8.52651282912120E-13, relative difference = -9.99999999999123E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.738E-25,  chi2_M: 3.175E+15
   max(B_n): 4.757E-14,  max(M): 3.532E+09,  min(M): 9.571E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.158E-30 (  9 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.05067560634217E+15, method 2 =  2.05067560634217E+15, relative difference =  5.48599689058906E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.71567426064868E-25, method 2 =  6.25277607468888E-13, relative difference = -9.99999999999131E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.716E-25,  chi2_M: 2.051E+15
   max(B_n): 4.671E-14,  max(M): 2.266E+09,  min(M): 9.392E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.976E-29 ( 10 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.46109416746097E+15, method 2 =  1.46109416746097E+15, relative difference = -1.71104645797362E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.78228733393535E-25, method 2 =  9.66338120633736E-13, relative difference = -9.99999999999631E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.782E-25,  chi2_M: 1.461E+15
   max(B_n): 3.886E-14,  max(M): 1.963E+09,  min(M): 1.818E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.438E-28 ( 11 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.00367930519800E+15, method 2 =  1.00367930519800E+15, relative difference =  1.24541772807940E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.71647507043014E-25, method 2 =  9.66338120633736E-13, relative difference = -9.99999999999645E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.716E-25,  chi2_M: 1.004E+15
   max(B_n): 3.703E-14,  max(M): 1.272E+09,  min(M): 5.126E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.952E-28 ( 12 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.17221546643586E+14, method 2 =  9.17221546643585E+14, relative difference =  6.81405710852607E-17
   2 methods of computing chi2_B that should agree: method 1 =  3.28219369433402E-24, method 2 =  9.66338120633736E-13, relative difference = -9.99999999993207E-01
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.282E-24,  chi2_M: 9.172E+14
   max(B_n): 8.915E-14,  max(M): 9.438E+08,  min(M): 8.034E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.360E-27 ( 13 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.13991630097167E+14, method 2 =  9.13991630097168E+14, relative difference = -3.41906850904945E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.44223389914434E-23, method 2 =  9.66338120633736E-13, relative difference = -9.99999999845970E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 7.442E-23,  chi2_M: 9.140E+14
   max(B_n): 3.623E-13,  max(M): 8.890E+08,  min(M): 8.595E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.624E-26 ( 14 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.13893233672680E+14, method 2 =  9.13893233672680E+14, relative difference = -2.05166197857150E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.72671744889682E-21, method 2 =  7.95807864051312E-13, relative difference = -9.99999995660467E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.727E-21,  chi2_M: 9.139E+14
   max(B_n): 1.688E-12,  max(M): 8.848E+08,  min(M): 8.693E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.848E-26 ( 15 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.13890868375840E+14, method 2 =  9.13890868375839E+14, relative difference =  6.83889096201110E-17
   2 methods of computing chi2_B that should agree: method 1 =  4.02255995403111E-20, method 2 =  9.66338120633736E-13, relative difference = -9.99999916746329E-01
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 4.023E-20,  chi2_M: 9.139E+14
   max(B_n): 8.087E-12,  max(M): 8.848E+08,  min(M): 8.707E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.793E-25 ( 16 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    2.00000009E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.13890264233450E+14, method 2 =  9.13890264233449E+14, relative difference =  3.41944774148697E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.32892437629569E-19, method 2 =  9.09494701772928E-13, relative difference = -9.99997948549939E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.329E-19,  chi2_M: 9.139E+14
   max(B_n): 3.891E-11,  max(M): 8.852E+08,  min(M): 8.714E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.833E-24 ( 17 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.13890228818101E+14, method 2 =  9.13890228818102E+14, relative difference = -6.83889574799687E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.10853251493349E-17, method 2 =  7.95807864051312E-13, relative difference = -9.99947010409295E-01
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 2.109E-17,  chi2_M: 9.139E+14
   max(B_n): 1.850E-10,  max(M): 8.853E+08,  min(M): 8.715E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.859E-24 ( 18 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    8.00000038E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.13890151287020E+14, method 2 =  9.13890151287019E+14, relative difference =  5.47111706254692E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.23302923254455E-16, method 2 =  8.52651282912120E-13, relative difference = -9.99007582853169E-01
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 4.233E-16,  chi2_M: 9.139E+14
   max(B_n): 8.309E-10,  max(M): 8.853E+08,  min(M): 8.715E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.281E-23 ( 19 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.13889937312209E+14, method 2 =  9.13889937312208E+14, relative difference =  2.05166937882527E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.40744472099670E-15, method 2 =  7.95807864051312E-13, relative difference = -9.86501893656320E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.407E-15,  chi2_M: 9.139E+14
   max(B_n): 3.000E-09,  max(M): 8.852E+08,  min(M): 8.716E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.069E-22 ( 20 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.13889729133432E+14, method 2 =  9.13889729133432E+14, relative difference = -6.83889948727881E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.48276361263166E-14, method 2 =  9.66338120633736E-13, relative difference = -9.49902151164960E-01
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 2.483E-14,  chi2_M: 9.139E+14
   max(B_n): 6.742E-09,  max(M): 8.852E+08,  min(M): 8.717E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-21 ( 21 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.13889678263045E+14, method 2 =  9.13889678263045E+14, relative difference =  3.41944993397828E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.48740712484417E-14, method 2 =  8.52651282912120E-13, relative difference = -9.00004894479195E-01
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 4.487E-14,  chi2_M: 9.139E+14
   max(B_n): 1.087E-08,  max(M): 8.851E+08,  min(M): 8.717E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.833E-21 ( 22 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.13889655084410E+14, method 2 =  9.13889655084410E+14, relative difference = -6.83890004140897E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.08169542288506E-13, method 2 =  1.08002495835535E-12, relative difference = -8.17926202772540E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.082E-13,  chi2_M: 9.139E+14
   max(B_n): 2.220E-08,  max(M): 8.851E+08,  min(M): 8.717E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.336E-20 ( 23 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.13889559904375E+14, method 2 =  9.13889559904376E+14, relative difference = -4.78723052756810E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.44869365338048E-12, method 2 =  2.27373675443232E-12, relative difference = -2.21640973950839E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.449E-12,  chi2_M: 9.139E+14
   max(B_n): 7.442E-08,  max(M): 8.851E+08,  min(M): 8.717E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.129E-19 ( 24 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.13889100733745E+14, method 2 =  9.13889100733745E+14, relative difference = -1.36778083795550E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.27274010662222E-11, method 2 =  3.37081473844592E-11, relative difference = -1.47623725717409E-02
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 3.273E-11,  chi2_M: 9.139E+14
   max(B_n): 3.262E-07,  max(M): 8.851E+08,  min(M): 8.717E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.456E-19 ( 25 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.13886881635504E+14, method 2 =  9.13886881635505E+14, relative difference = -8.20670495518865E-16
   2 methods of computing chi2_B that should agree: method 1 =  7.63302211715342E-10, method 2 =  7.64259766583564E-10, relative difference = -6.26851729635433E-04
   Diagnostics:    4.00000019E-03  sec.
   chi2_B: 7.633E-10,  chi2_M: 9.139E+14
   max(B_n): 1.543E-06,  max(M): 8.851E+08,  min(M): 8.717E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.637E-18 ( 26 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.13876157004010E+14, method 2 =  9.13876157004011E+14, relative difference = -1.36780021058643E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.78272861182908E-08, method 2 =  1.78281425178284E-08, relative difference = -2.40187699407047E-05
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 1.783E-08,  chi2_M: 9.139E+14
   max(B_n): 7.424E-06,  max(M): 8.851E+08,  min(M): 8.717E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.274E-17 ( 27 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.13824328269185E+14, method 2 =  9.13824328269186E+14, relative difference = -3.41969446788406E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.16371014421579E-07, method 2 =  4.16372017753019E-07, relative difference = -1.20485119820595E-06
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.164E-07,  chi2_M: 9.138E+14
   max(B_n): 3.591E-05,  max(M): 8.851E+08,  min(M): 8.717E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.158E-17 ( 28 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.13573905747549E+14, method 2 =  9.13573905747548E+14, relative difference =  4.10475822088142E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.72260083572291E-06, method 2 =  9.72260176013151E-06, relative difference = -4.75391600761016E-08
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.723E-06,  chi2_M: 9.136E+14
   max(B_n): 1.736E-04,  max(M): 8.850E+08,  min(M): 8.716E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.976E-16 ( 29 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.12365080943398E+14, method 2 =  9.12365080943399E+14, relative difference = -5.48026234720637E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.26792374213641E-04, method 2 =  2.26792374917295E-04, relative difference = -1.55131845621062E-09
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.268E-04,  chi2_M: 9.124E+14
   max(B_n): 8.383E-04,  max(M): 8.844E+08,  min(M): 8.710E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.438E-15 ( 30 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.06556597895829E+14, method 2 =  9.06556597895829E+14, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  5.26351221574128E-03, method 2 =  5.26351221657251E-03, relative difference = -7.89608798307327E-11
   Diagnostics:    0.00000000      sec.
   chi2_B: 5.264E-03,  chi2_M: 9.066E+14
   max(B_n): 4.038E-03,  max(M): 8.816E+08,  min(M): 8.682E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.952E-15 ( 31 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  8.79251534195424E+14, method 2 =  8.79251534195423E+14, relative difference =  7.10831856064850E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.19237890517518E-01, method 2 =  1.19237890518377E-01, relative difference = -3.60515784348142E-12
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.192E-01,  chi2_M: 8.793E+14
   max(B_n): 1.922E-02,  max(M): 8.683E+08,  min(M): 8.550E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.360E-14 ( 32 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  7.63201837364349E+14, method 2 =  7.63201837364348E+14, relative difference =  6.55134691141084E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.41745914587770E+00, method 2 =  2.41745914587830E+00, relative difference = -1.24089899002106E-13
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.417E+00,  chi2_M: 7.632E+14
   max(B_n): 8.649E-02,  max(M): 8.092E+08,  min(M): 7.963E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.624E-13 ( 33 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.31313358153423E+14, method 2 =  4.31313358153423E+14, relative difference =  4.34718740923633E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.19100649200428E+01, method 2 =  3.19100649200436E+01, relative difference = -1.13561910997566E-14
   Diagnostics:    0.00000000      sec.
   chi2_B: 3.191E+01,  chi2_M: 4.313E+14
   max(B_n): 3.138E-01,  max(M): 6.089E+08,  min(M): 5.980E+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.91358715512144E+13, method 2 =  8.91358715512143E+13, relative difference =  7.01176741892176E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.54029103675879E+02, method 2 =  1.54029103675880E+02, relative difference = -2.39878187808420E-15
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 1.540E+02,  chi2_M: 8.914E+13
   max(B_n): 6.883E-01,  max(M): 2.771E+08,  min(M): 2.716E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.793E-12 ( 35 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.74272518890711E+12, method 2 =  6.74272518890711E+12, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  2.72147750018532E+02, method 2 =  2.72147750018532E+02, relative difference = -8.35478799393891E-16
   Diagnostics:    4.00000019E-03  sec.
   chi2_B: 2.721E+02,  chi2_M: 6.743E+12
   max(B_n): 9.143E-01,  max(M): 7.624E+07,  min(M): 7.467E+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:    3.00000003E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.32428120000501E+11, method 2 =  3.32428120000501E+11, relative difference = -9.18020356549681E-17
   2 methods of computing chi2_B that should agree: method 1 =  3.13391783479996E+02, method 2 =  3.13391783479997E+02, relative difference = -4.53453327888830E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.134E+02,  chi2_M: 3.324E+11
   max(B_n): 9.810E-01,  max(M): 1.693E+07,  min(M): 1.658E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.859E-11 ( 37 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.46728845538188E+10, method 2 =  1.46728845538188E+10, relative difference = -1.29991388252018E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.23092071689763E+02, method 2 =  3.23092071689763E+02, relative difference = -3.51871332301771E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.231E+02,  chi2_M: 1.467E+10
   max(B_n): 9.960E-01,  max(M): 3.557E+06,  min(M): 3.483E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.281E-10 ( 38 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.32206732109650E+08, method 2 =  6.32206732109650E+08, relative difference =  3.77121227902728E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.25155382488755E+02, method 2 =  3.25155382488756E+02, relative difference = -8.74096230942352E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.252E+02,  chi2_M: 6.322E+08
   max(B_n): 9.992E-01,  max(M): 7.383E+05,  min(M): 7.230E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.71026160064948E+07, method 2 =  2.71026160064948E+07, relative difference =  6.87256591313768E-17
   2 methods of computing chi2_B that should agree: method 1 =  3.25584777467005E+02, method 2 =  3.25584777467006E+02, relative difference = -1.39670949736754E-15
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.256E+02,  chi2_M: 2.710E+07
   max(B_n): 9.998E-01,  max(M): 1.529E+05,  min(M): 1.497E+05
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-08 ( 40 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.16067010588398E+06, method 2 =  1.16067010588398E+06, relative difference =  3.00900284852960E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.25673731420176E+02, method 2 =  3.25673731420177E+02, relative difference = -8.72705001612028E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.257E+02,  chi2_M: 1.161E+06
   max(B_n): 1.000E+00,  max(M): 3.164E+04,  min(M): 3.098E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   4.2179999351501465      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
