Skip to content
Snippets Groups Projects
Commit b8a34ab1 authored by Tomáš Oberhuber's avatar Tomáš Oberhuber
Browse files

Apply suggestion to src/TNL/Solvers/Linear/TFQMR_impl.h

parent 21a5cf18
No related branches found
No related tags found
1 merge request!34Runge kutta
......@@ -89,7 +89,7 @@ bool TFQMR< Matrix >::solve( ConstVectorViewType b, VectorViewType x )
}
w.addVector( Au, -alpha );
d.addVector( u, 1.0, theta * theta * eta / alpha );
w_norm = lpNorm( w.getView(), 2.0 );
w_norm = lpNorm( w, 2.0 );
theta = w_norm / tau;
const RealType c = 1.0 / std::sqrt( 1.0 + theta * theta );
tau = tau * theta * c;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment