Peano
Loading...
Searching...
No Matches
compare_results Namespace Reference

Compare smoothing lengths that Peano yields and what you expect them to be. More...

Functions

 sml_search (i)
 

Variables

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

Detailed Description

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

Function Documentation

◆ sml_search()

compare_results.sml_search ( i)

Definition at line 170 of file compare_results.py.

Variable Documentation

◆ args

compare_results.args = parser.parse_args()

Definition at line 86 of file compare_results.py.

◆ choices

compare_results.choices

Definition at line 83 of file compare_results.py.

◆ comparison_tolerance

float compare_results.comparison_tolerance = 10.0 * h_tolerance

Definition at line 114 of file compare_results.py.

◆ coords

compare_results.coords = gas["Coordinates"][:]

Definition at line 136 of file compare_results.py.

◆ default

compare_results.default

Definition at line 49 of file compare_results.py.

◆ dest

compare_results.dest

Definition at line 40 of file compare_results.py.

◆ diff

float compare_results.diff = 1.0 - sml_peano[i] / sml_python[i]

Definition at line 216 of file compare_results.py.

◆ distances

compare_results.distances

Definition at line 163 of file compare_results.py.

◆ errors

int compare_results.errors = 0

Definition at line 206 of file compare_results.py.

◆ errorstr

tuple compare_results.errorstr
Initial value:
1= (
2 "Found difference: "
3 + "particle ID = {0:4d} "
4 + "h_peano = {1:.5e} "
5 + "h_python = {2:.5e} "
6 + "diff = {3:.5e}"
7)

Definition at line 207 of file compare_results.py.

◆ eta

compare_results.eta = args.eta

Definition at line 94 of file compare_results.py.

◆ gas

compare_results.gas = hfile["PartType0"]

Definition at line 135 of file compare_results.py.

◆ h_tolerance

compare_results.h_tolerance = args.h_tol

Definition at line 88 of file compare_results.py.

◆ help

compare_results.help

Definition at line 41 of file compare_results.py.

◆ hfile

compare_results.hfile = h5py.File(icfile, "r")

Definition at line 134 of file compare_results.py.

◆ icfile

compare_results.icfile = args.ic_filename

Definition at line 87 of file compare_results.py.

◆ ids

compare_results.ids = gas["ParticleIDs"][:]

Definition at line 138 of file compare_results.py.

◆ indexes

compare_results.indexes

Definition at line 163 of file compare_results.py.

◆ k

compare_results.k

Definition at line 163 of file compare_results.py.

◆ kernel

str compare_results.kernel = args.kernel

Definition at line 89 of file compare_results.py.

◆ ndim

compare_results.ndim = args.ndim

Definition at line 92 of file compare_results.py.

◆ neighbours

compare_results.neighbours = coords[indexes]

Definition at line 164 of file compare_results.py.

◆ nneigh

compare_results.nneigh = args.nneigh

Definition at line 95 of file compare_results.py.

◆ npart

compare_results.npart = coords.shape[0]

Definition at line 166 of file compare_results.py.

◆ parser

compare_results.parser
Initial value:
1= argparse.ArgumentParser(
2 description=
3)

Definition at line 18 of file compare_results.py.

◆ partData

compare_results.partData = reader.load()

Definition at line 146 of file compare_results.py.

◆ partIDs_peano

compare_results.partIDs_peano = partData.partid

Definition at line 193 of file compare_results.py.

◆ pool

compare_results.pool = multiprocessing.Pool()

Definition at line 183 of file compare_results.py.

◆ reader

compare_results.reader = ParticleVTUReader(vtufile=vtufile, verbose=False)

Definition at line 145 of file compare_results.py.

◆ required

compare_results.required

Definition at line 42 of file compare_results.py.

◆ sml_peano

compare_results.sml_peano = partData.smoothingLength

Definition at line 192 of file compare_results.py.

◆ sml_python

compare_results.sml_python = np.zeros((npart))

Definition at line 167 of file compare_results.py.

◆ sort_ic

compare_results.sort_ic = np.argsort(ids)

Definition at line 154 of file compare_results.py.

◆ sort_peano

compare_results.sort_peano = np.argsort(partIDs_peano)

Definition at line 196 of file compare_results.py.

◆ tree

compare_results.tree = KDTree(coords)

Definition at line 162 of file compare_results.py.

◆ type

compare_results.type

Definition at line 48 of file compare_results.py.

◆ vtufile

str compare_results.vtufile = "particles-1.vtu"

Definition at line 132 of file compare_results.py.