First import all the modules such as healpy and astropy needed for analyzing the structure
In [1]:
import healpix_util as hu
import astropy as ap
import numpy as np
from astropy.io import fits
from astropy.table import Table
import astropy.io.ascii as ascii
from astropy.io import fits
from astropy.constants import c
import matplotlib.pyplot as plt
import math as m
from math import pi
#from scipy.constants import c
import scipy.special as sp
from astroML.decorators import pickle_results
from scipy import integrate
import warnings
from sklearn.neighbors import BallTree
import pickle
import multiprocessing as mp
import time
from aptestmetricdt import *
from aptestmetricdz import *
from scipy.spatial import distance as d
from apcat import *
from progressbar import *
from tqdm import *
from functools import partial
import pymangle
from apdz import *
from apdt import *
from scipy.optimize import curve_fit
#from astroML.datasets import fetch_sdss_specgals
#from astroML.correlation import bootstrap_two_point_angular
%matplotlib inline
Read the data file (taken from http://cosmo.nyu.edu/~eak306/SDSS-LRG.html ) converted to ascii with comoving distance etc. in V01 reading from pkl files for faster read
In [16]:
# Getting back the objects:
with open('../output/randDR12Qbin1.pkl') as f: # Python 3: open(..., 'rb')
dat = pickle.load(f)
dat
Out[16]:
array([[ 1.360358, 4.269086, 0.514813],
[ 1.134544, 0.060924, 0.615164],
[ 0.835325, 5.941764, 0.184429],
...,
[ 0.786509, 2.980918, 0.480278],
[ 0.835873, 0.522808, -0.123566],
[ 1.197253, 3.449246, 0.619426]])
In [17]:
rr2d=np.zeros((20,20))
In [18]:
len(dat)
Out[18]:
48309
In [19]:
rng = np.array([[0, 0.02], [0, 0.02]])
In [20]:
rr2d
Out[20]:
array([[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.]])
In [21]:
%%time
dist0=d.cdist([dat[0],],dat,APdz)[0]
dist1=d.cdist([dat[0],],dat,APzdth)[0]
np.histogram2d(dist0, dist1, bins=20, range=rng)[0]
CPU times: user 459 ms, sys: 17.9 ms, total: 477 ms
Wall time: 465 ms
In [22]:
rr2d
Out[22]:
array([[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.]])
In [23]:
rr2d=np.zeros((20,20))
In [24]:
%%time
while(len(dat))>0:
dist0=d.cdist([dat[0],],dat,APdz)[0]
dist1=d.cdist([dat[0],],dat,APzdth)[0]
rr2d+=np.histogram2d(dist0, dist1, bins=20, range=rng)[0]
dat=np.delete(dat,0,axis=0)
if len(dat)%1000==0:
print len(dat)/1000
print rr2d
/Users/rohin/anaconda/lib/python2.7/site-packages/numpy/lib/function_base.py:804: RuntimeWarning: invalid value encountered in greater_equal
not_smaller_than_edge = (sample[:, i] >= edges[i][-1])
48
47
46
45
44
43
42
41
40
39
38
37
36
35
34
33
32
31
30
29
28
27
26
25
24
23
22
21
20
19
18
17
16
15
14
13
12
11
10
9
8
7
6
5
4
3
2
1
0
[[ 4.59620000e+04 1.10000000e+01 2.40000000e+01 4.30000000e+01
5.50000000e+01 6.60000000e+01 6.90000000e+01 7.70000000e+01
9.70000000e+01 1.04000000e+02 1.24000000e+02 1.45000000e+02
1.37000000e+02 1.46000000e+02 1.59000000e+02 1.80000000e+02
1.90000000e+02 2.11000000e+02 2.29000000e+02 2.09000000e+02]
[ 8.00000000e+00 1.40000000e+01 2.00000000e+01 3.60000000e+01
5.40000000e+01 5.20000000e+01 7.50000000e+01 9.30000000e+01
1.01000000e+02 1.14000000e+02 1.32000000e+02 1.28000000e+02
1.19000000e+02 1.60000000e+02 1.71000000e+02 1.78000000e+02
1.99000000e+02 1.76000000e+02 2.13000000e+02 2.23000000e+02]
[ 7.00000000e+00 1.90000000e+01 3.10000000e+01 4.20000000e+01
4.60000000e+01 7.20000000e+01 7.00000000e+01 8.40000000e+01
9.00000000e+01 1.26000000e+02 1.31000000e+02 1.28000000e+02
1.31000000e+02 1.70000000e+02 1.60000000e+02 1.67000000e+02
1.86000000e+02 2.09000000e+02 2.18000000e+02 2.44000000e+02]
[ 8.00000000e+00 2.40000000e+01 2.90000000e+01 2.80000000e+01
5.20000000e+01 5.20000000e+01 7.00000000e+01 8.60000000e+01
7.80000000e+01 1.17000000e+02 9.40000000e+01 1.39000000e+02
1.58000000e+02 1.47000000e+02 1.77000000e+02 1.67000000e+02
2.05000000e+02 1.93000000e+02 2.08000000e+02 2.32000000e+02]
[ 6.00000000e+00 1.90000000e+01 2.70000000e+01 4.50000000e+01
5.50000000e+01 5.10000000e+01 7.10000000e+01 9.60000000e+01
1.01000000e+02 1.26000000e+02 1.05000000e+02 1.32000000e+02
1.35000000e+02 1.43000000e+02 1.64000000e+02 1.81000000e+02
1.71000000e+02 2.18000000e+02 2.03000000e+02 2.22000000e+02]
[ 4.00000000e+00 2.40000000e+01 2.10000000e+01 3.90000000e+01
5.40000000e+01 5.90000000e+01 8.90000000e+01 8.20000000e+01
7.20000000e+01 9.20000000e+01 9.90000000e+01 1.39000000e+02
1.50000000e+02 1.61000000e+02 1.41000000e+02 1.77000000e+02
2.16000000e+02 1.95000000e+02 2.12000000e+02 2.02000000e+02]
[ 5.00000000e+00 1.20000000e+01 4.40000000e+01 5.30000000e+01
4.80000000e+01 5.40000000e+01 8.10000000e+01 8.40000000e+01
9.80000000e+01 8.10000000e+01 1.46000000e+02 1.31000000e+02
1.36000000e+02 1.48000000e+02 1.64000000e+02 1.53000000e+02
2.06000000e+02 2.26000000e+02 2.17000000e+02 2.06000000e+02]
[ 3.00000000e+00 1.40000000e+01 2.20000000e+01 5.20000000e+01
4.30000000e+01 5.30000000e+01 7.10000000e+01 9.00000000e+01
8.90000000e+01 9.90000000e+01 1.20000000e+02 1.18000000e+02
1.29000000e+02 1.50000000e+02 1.59000000e+02 1.56000000e+02
2.22000000e+02 2.04000000e+02 2.13000000e+02 2.43000000e+02]
[ 6.00000000e+00 1.70000000e+01 3.00000000e+01 3.10000000e+01
4.90000000e+01 6.60000000e+01 8.50000000e+01 9.10000000e+01
1.12000000e+02 1.11000000e+02 1.37000000e+02 1.32000000e+02
1.50000000e+02 1.57000000e+02 1.58000000e+02 1.76000000e+02
1.72000000e+02 1.90000000e+02 2.07000000e+02 2.18000000e+02]
[ 3.00000000e+00 1.30000000e+01 3.20000000e+01 4.50000000e+01
5.20000000e+01 5.20000000e+01 6.30000000e+01 7.60000000e+01
1.00000000e+02 9.80000000e+01 1.38000000e+02 1.24000000e+02
1.21000000e+02 1.45000000e+02 1.52000000e+02 1.73000000e+02
1.86000000e+02 1.74000000e+02 1.78000000e+02 1.93000000e+02]
[ 6.00000000e+00 2.50000000e+01 2.70000000e+01 4.00000000e+01
5.20000000e+01 7.40000000e+01 5.50000000e+01 9.30000000e+01
9.70000000e+01 9.50000000e+01 1.21000000e+02 1.06000000e+02
1.38000000e+02 1.60000000e+02 1.65000000e+02 1.77000000e+02
2.04000000e+02 1.93000000e+02 1.88000000e+02 1.97000000e+02]
[ 7.00000000e+00 2.30000000e+01 3.50000000e+01 4.40000000e+01
5.10000000e+01 4.90000000e+01 5.40000000e+01 8.60000000e+01
1.02000000e+02 1.06000000e+02 1.03000000e+02 1.31000000e+02
1.20000000e+02 1.50000000e+02 1.65000000e+02 1.73000000e+02
1.87000000e+02 2.11000000e+02 1.94000000e+02 2.15000000e+02]
[ 6.00000000e+00 1.10000000e+01 2.50000000e+01 5.60000000e+01
7.40000000e+01 6.10000000e+01 9.00000000e+01 8.30000000e+01
8.30000000e+01 1.13000000e+02 1.07000000e+02 1.30000000e+02
1.29000000e+02 1.64000000e+02 1.36000000e+02 1.65000000e+02
1.71000000e+02 1.88000000e+02 2.04000000e+02 2.14000000e+02]
[ 4.00000000e+00 1.30000000e+01 2.70000000e+01 3.70000000e+01
5.00000000e+01 5.30000000e+01 7.80000000e+01 8.80000000e+01
8.20000000e+01 1.02000000e+02 1.02000000e+02 1.21000000e+02
1.44000000e+02 1.43000000e+02 1.49000000e+02 1.90000000e+02
1.75000000e+02 1.98000000e+02 2.17000000e+02 2.37000000e+02]
[ 9.00000000e+00 1.00000000e+01 3.40000000e+01 3.40000000e+01
4.00000000e+01 7.00000000e+01 6.00000000e+01 8.90000000e+01
8.00000000e+01 1.11000000e+02 1.25000000e+02 1.39000000e+02
1.50000000e+02 1.66000000e+02 1.45000000e+02 1.67000000e+02
1.76000000e+02 2.01000000e+02 2.14000000e+02 2.27000000e+02]
[ 5.00000000e+00 1.80000000e+01 1.80000000e+01 4.90000000e+01
5.80000000e+01 5.40000000e+01 7.20000000e+01 7.50000000e+01
9.50000000e+01 1.14000000e+02 1.23000000e+02 1.19000000e+02
1.18000000e+02 1.42000000e+02 1.51000000e+02 1.74000000e+02
1.62000000e+02 1.96000000e+02 2.00000000e+02 2.34000000e+02]
[ 4.00000000e+00 2.00000000e+01 2.90000000e+01 3.90000000e+01
5.10000000e+01 5.80000000e+01 8.70000000e+01 8.80000000e+01
9.10000000e+01 8.50000000e+01 1.06000000e+02 1.22000000e+02
1.37000000e+02 1.37000000e+02 1.53000000e+02 1.77000000e+02
1.60000000e+02 2.05000000e+02 1.99000000e+02 1.90000000e+02]
[ 3.00000000e+00 2.30000000e+01 3.40000000e+01 3.80000000e+01
5.40000000e+01 5.20000000e+01 5.90000000e+01 9.30000000e+01
7.40000000e+01 9.80000000e+01 9.50000000e+01 1.19000000e+02
1.36000000e+02 1.48000000e+02 1.69000000e+02 1.64000000e+02
1.84000000e+02 1.76000000e+02 2.14000000e+02 2.00000000e+02]
[ 2.00000000e+00 1.50000000e+01 3.00000000e+01 3.60000000e+01
5.20000000e+01 6.30000000e+01 6.80000000e+01 9.60000000e+01
1.05000000e+02 9.30000000e+01 1.12000000e+02 1.40000000e+02
1.37000000e+02 1.34000000e+02 1.49000000e+02 1.47000000e+02
1.75000000e+02 1.88000000e+02 1.99000000e+02 2.02000000e+02]
[ 5.00000000e+00 1.60000000e+01 2.50000000e+01 3.40000000e+01
4.80000000e+01 7.00000000e+01 7.40000000e+01 8.20000000e+01
7.90000000e+01 8.90000000e+01 1.15000000e+02 1.21000000e+02
1.17000000e+02 1.51000000e+02 1.74000000e+02 1.52000000e+02
1.76000000e+02 1.81000000e+02 1.98000000e+02 2.01000000e+02]]
CPU times: user 3h 50min 9s, sys: 2min 23s, total: 3h 52min 32s
Wall time: 3h 52min 40s
In [25]:
rr2d
Out[25]:
array([[ 4.59620000e+04, 1.10000000e+01, 2.40000000e+01,
4.30000000e+01, 5.50000000e+01, 6.60000000e+01,
6.90000000e+01, 7.70000000e+01, 9.70000000e+01,
1.04000000e+02, 1.24000000e+02, 1.45000000e+02,
1.37000000e+02, 1.46000000e+02, 1.59000000e+02,
1.80000000e+02, 1.90000000e+02, 2.11000000e+02,
2.29000000e+02, 2.09000000e+02],
[ 8.00000000e+00, 1.40000000e+01, 2.00000000e+01,
3.60000000e+01, 5.40000000e+01, 5.20000000e+01,
7.50000000e+01, 9.30000000e+01, 1.01000000e+02,
1.14000000e+02, 1.32000000e+02, 1.28000000e+02,
1.19000000e+02, 1.60000000e+02, 1.71000000e+02,
1.78000000e+02, 1.99000000e+02, 1.76000000e+02,
2.13000000e+02, 2.23000000e+02],
[ 7.00000000e+00, 1.90000000e+01, 3.10000000e+01,
4.20000000e+01, 4.60000000e+01, 7.20000000e+01,
7.00000000e+01, 8.40000000e+01, 9.00000000e+01,
1.26000000e+02, 1.31000000e+02, 1.28000000e+02,
1.31000000e+02, 1.70000000e+02, 1.60000000e+02,
1.67000000e+02, 1.86000000e+02, 2.09000000e+02,
2.18000000e+02, 2.44000000e+02],
[ 8.00000000e+00, 2.40000000e+01, 2.90000000e+01,
2.80000000e+01, 5.20000000e+01, 5.20000000e+01,
7.00000000e+01, 8.60000000e+01, 7.80000000e+01,
1.17000000e+02, 9.40000000e+01, 1.39000000e+02,
1.58000000e+02, 1.47000000e+02, 1.77000000e+02,
1.67000000e+02, 2.05000000e+02, 1.93000000e+02,
2.08000000e+02, 2.32000000e+02],
[ 6.00000000e+00, 1.90000000e+01, 2.70000000e+01,
4.50000000e+01, 5.50000000e+01, 5.10000000e+01,
7.10000000e+01, 9.60000000e+01, 1.01000000e+02,
1.26000000e+02, 1.05000000e+02, 1.32000000e+02,
1.35000000e+02, 1.43000000e+02, 1.64000000e+02,
1.81000000e+02, 1.71000000e+02, 2.18000000e+02,
2.03000000e+02, 2.22000000e+02],
[ 4.00000000e+00, 2.40000000e+01, 2.10000000e+01,
3.90000000e+01, 5.40000000e+01, 5.90000000e+01,
8.90000000e+01, 8.20000000e+01, 7.20000000e+01,
9.20000000e+01, 9.90000000e+01, 1.39000000e+02,
1.50000000e+02, 1.61000000e+02, 1.41000000e+02,
1.77000000e+02, 2.16000000e+02, 1.95000000e+02,
2.12000000e+02, 2.02000000e+02],
[ 5.00000000e+00, 1.20000000e+01, 4.40000000e+01,
5.30000000e+01, 4.80000000e+01, 5.40000000e+01,
8.10000000e+01, 8.40000000e+01, 9.80000000e+01,
8.10000000e+01, 1.46000000e+02, 1.31000000e+02,
1.36000000e+02, 1.48000000e+02, 1.64000000e+02,
1.53000000e+02, 2.06000000e+02, 2.26000000e+02,
2.17000000e+02, 2.06000000e+02],
[ 3.00000000e+00, 1.40000000e+01, 2.20000000e+01,
5.20000000e+01, 4.30000000e+01, 5.30000000e+01,
7.10000000e+01, 9.00000000e+01, 8.90000000e+01,
9.90000000e+01, 1.20000000e+02, 1.18000000e+02,
1.29000000e+02, 1.50000000e+02, 1.59000000e+02,
1.56000000e+02, 2.22000000e+02, 2.04000000e+02,
2.13000000e+02, 2.43000000e+02],
[ 6.00000000e+00, 1.70000000e+01, 3.00000000e+01,
3.10000000e+01, 4.90000000e+01, 6.60000000e+01,
8.50000000e+01, 9.10000000e+01, 1.12000000e+02,
1.11000000e+02, 1.37000000e+02, 1.32000000e+02,
1.50000000e+02, 1.57000000e+02, 1.58000000e+02,
1.76000000e+02, 1.72000000e+02, 1.90000000e+02,
2.07000000e+02, 2.18000000e+02],
[ 3.00000000e+00, 1.30000000e+01, 3.20000000e+01,
4.50000000e+01, 5.20000000e+01, 5.20000000e+01,
6.30000000e+01, 7.60000000e+01, 1.00000000e+02,
9.80000000e+01, 1.38000000e+02, 1.24000000e+02,
1.21000000e+02, 1.45000000e+02, 1.52000000e+02,
1.73000000e+02, 1.86000000e+02, 1.74000000e+02,
1.78000000e+02, 1.93000000e+02],
[ 6.00000000e+00, 2.50000000e+01, 2.70000000e+01,
4.00000000e+01, 5.20000000e+01, 7.40000000e+01,
5.50000000e+01, 9.30000000e+01, 9.70000000e+01,
9.50000000e+01, 1.21000000e+02, 1.06000000e+02,
1.38000000e+02, 1.60000000e+02, 1.65000000e+02,
1.77000000e+02, 2.04000000e+02, 1.93000000e+02,
1.88000000e+02, 1.97000000e+02],
[ 7.00000000e+00, 2.30000000e+01, 3.50000000e+01,
4.40000000e+01, 5.10000000e+01, 4.90000000e+01,
5.40000000e+01, 8.60000000e+01, 1.02000000e+02,
1.06000000e+02, 1.03000000e+02, 1.31000000e+02,
1.20000000e+02, 1.50000000e+02, 1.65000000e+02,
1.73000000e+02, 1.87000000e+02, 2.11000000e+02,
1.94000000e+02, 2.15000000e+02],
[ 6.00000000e+00, 1.10000000e+01, 2.50000000e+01,
5.60000000e+01, 7.40000000e+01, 6.10000000e+01,
9.00000000e+01, 8.30000000e+01, 8.30000000e+01,
1.13000000e+02, 1.07000000e+02, 1.30000000e+02,
1.29000000e+02, 1.64000000e+02, 1.36000000e+02,
1.65000000e+02, 1.71000000e+02, 1.88000000e+02,
2.04000000e+02, 2.14000000e+02],
[ 4.00000000e+00, 1.30000000e+01, 2.70000000e+01,
3.70000000e+01, 5.00000000e+01, 5.30000000e+01,
7.80000000e+01, 8.80000000e+01, 8.20000000e+01,
1.02000000e+02, 1.02000000e+02, 1.21000000e+02,
1.44000000e+02, 1.43000000e+02, 1.49000000e+02,
1.90000000e+02, 1.75000000e+02, 1.98000000e+02,
2.17000000e+02, 2.37000000e+02],
[ 9.00000000e+00, 1.00000000e+01, 3.40000000e+01,
3.40000000e+01, 4.00000000e+01, 7.00000000e+01,
6.00000000e+01, 8.90000000e+01, 8.00000000e+01,
1.11000000e+02, 1.25000000e+02, 1.39000000e+02,
1.50000000e+02, 1.66000000e+02, 1.45000000e+02,
1.67000000e+02, 1.76000000e+02, 2.01000000e+02,
2.14000000e+02, 2.27000000e+02],
[ 5.00000000e+00, 1.80000000e+01, 1.80000000e+01,
4.90000000e+01, 5.80000000e+01, 5.40000000e+01,
7.20000000e+01, 7.50000000e+01, 9.50000000e+01,
1.14000000e+02, 1.23000000e+02, 1.19000000e+02,
1.18000000e+02, 1.42000000e+02, 1.51000000e+02,
1.74000000e+02, 1.62000000e+02, 1.96000000e+02,
2.00000000e+02, 2.34000000e+02],
[ 4.00000000e+00, 2.00000000e+01, 2.90000000e+01,
3.90000000e+01, 5.10000000e+01, 5.80000000e+01,
8.70000000e+01, 8.80000000e+01, 9.10000000e+01,
8.50000000e+01, 1.06000000e+02, 1.22000000e+02,
1.37000000e+02, 1.37000000e+02, 1.53000000e+02,
1.77000000e+02, 1.60000000e+02, 2.05000000e+02,
1.99000000e+02, 1.90000000e+02],
[ 3.00000000e+00, 2.30000000e+01, 3.40000000e+01,
3.80000000e+01, 5.40000000e+01, 5.20000000e+01,
5.90000000e+01, 9.30000000e+01, 7.40000000e+01,
9.80000000e+01, 9.50000000e+01, 1.19000000e+02,
1.36000000e+02, 1.48000000e+02, 1.69000000e+02,
1.64000000e+02, 1.84000000e+02, 1.76000000e+02,
2.14000000e+02, 2.00000000e+02],
[ 2.00000000e+00, 1.50000000e+01, 3.00000000e+01,
3.60000000e+01, 5.20000000e+01, 6.30000000e+01,
6.80000000e+01, 9.60000000e+01, 1.05000000e+02,
9.30000000e+01, 1.12000000e+02, 1.40000000e+02,
1.37000000e+02, 1.34000000e+02, 1.49000000e+02,
1.47000000e+02, 1.75000000e+02, 1.88000000e+02,
1.99000000e+02, 2.02000000e+02],
[ 5.00000000e+00, 1.60000000e+01, 2.50000000e+01,
3.40000000e+01, 4.80000000e+01, 7.00000000e+01,
7.40000000e+01, 8.20000000e+01, 7.90000000e+01,
8.90000000e+01, 1.15000000e+02, 1.21000000e+02,
1.17000000e+02, 1.51000000e+02, 1.74000000e+02,
1.52000000e+02, 1.76000000e+02, 1.81000000e+02,
1.98000000e+02, 2.01000000e+02]])
In [27]:
with open('randDR12QRRbin1.pkl','w') as f:
pickle.dump(rr2d,f)
rr2d
Out[27]:
array([[ 4.59620000e+04, 1.10000000e+01, 2.40000000e+01,
4.30000000e+01, 5.50000000e+01, 6.60000000e+01,
6.90000000e+01, 7.70000000e+01, 9.70000000e+01,
1.04000000e+02, 1.24000000e+02, 1.45000000e+02,
1.37000000e+02, 1.46000000e+02, 1.59000000e+02,
1.80000000e+02, 1.90000000e+02, 2.11000000e+02,
2.29000000e+02, 2.09000000e+02],
[ 8.00000000e+00, 1.40000000e+01, 2.00000000e+01,
3.60000000e+01, 5.40000000e+01, 5.20000000e+01,
7.50000000e+01, 9.30000000e+01, 1.01000000e+02,
1.14000000e+02, 1.32000000e+02, 1.28000000e+02,
1.19000000e+02, 1.60000000e+02, 1.71000000e+02,
1.78000000e+02, 1.99000000e+02, 1.76000000e+02,
2.13000000e+02, 2.23000000e+02],
[ 7.00000000e+00, 1.90000000e+01, 3.10000000e+01,
4.20000000e+01, 4.60000000e+01, 7.20000000e+01,
7.00000000e+01, 8.40000000e+01, 9.00000000e+01,
1.26000000e+02, 1.31000000e+02, 1.28000000e+02,
1.31000000e+02, 1.70000000e+02, 1.60000000e+02,
1.67000000e+02, 1.86000000e+02, 2.09000000e+02,
2.18000000e+02, 2.44000000e+02],
[ 8.00000000e+00, 2.40000000e+01, 2.90000000e+01,
2.80000000e+01, 5.20000000e+01, 5.20000000e+01,
7.00000000e+01, 8.60000000e+01, 7.80000000e+01,
1.17000000e+02, 9.40000000e+01, 1.39000000e+02,
1.58000000e+02, 1.47000000e+02, 1.77000000e+02,
1.67000000e+02, 2.05000000e+02, 1.93000000e+02,
2.08000000e+02, 2.32000000e+02],
[ 6.00000000e+00, 1.90000000e+01, 2.70000000e+01,
4.50000000e+01, 5.50000000e+01, 5.10000000e+01,
7.10000000e+01, 9.60000000e+01, 1.01000000e+02,
1.26000000e+02, 1.05000000e+02, 1.32000000e+02,
1.35000000e+02, 1.43000000e+02, 1.64000000e+02,
1.81000000e+02, 1.71000000e+02, 2.18000000e+02,
2.03000000e+02, 2.22000000e+02],
[ 4.00000000e+00, 2.40000000e+01, 2.10000000e+01,
3.90000000e+01, 5.40000000e+01, 5.90000000e+01,
8.90000000e+01, 8.20000000e+01, 7.20000000e+01,
9.20000000e+01, 9.90000000e+01, 1.39000000e+02,
1.50000000e+02, 1.61000000e+02, 1.41000000e+02,
1.77000000e+02, 2.16000000e+02, 1.95000000e+02,
2.12000000e+02, 2.02000000e+02],
[ 5.00000000e+00, 1.20000000e+01, 4.40000000e+01,
5.30000000e+01, 4.80000000e+01, 5.40000000e+01,
8.10000000e+01, 8.40000000e+01, 9.80000000e+01,
8.10000000e+01, 1.46000000e+02, 1.31000000e+02,
1.36000000e+02, 1.48000000e+02, 1.64000000e+02,
1.53000000e+02, 2.06000000e+02, 2.26000000e+02,
2.17000000e+02, 2.06000000e+02],
[ 3.00000000e+00, 1.40000000e+01, 2.20000000e+01,
5.20000000e+01, 4.30000000e+01, 5.30000000e+01,
7.10000000e+01, 9.00000000e+01, 8.90000000e+01,
9.90000000e+01, 1.20000000e+02, 1.18000000e+02,
1.29000000e+02, 1.50000000e+02, 1.59000000e+02,
1.56000000e+02, 2.22000000e+02, 2.04000000e+02,
2.13000000e+02, 2.43000000e+02],
[ 6.00000000e+00, 1.70000000e+01, 3.00000000e+01,
3.10000000e+01, 4.90000000e+01, 6.60000000e+01,
8.50000000e+01, 9.10000000e+01, 1.12000000e+02,
1.11000000e+02, 1.37000000e+02, 1.32000000e+02,
1.50000000e+02, 1.57000000e+02, 1.58000000e+02,
1.76000000e+02, 1.72000000e+02, 1.90000000e+02,
2.07000000e+02, 2.18000000e+02],
[ 3.00000000e+00, 1.30000000e+01, 3.20000000e+01,
4.50000000e+01, 5.20000000e+01, 5.20000000e+01,
6.30000000e+01, 7.60000000e+01, 1.00000000e+02,
9.80000000e+01, 1.38000000e+02, 1.24000000e+02,
1.21000000e+02, 1.45000000e+02, 1.52000000e+02,
1.73000000e+02, 1.86000000e+02, 1.74000000e+02,
1.78000000e+02, 1.93000000e+02],
[ 6.00000000e+00, 2.50000000e+01, 2.70000000e+01,
4.00000000e+01, 5.20000000e+01, 7.40000000e+01,
5.50000000e+01, 9.30000000e+01, 9.70000000e+01,
9.50000000e+01, 1.21000000e+02, 1.06000000e+02,
1.38000000e+02, 1.60000000e+02, 1.65000000e+02,
1.77000000e+02, 2.04000000e+02, 1.93000000e+02,
1.88000000e+02, 1.97000000e+02],
[ 7.00000000e+00, 2.30000000e+01, 3.50000000e+01,
4.40000000e+01, 5.10000000e+01, 4.90000000e+01,
5.40000000e+01, 8.60000000e+01, 1.02000000e+02,
1.06000000e+02, 1.03000000e+02, 1.31000000e+02,
1.20000000e+02, 1.50000000e+02, 1.65000000e+02,
1.73000000e+02, 1.87000000e+02, 2.11000000e+02,
1.94000000e+02, 2.15000000e+02],
[ 6.00000000e+00, 1.10000000e+01, 2.50000000e+01,
5.60000000e+01, 7.40000000e+01, 6.10000000e+01,
9.00000000e+01, 8.30000000e+01, 8.30000000e+01,
1.13000000e+02, 1.07000000e+02, 1.30000000e+02,
1.29000000e+02, 1.64000000e+02, 1.36000000e+02,
1.65000000e+02, 1.71000000e+02, 1.88000000e+02,
2.04000000e+02, 2.14000000e+02],
[ 4.00000000e+00, 1.30000000e+01, 2.70000000e+01,
3.70000000e+01, 5.00000000e+01, 5.30000000e+01,
7.80000000e+01, 8.80000000e+01, 8.20000000e+01,
1.02000000e+02, 1.02000000e+02, 1.21000000e+02,
1.44000000e+02, 1.43000000e+02, 1.49000000e+02,
1.90000000e+02, 1.75000000e+02, 1.98000000e+02,
2.17000000e+02, 2.37000000e+02],
[ 9.00000000e+00, 1.00000000e+01, 3.40000000e+01,
3.40000000e+01, 4.00000000e+01, 7.00000000e+01,
6.00000000e+01, 8.90000000e+01, 8.00000000e+01,
1.11000000e+02, 1.25000000e+02, 1.39000000e+02,
1.50000000e+02, 1.66000000e+02, 1.45000000e+02,
1.67000000e+02, 1.76000000e+02, 2.01000000e+02,
2.14000000e+02, 2.27000000e+02],
[ 5.00000000e+00, 1.80000000e+01, 1.80000000e+01,
4.90000000e+01, 5.80000000e+01, 5.40000000e+01,
7.20000000e+01, 7.50000000e+01, 9.50000000e+01,
1.14000000e+02, 1.23000000e+02, 1.19000000e+02,
1.18000000e+02, 1.42000000e+02, 1.51000000e+02,
1.74000000e+02, 1.62000000e+02, 1.96000000e+02,
2.00000000e+02, 2.34000000e+02],
[ 4.00000000e+00, 2.00000000e+01, 2.90000000e+01,
3.90000000e+01, 5.10000000e+01, 5.80000000e+01,
8.70000000e+01, 8.80000000e+01, 9.10000000e+01,
8.50000000e+01, 1.06000000e+02, 1.22000000e+02,
1.37000000e+02, 1.37000000e+02, 1.53000000e+02,
1.77000000e+02, 1.60000000e+02, 2.05000000e+02,
1.99000000e+02, 1.90000000e+02],
[ 3.00000000e+00, 2.30000000e+01, 3.40000000e+01,
3.80000000e+01, 5.40000000e+01, 5.20000000e+01,
5.90000000e+01, 9.30000000e+01, 7.40000000e+01,
9.80000000e+01, 9.50000000e+01, 1.19000000e+02,
1.36000000e+02, 1.48000000e+02, 1.69000000e+02,
1.64000000e+02, 1.84000000e+02, 1.76000000e+02,
2.14000000e+02, 2.00000000e+02],
[ 2.00000000e+00, 1.50000000e+01, 3.00000000e+01,
3.60000000e+01, 5.20000000e+01, 6.30000000e+01,
6.80000000e+01, 9.60000000e+01, 1.05000000e+02,
9.30000000e+01, 1.12000000e+02, 1.40000000e+02,
1.37000000e+02, 1.34000000e+02, 1.49000000e+02,
1.47000000e+02, 1.75000000e+02, 1.88000000e+02,
1.99000000e+02, 2.02000000e+02],
[ 5.00000000e+00, 1.60000000e+01, 2.50000000e+01,
3.40000000e+01, 4.80000000e+01, 7.00000000e+01,
7.40000000e+01, 8.20000000e+01, 7.90000000e+01,
8.90000000e+01, 1.15000000e+02, 1.21000000e+02,
1.17000000e+02, 1.51000000e+02, 1.74000000e+02,
1.52000000e+02, 1.76000000e+02, 1.81000000e+02,
1.98000000e+02, 2.01000000e+02]])
In [ ]:
In [28]:
# Getting back the objects:
with open('../output/randDR12Qbin2.pkl') as f: # Python 3: open(..., 'rb')
dat = pickle.load(f)
dat
Out[28]:
array([[ 1.618846, 3.273344, 0.29723 ],
[ 1.576134, 3.051385, 0.946603],
[ 1.710103, 3.953683, 0.110164],
...,
[ 1.621909, 3.488914, 0.172668],
[ 1.569409, 0.428843, 0.035387],
[ 1.849221, 2.659185, 0.276303]])
In [29]:
rr2d=np.zeros((20,20))
In [30]:
len(dat)
Out[30]:
42650
In [31]:
rr2d
Out[31]:
array([[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.]])
In [32]:
%%time
dist0=d.cdist([dat[0],],dat,APdz)[0]
dist1=d.cdist([dat[0],],dat,APzdth)[0]
print np.histogram2d(dist0, dist1, bins=20, range=rng)[0]
[[ 1. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]]
CPU times: user 457 ms, sys: 39.1 ms, total: 496 ms
Wall time: 499 ms
In [33]:
%%time
while(len(dat))>0:
dist0=d.cdist([dat[0],],dat,APdz)[0]
dist1=d.cdist([dat[0],],dat,APzdth)[0]
rr2d+=np.histogram2d(dist0, dist1, bins=20, range=rng)[0]
dat=np.delete(dat,0,axis=0)
if len(dat)%1000==0:
print len(dat)/1000
print rr2d
42
41
40
39
38
37
36
35
34
33
32
31
30
29
28
27
26
25
24
23
22
21
20
19
18
17
16
15
14
13
12
11
10
9
8
7
6
5
4
3
2
1
0
[[ 4.05300000e+04 5.00000000e+00 6.00000000e+00 7.00000000e+00
7.00000000e+00 1.10000000e+01 1.10000000e+01 1.50000000e+01
1.90000000e+01 1.80000000e+01 1.40000000e+01 2.00000000e+01
2.10000000e+01 2.40000000e+01 2.90000000e+01 3.20000000e+01
2.00000000e+01 3.60000000e+01 2.80000000e+01 3.80000000e+01]
[ 0.00000000e+00 5.00000000e+00 5.00000000e+00 2.00000000e+00
8.00000000e+00 1.20000000e+01 1.50000000e+01 1.80000000e+01
1.50000000e+01 1.20000000e+01 1.50000000e+01 2.70000000e+01
2.40000000e+01 1.70000000e+01 2.20000000e+01 2.50000000e+01
3.50000000e+01 4.20000000e+01 2.80000000e+01 4.50000000e+01]
[ 1.00000000e+00 4.00000000e+00 4.00000000e+00 4.00000000e+00
9.00000000e+00 1.20000000e+01 1.20000000e+01 1.40000000e+01
1.20000000e+01 2.20000000e+01 1.50000000e+01 2.10000000e+01
2.40000000e+01 2.30000000e+01 1.50000000e+01 2.40000000e+01
1.80000000e+01 2.50000000e+01 4.10000000e+01 2.70000000e+01]
[ 1.00000000e+00 4.00000000e+00 1.00000000e+00 1.20000000e+01
5.00000000e+00 1.00000000e+01 1.50000000e+01 1.30000000e+01
9.00000000e+00 1.70000000e+01 2.60000000e+01 1.90000000e+01
2.30000000e+01 2.90000000e+01 3.20000000e+01 1.30000000e+01
2.00000000e+01 3.60000000e+01 3.00000000e+01 3.40000000e+01]
[ 0.00000000e+00 1.00000000e+00 4.00000000e+00 9.00000000e+00
1.00000000e+01 9.00000000e+00 1.60000000e+01 1.10000000e+01
1.30000000e+01 1.90000000e+01 1.90000000e+01 2.00000000e+01
2.20000000e+01 2.50000000e+01 2.40000000e+01 2.10000000e+01
2.70000000e+01 2.50000000e+01 3.50000000e+01 3.30000000e+01]
[ 1.00000000e+00 1.00000000e+00 5.00000000e+00 4.00000000e+00
8.00000000e+00 5.00000000e+00 1.50000000e+01 1.00000000e+01
1.50000000e+01 2.80000000e+01 1.50000000e+01 1.60000000e+01
2.60000000e+01 2.10000000e+01 1.70000000e+01 3.30000000e+01
2.90000000e+01 2.50000000e+01 4.10000000e+01 2.70000000e+01]
[ 1.00000000e+00 0.00000000e+00 5.00000000e+00 5.00000000e+00
3.00000000e+00 8.00000000e+00 1.60000000e+01 1.20000000e+01
8.00000000e+00 1.70000000e+01 1.20000000e+01 2.10000000e+01
1.70000000e+01 2.70000000e+01 3.00000000e+01 2.30000000e+01
2.70000000e+01 2.50000000e+01 2.10000000e+01 3.10000000e+01]
[ 0.00000000e+00 2.00000000e+00 4.00000000e+00 4.00000000e+00
5.00000000e+00 8.00000000e+00 6.00000000e+00 1.90000000e+01
1.70000000e+01 1.60000000e+01 1.50000000e+01 2.00000000e+01
2.00000000e+01 2.80000000e+01 1.90000000e+01 1.70000000e+01
2.70000000e+01 3.00000000e+01 3.40000000e+01 3.50000000e+01]
[ 0.00000000e+00 1.00000000e+00 5.00000000e+00 3.00000000e+00
6.00000000e+00 8.00000000e+00 1.40000000e+01 8.00000000e+00
1.40000000e+01 1.50000000e+01 1.60000000e+01 1.80000000e+01
1.40000000e+01 2.30000000e+01 2.50000000e+01 2.00000000e+01
3.10000000e+01 3.40000000e+01 3.50000000e+01 3.10000000e+01]
[ 2.00000000e+00 3.00000000e+00 4.00000000e+00 3.00000000e+00
5.00000000e+00 1.30000000e+01 1.10000000e+01 1.40000000e+01
1.80000000e+01 1.30000000e+01 1.80000000e+01 2.20000000e+01
2.50000000e+01 2.70000000e+01 2.30000000e+01 2.70000000e+01
3.40000000e+01 3.50000000e+01 4.40000000e+01 3.00000000e+01]
[ 0.00000000e+00 5.00000000e+00 3.00000000e+00 7.00000000e+00
4.00000000e+00 7.00000000e+00 1.30000000e+01 1.50000000e+01
1.00000000e+01 1.70000000e+01 1.60000000e+01 2.90000000e+01
2.80000000e+01 3.30000000e+01 2.30000000e+01 1.80000000e+01
3.20000000e+01 2.30000000e+01 3.20000000e+01 3.50000000e+01]
[ 1.00000000e+00 4.00000000e+00 4.00000000e+00 3.00000000e+00
1.00000000e+01 1.10000000e+01 9.00000000e+00 9.00000000e+00
1.70000000e+01 1.60000000e+01 3.00000000e+01 1.60000000e+01
1.70000000e+01 1.70000000e+01 2.90000000e+01 3.10000000e+01
1.60000000e+01 2.60000000e+01 4.50000000e+01 2.60000000e+01]
[ 1.00000000e+00 2.00000000e+00 4.00000000e+00 6.00000000e+00
3.00000000e+00 8.00000000e+00 1.20000000e+01 1.00000000e+01
1.60000000e+01 1.80000000e+01 2.30000000e+01 1.40000000e+01
1.80000000e+01 1.60000000e+01 2.10000000e+01 1.60000000e+01
2.60000000e+01 3.30000000e+01 3.20000000e+01 2.60000000e+01]
[ 1.00000000e+00 4.00000000e+00 5.00000000e+00 6.00000000e+00
9.00000000e+00 1.00000000e+01 1.40000000e+01 1.30000000e+01
1.60000000e+01 1.80000000e+01 9.00000000e+00 2.50000000e+01
1.40000000e+01 2.40000000e+01 2.90000000e+01 2.90000000e+01
2.60000000e+01 3.00000000e+01 3.70000000e+01 2.50000000e+01]
[ 1.00000000e+00 1.00000000e+00 5.00000000e+00 1.10000000e+01
4.00000000e+00 9.00000000e+00 1.10000000e+01 1.40000000e+01
1.30000000e+01 2.30000000e+01 2.20000000e+01 1.40000000e+01
1.80000000e+01 2.60000000e+01 3.70000000e+01 2.30000000e+01
3.20000000e+01 2.40000000e+01 2.80000000e+01 2.60000000e+01]
[ 1.00000000e+00 2.00000000e+00 6.00000000e+00 7.00000000e+00
9.00000000e+00 1.60000000e+01 8.00000000e+00 1.00000000e+01
2.00000000e+01 1.60000000e+01 1.40000000e+01 2.10000000e+01
2.40000000e+01 1.50000000e+01 1.80000000e+01 2.30000000e+01
2.80000000e+01 3.50000000e+01 3.50000000e+01 3.90000000e+01]
[ 3.00000000e+00 3.00000000e+00 4.00000000e+00 7.00000000e+00
8.00000000e+00 2.00000000e+00 1.40000000e+01 1.40000000e+01
1.50000000e+01 1.10000000e+01 2.70000000e+01 2.20000000e+01
1.20000000e+01 1.90000000e+01 2.40000000e+01 2.30000000e+01
2.90000000e+01 2.50000000e+01 2.50000000e+01 2.90000000e+01]
[ 0.00000000e+00 2.00000000e+00 1.00000000e+00 8.00000000e+00
1.20000000e+01 7.00000000e+00 1.10000000e+01 1.40000000e+01
1.50000000e+01 1.70000000e+01 2.20000000e+01 1.30000000e+01
2.40000000e+01 2.20000000e+01 2.90000000e+01 2.40000000e+01
3.60000000e+01 2.10000000e+01 2.70000000e+01 3.50000000e+01]
[ 0.00000000e+00 4.00000000e+00 3.00000000e+00 7.00000000e+00
8.00000000e+00 8.00000000e+00 9.00000000e+00 1.60000000e+01
1.40000000e+01 1.70000000e+01 1.90000000e+01 1.20000000e+01
1.80000000e+01 2.70000000e+01 2.50000000e+01 1.70000000e+01
2.50000000e+01 2.90000000e+01 3.60000000e+01 3.30000000e+01]
[ 2.00000000e+00 1.00000000e+00 7.00000000e+00 7.00000000e+00
9.00000000e+00 9.00000000e+00 1.30000000e+01 8.00000000e+00
1.10000000e+01 1.10000000e+01 1.40000000e+01 1.30000000e+01
2.20000000e+01 1.70000000e+01 1.50000000e+01 3.10000000e+01
2.30000000e+01 3.80000000e+01 2.60000000e+01 3.40000000e+01]]
CPU times: user 2h 51min 31s, sys: 1min 34s, total: 2h 53min 6s
Wall time: 2h 52min 45s
In [34]:
rr2d
Out[34]:
array([[ 4.05300000e+04, 5.00000000e+00, 6.00000000e+00,
7.00000000e+00, 7.00000000e+00, 1.10000000e+01,
1.10000000e+01, 1.50000000e+01, 1.90000000e+01,
1.80000000e+01, 1.40000000e+01, 2.00000000e+01,
2.10000000e+01, 2.40000000e+01, 2.90000000e+01,
3.20000000e+01, 2.00000000e+01, 3.60000000e+01,
2.80000000e+01, 3.80000000e+01],
[ 0.00000000e+00, 5.00000000e+00, 5.00000000e+00,
2.00000000e+00, 8.00000000e+00, 1.20000000e+01,
1.50000000e+01, 1.80000000e+01, 1.50000000e+01,
1.20000000e+01, 1.50000000e+01, 2.70000000e+01,
2.40000000e+01, 1.70000000e+01, 2.20000000e+01,
2.50000000e+01, 3.50000000e+01, 4.20000000e+01,
2.80000000e+01, 4.50000000e+01],
[ 1.00000000e+00, 4.00000000e+00, 4.00000000e+00,
4.00000000e+00, 9.00000000e+00, 1.20000000e+01,
1.20000000e+01, 1.40000000e+01, 1.20000000e+01,
2.20000000e+01, 1.50000000e+01, 2.10000000e+01,
2.40000000e+01, 2.30000000e+01, 1.50000000e+01,
2.40000000e+01, 1.80000000e+01, 2.50000000e+01,
4.10000000e+01, 2.70000000e+01],
[ 1.00000000e+00, 4.00000000e+00, 1.00000000e+00,
1.20000000e+01, 5.00000000e+00, 1.00000000e+01,
1.50000000e+01, 1.30000000e+01, 9.00000000e+00,
1.70000000e+01, 2.60000000e+01, 1.90000000e+01,
2.30000000e+01, 2.90000000e+01, 3.20000000e+01,
1.30000000e+01, 2.00000000e+01, 3.60000000e+01,
3.00000000e+01, 3.40000000e+01],
[ 0.00000000e+00, 1.00000000e+00, 4.00000000e+00,
9.00000000e+00, 1.00000000e+01, 9.00000000e+00,
1.60000000e+01, 1.10000000e+01, 1.30000000e+01,
1.90000000e+01, 1.90000000e+01, 2.00000000e+01,
2.20000000e+01, 2.50000000e+01, 2.40000000e+01,
2.10000000e+01, 2.70000000e+01, 2.50000000e+01,
3.50000000e+01, 3.30000000e+01],
[ 1.00000000e+00, 1.00000000e+00, 5.00000000e+00,
4.00000000e+00, 8.00000000e+00, 5.00000000e+00,
1.50000000e+01, 1.00000000e+01, 1.50000000e+01,
2.80000000e+01, 1.50000000e+01, 1.60000000e+01,
2.60000000e+01, 2.10000000e+01, 1.70000000e+01,
3.30000000e+01, 2.90000000e+01, 2.50000000e+01,
4.10000000e+01, 2.70000000e+01],
[ 1.00000000e+00, 0.00000000e+00, 5.00000000e+00,
5.00000000e+00, 3.00000000e+00, 8.00000000e+00,
1.60000000e+01, 1.20000000e+01, 8.00000000e+00,
1.70000000e+01, 1.20000000e+01, 2.10000000e+01,
1.70000000e+01, 2.70000000e+01, 3.00000000e+01,
2.30000000e+01, 2.70000000e+01, 2.50000000e+01,
2.10000000e+01, 3.10000000e+01],
[ 0.00000000e+00, 2.00000000e+00, 4.00000000e+00,
4.00000000e+00, 5.00000000e+00, 8.00000000e+00,
6.00000000e+00, 1.90000000e+01, 1.70000000e+01,
1.60000000e+01, 1.50000000e+01, 2.00000000e+01,
2.00000000e+01, 2.80000000e+01, 1.90000000e+01,
1.70000000e+01, 2.70000000e+01, 3.00000000e+01,
3.40000000e+01, 3.50000000e+01],
[ 0.00000000e+00, 1.00000000e+00, 5.00000000e+00,
3.00000000e+00, 6.00000000e+00, 8.00000000e+00,
1.40000000e+01, 8.00000000e+00, 1.40000000e+01,
1.50000000e+01, 1.60000000e+01, 1.80000000e+01,
1.40000000e+01, 2.30000000e+01, 2.50000000e+01,
2.00000000e+01, 3.10000000e+01, 3.40000000e+01,
3.50000000e+01, 3.10000000e+01],
[ 2.00000000e+00, 3.00000000e+00, 4.00000000e+00,
3.00000000e+00, 5.00000000e+00, 1.30000000e+01,
1.10000000e+01, 1.40000000e+01, 1.80000000e+01,
1.30000000e+01, 1.80000000e+01, 2.20000000e+01,
2.50000000e+01, 2.70000000e+01, 2.30000000e+01,
2.70000000e+01, 3.40000000e+01, 3.50000000e+01,
4.40000000e+01, 3.00000000e+01],
[ 0.00000000e+00, 5.00000000e+00, 3.00000000e+00,
7.00000000e+00, 4.00000000e+00, 7.00000000e+00,
1.30000000e+01, 1.50000000e+01, 1.00000000e+01,
1.70000000e+01, 1.60000000e+01, 2.90000000e+01,
2.80000000e+01, 3.30000000e+01, 2.30000000e+01,
1.80000000e+01, 3.20000000e+01, 2.30000000e+01,
3.20000000e+01, 3.50000000e+01],
[ 1.00000000e+00, 4.00000000e+00, 4.00000000e+00,
3.00000000e+00, 1.00000000e+01, 1.10000000e+01,
9.00000000e+00, 9.00000000e+00, 1.70000000e+01,
1.60000000e+01, 3.00000000e+01, 1.60000000e+01,
1.70000000e+01, 1.70000000e+01, 2.90000000e+01,
3.10000000e+01, 1.60000000e+01, 2.60000000e+01,
4.50000000e+01, 2.60000000e+01],
[ 1.00000000e+00, 2.00000000e+00, 4.00000000e+00,
6.00000000e+00, 3.00000000e+00, 8.00000000e+00,
1.20000000e+01, 1.00000000e+01, 1.60000000e+01,
1.80000000e+01, 2.30000000e+01, 1.40000000e+01,
1.80000000e+01, 1.60000000e+01, 2.10000000e+01,
1.60000000e+01, 2.60000000e+01, 3.30000000e+01,
3.20000000e+01, 2.60000000e+01],
[ 1.00000000e+00, 4.00000000e+00, 5.00000000e+00,
6.00000000e+00, 9.00000000e+00, 1.00000000e+01,
1.40000000e+01, 1.30000000e+01, 1.60000000e+01,
1.80000000e+01, 9.00000000e+00, 2.50000000e+01,
1.40000000e+01, 2.40000000e+01, 2.90000000e+01,
2.90000000e+01, 2.60000000e+01, 3.00000000e+01,
3.70000000e+01, 2.50000000e+01],
[ 1.00000000e+00, 1.00000000e+00, 5.00000000e+00,
1.10000000e+01, 4.00000000e+00, 9.00000000e+00,
1.10000000e+01, 1.40000000e+01, 1.30000000e+01,
2.30000000e+01, 2.20000000e+01, 1.40000000e+01,
1.80000000e+01, 2.60000000e+01, 3.70000000e+01,
2.30000000e+01, 3.20000000e+01, 2.40000000e+01,
2.80000000e+01, 2.60000000e+01],
[ 1.00000000e+00, 2.00000000e+00, 6.00000000e+00,
7.00000000e+00, 9.00000000e+00, 1.60000000e+01,
8.00000000e+00, 1.00000000e+01, 2.00000000e+01,
1.60000000e+01, 1.40000000e+01, 2.10000000e+01,
2.40000000e+01, 1.50000000e+01, 1.80000000e+01,
2.30000000e+01, 2.80000000e+01, 3.50000000e+01,
3.50000000e+01, 3.90000000e+01],
[ 3.00000000e+00, 3.00000000e+00, 4.00000000e+00,
7.00000000e+00, 8.00000000e+00, 2.00000000e+00,
1.40000000e+01, 1.40000000e+01, 1.50000000e+01,
1.10000000e+01, 2.70000000e+01, 2.20000000e+01,
1.20000000e+01, 1.90000000e+01, 2.40000000e+01,
2.30000000e+01, 2.90000000e+01, 2.50000000e+01,
2.50000000e+01, 2.90000000e+01],
[ 0.00000000e+00, 2.00000000e+00, 1.00000000e+00,
8.00000000e+00, 1.20000000e+01, 7.00000000e+00,
1.10000000e+01, 1.40000000e+01, 1.50000000e+01,
1.70000000e+01, 2.20000000e+01, 1.30000000e+01,
2.40000000e+01, 2.20000000e+01, 2.90000000e+01,
2.40000000e+01, 3.60000000e+01, 2.10000000e+01,
2.70000000e+01, 3.50000000e+01],
[ 0.00000000e+00, 4.00000000e+00, 3.00000000e+00,
7.00000000e+00, 8.00000000e+00, 8.00000000e+00,
9.00000000e+00, 1.60000000e+01, 1.40000000e+01,
1.70000000e+01, 1.90000000e+01, 1.20000000e+01,
1.80000000e+01, 2.70000000e+01, 2.50000000e+01,
1.70000000e+01, 2.50000000e+01, 2.90000000e+01,
3.60000000e+01, 3.30000000e+01],
[ 2.00000000e+00, 1.00000000e+00, 7.00000000e+00,
7.00000000e+00, 9.00000000e+00, 9.00000000e+00,
1.30000000e+01, 8.00000000e+00, 1.10000000e+01,
1.10000000e+01, 1.40000000e+01, 1.30000000e+01,
2.20000000e+01, 1.70000000e+01, 1.50000000e+01,
3.10000000e+01, 2.30000000e+01, 3.80000000e+01,
2.60000000e+01, 3.40000000e+01]])
In [35]:
with open('randDR12QRRbin2.pkl','w') as f:
pickle.dump(rr2d,f)
rr2d
Out[35]:
array([[ 4.05300000e+04, 5.00000000e+00, 6.00000000e+00,
7.00000000e+00, 7.00000000e+00, 1.10000000e+01,
1.10000000e+01, 1.50000000e+01, 1.90000000e+01,
1.80000000e+01, 1.40000000e+01, 2.00000000e+01,
2.10000000e+01, 2.40000000e+01, 2.90000000e+01,
3.20000000e+01, 2.00000000e+01, 3.60000000e+01,
2.80000000e+01, 3.80000000e+01],
[ 0.00000000e+00, 5.00000000e+00, 5.00000000e+00,
2.00000000e+00, 8.00000000e+00, 1.20000000e+01,
1.50000000e+01, 1.80000000e+01, 1.50000000e+01,
1.20000000e+01, 1.50000000e+01, 2.70000000e+01,
2.40000000e+01, 1.70000000e+01, 2.20000000e+01,
2.50000000e+01, 3.50000000e+01, 4.20000000e+01,
2.80000000e+01, 4.50000000e+01],
[ 1.00000000e+00, 4.00000000e+00, 4.00000000e+00,
4.00000000e+00, 9.00000000e+00, 1.20000000e+01,
1.20000000e+01, 1.40000000e+01, 1.20000000e+01,
2.20000000e+01, 1.50000000e+01, 2.10000000e+01,
2.40000000e+01, 2.30000000e+01, 1.50000000e+01,
2.40000000e+01, 1.80000000e+01, 2.50000000e+01,
4.10000000e+01, 2.70000000e+01],
[ 1.00000000e+00, 4.00000000e+00, 1.00000000e+00,
1.20000000e+01, 5.00000000e+00, 1.00000000e+01,
1.50000000e+01, 1.30000000e+01, 9.00000000e+00,
1.70000000e+01, 2.60000000e+01, 1.90000000e+01,
2.30000000e+01, 2.90000000e+01, 3.20000000e+01,
1.30000000e+01, 2.00000000e+01, 3.60000000e+01,
3.00000000e+01, 3.40000000e+01],
[ 0.00000000e+00, 1.00000000e+00, 4.00000000e+00,
9.00000000e+00, 1.00000000e+01, 9.00000000e+00,
1.60000000e+01, 1.10000000e+01, 1.30000000e+01,
1.90000000e+01, 1.90000000e+01, 2.00000000e+01,
2.20000000e+01, 2.50000000e+01, 2.40000000e+01,
2.10000000e+01, 2.70000000e+01, 2.50000000e+01,
3.50000000e+01, 3.30000000e+01],
[ 1.00000000e+00, 1.00000000e+00, 5.00000000e+00,
4.00000000e+00, 8.00000000e+00, 5.00000000e+00,
1.50000000e+01, 1.00000000e+01, 1.50000000e+01,
2.80000000e+01, 1.50000000e+01, 1.60000000e+01,
2.60000000e+01, 2.10000000e+01, 1.70000000e+01,
3.30000000e+01, 2.90000000e+01, 2.50000000e+01,
4.10000000e+01, 2.70000000e+01],
[ 1.00000000e+00, 0.00000000e+00, 5.00000000e+00,
5.00000000e+00, 3.00000000e+00, 8.00000000e+00,
1.60000000e+01, 1.20000000e+01, 8.00000000e+00,
1.70000000e+01, 1.20000000e+01, 2.10000000e+01,
1.70000000e+01, 2.70000000e+01, 3.00000000e+01,
2.30000000e+01, 2.70000000e+01, 2.50000000e+01,
2.10000000e+01, 3.10000000e+01],
[ 0.00000000e+00, 2.00000000e+00, 4.00000000e+00,
4.00000000e+00, 5.00000000e+00, 8.00000000e+00,
6.00000000e+00, 1.90000000e+01, 1.70000000e+01,
1.60000000e+01, 1.50000000e+01, 2.00000000e+01,
2.00000000e+01, 2.80000000e+01, 1.90000000e+01,
1.70000000e+01, 2.70000000e+01, 3.00000000e+01,
3.40000000e+01, 3.50000000e+01],
[ 0.00000000e+00, 1.00000000e+00, 5.00000000e+00,
3.00000000e+00, 6.00000000e+00, 8.00000000e+00,
1.40000000e+01, 8.00000000e+00, 1.40000000e+01,
1.50000000e+01, 1.60000000e+01, 1.80000000e+01,
1.40000000e+01, 2.30000000e+01, 2.50000000e+01,
2.00000000e+01, 3.10000000e+01, 3.40000000e+01,
3.50000000e+01, 3.10000000e+01],
[ 2.00000000e+00, 3.00000000e+00, 4.00000000e+00,
3.00000000e+00, 5.00000000e+00, 1.30000000e+01,
1.10000000e+01, 1.40000000e+01, 1.80000000e+01,
1.30000000e+01, 1.80000000e+01, 2.20000000e+01,
2.50000000e+01, 2.70000000e+01, 2.30000000e+01,
2.70000000e+01, 3.40000000e+01, 3.50000000e+01,
4.40000000e+01, 3.00000000e+01],
[ 0.00000000e+00, 5.00000000e+00, 3.00000000e+00,
7.00000000e+00, 4.00000000e+00, 7.00000000e+00,
1.30000000e+01, 1.50000000e+01, 1.00000000e+01,
1.70000000e+01, 1.60000000e+01, 2.90000000e+01,
2.80000000e+01, 3.30000000e+01, 2.30000000e+01,
1.80000000e+01, 3.20000000e+01, 2.30000000e+01,
3.20000000e+01, 3.50000000e+01],
[ 1.00000000e+00, 4.00000000e+00, 4.00000000e+00,
3.00000000e+00, 1.00000000e+01, 1.10000000e+01,
9.00000000e+00, 9.00000000e+00, 1.70000000e+01,
1.60000000e+01, 3.00000000e+01, 1.60000000e+01,
1.70000000e+01, 1.70000000e+01, 2.90000000e+01,
3.10000000e+01, 1.60000000e+01, 2.60000000e+01,
4.50000000e+01, 2.60000000e+01],
[ 1.00000000e+00, 2.00000000e+00, 4.00000000e+00,
6.00000000e+00, 3.00000000e+00, 8.00000000e+00,
1.20000000e+01, 1.00000000e+01, 1.60000000e+01,
1.80000000e+01, 2.30000000e+01, 1.40000000e+01,
1.80000000e+01, 1.60000000e+01, 2.10000000e+01,
1.60000000e+01, 2.60000000e+01, 3.30000000e+01,
3.20000000e+01, 2.60000000e+01],
[ 1.00000000e+00, 4.00000000e+00, 5.00000000e+00,
6.00000000e+00, 9.00000000e+00, 1.00000000e+01,
1.40000000e+01, 1.30000000e+01, 1.60000000e+01,
1.80000000e+01, 9.00000000e+00, 2.50000000e+01,
1.40000000e+01, 2.40000000e+01, 2.90000000e+01,
2.90000000e+01, 2.60000000e+01, 3.00000000e+01,
3.70000000e+01, 2.50000000e+01],
[ 1.00000000e+00, 1.00000000e+00, 5.00000000e+00,
1.10000000e+01, 4.00000000e+00, 9.00000000e+00,
1.10000000e+01, 1.40000000e+01, 1.30000000e+01,
2.30000000e+01, 2.20000000e+01, 1.40000000e+01,
1.80000000e+01, 2.60000000e+01, 3.70000000e+01,
2.30000000e+01, 3.20000000e+01, 2.40000000e+01,
2.80000000e+01, 2.60000000e+01],
[ 1.00000000e+00, 2.00000000e+00, 6.00000000e+00,
7.00000000e+00, 9.00000000e+00, 1.60000000e+01,
8.00000000e+00, 1.00000000e+01, 2.00000000e+01,
1.60000000e+01, 1.40000000e+01, 2.10000000e+01,
2.40000000e+01, 1.50000000e+01, 1.80000000e+01,
2.30000000e+01, 2.80000000e+01, 3.50000000e+01,
3.50000000e+01, 3.90000000e+01],
[ 3.00000000e+00, 3.00000000e+00, 4.00000000e+00,
7.00000000e+00, 8.00000000e+00, 2.00000000e+00,
1.40000000e+01, 1.40000000e+01, 1.50000000e+01,
1.10000000e+01, 2.70000000e+01, 2.20000000e+01,
1.20000000e+01, 1.90000000e+01, 2.40000000e+01,
2.30000000e+01, 2.90000000e+01, 2.50000000e+01,
2.50000000e+01, 2.90000000e+01],
[ 0.00000000e+00, 2.00000000e+00, 1.00000000e+00,
8.00000000e+00, 1.20000000e+01, 7.00000000e+00,
1.10000000e+01, 1.40000000e+01, 1.50000000e+01,
1.70000000e+01, 2.20000000e+01, 1.30000000e+01,
2.40000000e+01, 2.20000000e+01, 2.90000000e+01,
2.40000000e+01, 3.60000000e+01, 2.10000000e+01,
2.70000000e+01, 3.50000000e+01],
[ 0.00000000e+00, 4.00000000e+00, 3.00000000e+00,
7.00000000e+00, 8.00000000e+00, 8.00000000e+00,
9.00000000e+00, 1.60000000e+01, 1.40000000e+01,
1.70000000e+01, 1.90000000e+01, 1.20000000e+01,
1.80000000e+01, 2.70000000e+01, 2.50000000e+01,
1.70000000e+01, 2.50000000e+01, 2.90000000e+01,
3.60000000e+01, 3.30000000e+01],
[ 2.00000000e+00, 1.00000000e+00, 7.00000000e+00,
7.00000000e+00, 9.00000000e+00, 9.00000000e+00,
1.30000000e+01, 8.00000000e+00, 1.10000000e+01,
1.10000000e+01, 1.40000000e+01, 1.30000000e+01,
2.20000000e+01, 1.70000000e+01, 1.50000000e+01,
3.10000000e+01, 2.30000000e+01, 3.80000000e+01,
2.60000000e+01, 3.40000000e+01]])
In [ ]:
In [36]:
# Getting back the objects:
with open('../output/randDR12Qbin3.pkl') as f: # Python 3: open(..., 'rb')
dat = pickle.load(f)
dat
Out[36]:
array([[ 2.30909700e+00, 3.30000000e-05, 3.10210000e-01],
[ 2.49794100e+00, 4.80000000e-05, 2.61357000e-01],
[ 2.33265500e+00, 1.00000000e-04, -2.31260000e-02],
...,
[ 2.41550800e+00, 6.28316600e+00, 1.76571000e-01],
[ 2.45101500e+00, 6.28317000e+00, 5.05355000e-01],
[ 2.39766700e+00, 6.28318400e+00, 6.06452000e-01]])
In [37]:
rr2d=np.zeros((20,20))
In [38]:
len(dat)
Out[38]:
139641
In [39]:
rr2d
Out[39]:
array([[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.]])
In [40]:
%%time
dist0=d.cdist([dat[0],],dat,APdz)[0]
dist1=d.cdist([dat[0],],dat,APzdth)[0]
print np.histogram2d(dist0, dist1, bins=20, range=rng)[0]
[[ 1. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 1. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]]
CPU times: user 1.31 s, sys: 116 ms, total: 1.43 s
Wall time: 1.47 s
In [41]:
%%time
while(len(dat))>0:
dist0=d.cdist([dat[0],],dat,APdz)[0]
dist1=d.cdist([dat[0],],dat,APzdth)[0]
rr2d+=np.histogram2d(dist0, dist1, bins=20, range=rng)[0]
dat=np.delete(dat,0,axis=0)
if len(dat)%1000==0:
print len(dat)/1000
print rr2d
139
138
137
136
135
134
133
132
131
130
129
128
127
126
125
124
123
122
121
120
119
118
117
116
115
114
113
112
111
110
109
108
107
106
105
104
103
102
101
100
99
98
97
96
95
94
93
92
91
90
89
88
87
86
85
84
83
82
81
80
79
78
77
76
75
74
73
72
71
70
69
68
67
66
65
64
63
62
61
60
59
58
57
56
55
54
53
52
51
50
49
48
47
46
45
44
43
42
41
40
39
38
37
36
35
34
33
32
31
30
29
28
27
26
25
24
23
22
21
20
19
18
17
16
15
14
13
12
11
10
9
8
7
6
5
4
3
2
1
0
[[ 1.32956000e+05 5.40000000e+01 6.30000000e+01 8.40000000e+01
8.20000000e+01 1.08000000e+02 1.04000000e+02 1.04000000e+02
1.20000000e+02 1.65000000e+02 1.63000000e+02 1.58000000e+02
1.78000000e+02 1.63000000e+02 1.58000000e+02 1.96000000e+02
1.93000000e+02 2.31000000e+02 2.42000000e+02 2.71000000e+02]
[ 2.00000000e+01 6.00000000e+01 8.10000000e+01 7.10000000e+01
1.07000000e+02 9.70000000e+01 9.40000000e+01 1.07000000e+02
1.54000000e+02 1.51000000e+02 1.29000000e+02 1.75000000e+02
1.57000000e+02 1.91000000e+02 1.85000000e+02 2.09000000e+02
2.07000000e+02 2.22000000e+02 2.60000000e+02 2.59000000e+02]
[ 1.70000000e+01 5.80000000e+01 6.80000000e+01 8.10000000e+01
9.80000000e+01 1.21000000e+02 8.90000000e+01 1.11000000e+02
1.42000000e+02 1.24000000e+02 1.54000000e+02 1.46000000e+02
1.93000000e+02 2.01000000e+02 1.79000000e+02 2.01000000e+02
2.11000000e+02 2.23000000e+02 2.36000000e+02 2.65000000e+02]
[ 2.40000000e+01 5.90000000e+01 6.90000000e+01 7.40000000e+01
7.90000000e+01 1.06000000e+02 1.07000000e+02 1.16000000e+02
1.25000000e+02 1.33000000e+02 1.60000000e+02 1.76000000e+02
1.75000000e+02 1.37000000e+02 1.79000000e+02 1.88000000e+02
2.04000000e+02 2.29000000e+02 2.46000000e+02 2.59000000e+02]
[ 1.60000000e+01 5.40000000e+01 6.00000000e+01 8.10000000e+01
8.60000000e+01 8.10000000e+01 1.01000000e+02 1.30000000e+02
1.15000000e+02 1.48000000e+02 1.58000000e+02 1.58000000e+02
1.63000000e+02 1.81000000e+02 1.85000000e+02 2.09000000e+02
2.06000000e+02 2.24000000e+02 2.14000000e+02 2.57000000e+02]
[ 2.70000000e+01 5.00000000e+01 6.30000000e+01 8.00000000e+01
7.90000000e+01 9.70000000e+01 1.11000000e+02 1.35000000e+02
1.26000000e+02 1.58000000e+02 1.46000000e+02 1.79000000e+02
1.69000000e+02 1.86000000e+02 1.78000000e+02 2.01000000e+02
2.22000000e+02 2.56000000e+02 2.29000000e+02 2.36000000e+02]
[ 2.20000000e+01 5.70000000e+01 7.50000000e+01 7.60000000e+01
1.00000000e+02 9.30000000e+01 1.05000000e+02 9.50000000e+01
1.14000000e+02 1.36000000e+02 1.42000000e+02 1.56000000e+02
1.65000000e+02 2.02000000e+02 1.63000000e+02 1.98000000e+02
2.25000000e+02 2.29000000e+02 2.16000000e+02 2.67000000e+02]
[ 1.80000000e+01 3.40000000e+01 6.20000000e+01 8.10000000e+01
9.40000000e+01 8.90000000e+01 1.13000000e+02 1.26000000e+02
1.22000000e+02 1.47000000e+02 1.45000000e+02 1.52000000e+02
1.65000000e+02 1.88000000e+02 1.94000000e+02 1.94000000e+02
2.17000000e+02 2.17000000e+02 2.12000000e+02 2.44000000e+02]
[ 1.40000000e+01 5.30000000e+01 5.60000000e+01 6.00000000e+01
7.50000000e+01 9.40000000e+01 9.40000000e+01 1.04000000e+02
1.24000000e+02 1.48000000e+02 1.58000000e+02 1.75000000e+02
1.68000000e+02 1.75000000e+02 1.89000000e+02 2.26000000e+02
2.09000000e+02 2.61000000e+02 2.15000000e+02 2.64000000e+02]
[ 2.10000000e+01 3.80000000e+01 6.00000000e+01 7.70000000e+01
7.50000000e+01 8.30000000e+01 1.19000000e+02 1.10000000e+02
1.14000000e+02 1.39000000e+02 1.53000000e+02 1.65000000e+02
1.59000000e+02 1.81000000e+02 1.98000000e+02 2.28000000e+02
2.09000000e+02 2.12000000e+02 2.38000000e+02 2.34000000e+02]
[ 2.10000000e+01 2.70000000e+01 7.00000000e+01 7.20000000e+01
8.60000000e+01 8.10000000e+01 1.09000000e+02 1.08000000e+02
1.19000000e+02 1.69000000e+02 1.59000000e+02 1.57000000e+02
1.63000000e+02 1.84000000e+02 2.00000000e+02 1.91000000e+02
2.08000000e+02 2.03000000e+02 2.17000000e+02 2.32000000e+02]
[ 1.80000000e+01 4.80000000e+01 4.50000000e+01 6.50000000e+01
8.70000000e+01 9.00000000e+01 1.14000000e+02 1.11000000e+02
1.34000000e+02 1.19000000e+02 1.43000000e+02 1.54000000e+02
1.76000000e+02 1.64000000e+02 1.88000000e+02 2.13000000e+02
1.99000000e+02 2.14000000e+02 2.49000000e+02 2.46000000e+02]
[ 1.00000000e+01 4.80000000e+01 4.80000000e+01 7.20000000e+01
8.10000000e+01 1.06000000e+02 1.14000000e+02 1.05000000e+02
1.32000000e+02 1.30000000e+02 1.61000000e+02 1.61000000e+02
1.87000000e+02 1.77000000e+02 1.91000000e+02 2.19000000e+02
2.10000000e+02 2.34000000e+02 2.15000000e+02 2.21000000e+02]
[ 1.40000000e+01 3.50000000e+01 4.50000000e+01 6.30000000e+01
7.70000000e+01 8.90000000e+01 9.10000000e+01 8.80000000e+01
1.23000000e+02 1.48000000e+02 1.56000000e+02 1.51000000e+02
1.54000000e+02 1.77000000e+02 1.88000000e+02 2.09000000e+02
1.87000000e+02 2.13000000e+02 2.28000000e+02 2.24000000e+02]
[ 1.30000000e+01 4.00000000e+01 5.30000000e+01 6.40000000e+01
7.50000000e+01 8.80000000e+01 1.04000000e+02 1.31000000e+02
1.16000000e+02 1.43000000e+02 1.39000000e+02 1.48000000e+02
1.78000000e+02 1.69000000e+02 1.88000000e+02 1.84000000e+02
2.02000000e+02 2.18000000e+02 2.21000000e+02 2.04000000e+02]
[ 1.30000000e+01 2.90000000e+01 5.10000000e+01 5.30000000e+01
6.10000000e+01 8.70000000e+01 1.07000000e+02 1.18000000e+02
1.20000000e+02 1.21000000e+02 1.39000000e+02 1.30000000e+02
1.71000000e+02 1.57000000e+02 1.86000000e+02 1.87000000e+02
2.06000000e+02 2.12000000e+02 1.96000000e+02 2.44000000e+02]
[ 1.00000000e+01 3.70000000e+01 5.10000000e+01 6.80000000e+01
8.80000000e+01 7.30000000e+01 6.90000000e+01 1.02000000e+02
1.09000000e+02 1.25000000e+02 1.32000000e+02 1.47000000e+02
1.93000000e+02 1.91000000e+02 1.72000000e+02 1.88000000e+02
1.96000000e+02 2.39000000e+02 2.39000000e+02 2.42000000e+02]
[ 1.50000000e+01 2.50000000e+01 5.20000000e+01 5.90000000e+01
7.30000000e+01 6.90000000e+01 8.30000000e+01 9.80000000e+01
1.03000000e+02 1.37000000e+02 1.78000000e+02 1.51000000e+02
1.73000000e+02 1.65000000e+02 1.88000000e+02 2.03000000e+02
1.97000000e+02 2.07000000e+02 2.31000000e+02 2.41000000e+02]
[ 1.60000000e+01 2.60000000e+01 5.40000000e+01 4.70000000e+01
6.00000000e+01 7.30000000e+01 1.06000000e+02 1.05000000e+02
1.14000000e+02 1.27000000e+02 1.45000000e+02 1.49000000e+02
1.59000000e+02 1.94000000e+02 1.68000000e+02 2.08000000e+02
2.32000000e+02 2.26000000e+02 2.04000000e+02 2.33000000e+02]
[ 3.00000000e+00 2.90000000e+01 4.20000000e+01 6.20000000e+01
6.50000000e+01 7.40000000e+01 9.00000000e+01 8.30000000e+01
1.28000000e+02 1.19000000e+02 1.50000000e+02 1.18000000e+02
1.71000000e+02 1.75000000e+02 1.76000000e+02 1.60000000e+02
2.24000000e+02 2.46000000e+02 2.28000000e+02 2.19000000e+02]]
CPU times: user 1d 3h 9min 39s, sys: 29min 10s, total: 1d 3h 38min 50s
Wall time: 1d 10h 53min 25s
In [42]:
rr2d
Out[42]:
array([[ 1.32956000e+05, 5.40000000e+01, 6.30000000e+01,
8.40000000e+01, 8.20000000e+01, 1.08000000e+02,
1.04000000e+02, 1.04000000e+02, 1.20000000e+02,
1.65000000e+02, 1.63000000e+02, 1.58000000e+02,
1.78000000e+02, 1.63000000e+02, 1.58000000e+02,
1.96000000e+02, 1.93000000e+02, 2.31000000e+02,
2.42000000e+02, 2.71000000e+02],
[ 2.00000000e+01, 6.00000000e+01, 8.10000000e+01,
7.10000000e+01, 1.07000000e+02, 9.70000000e+01,
9.40000000e+01, 1.07000000e+02, 1.54000000e+02,
1.51000000e+02, 1.29000000e+02, 1.75000000e+02,
1.57000000e+02, 1.91000000e+02, 1.85000000e+02,
2.09000000e+02, 2.07000000e+02, 2.22000000e+02,
2.60000000e+02, 2.59000000e+02],
[ 1.70000000e+01, 5.80000000e+01, 6.80000000e+01,
8.10000000e+01, 9.80000000e+01, 1.21000000e+02,
8.90000000e+01, 1.11000000e+02, 1.42000000e+02,
1.24000000e+02, 1.54000000e+02, 1.46000000e+02,
1.93000000e+02, 2.01000000e+02, 1.79000000e+02,
2.01000000e+02, 2.11000000e+02, 2.23000000e+02,
2.36000000e+02, 2.65000000e+02],
[ 2.40000000e+01, 5.90000000e+01, 6.90000000e+01,
7.40000000e+01, 7.90000000e+01, 1.06000000e+02,
1.07000000e+02, 1.16000000e+02, 1.25000000e+02,
1.33000000e+02, 1.60000000e+02, 1.76000000e+02,
1.75000000e+02, 1.37000000e+02, 1.79000000e+02,
1.88000000e+02, 2.04000000e+02, 2.29000000e+02,
2.46000000e+02, 2.59000000e+02],
[ 1.60000000e+01, 5.40000000e+01, 6.00000000e+01,
8.10000000e+01, 8.60000000e+01, 8.10000000e+01,
1.01000000e+02, 1.30000000e+02, 1.15000000e+02,
1.48000000e+02, 1.58000000e+02, 1.58000000e+02,
1.63000000e+02, 1.81000000e+02, 1.85000000e+02,
2.09000000e+02, 2.06000000e+02, 2.24000000e+02,
2.14000000e+02, 2.57000000e+02],
[ 2.70000000e+01, 5.00000000e+01, 6.30000000e+01,
8.00000000e+01, 7.90000000e+01, 9.70000000e+01,
1.11000000e+02, 1.35000000e+02, 1.26000000e+02,
1.58000000e+02, 1.46000000e+02, 1.79000000e+02,
1.69000000e+02, 1.86000000e+02, 1.78000000e+02,
2.01000000e+02, 2.22000000e+02, 2.56000000e+02,
2.29000000e+02, 2.36000000e+02],
[ 2.20000000e+01, 5.70000000e+01, 7.50000000e+01,
7.60000000e+01, 1.00000000e+02, 9.30000000e+01,
1.05000000e+02, 9.50000000e+01, 1.14000000e+02,
1.36000000e+02, 1.42000000e+02, 1.56000000e+02,
1.65000000e+02, 2.02000000e+02, 1.63000000e+02,
1.98000000e+02, 2.25000000e+02, 2.29000000e+02,
2.16000000e+02, 2.67000000e+02],
[ 1.80000000e+01, 3.40000000e+01, 6.20000000e+01,
8.10000000e+01, 9.40000000e+01, 8.90000000e+01,
1.13000000e+02, 1.26000000e+02, 1.22000000e+02,
1.47000000e+02, 1.45000000e+02, 1.52000000e+02,
1.65000000e+02, 1.88000000e+02, 1.94000000e+02,
1.94000000e+02, 2.17000000e+02, 2.17000000e+02,
2.12000000e+02, 2.44000000e+02],
[ 1.40000000e+01, 5.30000000e+01, 5.60000000e+01,
6.00000000e+01, 7.50000000e+01, 9.40000000e+01,
9.40000000e+01, 1.04000000e+02, 1.24000000e+02,
1.48000000e+02, 1.58000000e+02, 1.75000000e+02,
1.68000000e+02, 1.75000000e+02, 1.89000000e+02,
2.26000000e+02, 2.09000000e+02, 2.61000000e+02,
2.15000000e+02, 2.64000000e+02],
[ 2.10000000e+01, 3.80000000e+01, 6.00000000e+01,
7.70000000e+01, 7.50000000e+01, 8.30000000e+01,
1.19000000e+02, 1.10000000e+02, 1.14000000e+02,
1.39000000e+02, 1.53000000e+02, 1.65000000e+02,
1.59000000e+02, 1.81000000e+02, 1.98000000e+02,
2.28000000e+02, 2.09000000e+02, 2.12000000e+02,
2.38000000e+02, 2.34000000e+02],
[ 2.10000000e+01, 2.70000000e+01, 7.00000000e+01,
7.20000000e+01, 8.60000000e+01, 8.10000000e+01,
1.09000000e+02, 1.08000000e+02, 1.19000000e+02,
1.69000000e+02, 1.59000000e+02, 1.57000000e+02,
1.63000000e+02, 1.84000000e+02, 2.00000000e+02,
1.91000000e+02, 2.08000000e+02, 2.03000000e+02,
2.17000000e+02, 2.32000000e+02],
[ 1.80000000e+01, 4.80000000e+01, 4.50000000e+01,
6.50000000e+01, 8.70000000e+01, 9.00000000e+01,
1.14000000e+02, 1.11000000e+02, 1.34000000e+02,
1.19000000e+02, 1.43000000e+02, 1.54000000e+02,
1.76000000e+02, 1.64000000e+02, 1.88000000e+02,
2.13000000e+02, 1.99000000e+02, 2.14000000e+02,
2.49000000e+02, 2.46000000e+02],
[ 1.00000000e+01, 4.80000000e+01, 4.80000000e+01,
7.20000000e+01, 8.10000000e+01, 1.06000000e+02,
1.14000000e+02, 1.05000000e+02, 1.32000000e+02,
1.30000000e+02, 1.61000000e+02, 1.61000000e+02,
1.87000000e+02, 1.77000000e+02, 1.91000000e+02,
2.19000000e+02, 2.10000000e+02, 2.34000000e+02,
2.15000000e+02, 2.21000000e+02],
[ 1.40000000e+01, 3.50000000e+01, 4.50000000e+01,
6.30000000e+01, 7.70000000e+01, 8.90000000e+01,
9.10000000e+01, 8.80000000e+01, 1.23000000e+02,
1.48000000e+02, 1.56000000e+02, 1.51000000e+02,
1.54000000e+02, 1.77000000e+02, 1.88000000e+02,
2.09000000e+02, 1.87000000e+02, 2.13000000e+02,
2.28000000e+02, 2.24000000e+02],
[ 1.30000000e+01, 4.00000000e+01, 5.30000000e+01,
6.40000000e+01, 7.50000000e+01, 8.80000000e+01,
1.04000000e+02, 1.31000000e+02, 1.16000000e+02,
1.43000000e+02, 1.39000000e+02, 1.48000000e+02,
1.78000000e+02, 1.69000000e+02, 1.88000000e+02,
1.84000000e+02, 2.02000000e+02, 2.18000000e+02,
2.21000000e+02, 2.04000000e+02],
[ 1.30000000e+01, 2.90000000e+01, 5.10000000e+01,
5.30000000e+01, 6.10000000e+01, 8.70000000e+01,
1.07000000e+02, 1.18000000e+02, 1.20000000e+02,
1.21000000e+02, 1.39000000e+02, 1.30000000e+02,
1.71000000e+02, 1.57000000e+02, 1.86000000e+02,
1.87000000e+02, 2.06000000e+02, 2.12000000e+02,
1.96000000e+02, 2.44000000e+02],
[ 1.00000000e+01, 3.70000000e+01, 5.10000000e+01,
6.80000000e+01, 8.80000000e+01, 7.30000000e+01,
6.90000000e+01, 1.02000000e+02, 1.09000000e+02,
1.25000000e+02, 1.32000000e+02, 1.47000000e+02,
1.93000000e+02, 1.91000000e+02, 1.72000000e+02,
1.88000000e+02, 1.96000000e+02, 2.39000000e+02,
2.39000000e+02, 2.42000000e+02],
[ 1.50000000e+01, 2.50000000e+01, 5.20000000e+01,
5.90000000e+01, 7.30000000e+01, 6.90000000e+01,
8.30000000e+01, 9.80000000e+01, 1.03000000e+02,
1.37000000e+02, 1.78000000e+02, 1.51000000e+02,
1.73000000e+02, 1.65000000e+02, 1.88000000e+02,
2.03000000e+02, 1.97000000e+02, 2.07000000e+02,
2.31000000e+02, 2.41000000e+02],
[ 1.60000000e+01, 2.60000000e+01, 5.40000000e+01,
4.70000000e+01, 6.00000000e+01, 7.30000000e+01,
1.06000000e+02, 1.05000000e+02, 1.14000000e+02,
1.27000000e+02, 1.45000000e+02, 1.49000000e+02,
1.59000000e+02, 1.94000000e+02, 1.68000000e+02,
2.08000000e+02, 2.32000000e+02, 2.26000000e+02,
2.04000000e+02, 2.33000000e+02],
[ 3.00000000e+00, 2.90000000e+01, 4.20000000e+01,
6.20000000e+01, 6.50000000e+01, 7.40000000e+01,
9.00000000e+01, 8.30000000e+01, 1.28000000e+02,
1.19000000e+02, 1.50000000e+02, 1.18000000e+02,
1.71000000e+02, 1.75000000e+02, 1.76000000e+02,
1.60000000e+02, 2.24000000e+02, 2.46000000e+02,
2.28000000e+02, 2.19000000e+02]])
In [43]:
with open('randDR12QRRbin3.pkl','w') as f:
pickle.dump(rr2d,f)
rr2d
Out[43]:
array([[ 1.32956000e+05, 5.40000000e+01, 6.30000000e+01,
8.40000000e+01, 8.20000000e+01, 1.08000000e+02,
1.04000000e+02, 1.04000000e+02, 1.20000000e+02,
1.65000000e+02, 1.63000000e+02, 1.58000000e+02,
1.78000000e+02, 1.63000000e+02, 1.58000000e+02,
1.96000000e+02, 1.93000000e+02, 2.31000000e+02,
2.42000000e+02, 2.71000000e+02],
[ 2.00000000e+01, 6.00000000e+01, 8.10000000e+01,
7.10000000e+01, 1.07000000e+02, 9.70000000e+01,
9.40000000e+01, 1.07000000e+02, 1.54000000e+02,
1.51000000e+02, 1.29000000e+02, 1.75000000e+02,
1.57000000e+02, 1.91000000e+02, 1.85000000e+02,
2.09000000e+02, 2.07000000e+02, 2.22000000e+02,
2.60000000e+02, 2.59000000e+02],
[ 1.70000000e+01, 5.80000000e+01, 6.80000000e+01,
8.10000000e+01, 9.80000000e+01, 1.21000000e+02,
8.90000000e+01, 1.11000000e+02, 1.42000000e+02,
1.24000000e+02, 1.54000000e+02, 1.46000000e+02,
1.93000000e+02, 2.01000000e+02, 1.79000000e+02,
2.01000000e+02, 2.11000000e+02, 2.23000000e+02,
2.36000000e+02, 2.65000000e+02],
[ 2.40000000e+01, 5.90000000e+01, 6.90000000e+01,
7.40000000e+01, 7.90000000e+01, 1.06000000e+02,
1.07000000e+02, 1.16000000e+02, 1.25000000e+02,
1.33000000e+02, 1.60000000e+02, 1.76000000e+02,
1.75000000e+02, 1.37000000e+02, 1.79000000e+02,
1.88000000e+02, 2.04000000e+02, 2.29000000e+02,
2.46000000e+02, 2.59000000e+02],
[ 1.60000000e+01, 5.40000000e+01, 6.00000000e+01,
8.10000000e+01, 8.60000000e+01, 8.10000000e+01,
1.01000000e+02, 1.30000000e+02, 1.15000000e+02,
1.48000000e+02, 1.58000000e+02, 1.58000000e+02,
1.63000000e+02, 1.81000000e+02, 1.85000000e+02,
2.09000000e+02, 2.06000000e+02, 2.24000000e+02,
2.14000000e+02, 2.57000000e+02],
[ 2.70000000e+01, 5.00000000e+01, 6.30000000e+01,
8.00000000e+01, 7.90000000e+01, 9.70000000e+01,
1.11000000e+02, 1.35000000e+02, 1.26000000e+02,
1.58000000e+02, 1.46000000e+02, 1.79000000e+02,
1.69000000e+02, 1.86000000e+02, 1.78000000e+02,
2.01000000e+02, 2.22000000e+02, 2.56000000e+02,
2.29000000e+02, 2.36000000e+02],
[ 2.20000000e+01, 5.70000000e+01, 7.50000000e+01,
7.60000000e+01, 1.00000000e+02, 9.30000000e+01,
1.05000000e+02, 9.50000000e+01, 1.14000000e+02,
1.36000000e+02, 1.42000000e+02, 1.56000000e+02,
1.65000000e+02, 2.02000000e+02, 1.63000000e+02,
1.98000000e+02, 2.25000000e+02, 2.29000000e+02,
2.16000000e+02, 2.67000000e+02],
[ 1.80000000e+01, 3.40000000e+01, 6.20000000e+01,
8.10000000e+01, 9.40000000e+01, 8.90000000e+01,
1.13000000e+02, 1.26000000e+02, 1.22000000e+02,
1.47000000e+02, 1.45000000e+02, 1.52000000e+02,
1.65000000e+02, 1.88000000e+02, 1.94000000e+02,
1.94000000e+02, 2.17000000e+02, 2.17000000e+02,
2.12000000e+02, 2.44000000e+02],
[ 1.40000000e+01, 5.30000000e+01, 5.60000000e+01,
6.00000000e+01, 7.50000000e+01, 9.40000000e+01,
9.40000000e+01, 1.04000000e+02, 1.24000000e+02,
1.48000000e+02, 1.58000000e+02, 1.75000000e+02,
1.68000000e+02, 1.75000000e+02, 1.89000000e+02,
2.26000000e+02, 2.09000000e+02, 2.61000000e+02,
2.15000000e+02, 2.64000000e+02],
[ 2.10000000e+01, 3.80000000e+01, 6.00000000e+01,
7.70000000e+01, 7.50000000e+01, 8.30000000e+01,
1.19000000e+02, 1.10000000e+02, 1.14000000e+02,
1.39000000e+02, 1.53000000e+02, 1.65000000e+02,
1.59000000e+02, 1.81000000e+02, 1.98000000e+02,
2.28000000e+02, 2.09000000e+02, 2.12000000e+02,
2.38000000e+02, 2.34000000e+02],
[ 2.10000000e+01, 2.70000000e+01, 7.00000000e+01,
7.20000000e+01, 8.60000000e+01, 8.10000000e+01,
1.09000000e+02, 1.08000000e+02, 1.19000000e+02,
1.69000000e+02, 1.59000000e+02, 1.57000000e+02,
1.63000000e+02, 1.84000000e+02, 2.00000000e+02,
1.91000000e+02, 2.08000000e+02, 2.03000000e+02,
2.17000000e+02, 2.32000000e+02],
[ 1.80000000e+01, 4.80000000e+01, 4.50000000e+01,
6.50000000e+01, 8.70000000e+01, 9.00000000e+01,
1.14000000e+02, 1.11000000e+02, 1.34000000e+02,
1.19000000e+02, 1.43000000e+02, 1.54000000e+02,
1.76000000e+02, 1.64000000e+02, 1.88000000e+02,
2.13000000e+02, 1.99000000e+02, 2.14000000e+02,
2.49000000e+02, 2.46000000e+02],
[ 1.00000000e+01, 4.80000000e+01, 4.80000000e+01,
7.20000000e+01, 8.10000000e+01, 1.06000000e+02,
1.14000000e+02, 1.05000000e+02, 1.32000000e+02,
1.30000000e+02, 1.61000000e+02, 1.61000000e+02,
1.87000000e+02, 1.77000000e+02, 1.91000000e+02,
2.19000000e+02, 2.10000000e+02, 2.34000000e+02,
2.15000000e+02, 2.21000000e+02],
[ 1.40000000e+01, 3.50000000e+01, 4.50000000e+01,
6.30000000e+01, 7.70000000e+01, 8.90000000e+01,
9.10000000e+01, 8.80000000e+01, 1.23000000e+02,
1.48000000e+02, 1.56000000e+02, 1.51000000e+02,
1.54000000e+02, 1.77000000e+02, 1.88000000e+02,
2.09000000e+02, 1.87000000e+02, 2.13000000e+02,
2.28000000e+02, 2.24000000e+02],
[ 1.30000000e+01, 4.00000000e+01, 5.30000000e+01,
6.40000000e+01, 7.50000000e+01, 8.80000000e+01,
1.04000000e+02, 1.31000000e+02, 1.16000000e+02,
1.43000000e+02, 1.39000000e+02, 1.48000000e+02,
1.78000000e+02, 1.69000000e+02, 1.88000000e+02,
1.84000000e+02, 2.02000000e+02, 2.18000000e+02,
2.21000000e+02, 2.04000000e+02],
[ 1.30000000e+01, 2.90000000e+01, 5.10000000e+01,
5.30000000e+01, 6.10000000e+01, 8.70000000e+01,
1.07000000e+02, 1.18000000e+02, 1.20000000e+02,
1.21000000e+02, 1.39000000e+02, 1.30000000e+02,
1.71000000e+02, 1.57000000e+02, 1.86000000e+02,
1.87000000e+02, 2.06000000e+02, 2.12000000e+02,
1.96000000e+02, 2.44000000e+02],
[ 1.00000000e+01, 3.70000000e+01, 5.10000000e+01,
6.80000000e+01, 8.80000000e+01, 7.30000000e+01,
6.90000000e+01, 1.02000000e+02, 1.09000000e+02,
1.25000000e+02, 1.32000000e+02, 1.47000000e+02,
1.93000000e+02, 1.91000000e+02, 1.72000000e+02,
1.88000000e+02, 1.96000000e+02, 2.39000000e+02,
2.39000000e+02, 2.42000000e+02],
[ 1.50000000e+01, 2.50000000e+01, 5.20000000e+01,
5.90000000e+01, 7.30000000e+01, 6.90000000e+01,
8.30000000e+01, 9.80000000e+01, 1.03000000e+02,
1.37000000e+02, 1.78000000e+02, 1.51000000e+02,
1.73000000e+02, 1.65000000e+02, 1.88000000e+02,
2.03000000e+02, 1.97000000e+02, 2.07000000e+02,
2.31000000e+02, 2.41000000e+02],
[ 1.60000000e+01, 2.60000000e+01, 5.40000000e+01,
4.70000000e+01, 6.00000000e+01, 7.30000000e+01,
1.06000000e+02, 1.05000000e+02, 1.14000000e+02,
1.27000000e+02, 1.45000000e+02, 1.49000000e+02,
1.59000000e+02, 1.94000000e+02, 1.68000000e+02,
2.08000000e+02, 2.32000000e+02, 2.26000000e+02,
2.04000000e+02, 2.33000000e+02],
[ 3.00000000e+00, 2.90000000e+01, 4.20000000e+01,
6.20000000e+01, 6.50000000e+01, 7.40000000e+01,
9.00000000e+01, 8.30000000e+01, 1.28000000e+02,
1.19000000e+02, 1.50000000e+02, 1.18000000e+02,
1.71000000e+02, 1.75000000e+02, 1.76000000e+02,
1.60000000e+02, 2.24000000e+02, 2.46000000e+02,
2.28000000e+02, 2.19000000e+02]])
In [ ]:
In [44]:
# Getting back the objects:
with open('../output/randDR12Qbin4.pkl') as f: # Python 3: open(..., 'rb')
dat = pickle.load(f)
dat
Out[44]:
array([[ 3.088839, 3.509669, 0.526771],
[ 2.906448, 5.746496, 0.024581],
[ 2.879372, 2.298062, 0.193726],
...,
[ 3.166911, 0.098733, 0.085691],
[ 3.064749, 6.228745, 0.393404],
[ 3.113391, 5.734536, 0.043476]])
In [45]:
rr2d=np.zeros((20,20))
In [46]:
len(dat)
Out[46]:
42945
In [47]:
rr2d
Out[47]:
array([[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.],
[ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0.]])
In [48]:
%%time
dist0=d.cdist([dat[0],],dat,APdz)[0]
dist1=d.cdist([dat[0],],dat,APzdth)[0]
print np.histogram2d(dist0, dist1, bins=20, range=rng)[0]
[[ 1. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]
[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.
0. 0.]]
CPU times: user 410 ms, sys: 51.4 ms, total: 461 ms
Wall time: 422 ms
In [49]:
%%time
while(len(dat))>0:
dist0=d.cdist([dat[0],],dat,APdz)[0]
dist1=d.cdist([dat[0],],dat,APzdth)[0]
rr2d+=np.histogram2d(dist0, dist1, bins=20, range=rng)[0]
dat=np.delete(dat,0,axis=0)
if len(dat)%1000==0:
print len(dat)/1000
print rr2d
42
41
40
39
38
37
36
35
34
33
32
31
30
29
28
27
26
25
24
23
22
21
20
19
18
17
16
15
14
13
12
11
10
9
8
7
6
5
4
3
2
1
0
[[ 4.08440000e+04 1.00000000e+00 3.00000000e+00 1.00000000e+00
4.00000000e+00 4.00000000e+00 1.00000000e+00 4.00000000e+00
3.00000000e+00 1.00000000e+01 1.00000000e+01 8.00000000e+00
1.10000000e+01 6.00000000e+00 7.00000000e+00 8.00000000e+00
1.70000000e+01 1.20000000e+01 1.80000000e+01 1.70000000e+01]
[ 0.00000000e+00 2.00000000e+00 0.00000000e+00 3.00000000e+00
2.00000000e+00 5.00000000e+00 2.00000000e+00 6.00000000e+00
9.00000000e+00 7.00000000e+00 8.00000000e+00 1.10000000e+01
8.00000000e+00 1.50000000e+01 8.00000000e+00 8.00000000e+00
1.70000000e+01 8.00000000e+00 8.00000000e+00 1.40000000e+01]
[ 0.00000000e+00 1.00000000e+00 1.00000000e+00 0.00000000e+00
2.00000000e+00 4.00000000e+00 6.00000000e+00 6.00000000e+00
6.00000000e+00 4.00000000e+00 1.20000000e+01 7.00000000e+00
1.30000000e+01 1.10000000e+01 6.00000000e+00 1.40000000e+01
9.00000000e+00 1.30000000e+01 1.00000000e+01 8.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 2.00000000e+00 4.00000000e+00
4.00000000e+00 4.00000000e+00 6.00000000e+00 5.00000000e+00
1.00000000e+01 5.00000000e+00 8.00000000e+00 1.30000000e+01
7.00000000e+00 1.20000000e+01 1.10000000e+01 1.00000000e+01
1.00000000e+01 9.00000000e+00 1.00000000e+01 6.00000000e+00]
[ 0.00000000e+00 5.00000000e+00 1.00000000e+00 0.00000000e+00
5.00000000e+00 3.00000000e+00 6.00000000e+00 2.00000000e+00
3.00000000e+00 4.00000000e+00 4.00000000e+00 8.00000000e+00
9.00000000e+00 7.00000000e+00 8.00000000e+00 1.40000000e+01
7.00000000e+00 7.00000000e+00 1.20000000e+01 1.70000000e+01]
[ 0.00000000e+00 0.00000000e+00 3.00000000e+00 1.00000000e+00
4.00000000e+00 5.00000000e+00 3.00000000e+00 6.00000000e+00
9.00000000e+00 5.00000000e+00 3.00000000e+00 8.00000000e+00
6.00000000e+00 1.10000000e+01 6.00000000e+00 1.10000000e+01
1.40000000e+01 1.40000000e+01 1.10000000e+01 9.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 2.00000000e+00 2.00000000e+00
1.00000000e+00 2.00000000e+00 3.00000000e+00 3.00000000e+00
4.00000000e+00 4.00000000e+00 3.00000000e+00 6.00000000e+00
9.00000000e+00 6.00000000e+00 1.40000000e+01 1.30000000e+01
6.00000000e+00 7.00000000e+00 1.40000000e+01 1.60000000e+01]
[ 1.00000000e+00 1.00000000e+00 5.00000000e+00 2.00000000e+00
3.00000000e+00 4.00000000e+00 5.00000000e+00 6.00000000e+00
8.00000000e+00 2.00000000e+00 1.10000000e+01 9.00000000e+00
9.00000000e+00 1.50000000e+01 8.00000000e+00 1.00000000e+01
1.20000000e+01 1.20000000e+01 1.30000000e+01 7.00000000e+00]
[ 1.00000000e+00 2.00000000e+00 2.00000000e+00 1.00000000e+00
3.00000000e+00 2.00000000e+00 3.00000000e+00 6.00000000e+00
1.20000000e+01 9.00000000e+00 7.00000000e+00 5.00000000e+00
1.00000000e+01 1.50000000e+01 1.50000000e+01 1.40000000e+01
1.00000000e+01 1.50000000e+01 1.40000000e+01 1.10000000e+01]
[ 0.00000000e+00 1.00000000e+00 2.00000000e+00 5.00000000e+00
6.00000000e+00 7.00000000e+00 0.00000000e+00 5.00000000e+00
4.00000000e+00 1.20000000e+01 3.00000000e+00 9.00000000e+00
1.10000000e+01 1.10000000e+01 1.00000000e+01 5.00000000e+00
1.30000000e+01 1.20000000e+01 9.00000000e+00 1.70000000e+01]
[ 0.00000000e+00 3.00000000e+00 1.00000000e+00 2.00000000e+00
2.00000000e+00 4.00000000e+00 5.00000000e+00 3.00000000e+00
1.50000000e+01 7.00000000e+00 8.00000000e+00 8.00000000e+00
8.00000000e+00 8.00000000e+00 1.00000000e+01 1.00000000e+01
1.10000000e+01 1.00000000e+01 8.00000000e+00 1.00000000e+01]
[ 0.00000000e+00 0.00000000e+00 0.00000000e+00 2.00000000e+00
2.00000000e+00 7.00000000e+00 5.00000000e+00 2.00000000e+00
9.00000000e+00 6.00000000e+00 8.00000000e+00 4.00000000e+00
8.00000000e+00 9.00000000e+00 7.00000000e+00 9.00000000e+00
1.00000000e+01 8.00000000e+00 1.70000000e+01 6.00000000e+00]
[ 1.00000000e+00 2.00000000e+00 3.00000000e+00 1.00000000e+00
2.00000000e+00 4.00000000e+00 3.00000000e+00 5.00000000e+00
6.00000000e+00 6.00000000e+00 9.00000000e+00 8.00000000e+00
1.10000000e+01 3.00000000e+00 1.20000000e+01 1.50000000e+01
1.00000000e+01 1.20000000e+01 1.00000000e+01 1.40000000e+01]
[ 0.00000000e+00 4.00000000e+00 0.00000000e+00 2.00000000e+00
1.00000000e+00 3.00000000e+00 7.00000000e+00 6.00000000e+00
1.00000000e+01 6.00000000e+00 4.00000000e+00 1.40000000e+01
2.00000000e+00 1.00000000e+01 1.00000000e+01 5.00000000e+00
1.00000000e+01 1.40000000e+01 1.10000000e+01 1.30000000e+01]
[ 0.00000000e+00 3.00000000e+00 3.00000000e+00 3.00000000e+00
4.00000000e+00 3.00000000e+00 3.00000000e+00 6.00000000e+00
3.00000000e+00 4.00000000e+00 7.00000000e+00 1.20000000e+01
1.10000000e+01 1.10000000e+01 1.00000000e+01 1.30000000e+01
1.30000000e+01 9.00000000e+00 1.10000000e+01 7.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 0.00000000e+00 1.00000000e+00
4.00000000e+00 3.00000000e+00 2.00000000e+00 6.00000000e+00
1.00000000e+00 6.00000000e+00 1.00000000e+01 8.00000000e+00
8.00000000e+00 5.00000000e+00 8.00000000e+00 1.40000000e+01
1.30000000e+01 1.30000000e+01 8.00000000e+00 8.00000000e+00]
[ 0.00000000e+00 1.00000000e+00 2.00000000e+00 1.00000000e+00
4.00000000e+00 3.00000000e+00 3.00000000e+00 4.00000000e+00
1.00000000e+01 1.10000000e+01 5.00000000e+00 4.00000000e+00
1.10000000e+01 1.00000000e+01 7.00000000e+00 1.20000000e+01
8.00000000e+00 1.00000000e+01 1.40000000e+01 1.50000000e+01]
[ 2.00000000e+00 1.00000000e+00 2.00000000e+00 3.00000000e+00
5.00000000e+00 6.00000000e+00 2.00000000e+00 2.00000000e+00
6.00000000e+00 5.00000000e+00 8.00000000e+00 6.00000000e+00
1.50000000e+01 5.00000000e+00 1.00000000e+01 1.00000000e+01
1.30000000e+01 1.60000000e+01 1.40000000e+01 1.50000000e+01]
[ 0.00000000e+00 1.00000000e+00 2.00000000e+00 4.00000000e+00
0.00000000e+00 2.00000000e+00 3.00000000e+00 4.00000000e+00
2.00000000e+00 3.00000000e+00 4.00000000e+00 6.00000000e+00
8.00000000e+00 1.30000000e+01 6.00000000e+00 1.10000000e+01
5.00000000e+00 1.30000000e+01 1.00000000e+01 1.30000000e+01]
[ 1.00000000e+00 0.00000000e+00 3.00000000e+00 1.00000000e+00
4.00000000e+00 5.00000000e+00 7.00000000e+00 8.00000000e+00
6.00000000e+00 8.00000000e+00 7.00000000e+00 6.00000000e+00
1.30000000e+01 7.00000000e+00 1.50000000e+01 1.30000000e+01
6.00000000e+00 1.20000000e+01 1.40000000e+01 1.20000000e+01]]
CPU times: user 2h 4min 44s, sys: 53.6 s, total: 2h 5min 38s
Wall time: 2h 5min 12s
In [50]:
rr2d
Out[50]:
array([[ 4.08440000e+04, 1.00000000e+00, 3.00000000e+00,
1.00000000e+00, 4.00000000e+00, 4.00000000e+00,
1.00000000e+00, 4.00000000e+00, 3.00000000e+00,
1.00000000e+01, 1.00000000e+01, 8.00000000e+00,
1.10000000e+01, 6.00000000e+00, 7.00000000e+00,
8.00000000e+00, 1.70000000e+01, 1.20000000e+01,
1.80000000e+01, 1.70000000e+01],
[ 0.00000000e+00, 2.00000000e+00, 0.00000000e+00,
3.00000000e+00, 2.00000000e+00, 5.00000000e+00,
2.00000000e+00, 6.00000000e+00, 9.00000000e+00,
7.00000000e+00, 8.00000000e+00, 1.10000000e+01,
8.00000000e+00, 1.50000000e+01, 8.00000000e+00,
8.00000000e+00, 1.70000000e+01, 8.00000000e+00,
8.00000000e+00, 1.40000000e+01],
[ 0.00000000e+00, 1.00000000e+00, 1.00000000e+00,
0.00000000e+00, 2.00000000e+00, 4.00000000e+00,
6.00000000e+00, 6.00000000e+00, 6.00000000e+00,
4.00000000e+00, 1.20000000e+01, 7.00000000e+00,
1.30000000e+01, 1.10000000e+01, 6.00000000e+00,
1.40000000e+01, 9.00000000e+00, 1.30000000e+01,
1.00000000e+01, 8.00000000e+00],
[ 0.00000000e+00, 0.00000000e+00, 2.00000000e+00,
4.00000000e+00, 4.00000000e+00, 4.00000000e+00,
6.00000000e+00, 5.00000000e+00, 1.00000000e+01,
5.00000000e+00, 8.00000000e+00, 1.30000000e+01,
7.00000000e+00, 1.20000000e+01, 1.10000000e+01,
1.00000000e+01, 1.00000000e+01, 9.00000000e+00,
1.00000000e+01, 6.00000000e+00],
[ 0.00000000e+00, 5.00000000e+00, 1.00000000e+00,
0.00000000e+00, 5.00000000e+00, 3.00000000e+00,
6.00000000e+00, 2.00000000e+00, 3.00000000e+00,
4.00000000e+00, 4.00000000e+00, 8.00000000e+00,
9.00000000e+00, 7.00000000e+00, 8.00000000e+00,
1.40000000e+01, 7.00000000e+00, 7.00000000e+00,
1.20000000e+01, 1.70000000e+01],
[ 0.00000000e+00, 0.00000000e+00, 3.00000000e+00,
1.00000000e+00, 4.00000000e+00, 5.00000000e+00,
3.00000000e+00, 6.00000000e+00, 9.00000000e+00,
5.00000000e+00, 3.00000000e+00, 8.00000000e+00,
6.00000000e+00, 1.10000000e+01, 6.00000000e+00,
1.10000000e+01, 1.40000000e+01, 1.40000000e+01,
1.10000000e+01, 9.00000000e+00],
[ 0.00000000e+00, 0.00000000e+00, 2.00000000e+00,
2.00000000e+00, 1.00000000e+00, 2.00000000e+00,
3.00000000e+00, 3.00000000e+00, 4.00000000e+00,
4.00000000e+00, 3.00000000e+00, 6.00000000e+00,
9.00000000e+00, 6.00000000e+00, 1.40000000e+01,
1.30000000e+01, 6.00000000e+00, 7.00000000e+00,
1.40000000e+01, 1.60000000e+01],
[ 1.00000000e+00, 1.00000000e+00, 5.00000000e+00,
2.00000000e+00, 3.00000000e+00, 4.00000000e+00,
5.00000000e+00, 6.00000000e+00, 8.00000000e+00,
2.00000000e+00, 1.10000000e+01, 9.00000000e+00,
9.00000000e+00, 1.50000000e+01, 8.00000000e+00,
1.00000000e+01, 1.20000000e+01, 1.20000000e+01,
1.30000000e+01, 7.00000000e+00],
[ 1.00000000e+00, 2.00000000e+00, 2.00000000e+00,
1.00000000e+00, 3.00000000e+00, 2.00000000e+00,
3.00000000e+00, 6.00000000e+00, 1.20000000e+01,
9.00000000e+00, 7.00000000e+00, 5.00000000e+00,
1.00000000e+01, 1.50000000e+01, 1.50000000e+01,
1.40000000e+01, 1.00000000e+01, 1.50000000e+01,
1.40000000e+01, 1.10000000e+01],
[ 0.00000000e+00, 1.00000000e+00, 2.00000000e+00,
5.00000000e+00, 6.00000000e+00, 7.00000000e+00,
0.00000000e+00, 5.00000000e+00, 4.00000000e+00,
1.20000000e+01, 3.00000000e+00, 9.00000000e+00,
1.10000000e+01, 1.10000000e+01, 1.00000000e+01,
5.00000000e+00, 1.30000000e+01, 1.20000000e+01,
9.00000000e+00, 1.70000000e+01],
[ 0.00000000e+00, 3.00000000e+00, 1.00000000e+00,
2.00000000e+00, 2.00000000e+00, 4.00000000e+00,
5.00000000e+00, 3.00000000e+00, 1.50000000e+01,
7.00000000e+00, 8.00000000e+00, 8.00000000e+00,
8.00000000e+00, 8.00000000e+00, 1.00000000e+01,
1.00000000e+01, 1.10000000e+01, 1.00000000e+01,
8.00000000e+00, 1.00000000e+01],
[ 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,
2.00000000e+00, 2.00000000e+00, 7.00000000e+00,
5.00000000e+00, 2.00000000e+00, 9.00000000e+00,
6.00000000e+00, 8.00000000e+00, 4.00000000e+00,
8.00000000e+00, 9.00000000e+00, 7.00000000e+00,
9.00000000e+00, 1.00000000e+01, 8.00000000e+00,
1.70000000e+01, 6.00000000e+00],
[ 1.00000000e+00, 2.00000000e+00, 3.00000000e+00,
1.00000000e+00, 2.00000000e+00, 4.00000000e+00,
3.00000000e+00, 5.00000000e+00, 6.00000000e+00,
6.00000000e+00, 9.00000000e+00, 8.00000000e+00,
1.10000000e+01, 3.00000000e+00, 1.20000000e+01,
1.50000000e+01, 1.00000000e+01, 1.20000000e+01,
1.00000000e+01, 1.40000000e+01],
[ 0.00000000e+00, 4.00000000e+00, 0.00000000e+00,
2.00000000e+00, 1.00000000e+00, 3.00000000e+00,
7.00000000e+00, 6.00000000e+00, 1.00000000e+01,
6.00000000e+00, 4.00000000e+00, 1.40000000e+01,
2.00000000e+00, 1.00000000e+01, 1.00000000e+01,
5.00000000e+00, 1.00000000e+01, 1.40000000e+01,
1.10000000e+01, 1.30000000e+01],
[ 0.00000000e+00, 3.00000000e+00, 3.00000000e+00,
3.00000000e+00, 4.00000000e+00, 3.00000000e+00,
3.00000000e+00, 6.00000000e+00, 3.00000000e+00,
4.00000000e+00, 7.00000000e+00, 1.20000000e+01,
1.10000000e+01, 1.10000000e+01, 1.00000000e+01,
1.30000000e+01, 1.30000000e+01, 9.00000000e+00,
1.10000000e+01, 7.00000000e+00],
[ 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,
1.00000000e+00, 4.00000000e+00, 3.00000000e+00,
2.00000000e+00, 6.00000000e+00, 1.00000000e+00,
6.00000000e+00, 1.00000000e+01, 8.00000000e+00,
8.00000000e+00, 5.00000000e+00, 8.00000000e+00,
1.40000000e+01, 1.30000000e+01, 1.30000000e+01,
8.00000000e+00, 8.00000000e+00],
[ 0.00000000e+00, 1.00000000e+00, 2.00000000e+00,
1.00000000e+00, 4.00000000e+00, 3.00000000e+00,
3.00000000e+00, 4.00000000e+00, 1.00000000e+01,
1.10000000e+01, 5.00000000e+00, 4.00000000e+00,
1.10000000e+01, 1.00000000e+01, 7.00000000e+00,
1.20000000e+01, 8.00000000e+00, 1.00000000e+01,
1.40000000e+01, 1.50000000e+01],
[ 2.00000000e+00, 1.00000000e+00, 2.00000000e+00,
3.00000000e+00, 5.00000000e+00, 6.00000000e+00,
2.00000000e+00, 2.00000000e+00, 6.00000000e+00,
5.00000000e+00, 8.00000000e+00, 6.00000000e+00,
1.50000000e+01, 5.00000000e+00, 1.00000000e+01,
1.00000000e+01, 1.30000000e+01, 1.60000000e+01,
1.40000000e+01, 1.50000000e+01],
[ 0.00000000e+00, 1.00000000e+00, 2.00000000e+00,
4.00000000e+00, 0.00000000e+00, 2.00000000e+00,
3.00000000e+00, 4.00000000e+00, 2.00000000e+00,
3.00000000e+00, 4.00000000e+00, 6.00000000e+00,
8.00000000e+00, 1.30000000e+01, 6.00000000e+00,
1.10000000e+01, 5.00000000e+00, 1.30000000e+01,
1.00000000e+01, 1.30000000e+01],
[ 1.00000000e+00, 0.00000000e+00, 3.00000000e+00,
1.00000000e+00, 4.00000000e+00, 5.00000000e+00,
7.00000000e+00, 8.00000000e+00, 6.00000000e+00,
8.00000000e+00, 7.00000000e+00, 6.00000000e+00,
1.30000000e+01, 7.00000000e+00, 1.50000000e+01,
1.30000000e+01, 6.00000000e+00, 1.20000000e+01,
1.40000000e+01, 1.20000000e+01]])
In [51]:
with open('randDR12QRRbin4.pkl','w') as f:
pickle.dump(rr2d,f)
rr2d
Out[51]:
array([[ 4.08440000e+04, 1.00000000e+00, 3.00000000e+00,
1.00000000e+00, 4.00000000e+00, 4.00000000e+00,
1.00000000e+00, 4.00000000e+00, 3.00000000e+00,
1.00000000e+01, 1.00000000e+01, 8.00000000e+00,
1.10000000e+01, 6.00000000e+00, 7.00000000e+00,
8.00000000e+00, 1.70000000e+01, 1.20000000e+01,
1.80000000e+01, 1.70000000e+01],
[ 0.00000000e+00, 2.00000000e+00, 0.00000000e+00,
3.00000000e+00, 2.00000000e+00, 5.00000000e+00,
2.00000000e+00, 6.00000000e+00, 9.00000000e+00,
7.00000000e+00, 8.00000000e+00, 1.10000000e+01,
8.00000000e+00, 1.50000000e+01, 8.00000000e+00,
8.00000000e+00, 1.70000000e+01, 8.00000000e+00,
8.00000000e+00, 1.40000000e+01],
[ 0.00000000e+00, 1.00000000e+00, 1.00000000e+00,
0.00000000e+00, 2.00000000e+00, 4.00000000e+00,
6.00000000e+00, 6.00000000e+00, 6.00000000e+00,
4.00000000e+00, 1.20000000e+01, 7.00000000e+00,
1.30000000e+01, 1.10000000e+01, 6.00000000e+00,
1.40000000e+01, 9.00000000e+00, 1.30000000e+01,
1.00000000e+01, 8.00000000e+00],
[ 0.00000000e+00, 0.00000000e+00, 2.00000000e+00,
4.00000000e+00, 4.00000000e+00, 4.00000000e+00,
6.00000000e+00, 5.00000000e+00, 1.00000000e+01,
5.00000000e+00, 8.00000000e+00, 1.30000000e+01,
7.00000000e+00, 1.20000000e+01, 1.10000000e+01,
1.00000000e+01, 1.00000000e+01, 9.00000000e+00,
1.00000000e+01, 6.00000000e+00],
[ 0.00000000e+00, 5.00000000e+00, 1.00000000e+00,
0.00000000e+00, 5.00000000e+00, 3.00000000e+00,
6.00000000e+00, 2.00000000e+00, 3.00000000e+00,
4.00000000e+00, 4.00000000e+00, 8.00000000e+00,
9.00000000e+00, 7.00000000e+00, 8.00000000e+00,
1.40000000e+01, 7.00000000e+00, 7.00000000e+00,
1.20000000e+01, 1.70000000e+01],
[ 0.00000000e+00, 0.00000000e+00, 3.00000000e+00,
1.00000000e+00, 4.00000000e+00, 5.00000000e+00,
3.00000000e+00, 6.00000000e+00, 9.00000000e+00,
5.00000000e+00, 3.00000000e+00, 8.00000000e+00,
6.00000000e+00, 1.10000000e+01, 6.00000000e+00,
1.10000000e+01, 1.40000000e+01, 1.40000000e+01,
1.10000000e+01, 9.00000000e+00],
[ 0.00000000e+00, 0.00000000e+00, 2.00000000e+00,
2.00000000e+00, 1.00000000e+00, 2.00000000e+00,
3.00000000e+00, 3.00000000e+00, 4.00000000e+00,
4.00000000e+00, 3.00000000e+00, 6.00000000e+00,
9.00000000e+00, 6.00000000e+00, 1.40000000e+01,
1.30000000e+01, 6.00000000e+00, 7.00000000e+00,
1.40000000e+01, 1.60000000e+01],
[ 1.00000000e+00, 1.00000000e+00, 5.00000000e+00,
2.00000000e+00, 3.00000000e+00, 4.00000000e+00,
5.00000000e+00, 6.00000000e+00, 8.00000000e+00,
2.00000000e+00, 1.10000000e+01, 9.00000000e+00,
9.00000000e+00, 1.50000000e+01, 8.00000000e+00,
1.00000000e+01, 1.20000000e+01, 1.20000000e+01,
1.30000000e+01, 7.00000000e+00],
[ 1.00000000e+00, 2.00000000e+00, 2.00000000e+00,
1.00000000e+00, 3.00000000e+00, 2.00000000e+00,
3.00000000e+00, 6.00000000e+00, 1.20000000e+01,
9.00000000e+00, 7.00000000e+00, 5.00000000e+00,
1.00000000e+01, 1.50000000e+01, 1.50000000e+01,
1.40000000e+01, 1.00000000e+01, 1.50000000e+01,
1.40000000e+01, 1.10000000e+01],
[ 0.00000000e+00, 1.00000000e+00, 2.00000000e+00,
5.00000000e+00, 6.00000000e+00, 7.00000000e+00,
0.00000000e+00, 5.00000000e+00, 4.00000000e+00,
1.20000000e+01, 3.00000000e+00, 9.00000000e+00,
1.10000000e+01, 1.10000000e+01, 1.00000000e+01,
5.00000000e+00, 1.30000000e+01, 1.20000000e+01,
9.00000000e+00, 1.70000000e+01],
[ 0.00000000e+00, 3.00000000e+00, 1.00000000e+00,
2.00000000e+00, 2.00000000e+00, 4.00000000e+00,
5.00000000e+00, 3.00000000e+00, 1.50000000e+01,
7.00000000e+00, 8.00000000e+00, 8.00000000e+00,
8.00000000e+00, 8.00000000e+00, 1.00000000e+01,
1.00000000e+01, 1.10000000e+01, 1.00000000e+01,
8.00000000e+00, 1.00000000e+01],
[ 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,
2.00000000e+00, 2.00000000e+00, 7.00000000e+00,
5.00000000e+00, 2.00000000e+00, 9.00000000e+00,
6.00000000e+00, 8.00000000e+00, 4.00000000e+00,
8.00000000e+00, 9.00000000e+00, 7.00000000e+00,
9.00000000e+00, 1.00000000e+01, 8.00000000e+00,
1.70000000e+01, 6.00000000e+00],
[ 1.00000000e+00, 2.00000000e+00, 3.00000000e+00,
1.00000000e+00, 2.00000000e+00, 4.00000000e+00,
3.00000000e+00, 5.00000000e+00, 6.00000000e+00,
6.00000000e+00, 9.00000000e+00, 8.00000000e+00,
1.10000000e+01, 3.00000000e+00, 1.20000000e+01,
1.50000000e+01, 1.00000000e+01, 1.20000000e+01,
1.00000000e+01, 1.40000000e+01],
[ 0.00000000e+00, 4.00000000e+00, 0.00000000e+00,
2.00000000e+00, 1.00000000e+00, 3.00000000e+00,
7.00000000e+00, 6.00000000e+00, 1.00000000e+01,
6.00000000e+00, 4.00000000e+00, 1.40000000e+01,
2.00000000e+00, 1.00000000e+01, 1.00000000e+01,
5.00000000e+00, 1.00000000e+01, 1.40000000e+01,
1.10000000e+01, 1.30000000e+01],
[ 0.00000000e+00, 3.00000000e+00, 3.00000000e+00,
3.00000000e+00, 4.00000000e+00, 3.00000000e+00,
3.00000000e+00, 6.00000000e+00, 3.00000000e+00,
4.00000000e+00, 7.00000000e+00, 1.20000000e+01,
1.10000000e+01, 1.10000000e+01, 1.00000000e+01,
1.30000000e+01, 1.30000000e+01, 9.00000000e+00,
1.10000000e+01, 7.00000000e+00],
[ 0.00000000e+00, 0.00000000e+00, 0.00000000e+00,
1.00000000e+00, 4.00000000e+00, 3.00000000e+00,
2.00000000e+00, 6.00000000e+00, 1.00000000e+00,
6.00000000e+00, 1.00000000e+01, 8.00000000e+00,
8.00000000e+00, 5.00000000e+00, 8.00000000e+00,
1.40000000e+01, 1.30000000e+01, 1.30000000e+01,
8.00000000e+00, 8.00000000e+00],
[ 0.00000000e+00, 1.00000000e+00, 2.00000000e+00,
1.00000000e+00, 4.00000000e+00, 3.00000000e+00,
3.00000000e+00, 4.00000000e+00, 1.00000000e+01,
1.10000000e+01, 5.00000000e+00, 4.00000000e+00,
1.10000000e+01, 1.00000000e+01, 7.00000000e+00,
1.20000000e+01, 8.00000000e+00, 1.00000000e+01,
1.40000000e+01, 1.50000000e+01],
[ 2.00000000e+00, 1.00000000e+00, 2.00000000e+00,
3.00000000e+00, 5.00000000e+00, 6.00000000e+00,
2.00000000e+00, 2.00000000e+00, 6.00000000e+00,
5.00000000e+00, 8.00000000e+00, 6.00000000e+00,
1.50000000e+01, 5.00000000e+00, 1.00000000e+01,
1.00000000e+01, 1.30000000e+01, 1.60000000e+01,
1.40000000e+01, 1.50000000e+01],
[ 0.00000000e+00, 1.00000000e+00, 2.00000000e+00,
4.00000000e+00, 0.00000000e+00, 2.00000000e+00,
3.00000000e+00, 4.00000000e+00, 2.00000000e+00,
3.00000000e+00, 4.00000000e+00, 6.00000000e+00,
8.00000000e+00, 1.30000000e+01, 6.00000000e+00,
1.10000000e+01, 5.00000000e+00, 1.30000000e+01,
1.00000000e+01, 1.30000000e+01],
[ 1.00000000e+00, 0.00000000e+00, 3.00000000e+00,
1.00000000e+00, 4.00000000e+00, 5.00000000e+00,
7.00000000e+00, 8.00000000e+00, 6.00000000e+00,
8.00000000e+00, 7.00000000e+00, 6.00000000e+00,
1.30000000e+01, 7.00000000e+00, 1.50000000e+01,
1.30000000e+01, 6.00000000e+00, 1.20000000e+01,
1.40000000e+01, 1.20000000e+01]])
In [ ]:
Content source: rohinkumar/galsurveystudy
Similar notebooks: