test_delaunay.py
No OneTemporary
Actions

Subscribers

None

File Metadata

	#!/usr/bin/env python

	from pNbody import *
	from pNbody import ic
	from pNbody import libutil

	from numpy import *


	import sys
	import time

	import copy

	from PyGear import gadget


	nb = Nbody("snap.dat",ftype="gadget")


	###################################
	# init PyGear

	gadget.InitMPI() # init MPI
	gadget.InitDefaultParameters() # init default parameters
	gadget.Info()

	params = {}
	params['ErrTolTheta'] = 0.7
	params['DesNumNgb'] = 15
	params['MaxNumNgbDeviation'] = 3
	params['UnitLength_in_cm'] = 3.085e+21
	params['UnitMass_in_g'] = 4.435693e+44
	params['UnitVelocity_in_cm_per_s'] = 97824708.2699
	params['BoxSize'] = 1.0
	params['PartAllocFactor'] = 10

	gadget.SetParameters(params)


	###################################
	# load particles

	gadget.LoadParticles(array(nb.npart),nb.pos,nb.vel,nb.mass,nb.num,nb.tpe)

	###################################
	# compute Delaunay

	t1 = time.time()
	gadget.ConstructDelaunay()
	t2 = time.time()

	print "Delaunay","in",(t2-t1),"s"