complex mul(complex alpha, complex beta)
{
complex gamma;

gamma.re = alpha.re*beta.re - alpha.im*beta.im;
gamma.im = alpha.im*beta.re + alpha.re*beta.im;

return gamma;
}
これで実行すると、mat_mul呼び出す所でとまってしまいます。