using Test
using PotentialCalculation.Atoms

@testset "atoms" begin

    a = AtomOnlySymbol("H")

    am = convert(AtomWithMass, a)
    aa = convert(AtomOnlySymbol, "H")

    @test a.id == am.id

    x = AtomWithMass("x", 3.2)

    O = AtomWithMass("O")

    @test am.mass == masses["H"]
    @test O.mass == masses["O"]
end
