blas test for bug #4017

test_lock.c

#include <cblas.h>
#include <omp.h>
#include <stdio.h>

#define DATASIZE 100 

double A[DATASIZE*DATASIZE];
double B[DATASIZE*DATASIZE];
double C[DATASIZE*DATASIZE];

void main()
{
    printf("%d\n", omp_get_max_threads());
    cblas_dgemm(CblasRowMajor, CblasNoTrans, CblasNoTrans, 
                DATASIZE, DATASIZE, DATASIZE, 1, A, DATASIZE, B, DATASIZE, 0, C, DATASIZE);
    printf("%d\n", omp_get_max_threads());
    omp_set_num_threads(4);
    cblas_dgemm(CblasRowMajor, CblasNoTrans, CblasNoTrans, 
                DATASIZE, DATASIZE, DATASIZE, 1, A, DATASIZE, B, DATASIZE, 0, C, DATASIZE);
}