diff --git a/MathLib/LinAlg/Solvers/BiCGStab.cpp b/MathLib/LinAlg/Solvers/BiCGStab.cpp
index b81eac29e978e1094e9e0d4e634095ae7117426b..a63da21968505448aa71a3cecc8389e37f8d6f0e 100644
--- a/MathLib/LinAlg/Solvers/BiCGStab.cpp
+++ b/MathLib/LinAlg/Solvers/BiCGStab.cpp
@@ -59,7 +59,7 @@ unsigned BiCGStab(CRSMatrix<double> const& A, double* const b, double* const x,
 			blas::copy(N, r, p); // p = r
 		else {
 //			blas::axpy(N, -omega, v, p); // p = (p-omega v)*beta+r
-//			const double beta = rho1 * alpha / (rho2 * omega);
+			const double beta = rho1 * alpha / (rho2 * omega);
 //			blas::scal(N, beta, p);
 //			blas::axpy(N, D_ONE, r, p);
 			// p = (p-omega v)*beta+r