parallel interface working
simulate_atomic working, simulate_density needs help
simulate_atomic working
gpu code draft in