Commit 1ea676ef authored by Kirill Terekhov's avatar Kirill Terekhov

Prevent breakdown in bcgs

parent b5eb3ba5
......@@ -1685,7 +1685,8 @@ perturbate:
#pragma omp single
#endif
{
tempa = tempb = 0;
tempa = 0.0; //can be negative
tempb = 1.0e-100; //tempb is positive since squared
}
#if defined(USE_OMP)
#pragma omp for reduction(+:tempa,tempb)
......
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