diff options
Diffstat (limited to 'src/lib/netlist/plib/gmres.h')
-rw-r--r-- | src/lib/netlist/plib/gmres.h | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/src/lib/netlist/plib/gmres.h b/src/lib/netlist/plib/gmres.h index 840da0cfb45..686e1444d36 100644 --- a/src/lib/netlist/plib/gmres.h +++ b/src/lib/netlist/plib/gmres.h @@ -259,7 +259,7 @@ namespace plib m_ht[j][k] = vec_mult<FT>(n, m_v[kp1], m_v[j]); vec_add_mult_scalar(n, m_v[kp1], m_v[j], -m_ht[j][k]); } - m_ht[kp1][k] = std::sqrt(vec_mult2<FT>(n, m_v[kp1])); + m_ht[kp1][k] = plib::sqrt(vec_mult2<FT>(n, m_v[kp1])); // FIXME: comparison to zero if (m_ht[kp1][k] != plib::constants<FT>::zero()) @@ -268,7 +268,7 @@ namespace plib for (std::size_t j = 0; j < k; j++) givens_mult(m_c[j], m_s[j], m_ht[j][k], m_ht[j+1][k]); - const float_type mu = reciprocal(std::hypot(m_ht[k][k], m_ht[kp1][k])); + const float_type mu = reciprocal(plib::hypot(m_ht[k][k], m_ht[kp1][k])); m_c[k] = m_ht[k][k] * mu; m_s[k] = -m_ht[kp1][k] * mu; @@ -277,7 +277,7 @@ namespace plib givens_mult(m_c[k], m_s[k], m_g[k], m_g[kp1]); - FT rho = std::abs(m_g[kp1]); + FT rho = plib::abs(m_g[kp1]); // FIXME .. itr_used = itr_used + 1; @@ -359,12 +359,12 @@ namespace plib ops.solve_inplace(Ax); - const float_type rho_to_accuracy = std::sqrt(vec_mult2<FT>(n, Ax)) / accuracy; + const float_type rho_to_accuracy = plib::sqrt(vec_mult2<FT>(n, Ax)) / accuracy; rho_delta = accuracy * rho_to_accuracy; } else - rho_delta = accuracy * std::sqrt(static_cast<FT>(n)); + rho_delta = accuracy * plib::sqrt(static_cast<FT>(n)); /* * Using @@ -389,7 +389,7 @@ namespace plib ops.solve_inplace(residual); - rho = std::sqrt(vec_mult2<FT>(n, residual)); + rho = plib::sqrt(vec_mult2<FT>(n, residual)); if (rho < rho_delta) return itr_used + 1; |