compare_particle_positions.py File Reference

Go to the source code of this file.

Namespaces
namespace	compare_particle_positions
	Compare whether initial states and states after 1 loop through box are identical.

Variables
	compare_particle_positions.parser
	compare_particle_positions.dest
	compare_particle_positions.help
	compare_particle_positions.default
	compare_particle_positions.type
	compare_particle_positions.required
	compare_particle_positions.args = parser.parse_args()
	compare_particle_positions.ndim = args.ndim
	compare_particle_positions.comparison_tolerance = args.tolerance
	compare_particle_positions.boxsize = args.boxsize
	compare_particle_positions.end_time = args.t_end
	compare_particle_positions.pvdfile = args.pvdfile
	compare_particle_positions.reader_0 = ParticleVTUReader(pvdfile=pvdfile, snapshot_time=0.0)
	compare_particle_positions.partData_0 = reader_0.load()
	compare_particle_positions.reader_1 = ParticleVTUReader(pvdfile=pvdfile, snapshot_time=end_time)
	compare_particle_positions.partData_1 = reader_1.load()
	compare_particle_positions.actual_end_time = reader_1.get_snapshot_time()
	compare_particle_positions.coords_0 = partData_0.x[:, :ndim]
	compare_particle_positions.vels_0 = partData_0.v[:, :ndim]
	compare_particle_positions.ids_0 = partData_0.partid
	compare_particle_positions.sort0 = np.argsort(ids_0)
	compare_particle_positions.coords_expect = coords_0 + actual_end_time * vels_0
	compare_particle_positions.mask = coords_expect[:, dim] > boxsize
	compare_particle_positions.coords_1 = partData_1.x[:, :ndim]
	compare_particle_positions.ids_1 = partData_1.partid
	compare_particle_positions.sort1 = np.argsort(ids_1)
int	compare_particle_positions.errors = 0
float	compare_particle_positions.diff = coords_expect[i, n] / coords_1[i, n] - 1.0