#ifdef _WIN32
#include <windows.h>
#else
#include <sys/time.h>
#include <linux/limits.h>
#endif

#include <malloc.h>
#include <stdio.h>
#include <stdint.h>
#include <cblas.h>
#define ALIGNMET_SIZE 16

#ifdef __MINGW32__
#define _aligned_malloc __mingw_aligned_malloc
#define _aligned_free  __mingw_aligned_free
#endif //MINGW

typedef struct Matrix {
   float *matrix;
   uint32_t rows;
   uint32_t columns;

} Matrix;

Matrix makeInputMatrix(int rows, int cols) {
    Matrix matrix;
    matrix.rows = rows;
    matrix.columns = cols;
    int64_t expectedSize = (matrix.rows*matrix.columns * sizeof(float) + ALIGNMET_SIZE)&(~(ALIGNMET_SIZE - 1));
    //matrix.matrix = (float*)blas_memory_alloc(expectedSize);
    #if defined (_WIN32)
        matrix.matrix = (float*)_aligned_malloc(expectedSize, ALIGNMET_SIZE);
    #else
        matrix.matrix = (float*)memalign(ALIGNMET_SIZE,  expectedSize);
    #endif
    //srandom(getpid());
    for(int i = 0; i < rows; i++) {
        for(int j = 0; j < cols; j++) {
            matrix.matrix[i*cols + j] = abs((float)rand() / (float)RAND_MAX - 0.5f);
        }
    }
    return matrix;
}

Matrix multMatrixes(Matrix &a, Matrix &b, CBLAS_TRANSPOSE transA, CBLAS_TRANSPOSE transB, float alpha, float betta) {
    Matrix c  ;
    c.rows = a.rows;
    c.columns = b.columns;
    int64_t expectedSize = (c.rows * c.columns * sizeof(float) + ALIGNMET_SIZE)&(~(ALIGNMET_SIZE - 1));
    //c.matrix = (float*)blas_memory_alloc(expectedSize);
    #if defined (_WIN32)
        c.matrix = (float*)_aligned_malloc(expectedSize, ALIGNMET_SIZE);
    #else
        c.matrix = (float*)memalign(ALIGNMET_SIZE, expectedSize);
    #endif
    fprintf(stderr, "Start address for matrix C [%p]\n", c.matrix);
    cblas_sgemm(CblasRowMajor, transA, transB, a.rows, b.columns, a.columns, 
            alpha, a.matrix, a.columns, b.matrix,a.rows, betta, c.matrix, c.rows);
    return c;
}

int main(){ 
  //test
    Matrix A=makeInputMatrix(2048,1320);
    Matrix B=makeInputMatrix(1320,64);
    //segFault for bigger sizes 
    Matrix C=multMatrixes(A,B,CblasNoTrans,CblasNoTrans,1.0f,0);

    //...freeing mem  
    ///...///
    return 0;
}