diff options
author | alex <argriffi@ncsu.edu> | 2015-06-19 21:52:25 -0400 |
---|---|---|
committer | alex <argriffi@ncsu.edu> | 2015-06-19 21:52:25 -0400 |
commit | 1f79addc9d546690984896fc4a35298aab289ba5 (patch) | |
tree | 518a76ddea0521200e9c8db1ce2055c3e77e8eb2 /numpy/polynomial | |
parent | 0514d31299318bfa4c87abbfea6424eebaeec188 (diff) | |
download | numpy-1f79addc9d546690984896fc4a35298aab289ba5.tar.gz |
MAINT: improve laguerre and legendre efficiency
Diffstat (limited to 'numpy/polynomial')
-rw-r--r-- | numpy/polynomial/laguerre.py | 3 | ||||
-rw-r--r-- | numpy/polynomial/legendre.py | 1 |
2 files changed, 1 insertions, 3 deletions
diff --git a/numpy/polynomial/laguerre.py b/numpy/polynomial/laguerre.py index 8d2705d5d..fffe9e6b6 100644 --- a/numpy/polynomial/laguerre.py +++ b/numpy/polynomial/laguerre.py @@ -1688,8 +1688,7 @@ def laggauss(deg): # matrix is symmetric in this case in order to obtain better zeros. c = np.array([0]*deg + [1]) m = lagcompanion(c) - x = la.eigvals(m) - x.sort() + x = la.eigvalsh(m) # improve roots by one application of Newton dy = lagval(x, c) diff --git a/numpy/polynomial/legendre.py b/numpy/polynomial/legendre.py index e97de703f..c91cb72ec 100644 --- a/numpy/polynomial/legendre.py +++ b/numpy/polynomial/legendre.py @@ -1712,7 +1712,6 @@ def leggauss(deg): c = np.array([0]*deg + [1]) m = legcompanion(c) x = la.eigvalsh(m) - x.sort() # improve roots by one application of Newton dy = legval(x, c) |