#include <math.h>
#include <stdio.h>
#include "ctest2.h"

#define N 2
#define N4 8
/* 4 x N matrix */
double A[N4]  ;
/* N x 1 matrix */
double B[N]  ;
double C[4] = {0.0,0.0,0.0,0.0} ;


void iprod(void)
{
    int k ;
    double sum ;
    for (k = 0 ; k < N4 ; k++) {
	A[k] = 1.0 ;
    }
    for (k = 0 ; k < N ; k++) {
	B[k] = 1.0 ;
    }
    cblas_dgemm(CblasColMajor, CblasNoTrans, CblasNoTrans, 4, 1, N,  1.0, A, 4, B, N, 0.0, C, 4) ; 
    for (k = 0 ; k < 4 ; k++) {
    printf("C[%d] = %g\n", k, C[k]) ;
    }
}