From 46d841aabf9d30487eecb796fb5b2fa27597b3f0 Mon Sep 17 00:00:00 2001
From: jiri <jiri>
Date: Mon, 31 Jan 2005 13:43:38 +0000
Subject: [PATCH] *** empty log message ***

---
 davidson.h | 9 +++++++++
 1 file changed, 9 insertions(+)
 create mode 100644 davidson.h

diff --git a/davidson.h b/davidson.h
new file mode 100644
index 0000000..9adfbc8
--- /dev/null
+++ b/davidson.h
@@ -0,0 +1,9 @@
+
+//Davidson diagonalization of real symmetric matrix
+
+//matrix can be any class which provides nrows(), diagonalof() and operator*(const  NRVec<T>) methods
+//does not even have to be explicitly stored - direct CI
+//n<0 highest eigenvalues, n>0 lowest eigenvalues
+
+template <class T, class matrix>
+extern void davidson(const matrix &a, NRVec<T> *vecs /*input-output*/, T *vals, const int n=1, const double eps=1e-10);