summaryrefslogtreecommitdiffstatshomepage
path: root/src/lib/netlist/plib/gmres.h
diff options
context:
space:
mode:
Diffstat (limited to 'src/lib/netlist/plib/gmres.h')
-rw-r--r--src/lib/netlist/plib/gmres.h8
1 files changed, 4 insertions, 4 deletions
diff --git a/src/lib/netlist/plib/gmres.h b/src/lib/netlist/plib/gmres.h
index edda011dc19..352ce721722 100644
--- a/src/lib/netlist/plib/gmres.h
+++ b/src/lib/netlist/plib/gmres.h
@@ -8,9 +8,9 @@
#ifndef PLIB_GMRES_H_
#define PLIB_GMRES_H_
-#include "pconfig.h"
#include "mat_cr.h"
#include "parray.h"
+#include "pconfig.h"
#include "vector_ops.h"
#include <algorithm>
@@ -246,13 +246,13 @@ namespace plib
if (rho < rho_delta)
return itr_used + 1;
- vec_set_scalar(RESTART+1, m_g, NL_FCONST(0.0));
+ vec_set_scalar(RESTART+1, m_g, constants<FT>::zero);
m_g[0] = rho;
//for (std::size_t i = 0; i < mr + 1; i++)
// vec_set_scalar(mr, m_ht[i], NL_FCONST(0.0));
- vec_mult_scalar(n, residual, NL_FCONST(1.0) / rho, m_v[0]);
+ vec_mult_scalar(n, residual, constants<FT>::one / rho, m_v[0]);
for (std::size_t k = 0; k < RESTART; k++)
{
@@ -269,7 +269,7 @@ namespace plib
m_ht[kp1][k] = std::sqrt(vec_mult2<FT>(n, m_v[kp1]));
if (m_ht[kp1][k] != 0.0)
- vec_scale(n, m_v[kp1], NL_FCONST(1.0) / m_ht[kp1][k]);
+ vec_scale(n, m_v[kp1], constants<FT>::one / m_ht[kp1][k]);
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]);