python - numpy ravel versus flat in slice assignment -


according documentation, ndarray.flat iterator on array while ndarray.ravel returns flattened array (when possible). question is, when should use 1 or other? 1 preferred rvalue in assignment 1 in code below?

import numpy np  x = np.arange(2).reshape((2,1,1)) y = np.arange(3).reshape((1,3,1)) z = np.arange(5).reshape((1,1,5))  mask = np.random.choice([true, false], size=(2,3,5)) # netcdf4 module wants kind of boolean indexing: nc4slice = tuple(mask.any(axis=axis) axis in ((1,2),(2,0),(0,1))) indices = np.ix_(*nc4slice)  ncrds = 3 npnts = (np.broadcast(*indices)).size points = np.empty((npnts, ncrds)) i,crd in enumerate(np.broadcast_arrays(x,y,z)):     # should use ndarray.flat ...     points[:,i] = crd[indices].flat     # ... or ndarray.ravel():     points[:,i] = crd[indices].ravel()

you don't need either. crd[mask] 1-d. if did, numpy calls np.asarray(rhs) first, same if no copy needed ravel. when copy needed, guess ravel may faster (i did not time it).

if knew copy might needed, , here know nothing needed, reshaping points fastest. since don't need fastest, more matter of taste, , use ravel.


Comments

Post a Comment

Popular posts from this blog

iphone - Three second countdown in cocos2d -

i trying have 3 second countdown display when level beat in game. did: -(void) gamesegmentbeat { [self pauseschedulerandactions]; // overlay game opaque background ccsprite *opaquebg = [ccsprite spritewithfile:@"background1.png"]; opaquebg.position = screencenter; [self addchild:opaquebg z:10000]; // these 3 labels 3 seconds cclabelttf *three = [cclabelttf labelwithstring:@"3" fontname:@"helvetica" fontsize:100]; three.position = ccp(screencenter.x, screencenter.y); cclabelttf *two = [cclabelttf labelwithstring:@"2" fontname:@"helvetica" fontsize:100]; two.position = ccp(screencenter.x, screencenter.y); cclabelttf *one = [cclabelttf labelwithstring:@"1" fontname:@"helvetica" fontsize:100]; one.position = ccp(screencenter.x, screencenter.y); // secondspast specified in update method secondspast = 0; if (secondspast == 1) { [self addchild:thr...
Read more

hyperlink - how to do url routing in php -

sorry, don't know should exact technical term problem statement title of thread might misleading/confusing. that's reason why not able search problem solution on google. okay, i'm developing simple social network users. want each user have own profile link. for example, if username of user 'xyz' profile link should ' http://example.com/user/xyz now how create "virtual" link? edit: okay, learned it's mod_rewrite , got similar question how dynamically rewrite url facebook routes php delicate matter. mod_rewrite can in beginning, after becoming mess , f*** rather help, when "big" software needs scalable. remembering "routes" not part of specific application, nomenclature "server" interprets request "client" , returns specific value "client" configuration of request. ie, "router" server. what recommend study http protocol, rest , on. this surely in many things, knowi...
Read more

c - Avoiding Extra Malloc in Linked List (node->next = NULL) -

in linked list, i'm trying avoid malloc ing node without adding bunch of if statements , such. have following: polynomial create() { polynomial head = null; polynomial temp = null; int numterms, coefficient, exponent; int counter = 0; printf("enter number of terms in polynomial: "); scanf ("%d", &numterms); printf("\n"); if (numterms == 0) { head = malloc(sizeof(term)); head->coef = 0; head->exp = 0; head->next = null; return head; } while (numterms != counter) { // ask input printf("enter coefficient , exponent of term %d: ", (counter + 1)); scanf("%d %d", &coefficient, &exponent); printf("\n"); // create term if (temp == null) temp = malloc(sizeof(term)); temp->coef = coefficient; temp->exp = exponent; temp->next = null; //if((numterms - 1) != counter) tem...
Read more