This is REGCOIL_PM,
a regularized least-squares method for designing permanent magnets for stellarators.
 Successfully read parameters from regcoil_nml namelist in regcoil_in.benchmark.
 Resolution parameters:
   ntheta_plasma  =   96
   ntheta_coil    =   96
   nzeta_plasma   =    4
   nzeta_coil     =    4
   mpol_magnetization =   24
   ntor_magnetization =    0
   ns_magnetization =    1
   ns_integration   =    2
 Not assuming stellarator symmetry
 We will use the following values of the regularization weight lambda:
 0.000E+00 1.000E-34 4.833E-34 2.336E-33 1.129E-32 5.456E-32 2.637E-31 1.274E-30 6.158E-30 2.976E-29 1.438E-28 6.952E-28 3.360E-27 1.624E-26 7.848E-26 3.793E-25 1.833E-24 8.859E-24 4.281E-23 2.069E-22 1.000E-21 4.833E-21 2.336E-20 1.129E-19 5.456E-19 2.637E-18 1.274E-17 6.158E-17 2.976E-16 1.438E-15 6.952E-15 3.360E-14 1.624E-13 7.848E-13 3.793E-12 1.833E-11 8.859E-11 4.281E-10 2.069E-09 1.000E-08
 Initializing plasma surface.
 Reading FOCUS format data from file:plasma
 Number of modes for Bnormal read from FOCUS file:           1
 Plasma surface area: 4.145E+02 m^2. Volume: 7.249E+01 m^3.
 G =    0.0000000000000000       ; curpol =    1.0000000000000000     
 Done initializing plasma surface. Took    1.60000008E-02  sec.
 Initializing coil surface.
   Building a plain circular torus.
   Evaluating coil surface & derivatives:   1.09999999E-02  sec.
 s_integration:  0.21132486540518713       0.78867513459481287     
 s_weights:  0.50000000000000011       0.50000000000000011     
 s_magnetization  0.50000000000000000     
 s_magnetization_weights   1.0000000000000000     
 interpolate_magnetization_to_integration:
   1.0000000000000000     
   1.0000000000000000     
 Maximum allowed uniform d before singularity:  2.9000E+01
 Coil surface area: 1.184E+03 m^2. Volume: 5.918E+02 m^3.
 Done initializing coil surface. Took    1.30000003E-02  sec.
 Using Bn coefficients in FOCUS file.
 Initializing basis functions
 Done. Took   1.00000005E-03 sec.
 Building inductance matrix.
   Number of OpenMP threads:           4
 Done. Took   1.64400005     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  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  1.000E-03 sec.
  js=   1, ks=   1, j_RZetaZ= 3, k_RZetaZ= 3 dgemm took  0.000E+00 sec.
 Total time for matrix_B:   1.00000005E-03 sec.
 Regularization term:   0.00000000     sec.
 Optimal LWORK:        9408
 Solving system for lambda= 0.000E+00 (  1 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.55922605675235E+30, method 2 =  1.55922605675235E+30, relative difference =  2.70783356420681E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.24787342462598E-15, method 2 = -3.41060513164848E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.248E-15,  chi2_M: 1.559E+30
   max(B_n): 5.378E-09,  max(M): 9.986E+16,  min(M): 8.083E+15
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-34 (  2 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.14349600669852E+21, method 2 =  2.14349600669851E+21, relative difference =  7.33784432107517E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.02183272806864E-21, method 2 = -1.67688085639384E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.022E-21,  chi2_M: 2.143E+21
   max(B_n): 7.978E-12,  max(M): 1.978E+12,  min(M): 7.009E+11
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.833E-34 (  3 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.65732770301484E+22, method 2 =  5.65732770301484E+22, relative difference = -2.96557259553115E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.51506760663680E-22, method 2 = -6.25277607468888E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.515E-22,  chi2_M: 5.657E+22
   max(B_n): 4.158E-12,  max(M): 7.834E+12,  min(M): 5.725E+12
 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 =  1.24834995744182E+19, method 2 =  1.24834995744182E+19, relative difference =  6.56226240980329E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.15680451880741E-25, method 2 =  3.41060513164848E-13, relative difference = -9.99999999997563E-01
   Diagnostics:    2.00000009E-03  sec.
   chi2_B: 4.157E-25,  chi2_M: 1.248E+19
   max(B_n): 1.247E-13,  max(M): 1.654E+11,  min(M): 2.915E+10
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.129E-32 (  5 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.95862341424933E+20, method 2 =  2.95862341424932E+20, relative difference =  1.10754210360747E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.80227555737128E-24, method 2 =  2.50111042987555E-12, relative difference = -9.99999999997759E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.802E-24,  chi2_M: 2.959E+20
   max(B_n): 1.950E-13,  max(M): 6.691E+11,  min(M): 3.060E+11
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.456E-32 (  6 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.86248629223507E+17, method 2 =  3.86248629223507E+17, relative difference = -8.28481904630472E-17
   2 methods of computing chi2_B that should agree: method 1 =  1.00095987564196E-25, method 2 =  1.42108547152020E-13, relative difference = -9.99999999998591E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.001E-25,  chi2_M: 3.862E+17
   max(B_n): 6.073E-14,  max(M): 3.218E+10,  min(M): 2.048E+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:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.92899199968552E+17, method 2 =  5.92899199968551E+17, relative difference =  1.07944149702672E-15
   2 methods of computing chi2_B that should agree: method 1 =  6.01482705972452E-26, method 2 = -7.38964445190504E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.015E-26,  chi2_M: 5.929E+17
   max(B_n): 4.485E-14,  max(M): 3.475E+10,  min(M): 5.258E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.274E-30 (  8 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.06668803955032E+19, method 2 =  5.06668803955033E+19, relative difference = -2.42525292737197E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.17593994139967E-25, method 2 = -5.28643795405515E-12, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.176E-25,  chi2_M: 5.067E+19
   max(B_n): 9.259E-14,  max(M): 2.470E+11,  min(M): 1.417E+11
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.158E-30 (  9 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.11642217501624E+16, method 2 =  6.11642217501623E+16, relative difference =  7.84772512859981E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.11046168428296E-26, method 2 = -4.54747350886464E-13, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.110E-26,  chi2_M: 6.116E+16
   max(B_n): 2.331E-14,  max(M): 1.012E+10,  min(M): 3.581E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.976E-29 ( 10 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  6.19291256471138E+15, method 2 =  6.19291256471138E+15, relative difference =  5.65162185551237E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.12409408421177E-26, method 2 =  2.84217094304040E-14, relative difference = -9.99999999998505E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.124E-26,  chi2_M: 6.193E+15
   max(B_n): 2.465E-14,  max(M): 3.079E+09,  min(M): 1.284E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.438E-28 ( 11 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.58697798196042E+15, method 2 =  5.58697798196041E+15, relative difference =  3.57975278667238E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.40151003843597E-26, method 2 =  2.84217094304040E-14, relative difference = -9.99999999998310E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.402E-26,  chi2_M: 5.587E+15
   max(B_n): 2.809E-14,  max(M): 2.333E+09,  min(M): 1.984E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.952E-28 ( 12 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.56408748896554E+15, method 2 =  5.56408748896553E+15, relative difference =  4.49309973101231E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.33390053784398E-26, method 2 =  2.84217094304040E-14, relative difference = -9.99999999993432E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.334E-26,  chi2_M: 5.564E+15
   max(B_n): 4.274E-14,  max(M): 2.209E+09,  min(M): 2.108E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.360E-27 ( 13 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.56302438051923E+15, method 2 =  5.56302438051923E+15, relative difference =  3.59516669925744E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.72049659661469E-24, method 2 =  0.00000000000000E+00, relative difference =  1.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.720E-24,  chi2_M: 5.563E+15
   max(B_n): 1.129E-13,  max(M): 2.188E+09,  min(M): 2.135E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.624E-26 ( 14 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.56295414421582E+15, method 2 =  5.56295414421582E+15, relative difference = -5.39281813624026E-16
   2 methods of computing chi2_B that should agree: method 1 =  3.97174947228049E-23, method 2 =  2.84217094304040E-14, relative difference = -9.99999997205130E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.972E-23,  chi2_M: 5.563E+15
   max(B_n): 4.707E-13,  max(M): 2.188E+09,  min(M): 2.142E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.848E-26 ( 15 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.56295163292266E+15, method 2 =  5.56295163292266E+15, relative difference = -7.19042742763967E-16
   2 methods of computing chi2_B that should agree: method 1 =  9.26202850740405E-22, method 2 =  8.52651282912120E-14, relative difference = -9.99999978274756E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 9.262E-22,  chi2_M: 5.563E+15
   max(B_n): 2.249E-12,  max(M): 2.189E+09,  min(M): 2.143E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.793E-25 ( 16 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.56295152507093E+15, method 2 =  5.56295152507092E+15, relative difference =  8.98803445880512E-17
   2 methods of computing chi2_B that should agree: method 1 =  2.16348435791659E-20, method 2 =  1.42108547152020E-13, relative difference = -9.99999695516692E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.163E-20,  chi2_M: 5.563E+15
   max(B_n): 1.084E-11,  max(M): 2.189E+09,  min(M): 2.143E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.833E-24 ( 17 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.56295152018427E+15, method 2 =  5.56295152018427E+15, relative difference =  8.98803446670047E-17
   2 methods of computing chi2_B that should agree: method 1 =  5.05310916122159E-19, method 2 =  8.52651282912120E-14, relative difference = -9.99988147372047E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.053E-19,  chi2_M: 5.563E+15
   max(B_n): 5.236E-11,  max(M): 2.190E+09,  min(M): 2.143E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.859E-24 ( 18 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.56295151784590E+15, method 2 =  5.56295151784590E+15, relative difference = -2.69641034114357E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.18025924737100E-17, method 2 =  5.68434188608080E-14, relative difference = -9.99584819400153E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.180E-17,  chi2_M: 5.563E+15
   max(B_n): 2.530E-10,  max(M): 2.190E+09,  min(M): 2.143E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.281E-23 ( 19 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.56295150762788E+15, method 2 =  5.56295150762789E+15, relative difference = -4.49401724349388E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.75674212684163E-16, method 2 =  1.42108547152020E-13, relative difference = -9.96127742104541E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.757E-16,  chi2_M: 5.563E+15
   max(B_n): 1.223E-09,  max(M): 2.190E+09,  min(M): 2.143E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.069E-22 ( 20 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.56295145826838E+15, method 2 =  5.56295145826837E+15, relative difference =  6.29162419671638E-16
   2 methods of computing chi2_B that should agree: method 1 =  6.43883008603086E-15, method 2 =  5.68434188608080E-14, relative difference = -7.96504384936133E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 6.439E-15,  chi2_M: 5.563E+15
   max(B_n): 5.910E-09,  max(M): 2.190E+09,  min(M): 2.143E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-21 ( 21 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.56295121974260E+15, method 2 =  5.56295121974260E+15, relative difference = -5.39282097127361E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.50376211206874E-13, method 2 =  2.27373675443232E-13, relative difference = -2.03831865891936E-01
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.504E-13,  chi2_M: 5.563E+15
   max(B_n): 2.855E-08,  max(M): 2.190E+09,  min(M): 2.143E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.833E-21 ( 22 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.56295006749775E+15, method 2 =  5.56295006749775E+15, relative difference =  8.98803681379983E-17
   2 methods of computing chi2_B that should agree: method 1 =  3.51058569881883E-12, method 2 =  3.66640051652212E-12, relative difference = -2.17103409464860E-02
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 3.511E-12,  chi2_M: 5.563E+15
   max(B_n): 1.377E-07,  max(M): 2.190E+09,  min(M): 2.143E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.336E-20 ( 23 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.56294450822327E+15, method 2 =  5.56294450822326E+15, relative difference =  6.29163205713489E-16
   2 methods of computing chi2_B that should agree: method 1 =  8.18484592022727E-11, method 2 =  8.19113665784243E-11, relative difference = -3.84144132126029E-04
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 8.185E-11,  chi2_M: 5.563E+15
   max(B_n): 6.600E-07,  max(M): 2.190E+09,  min(M): 2.143E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.129E-19 ( 24 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.56291772474458E+15, method 2 =  5.56291772474457E+15, relative difference =  3.59523562806572E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.90579857425191E-09, method 2 =  1.90593141269346E-09, relative difference = -3.48499085739387E-05
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.906E-09,  chi2_M: 5.563E+15
   max(B_n): 3.128E-06,  max(M): 2.189E+09,  min(M): 2.143E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 5.456E-19 ( 25 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.56278848634677E+15, method 2 =  5.56278848634676E+15, relative difference =  6.29180852119464E-16
   2 methods of computing chi2_B that should agree: method 1 =  4.44514816557278E-08, method 2 =  4.44515819708613E-08, relative difference = -1.12836531617732E-06
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 4.445E-08,  chi2_M: 5.563E+15
   max(B_n): 1.484E-05,  max(M): 2.189E+09,  min(M): 2.143E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.637E-18 ( 26 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.56216404587578E+15, method 2 =  5.56216404587577E+15, relative difference =  8.09037626881331E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.03798211899783E-06, method 2 =  1.03798225836726E-06, relative difference = -6.71347887883956E-08
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.038E-06,  chi2_M: 5.562E+15
   max(B_n): 7.123E-05,  max(M): 2.189E+09,  min(M): 2.143E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.274E-17 ( 27 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.55914766198925E+15, method 2 =  5.55914766198925E+15, relative difference =  2.69825536431829E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.42310019335367E-05, method 2 =  2.42310019586967E-05, relative difference = -5.19167638859698E-10
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.423E-05,  chi2_M: 5.559E+15
   max(B_n): 3.446E-04,  max(M): 2.189E+09,  min(M): 2.142E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.158E-17 ( 28 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.54460423518680E+15, method 2 =  5.54460423518680E+15, relative difference =  2.70533285402193E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.64487083167890E-04, method 2 =  5.64487083209997E-04, relative difference = -3.72967359277097E-11
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.645E-04,  chi2_M: 5.545E+15
   max(B_n): 1.664E-03,  max(M): 2.186E+09,  min(M): 2.139E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.976E-16 ( 29 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.47511529429921E+15, method 2 =  5.47511529429921E+15, relative difference =  4.56611389097696E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.30195232461707E-02, method 2 =  1.30195232462995E-02, relative difference = -4.94694098992385E-12
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 1.302E-02,  chi2_M: 5.475E+15
   max(B_n): 7.989E-03,  max(M): 2.172E+09,  min(M): 2.126E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.438E-15 ( 30 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.15701861162607E+15, method 2 =  5.15701861162607E+15, relative difference =  1.93910488852141E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.86425100252670E-01, method 2 =  2.86425100252615E-01, relative difference =  9.66127230622371E-14
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.864E-01,  chi2_M: 5.157E+15
   max(B_n): 3.743E-02,  max(M): 2.109E+09,  min(M): 2.062E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 6.952E-15 ( 31 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  3.95090348851853E+15, method 2 =  3.95090348851853E+15, relative difference =  3.79659995329943E-16
   2 methods of computing chi2_B that should agree: method 1 =  5.12514978235019E+00, method 2 =  5.12514978235035E+00, relative difference = -1.55101747146809E-14
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 5.125E+00,  chi2_M: 3.951E+15
   max(B_n): 1.578E-01,  max(M): 1.849E+09,  min(M): 1.803E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.360E-14 ( 32 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.53808471417410E+15, method 2 =  1.53808471417410E+15, relative difference = -8.12699059083497E-17
   2 methods of computing chi2_B that should agree: method 1 =  4.66008464359423E+01, method 2 =  4.66008464359423E+01, relative difference = -8.38608168212670E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 4.660E+01,  chi2_M: 1.538E+15
   max(B_n): 4.744E-01,  max(M): 1.155E+09,  min(M): 1.124E+09
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.624E-13 ( 33 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.93760435335390E+14, method 2 =  1.93760435335389E+14, relative difference =  7.25767361931166E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.37117528828683E+02, method 2 =  1.37117528828683E+02, relative difference = -4.14559825766910E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.371E+02,  chi2_M: 1.938E+14
   max(B_n): 8.137E-01,  max(M): 4.097E+08,  min(M): 3.990E+08
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 7.848E-13 ( 34 of  40)
   Additions:    1.00000005E-03  sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.14283632677513E+13, method 2 =  1.14283632677513E+13, relative difference =  7.69056976408954E-16
   2 methods of computing chi2_B that should agree: method 1 =  1.88898822472155E+02, method 2 =  1.88898822472155E+02, relative difference = -3.76149902080433E-16
   Diagnostics:    0.00000000      sec.
   chi2_B: 1.889E+02,  chi2_M: 1.143E+13
   max(B_n): 9.548E-01,  max(M): 9.940E+07,  min(M): 9.699E+07
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 3.793E-12 ( 35 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  5.26457024822474E+11, method 2 =  5.26457024822473E+11, relative difference =  1.10138901569511E-15
   2 methods of computing chi2_B that should agree: method 1 =  2.03248699406565E+02, method 2 =  2.03248699406565E+02, relative difference = -6.99185517875101E-17
   Diagnostics:    0.00000000      sec.
   chi2_B: 2.032E+02,  chi2_M: 5.265E+11
   max(B_n): 9.903E-01,  max(M): 2.133E+07,  min(M): 2.082E+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:    1.00000005E-03  sec.
   2 methods of computing chi2_M that should agree: method 1 =  2.28912602754667E+10, method 2 =  2.28912602754668E+10, relative difference = -9.99865595791594E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.06421654702505E+02, method 2 =  2.06421654702505E+02, relative difference =  0.00000000000000E+00
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.064E+02,  chi2_M: 2.289E+10
   max(B_n): 9.980E-01,  max(M): 4.447E+06,  min(M): 4.342E+06
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 8.859E-11 ( 37 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  9.83212216138827E+08, method 2 =  9.83212216138829E+08, relative difference = -7.88090678046121E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.07087456777410E+02, method 2 =  2.07087456777409E+02, relative difference =  7.54847272257749E-16
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.071E+02,  chi2_M: 9.832E+08
   max(B_n): 9.996E-01,  max(M): 9.216E+05,  min(M): 9.000E+05
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 4.281E-10 ( 38 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  4.21226696043964E+07, method 2 =  4.21226696043964E+07, relative difference =  2.65317250790280E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.07225622623326E+02, method 2 =  2.07225622623326E+02, relative difference =  6.85767258667288E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.072E+02,  chi2_M: 4.212E+07
   max(B_n): 9.999E-01,  max(M): 1.907E+05,  min(M): 1.863E+05
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 2.069E-09 ( 39 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  1.80366045977676E+06, method 2 =  1.80366045977676E+06, relative difference = -7.09983152957060E-16
   2 methods of computing chi2_B that should agree: method 1 =  2.07254228310585E+02, method 2 =  2.07254228310585E+02, relative difference = -6.85672607552597E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.073E+02,  chi2_M: 1.804E+06
   max(B_n): 1.000E+00,  max(M): 3.947E+04,  min(M): 3.855E+04
 Volume of magnetization region:  1.18E+00 meters^3.
 Solving system for lambda= 1.000E-08 ( 40 of  40)
   Additions:    0.00000000      sec.
   DSYSV:    0.00000000      sec.
   2 methods of computing chi2_M that should agree: method 1 =  7.72229053462690E+04, method 2 =  7.72229053462690E+04, relative difference =  0.00000000000000E+00
   2 methods of computing chi2_B that should agree: method 1 =  2.07260147961892E+02, method 2 =  2.07260147961892E+02, relative difference =  6.85653023745543E-17
   Diagnostics:    1.00000005E-03  sec.
   chi2_B: 2.073E+02,  chi2_M: 7.722E+04
   max(B_n): 1.000E+00,  max(M): 8.167E+03,  min(M): 7.976E+03
 Volume of magnetization region:  1.18E+00 meters^3.
 REGCOIL_PM complete. Total time=   1.7239999771118164      sec.
 You can run regcoilPlot regcoil_out.benchmark.nc to plot results.
