|
@@ -496,8 +496,8 @@ namespace Lemma {
|
496
|
496
|
CSolver = new Eigen::BiCGSTAB<Eigen::SparseMatrix<Complex, Eigen::RowMajor>, Eigen::IncompleteLUT<Complex> > [Omegas.size()];
|
497
|
497
|
for (int iw=0; iw<Omegas.size(); ++iw) {
|
498
|
498
|
Csym = Cvec[iw].selfadjointView<Eigen::Lower>();
|
499
|
|
- CSolver[iw].preconditioner().setDroptol(1e-3);
|
500
|
|
- //CSolver[iw].preconditioner().setFillfactor(1e3);
|
|
499
|
+ CSolver[iw].preconditioner().setDroptol(1e-12);
|
|
500
|
+ CSolver[iw].preconditioner().setFillfactor(1e2);
|
501
|
501
|
jsw_timer timer;
|
502
|
502
|
timer.begin();
|
503
|
503
|
/* Complex system */
|