uires removing mprotect (see
http://stackoverflow.com/questions/2663456/write-a-signal-handler-to-catch-sigsegv)
and it's impossible to know when to lock the region by mprotect again.
Thanks,
-Kibeom Kim
___
NumPy-Discussion mailing list
NumPy-
mexpr?
Thanks.
-Kibeom Kim
___
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
http://mail.scipy.org/mailman/listinfo/numpy-discussion
I wonder why numpy.norm doesn't call blas norm function by default.
is there any good reason? or it's just not coded?
-Kibeom Kim
On Sun, Mar 6, 2011 at 12:13 PM, Till Stensitzki wrote:
>
> >
> > Moreover, np.linalg.norm is slow compare to blas.
> >
&g
for the sum function part,
It's not a good way to fix but... if you want more accuracy
x=(np.array([1]*1 + [1e4], dtype=np.float32))
np.sum(x*x)
1.0001e+08
You can sort x from small numbers to bigger numbers before you call sum.
-Kibeom Kim
On Sat, Mar 5, 2011 at 6:27 PM, Xavier