TridiagOptions accepts printritz

TestSuite/inputs/input250.inp

@@ -18,9 +18,9 @@ Version=version
 OutputFile=data250
 InfiniteLoopKeptStates=100
 FiniteLoops 3
-  7 100 8
--14 100 8
- 14 100 8
+  7 100 0 
+-14 100 0
+ 14 100 0
 NumberOfTargetQns=2
 TargetElectronsUp0=8
 TargetElectronsDown0=8

src/Engine/CorrectionVectorSkeleton.h

@@ -314,6 +314,10 @@ private:
 
 		xr.resize(n);
 		psimag::BLAS::GEMV('N',n,n2,zone,&(V(0,0)),n,&(tmp[0]),1,zzero,&(xr[0]),1);
+
+		typename LanczosSolverType::ParametersSolverType params(ioIn_, "Tridiag");
+		bool b1 = (params.options.find("printritz") != PsimagLite::String::npos);
+		if (b1) printRitz(eigs, i0);
 	}
 
 	void triDiag(const VectorWithOffsetType& phi,
@@ -347,6 +351,15 @@ private:
 		return sum;
 	}
 
+	static void printRitz(const VectorRealType& eigs, SizeType i0)
+	{
+		std::cout<<"RitzEigenvalues: index="<<i0<<"|";
+		const SizeType n = eigs.size();
+		for (SizeType i = 0; i < n; ++i)
+			std::cout<<eigs[i]<<" ";
+		std::cout<<"\n";
+	}
+
 	InputValidatorType& ioIn_;
 	const TargetParamsType& tstStruct_;
 	const ModelType& model_;