12 Cpx z, *q0, *p, *s, *
t;
16 q0 = (
Cpx *) calloc(n,
sizeof(
Cpx));
17 for (i = 0; i < n; ++i) {
18 for (j = 0,
t =
b; j < n; ++j) {
20 for (k = 0, s = a + i * n; k < n; ++k, ++s, ++
t) {
21 z.
re +=
t->re * s->re +
t->im * s->im;
22 z.
im +=
t->im * s->re -
t->re * s->im;
26 for (j = 0, p = hm + i,
t = a; j <= i; ++j, p += n) {
28 for (k = 0, s = q0; k < n; ++k, ++
t, ++s) {
29 z.
re +=
t->re * s->
re -
t->im * s->
im;
30 z.
im +=
t->im * s->
re +
t->re * s->
im;
void utrnhm(Cpx *hm, Cpx *a, Cpx *b, int n)