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

int ISEED[4] = {9, 10, 11, 12};

void test_func(int n)
{
    int info = 0;
    int lda = n;
    double *A = (double*)malloc(lda * n * sizeof(double));
    double *work;

    const int idist = 2;
    const int ione = 1;

    for (int j = 0; j < n; j++)
         dlarnv_(&idist, ISEED, &n, &A[j * n]);

    double res = dlansy_("F", "U", &n, A, &lda, work);

    printf("n = %d nthr = %d res = %lf\n", n, omp_get_max_threads(), res);
}

int main()
{
    test_func(126);
    test_func(127);
    test_func(128);
    test_func(129);
    return 0;
}