<![CDATA[
#include <numeric.h>
int main() {
double complex a[2][3] = {complex(1,2),2,3,
4,5,6};
array double complex b[3][4] = {1,2,3,4,
5,6,7,8,
1,2,3,4};
array double complex cproductv1[2], cproductv2[3];
cproduct(a,cproductv1);
printf("cproduct(a, cproductv1) = %f\n", cproductv1);
cproduct(b,cproductv2);
printf("cproduct(b,cproductv2) = %f\n", cproductv2);
}
]]>
<![CDATA[
cproduct(a, cproductv1) = complex(6.000000,12.000000) complex(120.000000,0.000000)
cproduct(b,cproductv2) = complex(24.000000,0.000000) complex(1680.000000,0.000000) complex(24.000000,0.000000)
]]>