before abandon normalisation effort

[moebius2.git] / common.c

diff --git a/common.c b/common.c
index ba61341dc9bc45c457feba0ba87b112190df3c6a..34381e890494a489c95e0161856dccdc598a6ce6 100644 (file)
--- a/common.c
+++ b/common.c
@@ -43,6 +43,13 @@ void xprod(double r[D3], const double a[D3], const double b[D3]) {
   r[2]= a[0]*b[1] - a[1]*b[0];
 }
 
+double dotprod(const double a[D3], const double b[D3]) {
+  int k;
+  double result= 0;
+  K result += a[k] * b[k];
+  return result;
+}
+
 void libdie(const char *lib, int l, const char *str) {
   fprintf(stderr,"%s library call failed, line %d: %s\n", lib, l, str);
 }