Changing BiCGSTAB parameters initialization so that it works with custom types

This commit is contained in:
Alexandre Hoffmann 2022-11-29 19:37:46 +00:00 committed by Rasmus Munk Larsen
parent ab2b26fbc2
commit 23524ab6fc

View File

@ -51,9 +51,9 @@ bool bicgstab(const MatrixType& mat, const Rhs& rhs, Dest& x,
x.setZero();
return true;
}
Scalar rho = 1;
Scalar alpha = 1;
Scalar w = 1;
Scalar rho (1);
Scalar alpha (1);
Scalar w (1);
VectorType v = VectorType::Zero(n), p = VectorType::Zero(n);
VectorType y(n), z(n);