aboutsummaryrefslogtreecommitdiff
path: root/src
diff options
context:
space:
mode:
Diffstat (limited to 'src')
-rw-r--r--src/matrix_orthogonal.c1
-rw-r--r--src/matrix_pinv.c4
-rw-r--r--src/spherical_line.c3
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