Hi,<br><br> I am also interested in this. In my application there is a large 2d array, lets call it 'b' to keep the notation consistent in the thread.  b's  columns need to be recomputed often. Ideally this re-computation happens in a function. Lets call that function updater(b, col_index): The simplest example is where<br>
updater(b, col_index) is a matrix vector multiply, where the matrix or the vector changes.<br><br> Is there anyway apart from using ufuncs that I can make updater() write the result directly in b and not create a new temporary column that is then copied into b ?  Say for the matrix vector multiply example.<br>
I can write the matrix vector product in terms of ufuncs but will lose out in terms of speed.<br><br>In the best case scenario I would like to maintain 'b' in a csr sparse matrix form, as 'b' participates in a matrix vector multiply. I think csr would be asking for too much, but even ccs should help.  I dont want to clutter this thread with the sparsity issues though, any solution to the original question or pointers to solutions would be appreciated.<br>
<br>Thanks<br>  --srean<br><br><div class="gmail_quote">On Sat, Mar 26, 2011 at 12:10 PM, Hugo Gagnon <span dir="ltr"><<a href="mailto:sourceforge.numpy@user.fastmail.fm">sourceforge.numpy@user.fastmail.fm</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hello,<br>
<br>
Say I have a few 1d arrays and one 2d array which columns I want to be<br>
the 1d arrays.<br>
I also want all the a's arrays to share the *same data* with the b<br>
array.<br>
If I call my 1d arrays a1, a2, etc. and my 2d array b, then<br>
<br>
b[:,0] = a1[:]<br>
b[:,1] = a2[:]<br>
...<br>
<br>
won't work because apparently copying occurs.<br>
I tried it the other way around i.e.<br>
<br>
a1 = b[:,0]<br>
a2 = b[:,1]<br>
...<br>
<br>
and it works but that doesn't help me for my problem.<br>
Is there a way to reformulate the first code snippet above but with<br>
shallow copying?<br>
<br>
Thanks,<br>
--<br>
  Hugo Gagnon<br>
<font color="#888888">--<br>
  Hugo Gagnon<br>
<br>
_______________________________________________<br>
NumPy-Discussion mailing list<br>
<a href="mailto:NumPy-Discussion@scipy.org">NumPy-Discussion@scipy.org</a><br>
<a href="http://mail.scipy.org/mailman/listinfo/numpy-discussion" target="_blank">http://mail.scipy.org/mailman/listinfo/numpy-discussion</a><br>
</font></blockquote></div><br>