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

#define N 8

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

void doinit(){
    for(int i=0;i<N;i++) X[i]=1.0/(i+1);
    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);
    cmultAx(N,N,A,X,B);
    printf("B=\n");  cvecprint(N,B);
    cmultAx(N,N,AC,B,BB);
    printf("BB=\n"); cvecprint(N,BB);
    return 0;
}