* Generally useful stuff.
*/
-#include <gsl/gsl_vector.h>
-#include <gsl/gsl_blas.h>
-
#include "common.h"
double magnD(const double pq[D3]) {
r[1]= a[2]*b[0] - a[0]*b[2];
r[2]= a[0]*b[1] - a[1]*b[0];
}
+
+void diee(const char *what) { perror(what); exit(16); }
+void flushoutput(void) { if (fflush(stdout)||ferror(stdout)) diee("stdout"); }