diff --git a/src/svd.hh b/src/svd.hh index bb2ba97911cf8066cad6b659e9aa02c6d11aa348..8374d1ca318ee30281edaf815f74919d4e42db76 100644 --- a/src/svd.hh +++ b/src/svd.hh @@ -4,7 +4,7 @@ #ifndef SVD_HH #define SVD_HH -#include <cmath.h> +#include <cmath> template <class T> @@ -30,7 +30,7 @@ T pythag(T a, T b) /** Given a matrix a[1..m][1..n], this routine computes its singular value decomposition, A = - U ·W ·V^T . The matrix U replaces a on output. The diagonal matrix of singular values W is out- + U W V^T . The matrix U replaces a on output. The diagonal matrix of singular values W is out- put as a vector w[1..n]. The matrix V (not the transpose V T ) is output as v[1..n][1..n]. */ template <class T, int m, int n>