Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Added basic unit tests #123

Merged
merged 2 commits into from
Jan 12, 2025
Merged

Added basic unit tests #123

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
8614c03
added basic APCdef test
lllx125 Jan 12, 2025
4cd9357
added unit tests of charges
lllx125 Jan 12, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
157 changes: 36 additions & 121 deletions test/runtests.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -4,128 +4,43 @@ using Test


@testset "AtomicAndPhysicalConstants.jl" begin
# Write your tests here.
# Write your tests here.
end

@testset "set_units.jl" begin

# units should be available to users, including "amu" and "e"
# if not availble, an error would be thrown and the test would fail.
@test u"kg" == u"kg"
@test u"km" == u"km"
@test u"amu" == u"amu"
@test u"eV/c^2" == u"eV/c^2"
@test u"e" == u"e"
@test u"C" == u"C"
@test u"s" == u"s"


#setunits to default units
setunits()

#when setunits is runned all constants should be initialized to default units

@test m_electron ≈ 0.51099895000e6
@test m_proton ≈ 0.93827208816e9
@test m_neutron ≈ 0.93956542052e9
@test m_muon ≈ 105.6583755e6
@test m_helion ≈ 2.808391607035771e9
@test m_deuteron ≈ 1.87561294257e9

@test c_light ≈ 2.99792458e8
@test e_charge ≈ 1
@test r_e ≈ 2.8179403262e-15
@test h_planck ≈ 4.135667696e-15
@test mu_0_vac ≈ 1.25663706212e-6
@test eps_0_vac ≈ 8.854187812800385e-12


@test kg_per_amu ≈ 1.66053906660e-27
@test eV_per_amu ≈ 9.3149410242e8
@test N_avogadro ≈ 6.02214076e23
@test fine_structure ≈ 7.2973525693e-3

@test classical_radius_factor ≈ 1.4399645478508574e-9
@test r_p ≈ 1.5346982671888946e-18
@test h_bar_planck ≈ 6.582119568038699e-16
@test kg_per_eV ≈ 1.7826619216224324e-36

#set units to MKS
setunits(MKS)

#mass should have unit in kg instead of eV/c^2
@test m_electron ≈ 9.109383701540453e-31
@test m_proton ≈ 1.6726219236839978e-27
@test m_neutron ≈ 1.674927498034172e-27
@test m_muon ≈ 1.8835316270433453e-28
@test m_helion ≈ 5.006412778866699e-27
@test m_deuteron ≈ 3.3435837724217406e-27

#c_light shouldn't change
@test c_light ≈ 2.99792458e8

#e_charge should be in Coulomb
@test e_charge ≈ 1.602176634e-19

#unitless constants shouldn't change


@test kg_per_amu ≈ 1.66053906660e-27
@test eV_per_amu ≈ 9.3149410242e8
@test N_avogadro ≈ 6.02214076e23
@test fine_structure ≈ 7.2973525693e-3

#set units to MKS and mass to MeV, and mass should have unit in MeV/c^2
setunits(MKS, mass_unit="MeV/c^2")

@test m_electron ≈ 0.51099895000e6 / 10^6
@test m_proton ≈ 0.93827208816e9 / 10^6
@test m_neutron ≈ 0.93956542052e9 / 10^6
@test m_muon ≈ 105.6583755e6 / 10^6
@test m_helion ≈ 2.808391607035771e9 / 10^6
@test m_deuteron ≈ 1.87561294257e9 / 10^6

#set the unit of time to hour, the unit of speed should be m/hour
setunits(time_unit="hr")

@test c_light ≈ 2.99792458e8 * 3600

# set units to CGS, the speed should be in cm/s, and mass should be in gram
setunits(CGS)

@test c_light ≈ 2.99792458e10

@test m_electron ≈ 9.109383701540453e-28
@test m_proton ≈ 1.6726219236839978e-24
@test m_neutron ≈ 1.674927498034172e-24
@test m_muon ≈ 1.8835316270433453e-25
@test m_helion ≈ 5.006412778866699e-24
@test m_deuteron ≈ 3.3435837724217406e-24

#if unit has the wrong dimension exist a MethodError will be thrown
@test_throws ErrorException setunits(mass_unit="km")
@test_throws ErrorException setunits(charge_unit="eV/c^2")
@test_throws ErrorException setunits(time_unit="J/km")



#test massof and chargeof

#create some particles

H = Species("H")

#mass should be in amu and charge in elementary charge
setunits(mass_unit=u"eV/c^2")

# @test mass(H) ≈ 1.000 atol = 1e-4
# @test charge(H) ≈ 1
@test mass(H, "kg") ≈ 1.6605e-27 atol = 1e-27

# since Species("H") yields a neutral atom, the following test is wrong
# @test charge(H, "C") ≈ 1.6021e-19 atol = 1e-19

s = Species("electron")
@test mass(s, u"eV/c^2") ≈ 510998.95069 atol = 1
# test default APCdef settings
@APCdef

@testset "testAPCdef" begin
#constants should be of type float in the right unit
@test C_LIGHT ≈ 2.99792458e8
@test H_PLANCK ≈ 4.135667696e-15
@test H_BAR_PLANCK ≈ 6.582119568038699e-16
@test R_E ≈ 2.8179403205e-15
@test R_P ≈ 1.5346982640795807e-12
@test E_CHARGE ≈ 1
@test MU_0_VAC ≈ 1.25663706127e-6
@test EPS_0_VAC ≈ 8.8541878188e-12
@test CLASSICAL_RADIUS_FACTOR ≈ 1.4399645468825422e-15
@test FINE_STRUCTURE ≈ 0.0072973525643
@test N_AVOGADRO ≈ 6.02214076e23

#test massof() and chargeof()
H = Species("H")

@test massof(H) ≈ 9.394018682512946e8
@test chargeof(H) ≈ 0

e = Species("electron")
@test massof(e) ≈ 510998.95069
@test chargeof(e) ≈ -1
end

@testset "testSpecies" begin
#test '+' and '-' signs to charge
@test chargeof(Species("H+")) ≈ 1
@test chargeof(Species("Mg++")) ≈ 2
@test chargeof(Species("Al+3")) ≈ 3
@test chargeof(Species("Cl-")) ≈ -1
@test chargeof(Species("O--")) ≈ -2
@test chargeof(Species("N-3")) ≈ -3
end
Loading