Exemplo n.º 1
0
def test_nonrelat():
    a = 2.7e6
    rmin = 1e-7
    rmax = 50
    N = 1000
    Z = 92
    relat = 0
    R = mesh_exp(rmin, rmax, a, N)
    Rp = mesh_exp_deriv(rmin, rmax, a, N)
    V_tot = -Z / R
    energies = []
    for n in range(8):
        for l in range(n):
            E, P, Q = solve_radial_eigenproblem(0.0, n, l, -1, 1e-11, 100, V_tot, R, Rp, Z, relat, False, -10000, 0)
            energies.append(E)
def test_nonrelat():
    a = 2.7e6
    rmin = 1e-7
    rmax = 50
    N = 1000
    Z = 92
    relat = 0
    R = mesh_exp(rmin, rmax, a, N)
    Rp = mesh_exp_deriv(rmin, rmax, a, N)
    V_tot = -Z / R
    energies = []
    for n in range(8):
        for l in range(n):
            E, P, Q = solve_radial_eigenproblem(0.0, n, l, -1, 1e-11, 100,
                    V_tot, R, Rp, Z, relat, False, -10000, 0)
            energies.append(E)
Exemplo n.º 3
0
def test_relat():
    a = 6.2e7
    rmin = 1e-8
    rmax = 50
    N = 1000
    Z = 92
    relat = 0
    c = 137.0359895
    R = mesh_exp(rmin, rmax, a, N)
    Rp = mesh_exp_deriv(rmin, rmax, a, N)
    V_tot = -Z / R
    energies = []
    for n in range(8):
        for l in range(n):
            relats = [2]
            if l > 0:
                relats.append(3)
            for relat in relats:
                E, y = solve_radial_eigenproblem(c, n, l, relat, 1e-11, 100,
                        V_tot, R, Rp, Z, relat, False, -10000, 0)
                energies.append(E)
def test_relat():
    a = 6.2e7
    rmin = 1e-8
    rmax = 50
    N = 1000
    Z = 92
    relat = 0
    c = 137.0359895
    R = mesh_exp(rmin, rmax, a, N)
    Rp = mesh_exp_deriv(rmin, rmax, a, N)
    V_tot = -Z / R
    energies = []
    for n in range(8):
        for l in range(n):
            relats = [2]
            if l > 0:
                relats.append(3)
            for relat in relats:
                E, P, Q = solve_radial_eigenproblem(c, n, l, relat, 1e-11, 100,
                        V_tot, R, Rp, Z, relat, False, -10000, 0)
                energies.append(E)