X-Git-Url: http://www.chiark.greenend.org.uk/ucgi/~ian/git?p=moebius2.git;a=blobdiff_plain;f=common.h;h=6414303ad755c8234366e5990e450441f2b2dbd7;hp=7e0e86d7622da0698d992bdf94af0dfef3e4676e;hb=7f93208fa022d3c9f925f48e5469793bc3612b61;hpb=f1dc8f4f3623b15deddba8db85d181a0e1011c9e diff --git a/common.h b/common.h index 7e0e86d..6414303 100644 --- a/common.h +++ b/common.h @@ -8,6 +8,12 @@ #ifndef _GNU_SOURCE #define _GNU_SOURCE #endif +#ifndef _XOPEN_SOURCE +#define _XOPEN_SOURCE 1000 +#endif + +#include +#include #include #include @@ -18,11 +24,13 @@ #include #include #include +#include #include #include #include #include +#include #define D3 3 @@ -31,11 +39,16 @@ double hypotD2(const double p[D3], const double q[D3]); double hypotD2plus(const double p[D3], const double q[D3], double add); double magnD(const double pq[D3]); +double magnD2(const double pq[D3]); void xprod(double r[D3], const double a[D3], const double b[D3]); double dotprod(const double a[D3], const double b[D3]); +void xprod_norm(double r[D3], const double a[D3], const double b[D3], + double one, double absepsilon); +void normalise(double v[D3], double one, double absepsilon); void flushoutput(void); void diee(const char *what); +void fail(const char *emsg); void libdie(const char *lib, int l, const char *str); #define XA(w) ((w) ? (void)0 : libdie("X", __LINE__, #w))