#include <stdio.h>
#include <strings.h>
#include <complex.h>
#include <math.h>
#include "matrixlib.h"
//#include "tictoc.h"

#define N 16

complex A[N][N],AC[N][N];
complex X[N];

void doinit(){
    for(int i=0;i<N;i++) X[i]=i;
    for(int i=0;i<N;i++) for(int j=0;j<N;j++) {
        A[i][j]=cexp(-I*2*M_PI*i*j/N);
        AC[i][j]=cexp(I*2*M_PI*i*j/N)/N;
    }
}

int main(){
    complex B[N],BB[N];
    doinit();
    printf("A=\n");  cmatprint(N,N,A);
    printf("X=\n");  cvecprint(N,X);
//	tic();
    cmultAx(N,N,A,X,B);
//	double t=toc();
//	printf("The elapsed time was %g seconds.\n",t);
    printf("B=\n");  cvecprint(N,B);
//    cmultAx(N,N,AC,B,BB);
//    printf("BB=\n"); cvecprint(N,BB);
    return 0;
}