More efficient way of suming over a matrix in numpy/regular python -

-

good day, studies created cellular automato in 2 dimension. program running still try optimize it. piece of code bellow sums 8 neighbor cells of central cell in 2d array. after next cell gets defined function of sum. there faster way 2 for-loops?

before had 4 loops summation, 2 time slower it's now... 

n = len(mat1) m = len(mat1[0]) mat2 = np.zeros((n,m)) sumn = 0 start = time.time() in range(1,n-1):       j in range(1,m-1):         sumn = mat1[i-1,j-1] + mat1[i-1,j] + mat1[i-1,j+1] + mat1[i,j-1] +mat1[i,j+1] + mat1[i+1,j] + mat1[i+1,j+1]+ mat1[i+1,j-1]         if str(sumn) in b , mat1[i,j] == 0:             mat2[i,j] = 1         elif str(sumn) in s , mat1[i,j] == 1:             mat2[i,j] = 1         sumn = 0 end = time.time() print end - start

thanks xnx, included roll on matrix instead of looping on elements. after created boolean 2d numpy array use initialize next generation.

sumn = sum(np.roll(np.roll(mat1, i, 0), j, 1) in (-1, 0, 1) j in (-1, 0, 1)   if (i != 0 or j != 0)).flatten()  mat1 = mat1.flatten()  b = np.array(map(lambda x,l: ((int(x) == 0) , (str(int(l)) in b))     or ((int(x) == 1) , (str(int(l)) in s)), mat1,sumn)).reshape(n,m)                            mat2 = np.zeros((n,m))                                                       mat2[b] = 1                                                          mat2 = mat2.reshape(n,m)

a nice way given in this blog article:

nbrs_count = sum(np.roll(np.roll(mat1, i, 0), j, 1)                  in (-1, 0, 1) j in (-1, 0, 1)                  if (i != 0 or j != 0))

it works because numpy.roll(a, shift, axis) shifts elements along given axis, axis, specified number of places, shift wrapping, each of 8 cells can visited in turn if iterate on rows i=-1,0,1 , columns j=-1,0,1 take care not count centre cell itself.


Comments

Post a Comment

Popular posts from this blog

Hatching array of circles in AutoCAD using c# -

-
i beginner in autocad plugins .i have created array of circle using following code. now, have hatch circles in array .how that??how hatch circles @ once? moreover, how can return hatch method , use hatch objects.sorry language.thanks in advance. [assembly: commandclass(typeof(beamsection.class1))] namespace beamsection { public class class1 { //class1 obj = new class1(); [commandmethod("beamatsupport")] public void addlightweightpolyline() { // current document , database document acdoc = application.documentmanager.mdiactivedocument; database accurdb = acdoc.database; promptpointresult ppr = acdoc.editor.getpoint("\nselect starting point "); var ucs = acdoc.editor.currentusercoordinatesystem; point3d startingpt = ppr.value.transformby(ucs); list<double> radius = new list<double>(); radius.add(0.5);
Read more

ios - UITEXTFIELD InputView Uipicker not working in swift -

-
hello trying pop uipickerview view programmatically when user clicks on textfield. have tried doesn't doing anything. isn't working. nothing happening when click textfield class userprofiletableviewcontroller: uitableviewcontroller,uitextfielddelegate,uipickerviewdelegate,uipickerviewdatasource { var itempicker: uipickerview! = uipickerview() @iboutlet weak var gendertxtfield: uitextfield! var gender = ["male","female"] override func viewdidload() { super.viewdidload() gendertxtfield.delegate = self itempicker!.delegate = self itempicker!.datasource = self itempicker!.backgroundcolor = uicolor.blackcolor() self.gendertxtfield.inputview = itempicker } func numberofcomponentsinpickerview(pickerview: uipickerview) -> int{ return 1 } // returns # of rows in each component.. func pickerview(pickerview: uipickerview, numberofrowsincomponent component:
Read more

Python Pig Latin Translator -

-
so have been working on @ codecademy.com while trying learn python. don't know if blind , can't see why won't work or what. here's code won't let me use. function print variable if had typed in , print "empty" if didn't type anything. print "welcome english pig latin translator!" original = raw_input("what's name?") if len(original) > 0 print original else print "empty" if len(original) > 0 else the above lines contains syntax errors. try: if len(original) > 0: else: note added colons, ( : ).
Read more