From 76032eb1c3c6401474684451b839564adee16e4b Mon Sep 17 00:00:00 2001 From: Nuwan Yapa Date: Mon, 17 Jun 2024 18:18:34 -0400 Subject: [PATCH] Test against Pb208 potential --- Project.toml | 1 + test/Pb208.jl | 22 ++ test/Pb208FldsFSUGarnet.csv | 401 ++++++++++++++++++++++++++++++++++++ test/dirac.jl | 3 - 4 files changed, 424 insertions(+), 3 deletions(-) create mode 100644 test/Pb208.jl create mode 100644 test/Pb208FldsFSUGarnet.csv delete mode 100644 test/dirac.jl diff --git a/Project.toml b/Project.toml index 03f6800..0f0d325 100644 --- a/Project.toml +++ b/Project.toml @@ -1,2 +1,3 @@ [deps] DifferentialEquations = "0c46a032-eb83-5123-abaf-570d42b7fbaa" +Interpolations = "a98d9a8b-a2ab-59e6-89dd-64a1c18fca59" diff --git a/test/Pb208.jl b/test/Pb208.jl new file mode 100644 index 0000000..a948354 --- /dev/null +++ b/test/Pb208.jl @@ -0,0 +1,22 @@ +using DelimitedFiles, Interpolations, Plots +include("../dirac.jl") + +# test data generated from Hartree.f +# format: x S(x) V(x) R(x) A(x) +test_data = readdlm("test/Pb208FldsFSUGarnet.csv") +xs = test_data[:, 1] +Ss = test_data[:, 2] +Vs = test_data[:, 3] +Rs = test_data[:, 4] +As = test_data[:, 5] + +S_interp = linear_interpolation(xs, Ss) +V_interp = linear_interpolation(xs, Vs) +R_interp = linear_interpolation(xs, Rs) +A_interp = linear_interpolation(xs, As) + +Es = collect(840:0.5:940) +boundaryVals = [boundaryValue(-1, M_n, E, S_interp, V_interp, maximum(xs))^2 for E in Es] + +plot(Es, boundaryVals, yscale=:log10, label="g(r_max)^2") +xlabel!("E (MeV)") diff --git a/test/Pb208FldsFSUGarnet.csv b/test/Pb208FldsFSUGarnet.csv new file mode 100644 index 0000000..bf7f938 --- /dev/null +++ b/test/Pb208FldsFSUGarnet.csv @@ -0,0 +1,401 @@ + 0.000 0.41122E+03 0.33426E+03 -0.28674E+01 0.20447E-12 + 0.050 0.41121E+03 0.33425E+03 -0.28729E+01 0.25728E+02 + 0.100 0.41117E+03 0.33422E+03 -0.28897E+01 0.25727E+02 + 0.150 0.41111E+03 0.33417E+03 -0.29175E+01 0.25724E+02 + 0.200 0.41103E+03 0.33410E+03 -0.29560E+01 0.25720E+02 + 0.250 0.41092E+03 0.33401E+03 -0.30049E+01 0.25715E+02 + 0.300 0.41078E+03 0.33390E+03 -0.30638E+01 0.25709E+02 + 0.350 0.41062E+03 0.33376E+03 -0.31319E+01 0.25702E+02 + 0.400 0.41043E+03 0.33360E+03 -0.32088E+01 0.25694E+02 + 0.450 0.41021E+03 0.33341E+03 -0.32935E+01 0.25685E+02 + 0.500 0.40996E+03 0.33320E+03 -0.33854E+01 0.25675E+02 + 0.550 0.40968E+03 0.33296E+03 -0.34834E+01 0.25664E+02 + 0.600 0.40937E+03 0.33269E+03 -0.35868E+01 0.25652E+02 + 0.650 0.40903E+03 0.33240E+03 -0.36946E+01 0.25638E+02 + 0.700 0.40865E+03 0.33208E+03 -0.38057E+01 0.25624E+02 + 0.750 0.40824E+03 0.33172E+03 -0.39191E+01 0.25609E+02 + 0.800 0.40780E+03 0.33134E+03 -0.40338E+01 0.25593E+02 + 0.850 0.40733E+03 0.33093E+03 -0.41487E+01 0.25576E+02 + 0.900 0.40683E+03 0.33050E+03 -0.42630E+01 0.25557E+02 + 0.950 0.40629E+03 0.33003E+03 -0.43755E+01 0.25538E+02 + 1.000 0.40573E+03 0.32954E+03 -0.44854E+01 0.25519E+02 + 1.050 0.40514E+03 0.32903E+03 -0.45916E+01 0.25498E+02 + 1.100 0.40453E+03 0.32849E+03 -0.46934E+01 0.25476E+02 + 1.150 0.40389E+03 0.32794E+03 -0.47900E+01 0.25453E+02 + 1.200 0.40324E+03 0.32737E+03 -0.48805E+01 0.25430E+02 + 1.250 0.40257E+03 0.32678E+03 -0.49644E+01 0.25405E+02 + 1.300 0.40189E+03 0.32619E+03 -0.50411E+01 0.25380E+02 + 1.350 0.40120E+03 0.32558E+03 -0.51101E+01 0.25354E+02 + 1.400 0.40051E+03 0.32498E+03 -0.51710E+01 0.25327E+02 + 1.450 0.39982E+03 0.32437E+03 -0.52235E+01 0.25299E+02 + 1.500 0.39913E+03 0.32377E+03 -0.52675E+01 0.25271E+02 + 1.550 0.39846E+03 0.32318E+03 -0.53027E+01 0.25241E+02 + 1.600 0.39780E+03 0.32260E+03 -0.53293E+01 0.25211E+02 + 1.650 0.39716E+03 0.32204E+03 -0.53474E+01 0.25180E+02 + 1.700 0.39655E+03 0.32150E+03 -0.53572E+01 0.25148E+02 + 1.750 0.39596E+03 0.32099E+03 -0.53590E+01 0.25115E+02 + 1.800 0.39541E+03 0.32050E+03 -0.53531E+01 0.25082E+02 + 1.850 0.39489E+03 0.32005E+03 -0.53402E+01 0.25047E+02 + 1.900 0.39441E+03 0.31963E+03 -0.53207E+01 0.25012E+02 + 1.950 0.39398E+03 0.31925E+03 -0.52953E+01 0.24976E+02 + 2.000 0.39359E+03 0.31890E+03 -0.52646E+01 0.24939E+02 + 2.050 0.39325E+03 0.31860E+03 -0.52295E+01 0.24902E+02 + 2.100 0.39296E+03 0.31834E+03 -0.51906E+01 0.24863E+02 + 2.150 0.39272E+03 0.31813E+03 -0.51487E+01 0.24824E+02 + 2.200 0.39253E+03 0.31796E+03 -0.51047E+01 0.24784E+02 + 2.250 0.39239E+03 0.31783E+03 -0.50592E+01 0.24743E+02 + 2.300 0.39230E+03 0.31775E+03 -0.50132E+01 0.24701E+02 + 2.350 0.39225E+03 0.31771E+03 -0.49673E+01 0.24658E+02 + 2.400 0.39226E+03 0.31771E+03 -0.49222E+01 0.24614E+02 + 2.450 0.39231E+03 0.31774E+03 -0.48785E+01 0.24569E+02 + 2.500 0.39240E+03 0.31782E+03 -0.48369E+01 0.24524E+02 + 2.550 0.39254E+03 0.31792E+03 -0.47979E+01 0.24477E+02 + 2.600 0.39270E+03 0.31805E+03 -0.47619E+01 0.24430E+02 + 2.650 0.39290E+03 0.31822E+03 -0.47294E+01 0.24382E+02 + 2.700 0.39313E+03 0.31840E+03 -0.47005E+01 0.24333E+02 + 2.750 0.39338E+03 0.31860E+03 -0.46756E+01 0.24282E+02 + 2.800 0.39365E+03 0.31882E+03 -0.46547E+01 0.24231E+02 + 2.850 0.39394E+03 0.31905E+03 -0.46379E+01 0.24179E+02 + 2.900 0.39423E+03 0.31928E+03 -0.46253E+01 0.24126E+02 + 2.950 0.39454E+03 0.31952E+03 -0.46166E+01 0.24072E+02 + 3.000 0.39484E+03 0.31976E+03 -0.46118E+01 0.24017E+02 + 3.050 0.39515E+03 0.31999E+03 -0.46106E+01 0.23961E+02 + 3.100 0.39545E+03 0.32022E+03 -0.46128E+01 0.23905E+02 + 3.150 0.39575E+03 0.32045E+03 -0.46179E+01 0.23847E+02 + 3.200 0.39603E+03 0.32066E+03 -0.46257E+01 0.23788E+02 + 3.250 0.39631E+03 0.32085E+03 -0.46357E+01 0.23728E+02 + 3.300 0.39657E+03 0.32104E+03 -0.46475E+01 0.23668E+02 + 3.350 0.39681E+03 0.32120E+03 -0.46606E+01 0.23606E+02 + 3.400 0.39704E+03 0.32136E+03 -0.46746E+01 0.23543E+02 + 3.450 0.39726E+03 0.32149E+03 -0.46890E+01 0.23480E+02 + 3.500 0.39746E+03 0.32162E+03 -0.47033E+01 0.23415E+02 + 3.550 0.39764E+03 0.32173E+03 -0.47171E+01 0.23349E+02 + 3.600 0.39781E+03 0.32182E+03 -0.47300E+01 0.23283E+02 + 3.650 0.39797E+03 0.32190E+03 -0.47415E+01 0.23215E+02 + 3.700 0.39811E+03 0.32198E+03 -0.47513E+01 0.23147E+02 + 3.750 0.39825E+03 0.32204E+03 -0.47590E+01 0.23078E+02 + 3.800 0.39838E+03 0.32210E+03 -0.47644E+01 0.23007E+02 + 3.850 0.39849E+03 0.32215E+03 -0.47672E+01 0.22936E+02 + 3.900 0.39861E+03 0.32220E+03 -0.47672E+01 0.22864E+02 + 3.950 0.39871E+03 0.32225E+03 -0.47643E+01 0.22791E+02 + 4.000 0.39881E+03 0.32229E+03 -0.47585E+01 0.22717E+02 + 4.050 0.39891E+03 0.32234E+03 -0.47497E+01 0.22642E+02 + 4.100 0.39900E+03 0.32239E+03 -0.47379E+01 0.22566E+02 + 4.150 0.39909E+03 0.32243E+03 -0.47233E+01 0.22489E+02 + 4.200 0.39917E+03 0.32248E+03 -0.47060E+01 0.22411E+02 + 4.250 0.39923E+03 0.32252E+03 -0.46862E+01 0.22332E+02 + 4.300 0.39928E+03 0.32256E+03 -0.46641E+01 0.22252E+02 + 4.350 0.39931E+03 0.32259E+03 -0.46401E+01 0.22171E+02 + 4.400 0.39932E+03 0.32261E+03 -0.46144E+01 0.22090E+02 + 4.450 0.39929E+03 0.32261E+03 -0.45875E+01 0.22007E+02 + 4.500 0.39923E+03 0.32258E+03 -0.45597E+01 0.21923E+02 + 4.550 0.39911E+03 0.32252E+03 -0.45315E+01 0.21839E+02 + 4.600 0.39894E+03 0.32243E+03 -0.45032E+01 0.21753E+02 + 4.650 0.39869E+03 0.32229E+03 -0.44752E+01 0.21666E+02 + 4.700 0.39837E+03 0.32209E+03 -0.44480E+01 0.21579E+02 + 4.750 0.39795E+03 0.32183E+03 -0.44220E+01 0.21490E+02 + 4.800 0.39743E+03 0.32148E+03 -0.43975E+01 0.21401E+02 + 4.850 0.39679E+03 0.32105E+03 -0.43751E+01 0.21310E+02 + 4.900 0.39601E+03 0.32052E+03 -0.43549E+01 0.21218E+02 + 4.950 0.39509E+03 0.31987E+03 -0.43374E+01 0.21126E+02 + 5.000 0.39400E+03 0.31910E+03 -0.43229E+01 0.21032E+02 + 5.050 0.39274E+03 0.31819E+03 -0.43116E+01 0.20937E+02 + 5.100 0.39128E+03 0.31713E+03 -0.43039E+01 0.20842E+02 + 5.150 0.38962E+03 0.31591E+03 -0.42998E+01 0.20745E+02 + 5.200 0.38774E+03 0.31450E+03 -0.42997E+01 0.20647E+02 + 5.250 0.38563E+03 0.31292E+03 -0.43037E+01 0.20549E+02 + 5.300 0.38328E+03 0.31113E+03 -0.43119E+01 0.20449E+02 + 5.350 0.38067E+03 0.30913E+03 -0.43244E+01 0.20349E+02 + 5.400 0.37779E+03 0.30691E+03 -0.43413E+01 0.20247E+02 + 5.450 0.37464E+03 0.30446E+03 -0.43626E+01 0.20145E+02 + 5.500 0.37121E+03 0.30177E+03 -0.43885E+01 0.20042E+02 + 5.550 0.36750E+03 0.29884E+03 -0.44188E+01 0.19938E+02 + 5.600 0.36349E+03 0.29565E+03 -0.44536E+01 0.19833E+02 + 5.650 0.35919E+03 0.29221E+03 -0.44927E+01 0.19728E+02 + 5.700 0.35459E+03 0.28852E+03 -0.45362E+01 0.19621E+02 + 5.750 0.34969E+03 0.28455E+03 -0.45839E+01 0.19514E+02 + 5.800 0.34451E+03 0.28033E+03 -0.46356E+01 0.19406E+02 + 5.850 0.33903E+03 0.27585E+03 -0.46913E+01 0.19298E+02 + 5.900 0.33327E+03 0.27111E+03 -0.47506E+01 0.19189E+02 + 5.950 0.32724E+03 0.26612E+03 -0.48132E+01 0.19079E+02 + 6.000 0.32094E+03 0.26088E+03 -0.48790E+01 0.18969E+02 + 6.050 0.31439E+03 0.25540E+03 -0.49475E+01 0.18859E+02 + 6.100 0.30759E+03 0.24970E+03 -0.50182E+01 0.18748E+02 + 6.150 0.30056E+03 0.24377E+03 -0.50908E+01 0.18637E+02 + 6.200 0.29333E+03 0.23765E+03 -0.51645E+01 0.18525E+02 + 6.250 0.28589E+03 0.23133E+03 -0.52388E+01 0.18413E+02 + 6.300 0.27828E+03 0.22484E+03 -0.53130E+01 0.18301E+02 + 6.350 0.27051E+03 0.21819E+03 -0.53862E+01 0.18189E+02 + 6.400 0.26260E+03 0.21141E+03 -0.54576E+01 0.18077E+02 + 6.450 0.25457E+03 0.20450E+03 -0.55263E+01 0.17965E+02 + 6.500 0.24645E+03 0.19750E+03 -0.55912E+01 0.17853E+02 + 6.550 0.23825E+03 0.19043E+03 -0.56514E+01 0.17741E+02 + 6.600 0.23000E+03 0.18330E+03 -0.57056E+01 0.17629E+02 + 6.650 0.22172E+03 0.17613E+03 -0.57530E+01 0.17518E+02 + 6.700 0.21343E+03 0.16896E+03 -0.57923E+01 0.17406E+02 + 6.750 0.20515E+03 0.16180E+03 -0.58226E+01 0.17295E+02 + 6.800 0.19691E+03 0.15468E+03 -0.58429E+01 0.17185E+02 + 6.850 0.18873E+03 0.14761E+03 -0.58522E+01 0.17075E+02 + 6.900 0.18063E+03 0.14062E+03 -0.58498E+01 0.16965E+02 + 6.950 0.17262E+03 0.13372E+03 -0.58351E+01 0.16856E+02 + 7.000 0.16473E+03 0.12695E+03 -0.58074E+01 0.16747E+02 + 7.050 0.15698E+03 0.12030E+03 -0.57666E+01 0.16639E+02 + 7.100 0.14937E+03 0.11381E+03 -0.57124E+01 0.16532E+02 + 7.150 0.14193E+03 0.10749E+03 -0.56449E+01 0.16426E+02 + 7.200 0.13468E+03 0.10134E+03 -0.55643E+01 0.16320E+02 + 7.250 0.12761E+03 0.95389E+02 -0.54710E+01 0.16215E+02 + 7.300 0.12075E+03 0.89637E+02 -0.53655E+01 0.16111E+02 + 7.350 0.11410E+03 0.84095E+02 -0.52487E+01 0.16007E+02 + 7.400 0.10767E+03 0.78769E+02 -0.51212E+01 0.15905E+02 + 7.450 0.10147E+03 0.73666E+02 -0.49842E+01 0.15803E+02 + 7.500 0.95503E+02 0.68788E+02 -0.48387E+01 0.15702E+02 + 7.550 0.89772E+02 0.64137E+02 -0.46858E+01 0.15602E+02 + 7.600 0.84280E+02 0.59713E+02 -0.45267E+01 0.15503E+02 + 7.650 0.79027E+02 0.55517E+02 -0.43626E+01 0.15405E+02 + 7.700 0.74013E+02 0.51545E+02 -0.41946E+01 0.15308E+02 + 7.750 0.69238E+02 0.47793E+02 -0.40240E+01 0.15212E+02 + 7.800 0.64697E+02 0.44259E+02 -0.38519E+01 0.15117E+02 + 7.850 0.60388E+02 0.40935E+02 -0.36793E+01 0.15023E+02 + 7.900 0.56306E+02 0.37817E+02 -0.35072E+01 0.14930E+02 + 7.950 0.52446E+02 0.34897E+02 -0.33366E+01 0.14838E+02 + 8.000 0.48802E+02 0.32167E+02 -0.31682E+01 0.14747E+02 + 8.050 0.45368E+02 0.29621E+02 -0.30029E+01 0.14656E+02 + 8.100 0.42136E+02 0.27250E+02 -0.28412E+01 0.14567E+02 + 8.150 0.39100E+02 0.25046E+02 -0.26837E+01 0.14479E+02 + 8.200 0.36252E+02 0.23000E+02 -0.25310E+01 0.14392E+02 + 8.250 0.33584E+02 0.21104E+02 -0.23833E+01 0.14305E+02 + 8.300 0.31088E+02 0.19349E+02 -0.22409E+01 0.14220E+02 + 8.350 0.28755E+02 0.17728E+02 -0.21042E+01 0.14135E+02 + 8.400 0.26579E+02 0.16232E+02 -0.19733E+01 0.14052E+02 + 8.450 0.24551E+02 0.14853E+02 -0.18482E+01 0.13969E+02 + 8.500 0.22663E+02 0.13584E+02 -0.17291E+01 0.13888E+02 + 8.550 0.20907E+02 0.12416E+02 -0.16158E+01 0.13807E+02 + 8.600 0.19276E+02 0.11344E+02 -0.15084E+01 0.13727E+02 + 8.650 0.17762E+02 0.10360E+02 -0.14068E+01 0.13648E+02 + 8.700 0.16359E+02 0.94573E+01 -0.13108E+01 0.13570E+02 + 8.750 0.15059E+02 0.86306E+01 -0.12203E+01 0.13493E+02 + 8.800 0.13856E+02 0.78739E+01 -0.11351E+01 0.13416E+02 + 8.850 0.12743E+02 0.71817E+01 -0.10551E+01 0.13340E+02 + 8.900 0.11716E+02 0.65488E+01 -0.98000E+00 0.13266E+02 + 8.950 0.10767E+02 0.59707E+01 -0.90967E+00 0.13192E+02 + 9.000 0.98910E+01 0.54428E+01 -0.84387E+00 0.13119E+02 + 9.050 0.90836E+01 0.49610E+01 -0.78238E+00 0.13046E+02 + 9.100 0.83396E+01 0.45215E+01 -0.72498E+00 0.12975E+02 + 9.150 0.76544E+01 0.41208E+01 -0.67147E+00 0.12904E+02 + 9.200 0.70237E+01 0.37554E+01 -0.62163E+00 0.12834E+02 + 9.250 0.64434E+01 0.34225E+01 -0.57525E+00 0.12765E+02 + 9.300 0.59099E+01 0.31192E+01 -0.53212E+00 0.12696E+02 + 9.350 0.54194E+01 0.28429E+01 -0.49206E+00 0.12628E+02 + 9.400 0.49689E+01 0.25913E+01 -0.45487E+00 0.12561E+02 + 9.450 0.45551E+01 0.23623E+01 -0.42037E+00 0.12495E+02 + 9.500 0.41752E+01 0.21537E+01 -0.38838E+00 0.12429E+02 + 9.550 0.38265E+01 0.19638E+01 -0.35874E+00 0.12364E+02 + 9.600 0.35066E+01 0.17909E+01 -0.33130E+00 0.12300E+02 + 9.650 0.32132E+01 0.16336E+01 -0.30590E+00 0.12236E+02 + 9.700 0.29441E+01 0.14903E+01 -0.28239E+00 0.12173E+02 + 9.750 0.26975E+01 0.13599E+01 -0.26066E+00 0.12110E+02 + 9.800 0.24714E+01 0.12411E+01 -0.24057E+00 0.12049E+02 + 9.850 0.22642E+01 0.11330E+01 -0.22201E+00 0.11988E+02 + 9.900 0.20743E+01 0.10345E+01 -0.20486E+00 0.11927E+02 + 9.950 0.19004E+01 0.94486E+00 -0.18902E+00 0.11867E+02 + 10.000 0.17411E+01 0.86318E+00 -0.17439E+00 0.11808E+02 + 10.050 0.15952E+01 0.78876E+00 -0.16089E+00 0.11749E+02 + 10.100 0.14616E+01 0.72094E+00 -0.14843E+00 0.11691E+02 + 10.150 0.13392E+01 0.65914E+00 -0.13694E+00 0.11633E+02 + 10.200 0.12272E+01 0.60279E+00 -0.12633E+00 0.11576E+02 + 10.250 0.11246E+01 0.55142E+00 -0.11655E+00 0.11520E+02 + 10.300 0.10306E+01 0.50457E+00 -0.10752E+00 0.11464E+02 + 10.350 0.94461E+00 0.46183E+00 -0.99196E-01 0.11409E+02 + 10.400 0.86586E+00 0.42283E+00 -0.91520E-01 0.11354E+02 + 10.450 0.79375E+00 0.38723E+00 -0.84441E-01 0.11299E+02 + 10.500 0.72773E+00 0.35473E+00 -0.77913E-01 0.11246E+02 + 10.550 0.66728E+00 0.32505E+00 -0.71894E-01 0.11192E+02 + 10.600 0.61193E+00 0.29795E+00 -0.66344E-01 0.11139E+02 + 10.650 0.56124E+00 0.27318E+00 -0.61227E-01 0.11087E+02 + 10.700 0.51481E+00 0.25054E+00 -0.56508E-01 0.11035E+02 + 10.750 0.47230E+00 0.22984E+00 -0.52157E-01 0.10984E+02 + 10.800 0.43336E+00 0.21091E+00 -0.48145E-01 0.10933E+02 + 10.850 0.39768E+00 0.19359E+00 -0.44445E-01 0.10883E+02 + 10.900 0.36500E+00 0.17775E+00 -0.41034E-01 0.10833E+02 + 10.950 0.33506E+00 0.16325E+00 -0.37887E-01 0.10783E+02 + 11.000 0.30763E+00 0.14997E+00 -0.34986E-01 0.10734E+02 + 11.050 0.28249E+00 0.13781E+00 -0.32309E-01 0.10686E+02 + 11.100 0.25944E+00 0.12666E+00 -0.29841E-01 0.10638E+02 + 11.150 0.23832E+00 0.11645E+00 -0.27563E-01 0.10590E+02 + 11.200 0.21895E+00 0.10709E+00 -0.25462E-01 0.10543E+02 + 11.250 0.20120E+00 0.98512E-01 -0.23524E-01 0.10496E+02 + 11.300 0.18491E+00 0.90640E-01 -0.21736E-01 0.10449E+02 + 11.350 0.16998E+00 0.83418E-01 -0.20085E-01 0.10403E+02 + 11.400 0.15628E+00 0.76790E-01 -0.18562E-01 0.10358E+02 + 11.450 0.14370E+00 0.70706E-01 -0.17157E-01 0.10313E+02 + 11.500 0.13217E+00 0.65120E-01 -0.15859E-01 0.10268E+02 + 11.550 0.12158E+00 0.59989E-01 -0.14662E-01 0.10223E+02 + 11.600 0.11186E+00 0.55274E-01 -0.13556E-01 0.10179E+02 + 11.650 0.10294E+00 0.50942E-01 -0.12535E-01 0.10136E+02 + 11.700 0.94749E-01 0.46960E-01 -0.11592E-01 0.10092E+02 + 11.750 0.87224E-01 0.43298E-01 -0.10721E-01 0.10049E+02 + 11.800 0.80311E-01 0.39931E-01 -0.99169E-02 0.10007E+02 + 11.850 0.73960E-01 0.36833E-01 -0.91740E-02 0.99645E+01 + 11.900 0.68124E-01 0.33982E-01 -0.84877E-02 0.99226E+01 + 11.950 0.62760E-01 0.31359E-01 -0.78535E-02 0.98811E+01 + 12.000 0.57829E-01 0.28943E-01 -0.72676E-02 0.98399E+01 + 12.050 0.53295E-01 0.26719E-01 -0.67260E-02 0.97991E+01 + 12.100 0.49126E-01 0.24671E-01 -0.62255E-02 0.97586E+01 + 12.150 0.45291E-01 0.22784E-01 -0.57629E-02 0.97184E+01 + 12.200 0.41762E-01 0.21045E-01 -0.53352E-02 0.96786E+01 + 12.250 0.38516E-01 0.19442E-01 -0.49397E-02 0.96391E+01 + 12.300 0.35528E-01 0.17965E-01 -0.45741E-02 0.95999E+01 + 12.350 0.32778E-01 0.16603E-01 -0.42359E-02 0.95610E+01 + 12.400 0.30246E-01 0.15346E-01 -0.39232E-02 0.95225E+01 + 12.450 0.27915E-01 0.14187E-01 -0.36339E-02 0.94842E+01 + 12.500 0.25767E-01 0.13118E-01 -0.33663E-02 0.94463E+01 + 12.550 0.23789E-01 0.12132E-01 -0.31187E-02 0.94087E+01 + 12.600 0.21967E-01 0.11221E-01 -0.28896E-02 0.93713E+01 + 12.650 0.20287E-01 0.10380E-01 -0.26776E-02 0.93343E+01 + 12.700 0.18739E-01 0.96040E-02 -0.24814E-02 0.92976E+01 + 12.750 0.17312E-01 0.88871E-02 -0.22998E-02 0.92611E+01 + 12.800 0.15996E-01 0.82250E-02 -0.21317E-02 0.92249E+01 + 12.850 0.14783E-01 0.76132E-02 -0.19760E-02 0.91890E+01 + 12.900 0.13664E-01 0.70480E-02 -0.18319E-02 0.91534E+01 + 12.950 0.12631E-01 0.65257E-02 -0.16985E-02 0.91181E+01 + 13.000 0.11678E-01 0.60428E-02 -0.15749E-02 0.90830E+01 + 13.050 0.10799E-01 0.55965E-02 -0.14605E-02 0.90482E+01 + 13.100 0.99878E-02 0.51838E-02 -0.13545E-02 0.90137E+01 + 13.150 0.92386E-02 0.48022E-02 -0.12562E-02 0.89794E+01 + 13.200 0.85470E-02 0.44492E-02 -0.11653E-02 0.89454E+01 + 13.250 0.79082E-02 0.41227E-02 -0.10810E-02 0.89116E+01 + 13.300 0.73183E-02 0.38206E-02 -0.10029E-02 0.88781E+01 + 13.350 0.67733E-02 0.35410E-02 -0.93047E-03 0.88449E+01 + 13.400 0.62697E-02 0.32823E-02 -0.86337E-03 0.88119E+01 + 13.450 0.58044E-02 0.30429E-02 -0.80119E-03 0.87791E+01 + 13.500 0.53744E-02 0.28213E-02 -0.74354E-03 0.87466E+01 + 13.550 0.49769E-02 0.26160E-02 -0.69010E-03 0.87143E+01 + 13.600 0.46094E-02 0.24260E-02 -0.64055E-03 0.86823E+01 + 13.650 0.42695E-02 0.22500E-02 -0.59460E-03 0.86505E+01 + 13.700 0.39553E-02 0.20871E-02 -0.55199E-03 0.86189E+01 + 13.750 0.36646E-02 0.19361E-02 -0.51248E-03 0.85876E+01 + 13.800 0.33957E-02 0.17962E-02 -0.47583E-03 0.85564E+01 + 13.850 0.31470E-02 0.16666E-02 -0.44184E-03 0.85256E+01 + 13.900 0.29168E-02 0.15465E-02 -0.41031E-03 0.84949E+01 + 13.950 0.27037E-02 0.14352E-02 -0.38105E-03 0.84644E+01 + 14.000 0.25065E-02 0.13320E-02 -0.35391E-03 0.84342E+01 + 14.050 0.23240E-02 0.12364E-02 -0.32872E-03 0.84042E+01 + 14.100 0.21550E-02 0.11478E-02 -0.30535E-03 0.83744E+01 + 14.150 0.19985E-02 0.10656E-02 -0.28366E-03 0.83448E+01 + 14.200 0.18536E-02 0.98933E-03 -0.26353E-03 0.83154E+01 + 14.250 0.17194E-02 0.91864E-03 -0.24485E-03 0.82862E+01 + 14.300 0.15950E-02 0.85307E-03 -0.22750E-03 0.82573E+01 + 14.350 0.14798E-02 0.79225E-03 -0.21140E-03 0.82285E+01 + 14.400 0.13731E-02 0.73583E-03 -0.19645E-03 0.81999E+01 + 14.450 0.12742E-02 0.68348E-03 -0.18257E-03 0.81715E+01 + 14.500 0.11826E-02 0.63491E-03 -0.16968E-03 0.81434E+01 + 14.550 0.10976E-02 0.58984E-03 -0.15772E-03 0.81154E+01 + 14.600 0.10188E-02 0.54802E-03 -0.14660E-03 0.80876E+01 + 14.650 0.94582E-03 0.50919E-03 -0.13628E-03 0.80600E+01 + 14.700 0.87813E-03 0.47316E-03 -0.12669E-03 0.80326E+01 + 14.750 0.81535E-03 0.43971E-03 -0.11779E-03 0.80053E+01 + 14.800 0.75713E-03 0.40865E-03 -0.10952E-03 0.79783E+01 + 14.850 0.70313E-03 0.37982E-03 -0.10183E-03 0.79514E+01 + 14.900 0.65304E-03 0.35304E-03 -0.94690E-04 0.79248E+01 + 14.950 0.60657E-03 0.32818E-03 -0.88056E-04 0.78983E+01 + 15.000 0.56345E-03 0.30509E-03 -0.81891E-04 0.78719E+01 + 15.050 0.52345E-03 0.28365E-03 -0.76163E-04 0.78458E+01 + 15.100 0.48632E-03 0.26373E-03 -0.70839E-04 0.78198E+01 + 15.150 0.45187E-03 0.24522E-03 -0.65891E-04 0.77940E+01 + 15.200 0.41989E-03 0.22803E-03 -0.61292E-04 0.77683E+01 + 15.250 0.39020E-03 0.21206E-03 -0.57018E-04 0.77429E+01 + 15.300 0.36265E-03 0.19722E-03 -0.53044E-04 0.77176E+01 + 15.350 0.33706E-03 0.18343E-03 -0.49350E-04 0.76924E+01 + 15.400 0.31330E-03 0.17061E-03 -0.45916E-04 0.76675E+01 + 15.450 0.29124E-03 0.15870E-03 -0.42723E-04 0.76426E+01 + 15.500 0.27076E-03 0.14764E-03 -0.39755E-04 0.76180E+01 + 15.550 0.25173E-03 0.13735E-03 -0.36994E-04 0.75935E+01 + 15.600 0.23406E-03 0.12778E-03 -0.34427E-04 0.75692E+01 + 15.650 0.21764E-03 0.11889E-03 -0.32040E-04 0.75450E+01 + 15.700 0.20239E-03 0.11063E-03 -0.29820E-04 0.75209E+01 + 15.750 0.18822E-03 0.10294E-03 -0.27755E-04 0.74971E+01 + 15.800 0.17505E-03 0.95796E-04 -0.25835E-04 0.74733E+01 + 15.850 0.16282E-03 0.89151E-04 -0.24048E-04 0.74498E+01 + 15.900 0.15145E-03 0.82972E-04 -0.22386E-04 0.74263E+01 + 15.950 0.14089E-03 0.77226E-04 -0.20840E-04 0.74031E+01 + 16.000 0.13107E-03 0.71881E-04 -0.19402E-04 0.73799E+01 + 16.050 0.12194E-03 0.66910E-04 -0.18064E-04 0.73569E+01 + 16.100 0.11345E-03 0.62287E-04 -0.16819E-04 0.73341E+01 + 16.150 0.10556E-03 0.57985E-04 -0.15661E-04 0.73114E+01 + 16.200 0.98230E-04 0.53984E-04 -0.14583E-04 0.72888E+01 + 16.250 0.91412E-04 0.50261E-04 -0.13579E-04 0.72664E+01 + 16.300 0.85071E-04 0.46798E-04 -0.12646E-04 0.72441E+01 + 16.350 0.79176E-04 0.43575E-04 -0.11777E-04 0.72220E+01 + 16.400 0.73693E-04 0.40576E-04 -0.10968E-04 0.71999E+01 + 16.450 0.68593E-04 0.37786E-04 -0.10216E-04 0.71780E+01 + 16.500 0.63851E-04 0.35189E-04 -0.95151E-05 0.71563E+01 + 16.550 0.59439E-04 0.32772E-04 -0.88629E-05 0.71347E+01 + 16.600 0.55335E-04 0.30523E-04 -0.82558E-05 0.71132E+01 + 16.650 0.51518E-04 0.28429E-04 -0.76907E-05 0.70918E+01 + 16.700 0.47966E-04 0.26480E-04 -0.71645E-05 0.70706E+01 + 16.750 0.44662E-04 0.24666E-04 -0.66746E-05 0.70495E+01 + 16.800 0.41587E-04 0.22978E-04 -0.62185E-05 0.70285E+01 + 16.850 0.38726E-04 0.21406E-04 -0.57938E-05 0.70076E+01 + 16.900 0.36064E-04 0.19942E-04 -0.53983E-05 0.69869E+01 + 16.950 0.33587E-04 0.18579E-04 -0.50300E-05 0.69663E+01 + 17.000 0.31281E-04 0.17310E-04 -0.46870E-05 0.69458E+01 + 17.050 0.29135E-04 0.16129E-04 -0.43676E-05 0.69254E+01 + 17.100 0.27137E-04 0.15029E-04 -0.40702E-05 0.69052E+01 + 17.150 0.25278E-04 0.14004E-04 -0.37931E-05 0.68851E+01 + 17.200 0.23547E-04 0.13050E-04 -0.35351E-05 0.68651E+01 + 17.250 0.21936E-04 0.12162E-04 -0.32947E-05 0.68452E+01 + 17.300 0.20436E-04 0.11334E-04 -0.30708E-05 0.68254E+01 + 17.350 0.19040E-04 0.10563E-04 -0.28622E-05 0.68057E+01 + 17.400 0.17740E-04 0.98451E-05 -0.26679E-05 0.67861E+01 + 17.450 0.16529E-04 0.91762E-05 -0.24869E-05 0.67667E+01 + 17.500 0.15401E-04 0.85532E-05 -0.23182E-05 0.67474E+01 + 17.550 0.14352E-04 0.79728E-05 -0.21611E-05 0.67281E+01 + 17.600 0.13374E-04 0.74321E-05 -0.20147E-05 0.67090E+01 + 17.650 0.12463E-04 0.69283E-05 -0.18783E-05 0.66900E+01 + 17.700 0.11615E-04 0.64589E-05 -0.17512E-05 0.66711E+01 + 17.750 0.10825E-04 0.60216E-05 -0.16327E-05 0.66523E+01 + 17.800 0.10090E-04 0.56141E-05 -0.15224E-05 0.66336E+01 + 17.850 0.94042E-05 0.52344E-05 -0.14195E-05 0.66151E+01 + 17.900 0.87658E-05 0.48805E-05 -0.13236E-05 0.65966E+01 + 17.950 0.81710E-05 0.45508E-05 -0.12343E-05 0.65782E+01 + 18.000 0.76169E-05 0.42435E-05 -0.11510E-05 0.65599E+01 + 18.050 0.71006E-05 0.39571E-05 -0.10734E-05 0.65418E+01 + 18.100 0.66195E-05 0.36901E-05 -0.10011E-05 0.65237E+01 + 18.150 0.61713E-05 0.34413E-05 -0.93361E-06 0.65057E+01 + 18.200 0.57536E-05 0.32094E-05 -0.87075E-06 0.64879E+01 + 18.250 0.53644E-05 0.29932E-05 -0.81215E-06 0.64701E+01 + 18.300 0.50016E-05 0.27917E-05 -0.75751E-06 0.64524E+01 + 18.350 0.46635E-05 0.26039E-05 -0.70658E-06 0.64348E+01 + 18.400 0.43484E-05 0.24287E-05 -0.65909E-06 0.64173E+01 + 18.450 0.40546E-05 0.22655E-05 -0.61481E-06 0.63999E+01 + 18.500 0.37808E-05 0.21132E-05 -0.57352E-06 0.63826E+01 + 18.550 0.35255E-05 0.19713E-05 -0.53502E-06 0.63654E+01 + 18.600 0.32874E-05 0.18389E-05 -0.49912E-06 0.63484E+01 + 18.650 0.30655E-05 0.17155E-05 -0.46563E-06 0.63313E+01 + 18.700 0.28584E-05 0.16004E-05 -0.43440E-06 0.63144E+01 + 18.750 0.26654E-05 0.14930E-05 -0.40527E-06 0.62975E+01 + 18.800 0.24852E-05 0.13929E-05 -0.37809E-06 0.62808E+01 + 18.850 0.23172E-05 0.12994E-05 -0.35274E-06 0.62641E+01 + 18.900 0.21603E-05 0.12122E-05 -0.32907E-06 0.62476E+01 + 18.950 0.20139E-05 0.11308E-05 -0.30699E-06 0.62311E+01 + 19.000 0.18772E-05 0.10549E-05 -0.28637E-06 0.62147E+01 + 19.050 0.17496E-05 0.98393E-06 -0.26711E-06 0.61984E+01 + 19.100 0.16302E-05 0.91766E-06 -0.24911E-06 0.61821E+01 + 19.150 0.15187E-05 0.85571E-06 -0.23229E-06 0.61660E+01 + 19.200 0.14144E-05 0.79776E-06 -0.21655E-06 0.61499E+01 + 19.250 0.13167E-05 0.74351E-06 -0.20182E-06 0.61340E+01 + 19.300 0.12252E-05 0.69266E-06 -0.18801E-06 0.61181E+01 + 19.350 0.11394E-05 0.64494E-06 -0.17504E-06 0.61023E+01 + 19.400 0.10589E-05 0.60007E-06 -0.16285E-06 0.60865E+01 + 19.450 0.98318E-06 0.55778E-06 -0.15136E-06 0.60709E+01 + 19.500 0.91189E-06 0.51780E-06 -0.14050E-06 0.60553E+01 + 19.550 0.84464E-06 0.47986E-06 -0.13019E-06 0.60398E+01 + 19.600 0.78104E-06 0.44369E-06 -0.12037E-06 0.60244E+01 + 19.650 0.72075E-06 0.40901E-06 -0.11095E-06 0.60091E+01 + 19.700 0.66340E-06 0.37550E-06 -0.10186E-06 0.59939E+01 + 19.750 0.60866E-06 0.34284E-06 -0.93007E-07 0.59787E+01 + 19.800 0.55618E-06 0.31067E-06 -0.84303E-07 0.59636E+01 + 19.850 0.50564E-06 0.27860E-06 -0.75642E-07 0.59486E+01 + 19.900 0.45667E-06 0.24617E-06 -0.66908E-07 0.59336E+01 + 19.950 0.40895E-06 0.21286E-06 -0.57970E-07 0.59187E+01 + 20.000 0.36213E-06 0.17807E-06 -0.48673E-07 0.59039E+01 diff --git a/test/dirac.jl b/test/dirac.jl deleted file mode 100644 index 8b35112..0000000 --- a/test/dirac.jl +++ /dev/null @@ -1,3 +0,0 @@ -include("../dirac.jl") - -boundaryValue(-1, M_n, 940, r -> 0, r -> 0, 20) \ No newline at end of file