r[2]= a[0]*b[1] - a[1]*b[0];
}
+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 flushoutput(void) { if (fflush(stdout)||ferror(stdout)) diee("stdout"); }