Physics, Science, Programming, Deep Dives, Thoughts And Theory.
molecules at each timestep f self.x_dynamics=np.zeros(((self.n_molecules,self.dim,self.n_iters))) self.v_dynamics=np.zeros((self.n_molecules,self.n_iters)) #In each iteration for i in range(self.n_iters): #Save positions and velocities to the defined tensors self.x_dynamics[:,:,i]=self.positions