Test dgetrf. gcc -o test test_dgetrf.c -I OpenBLAS_install/include OpenBLAS_install/lib/libopenblas.a -lpthread -lgfortran

test_dgetrf.c

#include <stdio.h>
#include "cblas.h"
#include "lapacke.h"

int main(int argc, char* argv[])
{
  int* ipiv;
  int info;
  int i, j;
  int n;
  double * m, *x, *y;

   int LDB,LDA, N, NRHS;
 	char transp = 'N';

  if(argc<2){
    printf("Error\n");
    return 1;
  }
  n=atoi(argv[1]);
  printf("n=%d\n",n);
  
  m=(double*)malloc(sizeof(double)*n*n);
  x=(double*)malloc(sizeof(double)*n);
  y=(double*)malloc(sizeof(double)*n);
  ipiv=(int*)malloc(sizeof(int)*n);

  for (i=0; i<n; ++i) {
    x[i]=1.0;
    for (j=0; j<n; ++j) {
      m[i*n+j]=(rand()%100+1)/10.0;
      //      printf("m[%d,%d]=%lf\n",i,j, m[i*n+j]); 
    }
  }

  /* test cblas.h */
  cblas_dgemv(CblasColMajor, CblasNoTrans, n, n, 1.0, m, n,
	      x, 1, 0.0, y, 1);

  //  for (i=0; i<n; ++i)  printf("x[%d]=%lf\n",i, x[i]); 
  //for (i=0; i<n; ++i)  printf("y[%d]=%lf\n",i, y[i]); 

	LDB=n;
	LDA=n;
	N=n;
	NRHS=1;
	info=0;
	
	LAPACK_dgetrf(&N, &N, m, &LDA,ipiv, &info);
	
if (info != 0) fprintf(stderr, "dgetrf failure with error %d\n", info);

  LAPACK_dgetrs(&transp, &N, &NRHS, m, &LDA, ipiv, y, &LDB, &info);
  
  if (info != 0) fprintf(stderr, "failure with error %d\n", info);
  //  for (i=0; i<n; ++i) printf("%lf\n", y[i]);
  
  return 0;
}

Hi @vtjnash ,

I think you use Mingwin32 compiler. Could you try Mingw64?

Xianyi