def xtest8(self):
instance = PhiGRAPEInterface(**default_test_options)
instance.initialize_code()
instance.set_eps2(0.0**2)
instance.set_eta(0.01,0.02)
instance.new_particle(
[0.01,0.01],
[10.,-10.],
[0.0,0.0],
[0.0,0.0],
[-5.0,5.0],
[0.0,0.0],
[0.0,0.0],
[0.1,0.1])
instance.commit_particles()
#HAS NO RESULT...
result = instance.evolve_model(3.14159)
tnow=instance.get_time()['time']
print "after evolve(pi), tnow = %f" % (tnow)
#self.assertEqual( id1, 1)
"""
instance.evolve(instance.get_time(),1)
id2=instance.find_colliding_secondary(id1)
self.assertEqual( id2, 2)
self.assertAlmostEqual( tnow, 2.,2)
state1 = instance.get_state(id1)
state2 = instance.get_state(id2)
self.assertTrue( abs(state1['x'] - state2['x'])<0.2)
"""
instance.cleanup_code()
instance.stop()
def xtest8(self):
instance = PhiGRAPEInterface(**default_test_options)
instance.initialize_code()
instance.set_eps2(0.0**2)
instance.set_eta(0.01,0.02)
instance.new_particle(
[0.01,0.01],
[10.,-10.],
[0.0,0.0],
[0.0,0.0],
[-5.0,5.0],
[0.0,0.0],
[0.0,0.0],
[0.1,0.1])
instance.commit_particles()
#HAS NO RESULT...
result = instance.evolve_model(3.14159)
tnow=instance.get_time()['time']
print "after evolve(pi), tnow = %f" % (tnow)
#self.assertEqual( id1, 1)
"""
instance.evolve(instance.get_time(),1)
id2=instance.find_colliding_secondary(id1)
self.assertEqual( id2, 2)
self.assertAlmostEqual( tnow, 2.,2)
state1 = instance.get_state(id1)
state2 = instance.get_state(id2)
self.assertTrue( abs(state1['x'] - state2['x'])<0.2)
"""
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 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()
