Physics, Science, Programming, Deep Dives, Thoughts And Theory.
_pos=np.random.uniform(low=[0]*self.dim,high=self.box_dim,size=(n,self.dim)) #Store positions in temporary variable self._positions=_pos def _generate_initial_velocities(self,n,v_mean,v_std,dist="normal"): #Normal distribution with mean v_mean and std v_std if dist=="normal": self.v_mean=v_mean