* Generally useful stuff.
*/
-#include <gsl/gsl_vector.h>
-#include <gsl/gsl_blas.h>
-
#include "common.h"
double magnD(const double pq[D3]) {
return gsl_blas_dnrm2(&v);
}
+double magnD2(const double pq[D3]) {
+ int k;
+ double d2= 0;
+
+ K d2= ffsqa(pq[k], d2);
+ return d2;
+}
+
double hypotD(const double p[D3], const double q[D3]) {
int k;
double pq[D3];
r[2]= a[0]*b[1] - a[1]*b[0];
}
+void xprod_norm(double r[D3], const double a[D3], const double b[D3],
+ double absepsilon, double one) {
+ int k;
+ xprod(r,a,b);
+ double multby= one/(magnD(r) + absepsilon);
+ K r[k] *= multby;
+}
+
+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);
+}
+
+void gsldie(int l, const char *what, int status) {
+ fprintf(stderr,"gsl function failed, line %d: %s: %s\n",
+ l, what, gsl_strerror(status));
+ exit(-1);
+}
+
void diee(const char *what) { perror(what); exit(16); }
+void fail(const char *emsg) { fputs(emsg,stderr); exit(12); }
void flushoutput(void) { if (fflush(stdout)||ferror(stdout)) diee("stdout"); }