From 052df9b8795891c5ed692e3a9e8231db97d97562 Mon Sep 17 00:00:00 2001
From: jiri <jiri>
Date: Tue, 12 Sep 2006 20:27:41 +0000
Subject: [PATCH] *** empty log message ***

---
 nonclass.h | 22 ++++++++++++++++++++++
 1 file changed, 22 insertions(+)

diff --git a/nonclass.h b/nonclass.h
index 00c5cda..3e7c89e 100644
--- a/nonclass.h
+++ b/nonclass.h
@@ -233,4 +233,26 @@ return cblas_ddot(n,x,incx,y,incy);
 }
 
 
+
+//debugging aid: reconstruct an explicit matrix from the implicit version
+//which provides gemv only
+template<typename M, typename T>
+NRMat<T> reconstructmatrix(const M &implicitmat)
+{
+NRMat<T> r(implicitmat.nrows(),implicitmat.ncols());
+NRVec<T> rhs(0.,implicitmat.ncols());
+NRVec<T> tmp(implicitmat.nrows());
+for(int i=0; i<implicitmat.ncols(); ++i)
+	{
+	rhs[i]=1.;
+	implicitmat.gemv(0.,tmp,'n',1.,rhs);
+	for(int j=0; j<implicitmat.nrows(); ++j) r(j,i)=tmp[j];
+	rhs[i]=0.;
+	}
+return r;
+}
+
+
+
+
 #endif