compare_results.py File Reference

Go to the source code of this file.

Namespaces
namespace	compare_results
	Compare smoothing lengths that Peano yields and what you expect them to be.

Functions
	compare_results.sml_search (i)

Variables
	compare_results.parser
	compare_results.dest
	compare_results.help
	compare_results.required
	compare_results.type
	compare_results.default
	compare_results.choices
	compare_results.args = parser.parse_args()
	compare_results.icfile = args.ic_filename
	compare_results.h_tolerance = args.h_tol
	compare_results.kernel = args.kernel
	compare_results.ndim = args.ndim
	compare_results.eta = args.eta
	compare_results.nneigh = args.nneigh
float	compare_results.comparison_tolerance = 10.0 * h_tolerance
str	compare_results.vtufile = "particles-1.vtu"
	compare_results.hfile = h5py.File(icfile, "r")
	compare_results.gas = hfile["PartType0"]
	compare_results.coords = gas["Coordinates"][:]
	compare_results.ids = gas["ParticleIDs"][:]
	compare_results.reader = ParticleVTUReader(vtufile=vtufile, verbose=False)
	compare_results.partData = reader.load()
	compare_results.sort_ic = np.argsort(ids)
	compare_results.tree = KDTree(coords)
	compare_results.distances
	compare_results.indexes
	compare_results.k
	compare_results.neighbours = coords[indexes]
	compare_results.npart = coords.shape[0]
	compare_results.sml_python = np.zeros((npart))
	compare_results.pool = multiprocessing.Pool()
	compare_results.sml_peano = partData.smoothingLength
	compare_results.partIDs_peano = partData.partid
	compare_results.sort_peano = np.argsort(partIDs_peano)
int	compare_results.errors = 0
tuple	compare_results.errorstr
float	compare_results.diff = 1.0 - sml_peano[i] / sml_python[i]