python - Will results of numpy.as_strided depend on input dtype? -


will results of numpy.lib.stride_tricks.as_strided depend on dtype of numpy array?

this question arises definition of .strides, is

tuple of bytes step in each dimension when traversing array.

take following function i've used in other questions here. takes 1d or 2d array , creates overlapping windows of length window. result 1 dimension greater input.

def rwindows(a, window):     if a.ndim == 1:         = a.reshape(-1, 1)     shape = a.shape[0] - window + 1, window, a.shape[-1]     strides = (a.strides[0],) + a.strides     windows = np.lib.stride_tricks.as_strided(a, shape=shape, strides=strides)     return np.squeeze(windows)  # examples # rwindows(np.arange(5), window=2) # rwindows(np.arange(20).reshape((5,4)), window=2)

because of definition of strides , because, instance, otherwise equivalent arrays of dtype float32 , float64 have different strides, ever blow rwindows function above?

i've tried test it's been in non-exhaustive way , looking answer (1) explains whether disclaimer/warning function doc has i'm asking here , (2) explains why or why not otherwise-equivalent arrays different dtypes & strides yield different results in above.

no, warning as_strided 2 issues not related size of data , more result writing resulting view.

  1. first, there no protection assure view = as_strided(a . . . ) only points memory in a. why there deliberate preparation work done before calling as_strided. if algorithm off, can have view point memory not in a, , may indeed addressed garbage, other variables, or operating system. if write view, data can lost, misplaced, corrupted . . . or crash computer.

for specific example, how safe depends lot on inputs you're using. you've set strides a.strides dynamic. may want assert dtype of a isn't weird object.

if you're sure always have 2-d a larger window, fine algorithm, can assert make sure. if not, may want make sure as_strided output works n-d a arrays. instance:

shape = a.shape[0] - window + 1, window, a.shape[-1]

should be

shape = (a.shape[0] - window + 1, window) + a.shape[1:]

in order accept n-d input. probably never problem far referencing bad memory, current shape reference wrong data in a if had more dimensions.

  1. second, view created references same data blocks multiple times. if parallel write view (through view = foo or bar( . . ., out = view)), results can unpredictable , not expect.

that said, if afraid of problems , don't need write as_strided view (as don't convolution applications commonly used), can set writable = false, prevent both problems if strides and/or shape incorrect.

edit: pointed out @hpaulj, in addition 2 problems, if view makes copy (like .flatten() or fancy indexing large chunk of it), can cause memoryerror.


Comments

Post a Comment

Popular posts from this blog

python Tkinter Capturing keyboard events save as one single string -

is possible store keyboard events 1 string? code below stores 1 char , prints it. card reader or bar code reader, contains collection of character/ string not 1 character @ time. goal save char pressed text variable. from tkinter import * root = tk() def key(event): text= event.char text+= event.char print ("pressed", text) def callback(event): frame.focus_set() print ("clicked at", event.x, event.y) frame = frame(root, width=100, height=100) frame.bind("<key>", key) frame.bind("<button-1>", callback) frame.pack() root.mainloop() currently, creating text variable , everytime key() function gets called, , text stores last character typed. you can define text module level variable , , use module level text inside key function - text = '' def key(event): global text text+= event.char print("pressed", text)
Read more

android - InAppBilling registering BroadcastReceiver in AndroidManifest -

i'm implementing inapp billing in android app. works well, however, i'm trying decouple broadcast receiver activity manifest. suggestion in android's trivialdrive sample: // important: dynamically register broadcast messages updated purchases. // register receiver here instead of <receiver> in manifest // because call getpurchases() @ startup, therefore can ignore // broadcasts sent while app isn't running. // note: registering listener in activity bad idea, done here // because sample. regardless, receiver must registered after // iabhelper setup, before first call getpurchases(). currently there's class extends broadcastreceiver : public class iabbroadcastreceiver extends broadcastreceiver { /** * intent action receiver should filter for. */ public static final string action = "com.android.vending.billing.purchases_updated"; private final iabbroadcastlistener mlistener; public iabbroadcastreceiver(iabbroadcastlistener listener) { m
Read more

javascript - Z-index in d3.js -

i have map , when hover effect on region of map below yellow line , circle, last ones disappear. how can solve problem? http://jsfiddle.net/ot21n9qx/ var arc = g.append("path") .style("fill", "none") .style("stroke", "yellow") .style("stroke-width", 2) .attr("d", "m" + pathorigin[0] + "," + pathorigin[1] + " q" + svgpoint[0] + "," + pathorigin[1] + " " + svgpoint[0] + "," + svgpoint[1]); var circlesize = d3.scale.linear() .domain([0, 0.5, 1]) .range([4, 10, 4]); svg not have z order draws objects in order created. need create path , circle after create highlight regions. here forked fiddle: http://jsfiddle.net/o00fzgaf/1/
Read more