1 #ifndef _CONJUGATE_GRADIENT_H
2 #define _CONJUGATE_GRADIENT_H
4 #include <valarray>
6 double
7 inner(std::valarray<double> const &x,
8 std::valarray<double> const &y);
10 void
11 conjugate_gradient(double **A,
12 double *x,
13 double *b,
14 unsigned n,
15 double tol,
16 unsigned max_iterations);
17 void
18 conjugate_gradient(std::valarray<double> const &A,
19 std::valarray<double> &x,
20 std::valarray<double> const &b,
21 unsigned n, double tol,
22 unsigned max_iterations);
23 #endif // _CONJUGATE_GRADIENT_H