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