#define DBL double

typedef DBL VECTOR    [3];
typedef DBL VECTOR_4D [4];