2011-03-27 Brian Gough <bjg@network-theory.co.uk>
* jacobi.c (norm): compute norm of off-diagonal elements only as
per algorithm 8.4.3 in Golub and van Loan
2009-11-26 Brian Gough <bjg@network-theory.co.uk>
* qrstep.c (qrstep): handle the case where |mu| >> |d0| + |sd0| to
prevent infinite loop
2009-07-09 Brian Gough <bjg@network-theory.co.uk>
* symmv.c (gsl_eigen_symmv_free): handle NULL argument in free
* symm.c (gsl_eigen_symm_free): handle NULL argument in free
* nonsymmv.c (gsl_eigen_nonsymmv_free): handle NULL argument in free
* nonsymm.c (gsl_eigen_nonsymm_free): handle NULL argument in free
* hermv.c (gsl_eigen_hermv_free): handle NULL argument in free
* herm.c (gsl_eigen_herm_free): handle NULL argument in free
* genv.c (gsl_eigen_genv_free): handle NULL argument in free
* gensymmv.c (gsl_eigen_gensymmv_free): handle NULL argument in free
* gensymm.c (gsl_eigen_gensymm_free): handle NULL argument in free
* genhermv.c (gsl_eigen_genhermv_free): handle NULL argument in free
* genherm.c (gsl_eigen_genherm_free): handle NULL argument in free
* gen.c (gsl_eigen_gen_free): handle NULL argument in free
* francis.c (gsl_eigen_francis_free): handle NULL argument in free
2008-07-03 Brian Gough <bjg@hp2.network-theory.co.uk>
* Makefile.am (INCLUDES): use top_srcdir instead of top_builddir
2007-08-30 Brian Gough <bjg@network-theory.co.uk>
* test.c (test_eigen_symm): added new test case for underflow
* qrstep.c (trailing_eigenvalue): handle underflow for small tab,
also handle dt==0 directly.
2007-07-30 Brian Gough <bjg@network-theory.co.uk>
* testgen.c (main): use gsl_finite instead of finite
2006-09-14 Brian Gough <bjg@network-theory.co.uk>
* test.c (test_eigen_symm): fix duplicate VAL_DESC to ABS_DESC
(test_eigen_herm): fix duplicate VAL_DESC to ABS_DESC
2006-08-14 Brian Gough <bjg@network-theory.co.uk>
* unsymm.c (gsl_eigen_unsymm): support for unsymmetric matrices
(P.Alken)
2006-04-18 Brian Gough <bjg@network-theory.co.uk>
* test.c (test_eigenvalues): increase error bound
2006-03-26 Brian Gough <bjg@network-theory.co.uk>
* jacobi.c (gsl_eigen_invert_jacobi): use unsigned int for nrot
consistently
2005-08-22 Brian Gough <bjg@network-theory.co.uk>
* test.c (test_eigenvalues): increased test tolerance
2004-12-03 Brian Gough <bjg@network-theory.co.uk>
* jacobi.c (gsl_eigen_jacobi): use algorithm from Golub and Van
Loan
(gsl_eigen_invert_jacobi): new code
2003-01-02 Brian Gough <brian.gough@network-theory.co.uk>
* symmv.c (gsl_eigen_symmv): change gsl_isnan to isnan
* symm.c (gsl_eigen_symm): change gsl_isnan to isnan
* hermv.c (gsl_eigen_hermv): change gsl_isnan to isnan
* herm.c (gsl_eigen_herm): change gsl_isnan to isnan
2002-11-16 Brian Gough <bjg@network-theory.co.uk>
* symm.c (gsl_eigen_symm): prevent infinite loop for NaNs
* herm.c (gsl_eigen_herm): prevent infinite loop for NaNs
* symmv.c (gsl_eigen_symmv): prevent infinite loop for NaNs
* hermv.c (gsl_eigen_hermv): prevent infinite loop for NaNs
Fri Oct 18 17:58:35 2002 Brian Gough <bjg@network-theory.co.uk>
* test.c (main): added test cases to catch division by zero
(beta_r=0) in householdercomplex.c
Wed Aug 28 19:05:54 2002 Brian Gough <bjg@network-theory.co.uk>
* Makefile.am (test_LDADD): fix link order to have libgslsys near
end
Thu Aug 2 18:17:58 2001 Brian Gough <bjg@network-theory.co.uk>
* hermv.c (gsl_eigen_hermv): fixed index bug in accumulation of evec.
* symmv.c (gsl_eigen_symmv): fixed index bug in accumulation of evec.
* test.c (main): added two new test cases
* qrstep.c (trailing_eigenvalue): chose better value of mu when
dt=0, prevents infinite loop.
Sun Jul 1 22:43:45 2001 Brian Gough <bjg@network-theory.co.uk>
* modified to use new-style vector views
Wed Jun 20 12:30:38 2001 Brian Gough <bjg@network-theory.co.uk>
* hermv.c (gsl_eigen_hermv): additional workspace argument no
longer required
Mon Apr 23 10:31:01 2001 Brian Gough <bjg@network-theory.co.uk>
* unified error handling conventions to _e for error handling
functions and no suffix for plain functions, so _impl functions
are no longer needed.
* removed tests for EFAULT, since EFAULT should only apply to
invalid non-null pointers.
Fri Apr 13 20:33:18 2001 Brian Gough <bjg@network-theory.co.uk>
* eigen/test.c (test_invert_jacobi): removed matmult and replaced
by gsl_blas_dgemm