| Commit message (Collapse) | Author | Age | Files | Lines |
... | |
| |
|
|
|
|
|
|
|
| |
Move the dotblas_matrixproduct function from the _dotblas.c file to a
new cblasfuncs.c file in multiarray and rename it cblas_matrixproduct.
Modify it so that it can called directly from PyArray_MatrixProduct2 and
do so. Fix up numeric.py and core/__init__.py to reflect these changes.
|
|
|
|
|
| |
allows passing flags like --pdb to test files
also add call to files where its missing
|
|
|
|
| |
Allows changing the error message at some point.
|
| |
|
|
|
|
| |
https://github.com/numpy/numpy/issues/4733
|
| |
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
When a multi index is tracked and RemoveAxis can be called, the
size of the iterator may still change. This was causing failures
for example for the SVD, because the gufunc machinery requires
a temporarily larger iterator for output allocation.
Thanks to Jaime (jaime.frio@gmail.com) for noting that this is
plausible since the size check can be delayed pretty ok up
until GetIterNext (or similar functions).
Closes gh-4442
|
|
|
|
| |
residuals array
|
|
|
|
|
|
|
| |
The variable used was not defined, causing an error when the exception
was raised.
Closes #3146.
|
|
|
|
|
| |
The removal of CFLOAT for eig/eigvals is related to a bug, see c355550dd80.
Can be put back once that issue is fixed.
|
|
|
|
|
|
|
| |
Some symbols are not defined leading to:
ImportError Undefined symbol: PyGILState_Release
Closes gh-3255
|
| |
|
|\
| |
| | |
ENH: use dot in linalg.norm
|
| |
| |
| |
| | |
Also used ravel(order='K') to prevent copies.
|
| | |
|
|\ \
| | |
| | | |
DOC: matrix_rank
|
| | | |
|
| | |
| | |
| | |
| | |
| | |
| | | |
asfarray truncates the complex part, so it must be avoided for complex
types.
Closes gh-4156.
|
|/ /
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| | |
This fixes the following bug with longdouble:
>>> x = np.arange(10, dtype=np.longdouble)
>>> np.linalg.norm(x, ord=3)
Traceback (most recent call last):
File "<ipython-input-5-7ee53a8ac142>", line 1, in <module>
np.linalg.norm(x, ord=3)
File "/tmp/v/lib/python2.7/site-packages/numpy/linalg/linalg.py", line 2090, in norm
return add.reduce(absx**ord, axis=axis)**(1.0/ord)
UnboundLocalError: local variable 'absx' referenced before assignment
As well as the handling of minimal values for signed integers:
>>> x = np.array([-2**31], dtype=np.int32)
>>> np.linalg.norm(x, ord=3)
/tmp/v/lib/python2.7/site-packages/numpy/linalg/linalg.py:2090: RuntimeWarning: invalid value encountered in double_scalars
return add.reduce(absx**ord, axis=axis)**(1.0/ord)
nan
|
|\ \
| | |
| | | |
ENH: avoid expensive clears in fenv functions
|
| | |
| | |
| | |
| | |
| | |
| | |
| | |
| | |
| | |
| | | |
Clearing is 50-100 times more expensive than checking on x86, so check if there
is anything needs to be cleared first. This speeds up scalar operations
by 10%-20%.
Based on Arink Verma code in #3739.
Implement the functions as new C-API functions npy_get_floatstatus and
npy_clear_floatstatus in npy_math.
|
| | |
| | |
| | |
| | |
| | |
| | | |
Do not convert the UPLO argument to a bytestring, it is not necessary.
Distribute parts of the eigh tests into the appropriate TestEigvalsh or
TestEigh test class.
|
| | |
| | |
| | |
| | | |
to prevent unintentional use of wrong function. Restores 1.7 behavior.
|
| | |
| | |
| | |
| | |
| | |
| | | |
UPLO was cast to bytes and compared to a string which is always false in
python3.
closes gh-3977
|
|/ / |
|
| |
| |
| |
| |
| |
| | |
Strictly speaking, any non-standard name with "__" in it causes
undefined behavior in C. It's also completely unnecessary to use
any in this module.
|
| | |
|
|/
|
|
| |
At least OSX Accelerate fails for this case.
|
| |
|
| |
|
|
|
|
|
| |
Some of the functions removed were untested, and apparently did not
function correctly (cholesky_up, inner1d, maybe more).
|
| |
|
| |
|
|
|
|
|
|
| |
Also remove TestCase subclassing, so that generator tests work. Also
fix bugs in the existing generator tests that were not actually run
previously.
|
|
|
|
|
| |
This was previously correct, but the special case was missing
since the gufunc code cannot handle it.
|
|
|
|
|
|
| |
uninitialized use in _strided_to_strided_one_to_n_with_finish
wrong error NULL checks in PyArray_MapIterNew and einsum_sub_op_from_str
double free and memleak in umath_linalg.c.src init_@lapack_func@
|
|
|
|
|
|
|
|
| |
The qr factorization is not unique and the values returned by the raw
mode may differ between LAPACK versions. Consequently, the results
cannot be checked against known values.
Closes #3703.
|
|
|
|
|
|
|
|
|
|
|
| |
Combine the xerbla override tests with the previous _is_xerbla_safe
decorator. The xerbla checks are then made in a forked process where the
outputs are closed so that if xerbla prints a message it does not
appear.
It would be nice to figure out a Windows version of this test.
Closes #3704. Closes #3687.
|
|
|
|
|
|
|
| |
Run the 2to3 ws_comma fixer on *.py files. Some lines are now too long
and will need to be broken at some point. OTOH, some lines were already
too long and need to be broken at some point. Now seems as good a time
as any to do this with open PRs at a minimum.
|
|
|
|
| |
Now is as good a time as any with open PR's at a low.
|
|\
| |
| | |
Mistake in linalg.eig() documentation
|
| | |
|
| |
| |
| |
| |
| |
| |
| | |
This makes the inverse of a 0x0 array simply be 0x0 again. It
also modifies the no-empty array check in favor of a no-empty
*inner* array, since the gufuncs seem to handle the other case
fine.
|
| | |
|
| | |
|
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| | |
Now that Python < 2.6 is no longer supported we can use the errstate
context manager in places where constructs like
```
old = seterr(invalid='ignore')
try:
blah
finally:
seterr(**old)
```
were used.
|
| | |
|
| | |
|
| | |
|