From: Ian Jackson <ian@davenant.relativity.greenend.org.uk>
Date: Mon, 29 Sep 2008 00:33:57 +0000 (+0100)
Subject: do not use gsl_blas_dnrm2; instead just do addition and sqrt
X-Git-Url: https://www.chiark.greenend.org.uk/ucgi/~ian/git?a=commitdiff_plain;h=fb621741c0f6d16a6c465585bf423452aad069ab;p=moebius2.git

do not use gsl_blas_dnrm2; instead just do addition and sqrt
---

diff --git a/common.c b/common.c
index 3b036aa..bfe11b7 100644
--- a/common.c
+++ b/common.c
@@ -5,14 +5,7 @@
 #include "common.h"
 
 double magnD(const double pq[D3]) {
-  gsl_vector v;
-
-  v.size= D3;
-  v.stride= 1;
-  v.data= (double*)pq;
-  /* owner and block ought not to be used */
-
-  return gsl_blas_dnrm2(&v);
+  return sqrt(magnD2(pq));
 }
 
 double magnD2(const double pq[D3]) {