#include <stdio.h>
#include <strings.h>
#include <complex.h>
//#include <cilk/cilk.h>
#include "matrixlib.h"
#define cilk_for for

void matprint(int m,int n,double A[m][n]){
    for(int i=0;i<m;i++){
        for(int j=0;j<n;j++){
            printf("%g ",A[i][j]);
        }
        printf("\n");
    }
}
void cmatprint(int m,int n,complex A[m][n]){
    for(int i=0;i<m;i++){
        for(int j=0;j<n;j++){
            printf("(%g %g) ",A[i][j]);
        }
        printf("\n");
    }
}
void vecprint(int m,double X[m]){
    matprint(m,1,(double (*)[1])X);
}
void cvecprint(int m,complex X[m]){
    cmatprint(m,1,(complex (*)[1])X);
}
// Algorithm (a) from page 13 of Numerical Algorithms by Justin Solomon
void multAx(int m,int n,double A[m][n],
    double X[n],double B[m]){
    bzero(B,m*sizeof(double));
    cilk_for(int i=0;i<m;i++){
        for(int j=0;j<n;j++){
            B[i]=B[i]+A[i][j]*X[j];
        }
    }
}
void cmultAx(int m,int n,complex A[m][n],
    complex X[n],complex B[m]){
    bzero(B,m*sizeof(complex));
    cilk_for(int i=0;i<m;i++){
        for(int j=0;j<n;j++){
            B[i]=B[i]+A[i][j]*X[j];
        }
    }
}