15 for (i = 0; i < n; ++i, ++u) {
17 for (j = 0, q = v; j < n; ++j, ++a, ++q) {
18 u->
re += a->re * q->
re - a->im * q->
im;
19 u->
im += a->im * q->
re + a->re * q->
im;
31 for (k = 0; k < n; ++k, ++u, ++v) {
32 z.
re += u->re * v->
re + u->im * v->
im;
33 z.
im += u->re * v->
im - u->im * v->
re;
void cvmul(Cpx *u, Cpx *a, Cpx *v, int n)
Cpx cvnrm(Cpx *u, Cpx *v, int n)