compiles

[moebius2.git] / common.c

diff --git a/common.c b/common.c
index 47edc4ddbe0f6be0e5c275aa8192f190ecaf1967..900102513cf2c7305b31b85deb223a4c9b36f872 100644 (file)
--- a/common.c
+++ b/common.c
@@ -2,6 +2,9 @@
  * Generally useful stuff.
  */
 
+#include <gsl/gsl_vector.h>
+#include <gsl/gsl_blas.h>
+
 #include "common.h"
 
 double magnD(const double pq[D3]) {
@@ -9,10 +12,10 @@ double magnD(const double pq[D3]) {
 
   v.size= D3;
   v.stride= 1;
-  v.vector.data= pq;
+  v.data= (double*)pq;
   /* owner and block ought not to be used */
 
-  return gsl_blas_snrm2(&v);
+  return gsl_blas_dnrm2(&v);
 }
 
 double hypotD(const double p[D3], const double q[D3]) {
@@ -20,22 +23,25 @@ double hypotD(const double p[D3], const double q[D3]) {
   double pq[D3];
 
   K pq[k]= p[k] - q[k];
-  return hypotD1(pq);
+  return magnD(pq);
 }
 
 double hypotD2(const double p[D3], const double q[D3]) {
+  int k;
   double d2= 0;
+  
   K d2= ffsqa(p[k] - q[k], d2);
   return d2;
 }
 
 double hypotD2plus(const double p[D3], const double q[D3], double d2) {
+  int k;
   K d2= ffsqa(p[k] - q[k], d2);
   return d2;
 }
 
 void xprod(double r[D3], const double a[D3], const double b[D3]) {
-  r[0]= a[1]*b[2] - a[2]*b[1]);
-  r[1]= a[2]*b[0] - a[0]*b[2]);
-  r[2]= a[0]*b[1] - a[1]*b[0]);
+  r[0]= a[1]*b[2] - a[2]*b[1];
+  r[1]= a[2]*b[0] - a[0]*b[2];
+  r[2]= a[0]*b[1] - a[1]*b[0];
 }