static inline Fixpt dbl2fixpt(double x) { return x * int2fixpt(1); }
static inline double fixpt2dbl(Fixpt x) { return x / (1.0*int2fixpt(1)); }
static inline Fixpt fixpt_mul(Fixpt a, Fixpt b) {
static inline Fixpt dbl2fixpt(double x) { return x * int2fixpt(1); }
static inline double fixpt2dbl(Fixpt x) { return x / (1.0*int2fixpt(1)); }
static inline Fixpt fixpt_mul(Fixpt a, Fixpt b) {