Skip to content
Snippets Groups Projects
Commit 0fe2621d authored by Cristian Lalescu's avatar Cristian Lalescu
Browse files

Merge branch 'fix/read_stats' into 'master'

fixes read_moments iteration range

See merge request !104
parents 70c89fb4 6d68b177
Branches
Tags v1.0.2
1 merge request!104fixes read_moments iteration range
Pipeline #200549 passed
......@@ -282,7 +282,8 @@ def read_moments(
if type(iteration_range) == type(None):
iteration_range = [self.statistics['iter0'], self.statistics['iter1']]
mm = pp_file[base_group + '/moments/' + quantity][
iteration_range[0]//self.parameters['niter_stat']:iteration_range[1]//self.parameters['niter_stat']]
iteration_range[0]//self.parameters['niter_stat']:
iteration_range[1]//self.parameters['niter_stat'] + 1]
if return_full_history:
return mm
mean_mm = np.zeros(mm.shape[1:], dtype=mm.dtype)
......@@ -316,7 +317,8 @@ def read_histogram(
if type(nbins) == type(None):
nbins = self.parameters['histogram_bins']
hh = pp_file[base_group + '/histograms/' + quantity][
iteration_range[0]//niter_stat:iteration_range[1]//niter_stat]
iteration_range[0]//niter_stat:
iteration_range[1]//niter_stat + 1]
if (hh.shape[-1] == 4):
hh = hh[..., :3]
max_estimate /= 3**0.5
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment