In [1]:
import numpy as np
import os
import qp
In [2]:
QP_DIR = os.path.abspath(os.path.dirname(qp.__file__))
data_file = os.path.join(QP_DIR, 'data', 'test.hdf5')
In [3]:
ens = qp.read(data_file)
In [4]:
itr = qp.iterator(data_file, 10)
In [5]:
ens.npdf
Out[5]:
100
In [6]:
test_vals = np.linspace(0., 1., 11)
In [7]:
for start, end, ens_i in itr:
    print(start, end, np.max(ens[start:end].pdf(test_vals) - ens_i.pdf(test_vals)))
0 10 1.1102230246251565e-16
10 20 0.0
20 30 0.0
30 40 1.1102230246251565e-16
40 50 0.0
50 60 0.0
60 70 1.1102230246251565e-16
70 80 0.0
80 90 1.1102230246251565e-16
90 100 1.1102230246251565e-16
In [ ]: