Commit 1c791cc2 authored by Alvarez, Gonzalo's avatar Alvarez, Gonzalo
Browse files

printritz works when useComplex

parent 862e1c1a
......@@ -292,6 +292,10 @@ private:
if (T.n_col()!=T.n_row()) throw PsimagLite::RuntimeError("T is not square\n");
if (V.n_col()!=T.n_col()) throw PsimagLite::RuntimeError("V is not nxn2\n");
typename LanczosSolverType::ParametersSolverType params(ioIn_, "Tridiag");
bool b1 = (params.options.find("printritz") != PsimagLite::String::npos);
if (b1) printRitz(eigs, i0);
ComplexOrRealType zone = 1.0;
ComplexOrRealType zzero = 0.0;
......@@ -314,10 +318,6 @@ 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,
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment