#include <numeric.h> int main() { complex double a[3] = {complex(1,2), complex(2,3)}; complex double b[2][3] = {1,2,3,4,5,6}; int b1[2][3] = {1,2,3,4,5,6}; double c[2][3][5]; c[0][0][0] = 10; c[0][0][1] = 20; double complex sums; sums = csum(a); printf("csum(a) = %f\n", sums); sums = csum(b); printf("csum(b) = %f\n", sums); sums = csum(b1); printf("csum(b) = %f\n", sums); sums = csum(c); printf("csum(c) = %f\n", sums); }
csum(a) = complex(3.000000,5.000000) csum(b) = complex(21.000000,0.000000) csum(b) = complex(21.000000,0.000000) csum(c) = complex(30.000000,0.000000)