コード例 #1
0
ファイル: md_header.py プロジェクト: jensengroup/molstat
def initialize_particles(n_atoms, temperature, rho):

    U = initialize_positions(n_atoms, rho)
    box = initialize_box(n_atoms, rho)
    V = initialize_velocities(U, temperature)
    epot, F, Vir = fortran_lj(U, box)

    return U, V, F, box[0]
コード例 #2
0
def initialize_particles(n_atoms, temperature, rho):

    U = initialize_positions(n_atoms, rho)
    box = initialize_box(n_atoms, rho)
    V = initialize_velocities(U, temperature)
    epot, F, Vir = fortran_lj(U, box)

    return U, V, F, box[0]
コード例 #3
0
ファイル: md_header_3d.py プロジェクト: jensengroup/molstat
def lennard_jones(U, box):
    return fortran_lj(U, numpy.array([box,box,box]))
コード例 #4
0
ファイル: md_header_3d.py プロジェクト: larsbratholm/molstat
def lennard_jones(U, box):
    return fortran_lj(U, numpy.array([box, box, box]))