static inline flo vek_dot (Vek* x, Vek* y) {
int i;
flo d = 0.0;
for (i=0; i<ND; i++)
d += elt_at(x, i) * elt_at(y, i);
return d;
}
static inline Vek* vek_sub (Vek* d, Vek* s) {
int i;
for (i=0; i<ND; i++)
*elt_ati(d, i) -= elt_at(s, i);
return d;
}
static inline Vek* vek_cpy (Vek* x, Vek* y) {
int i;
for (i=0; i<ND; i++)
*elt_ati(x, i) = elt_at(y, i);
return x;
}
static void vek_out (Vek* x) {
int i;
for (i=0; i<ND; i++)
post(" %.2f", elt_at(x, i));
}
always_inline heap_elt_t * first (heap_header_t * h)
{ return elt_at (h, h->head); }
always_inline heap_elt_t * last (heap_header_t * h)
{ return elt_at (h, h->tail); }
