From fb621741c0f6d16a6c465585bf423452aad069ab Mon Sep 17 00:00:00 2001
From: Ian Jackson <ian@davenant.relativity.greenend.org.uk>
Date: Mon, 29 Sep 2008 01:33:57 +0100
Subject: [PATCH] do not use gsl_blas_dnrm2; instead just do addition and sqrt

---
 common.c | 9 +--------
 1 file changed, 1 insertion(+), 8 deletions(-)

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]) {
-- 
2.30.2