python - Count occurences of elements of a matrix fast -


let m , n d x d- , d-dimensonal numpy arrays of integers, respectively. want count number of triples of form (n(i), n(j), m(i,j)). result want numpy array such each entry counts number of occurences of such triple.

edit: m symmetric , don't want count triples i=j.

i'm using itertools.product (for loop on pairs) , numpy.bincount this, slow. there smarter way doing this, using numpy?

since arrays contain integers, can consider each triplet linearly index-able element. here's approach philosophy in mind , avoids loops, -

# form n(i) x n(j) array , append "columnar" m(i,j) array nn_arr = n[np.array(list(itertools.product(range(5), repeat=2)))]  nn_m_arr = np.concatenate((nn_arr,m.reshape(-1,1)),axis=1)  # linear indices version dims = nn_m_arr.max(0)+1 lidx = np.ravel_multi_index(nn_m_arr.t,dims)  # unique indices , counts _, idx, counts = np.unique(lidx,return_index=true,return_counts=true)  # corresponding unique triplets using unique indices , zip counts out = zip(map(tuple,nn_m_arr[idx]),counts)

sample run -

in [206]: m out[206]:  array([[1, 0, 0, 2, 0],        [1, 1, 2, 0, 2],        [0, 0, 2, 0, 1],        [2, 1, 2, 0, 2],        [1, 1, 1, 1, 0]])  in [207]: n out[207]: array([0, 1, 1, 1, 2])  in [208]: out out[208]:  [((0, 0, 1), 1),  ((0, 1, 0), 2),  ((0, 1, 2), 1),  ((0, 2, 0), 1),  ((1, 0, 0), 1),  ((1, 0, 1), 1),  ((1, 0, 2), 1),  ((1, 1, 0), 4),  ((1, 1, 1), 2),  ((1, 1, 2), 3),  ((1, 2, 1), 1),  ((1, 2, 2), 2),  ((2, 0, 1), 1),  ((2, 1, 1), 3),  ((2, 2, 0), 1)]

Comments

Post a Comment

Popular posts from this blog

routing - AngularJS State management ->load multiple states in one page -

Image
i developing site in angularjs ui-router. have situation on button click need change , load 3 state page , 3 looks image .please me how can load multiple pages. in 1 page you need configure router config this: $stateprovider .state("state1", { url: "/url1", controller: "firstcontroller", views: { view1: { templateurl: "templates/template1.html" }, view2: { templateurl: "templates/template2.html" }, view3: { templateurl: "templates/template3.html" } } }).state("state2", { url: "/url2", controller: "secondcontroller", views: { view1: { templateurl: "templates/template4.html" }, view2: { templateurl: "templates/template5.html" }, view3: { templateurl: "templates/template6.html" } }).state("state3", { ...
Read more

php - mysql multiple keyword search in any order -

this question exact duplicate of: php keyword search engine no results [duplicate] i have simple mysql database keyword search functional. results search not being returned if keywords not in same order entered in database. example searching "dog cat lion" return result, searching "dog lion cat" return no results. on how fix issue appreciated. here code. <!doctype html> <html> <head> <meta charset="utf-8"> <title>part search</title> </head> <body> <?php $username = "xxxx"; $password = "xxxx"; $hostname = "localhost"; //connection database mysql_connect($hostname, $username, $password); mysql_select_db("wesco"); $search = mysql_real_escape_string(trim($_post['searchterm'])); $find_parts = mysql_query("select * `parts` `keywords` '%$search%'"); while($row = mysql_fetch_as...
Read more

python - GRASS parser() error -

i'm trying use grass on python 2.7, i'm having problems @ setting script on idle, i'm getting error @ parser() function: here script: import os import sys gisbase = os.environ['gisbase'] = 'c:\program files (x86)\grass gis 7.0.1rc1' gisrc = 'c:\grassdata' gisdbase = 'c:\grassdata' location = 'newlocation' mapset = 'tc' ld_library_path = 'c:\program files (x86)\grass gis 7.0.1rc1\lib' path = 'c:\program files (x86)\grass gis 7.0.1rc1\etc';'c:\program files (x86)\grass gis 7.0.1rc1\etc\python';'c:\program files (x86)\grass gis 7.0.1rc1\lib';'c:\program files (x86)\grass gis 7.0.1rc1\bin';'c:\python27';'c:\program files (x86)\grass gis 7.0.1rc1\python27';'c:\program files (x86)\grass gis 7.0.1rc1\msys' pythonlib = 'c:\python27' pythonpath = 'c:\program files (x86)\grass gis 7.0.1rc1\etc\python' sys.path.append(os.path.join(os.environ['...
Read more