#include <stdio.h>

typedef struct {
double re;
double im;
} complex;

complex add(complex alpha, complex beta);
complex mul(complex alpha, complex beta);
void mat_mul(complex** a, complex** b, complex** c);
void print_mat(complex** a, char *str);

int main(void)
{
int n,m;
n = 2; m = 3;
a = allocate(n, m);
for (i = 0, i < n; i++) {
for (j = 0, j < m; j++) {
a[i][j].re = 1.0 * (i * n + j);
a[i][j].im = 1.0 * (i * n + j);
}
}

b = allocate(n, m);
for (i = 0, i < n; i++) {
for (j = 0, j < m; j++) {
b[i][j].re = 1.0 * (i * n + j);
b[i][j].im = 1.0 * (i * n + j);
}
}

complex c[2][2];

print_mat(a, "a");
print_mat(b, "b");

mat_mul(a, b, c);

print_mat(c, "a * b");
return 0;
}