def new_instance_of_capreole_and_nbody_code(self):
gridcode=Capreole(number_of_workers=self.number_of_workers)
gridcode.initialize_code()
nbodycode = PhiGRAPE(mode="gpu")
#nbodycode = Octgrav(mode="gpu")
nbodycode.initialize_code()
nbodycode.parameters.epsilon_squared = (self.size / (10000.0 * self.number_of_grid_points)) ** 2
nbodycode.commit_parameters()
result = HydroGridAndNbodyWithAccelerationTransfer(gridcode, nbodycode)
return result
def new_instance_of_capreole_and_nbody_code(self):
gridcode = Capreole(number_of_workers=self.number_of_workers)
gridcode.initialize_code()
nbodycode = PhiGRAPE(mode="gpu")
#nbodycode = Octgrav(mode="gpu")
nbodycode.initialize_code()
nbodycode.parameters.epsilon_squared = (
self.size / (10000.0 * self.number_of_grid_points))**2
nbodycode.commit_parameters()
result = HydroGridAndNbodyWithAccelerationTransfer(gridcode, nbodycode)
return result
def new_instance_of_athena_and_nbody_code(self):
gridcode=Athena(redirection="none", number_of_workers=1)
gridcode.initialize_code()
gridcode.parameters.gamma = self.gamma
gridcode.parameters.courant_number=0.3
gridcode.set_has_external_gravitational_potential(1)
nbodycode = PhiGRAPE(mode="gpu")
#nbodycode = Octgrav(mode="gpu")
nbodycode.initialize_code()
nbodycode.parameters.epsilon_squared = (self.size / (10000.0 * self.number_of_grid_points)) ** 2
nbodycode.commit_parameters()
result = HydroGridAndNbody(gridcode, nbodycode)
return result
def new_instance_of_athena_and_nbody_code(self):
gridcode = Athena(redirection="none", number_of_workers=1)
gridcode.initialize_code()
gridcode.parameters.gamma = self.gamma
gridcode.parameters.courant_number = 0.3
gridcode.set_has_external_gravitational_potential(1)
nbodycode = PhiGRAPE(mode="gpu")
#nbodycode = Octgrav(mode="gpu")
nbodycode.initialize_code()
nbodycode.parameters.epsilon_squared = (
self.size / (10000.0 * self.number_of_grid_points))**2
nbodycode.commit_parameters()
result = HydroGridAndNbody(gridcode, nbodycode)
return result
