int main() {
array double complex a[2][2] = {complex(2,3), -4,
3, -7};
array double complex a2[3][2] = {complex(2,3), -4,
3, -7,
1, 1};
array double b[2][2] = {2, -4,
3, -7};
double complex t;
t = ctrace(a);
printf("ctrace(a) = %f\n", t);
t = ctrace(a2);
printf("ctrace(a2) = %f\n", t);
t = ctrace(b);
printf("ctrace(b) = %f\n", t);
}
ctrace(a) = complex(-5.000000,3.000000)
ctrace(a2) = complex(-5.000000,3.000000)
ctrace(b) = complex(-5.000000,0.000000)