diff options
| -rw-r--r-- | src/matrix_orthogonal.c | 1 | ||||
| -rw-r--r-- | src/matrix_pinv.c | 4 | ||||
| -rw-r--r-- | src/spherical_line.c | 3 |
3 files changed, 7 insertions, 1 deletions
diff --git a/src/matrix_orthogonal.c b/src/matrix_orthogonal.c index 8acf650..54a38b8 100644 --- a/src/matrix_orthogonal.c +++ b/src/matrix_orthogonal.c @@ -3,6 +3,7 @@ iem_matrix written by Thomas Musil (c) IEM KUG Graz Austria 2002 - 2006 */ +#include "m_pd.h" /* -------------------------- matrix_orthogonal ------------------------------ */ diff --git a/src/matrix_pinv.c b/src/matrix_pinv.c index f1ede7b..c1ef278 100644 --- a/src/matrix_pinv.c +++ b/src/matrix_pinv.c @@ -3,6 +3,8 @@ iem_matrix written by Thomas Musil (c) IEM KUG Graz Austria 2002 - 2006 */ +#include "m_pd.h" + /* -------------------------- matrix_pinv ------------------------------ */ typedef struct _matrix_pinv @@ -400,7 +402,7 @@ static void matrix_pinv_mul2(t_matrix_pinv *x) { int n_ls=x->x_n_ls+x->x_n_phls; int n_ambi=x->x_n_ambi; - int n_ambi2=2*n_ambi; + int n_ambi2=2*x->x_n_ambi; int i, j, k; double *vec1, *beg1=x->x_transp; double *vec2, *beg2=x->x_inv_work2+n_ambi; diff --git a/src/spherical_line.c b/src/spherical_line.c index 026b5a8..2477dd8 100644 --- a/src/spherical_line.c +++ b/src/spherical_line.c @@ -3,6 +3,9 @@ iem_matrix written by Thomas Musil (c) IEM KUG Graz Austria 2002 - 2006 */ +#include "m_pd.h" +#include "math.h" + /* -------------------------- spherical_line ------------------------------ */ typedef struct _spherical_line |