sysrap/tests/njuffa_erfcinvf_test.cc

0001 #include "njuffa_erfcinvf.h"
0002 #include "NP.hh"
0003
0004 int main()
0005 {
0006     const int ni = 1000 ;
0007     NP* a = NP::Make<float>(ni, 2);
0008     float* aa = a->values<float>();
0009
0010     float SQRT2 = sqrtf(2.f) ;
0011
0012     for(int i=0 ; i < ni ; i++ )
0013     {
0014         float u = 2.f*float(i)/float(ni-1) ;
0015         float v = -SQRT2*njuffa_erfcinvf(u);
0016         aa[2*i+0] = u ;
0017         aa[2*i+1] = v ;
0018     }
0019     a->save("$FOLD/njuffa_erfcinvf_test.npy");
0020
0021     return 0 ;
0022 }