cblas domino check NaN

Author: bilam

jsrc/vd.c

@@ -216,6 +216,7 @@ static A jtlq(J jt,A w,D *det){A l;D c=inf,d=0,x;I n1,n,*s,wr;
  if(hascblas&&(wr==2)&&(s[0]==s[1])&&(AT(w)&FL)&&(s[0]>1)){
  int info; int m1=s[0]; int *ipiv;
  D worksize; D *work; int nn=-1;
+ NAN0;
  ipiv=MALLOC(m1*sizeof(int));
  w=cant1(w);
  jdgetrf_(&m1,&m1,DAV(w),&m1,ipiv,&info);
@@ -226,11 +227,13 @@ static A jtlq(J jt,A w,D *det){A l;D c=inf,d=0,x;I n1,n,*s,wr;
  jdgetri_(&m1,DAV(w),&m1,ipiv,work,&nn,&info);
  if(info){FREE(ipiv); FREE(work); ASSERT(!info,EVDOMAIN);}
  FREE(ipiv); FREE(work);
+ NAN1;
  *det=0.0;
  RETF(cant1(w));
  }else if(hascblas&&(wr==2)&&(s[0]==s[1])&&(AT(w)&CMPX)&&(s[0]>1)){
  int info; int m1=s[0]; int *ipiv;
  dcomplex worksize; dcomplex *work; int nn=-1;
+ NAN0;
  ipiv=MALLOC(m1*sizeof(int));
  w=cant1(w);
  jzgetrf_(&m1,&m1,(dcomplex*)ZAV(w),&m1,ipiv,&info);
@@ -241,6 +244,7 @@ static A jtlq(J jt,A w,D *det){A l;D c=inf,d=0,x;I n1,n,*s,wr;
  jzgetri_(&m1,(dcomplex*)ZAV(w),&m1,ipiv,work,&nn,&info);
  if(info){FREE(ipiv); FREE(work); ASSERT(!info,EVDOMAIN);}
  FREE(ipiv); FREE(work);
+ NAN1;
  *det=0.0;
  RETF(cant1(w));
 #if 0      // slow

test/g131.ijs

@@ -188,6 +188,8 @@ assert. 1e_8>e
 0 -: %.~$~0x
 0 -: %.~$~0r4
 
+'NaN error' -: %. etx 2 2 $ 0 1 2 _
+'NaN error' -: %. etx 2 2$0 1e_309 10 10
 
 'domain error' -:        %. etx 4 4$'abc'
 'domain error' -:        %. etx 4 4$u:'abc'