def test7(self): instance = PhiGRAPEInterface(**default_test_options)#(debugger="xterm") instance.initialize_code() instance.set_eps2(0.0**2) instance.set_eta(0.01,0.02) instance.commit_parameters() instance.new_particle( [1.0,1.0,1.0], [1.0,0.0,-1.0], [0.0,0.0,0.0], [0.0,0.0,0.0], [0.0,1.0,0.0], [0.0,0.0,0.0], [0.0,0.0,0.0] ) instance.recommit_particles() Ep=instance.get_potential_energy()['potential_energy'] Ek=instance.get_kinetic_energy()['kinetic_energy'] self.assertEqual( Ek, 0.5) self.assertEqual( Ep, -2.5) instance.delete_particle(2) instance.recommit_particles() n=instance.get_number_of_particles()['number_of_particles'] Ep=instance.get_potential_energy()['potential_energy'] Ek=instance.get_kinetic_energy()['kinetic_energy'] self.assertEqual( n, 2) self.assertEqual( Ek, 0.) self.assertEqual( Ep, -0.5) instance.cleanup_code() instance.stop()
def test7(self): instance = PhiGRAPEInterface(**default_test_options)#(debugger="xterm") instance.initialize_code() instance.set_eps2(0.0**2) instance.set_eta(0.01,0.02) instance.commit_parameters() instance.new_particle( [1.0,1.0,1.0], [1.0,0.0,-1.0], [0.0,0.0,0.0], [0.0,0.0,0.0], [0.0,1.0,0.0], [0.0,0.0,0.0], [0.0,0.0,0.0] ) instance.recommit_particles() Ep=instance.get_potential_energy()['potential_energy'] Ek=instance.get_kinetic_energy()['kinetic_energy'] self.assertEqual( Ek, 0.5) self.assertEqual( Ep, -2.5) instance.delete_particle(2) instance.recommit_particles() n=instance.get_number_of_particles()['number_of_particles'] Ep=instance.get_potential_energy()['potential_energy'] Ek=instance.get_kinetic_energy()['kinetic_energy'] self.assertEqual( n, 2) self.assertEqual( Ek, 0.) self.assertEqual( Ep, -0.5) instance.cleanup_code() instance.stop()
def test9(self): instance = PhiGRAPEInterface(**default_test_options) instance.initialize_code() instance.set_eps2(0) instance.commit_parameters() id1,errorcode = instance.new_particle(mass = 10.0, radius = 1.0, x = 0.0, y = 0.0, z = 0.0, vx = 0.0, vy = 0.0, vz = 0.0) id2,errorcode = instance.new_particle(mass = 1.0, radius = 1.0, x = 2.0, y = 0.0, z = 0.0, vx = 10.0, vy = 0.0, vz = 0.0) instance.commit_particles() potential, errorcode = instance.get_potential(id1) self.assertEquals(errorcode, 0) self.assertAlmostRelativeEquals(potential, -1.0 / numpy.sqrt(2.0**2), 8) total_potential, errorcode = instance.get_potential_energy() potentials, errorcode = instance.get_potential([id1, id2]) instance.cleanup_code() instance.stop() self.assertAlmostRelativeEquals(total_potential, numpy.sum(potentials * [10.0, 1.0]) / 2.0)
def test9(self): instance = PhiGRAPEInterface(**default_test_options) instance.initialize_code() instance.set_eps2(0) instance.commit_parameters() id1,errorcode = instance.new_particle(mass = 10.0, radius = 1.0, x = 0.0, y = 0.0, z = 0.0, vx = 0.0, vy = 0.0, vz = 0.0) id2,errorcode = instance.new_particle(mass = 1.0, radius = 1.0, x = 2.0, y = 0.0, z = 0.0, vx = 10.0, vy = 0.0, vz = 0.0) instance.commit_particles() potential, errorcode = instance.get_potential(id1) self.assertEquals(errorcode, 0) self.assertAlmostRelativeEquals(potential, -1.0 / numpy.sqrt(2.0**2), 8) total_potential, errorcode = instance.get_potential_energy() potentials, errorcode = instance.get_potential([id1, id2]) instance.cleanup_code() instance.stop() self.assertAlmostRelativeEquals(total_potential, numpy.sum(potentials * [10.0, 1.0]) / 2.0)