[Bug-gsl] GSL 1.9 test_cblas vs Intel MKL

[Albert Strasheim]

Hello

I just compiled GSL 1.9 and linked against Intel MKL for the BLAS functions. 
All the tests pass, with the exception of some tests in test_cblas. The list 
of failures is included at the end of this message.

Choosing one at random:

FAIL: sscal(case 112) (0 observed vs 0.65100002288818359 expected) [148]

the code is:

int N = 1;
float alpha = 0.0f;
float X[] = { 0.651f };
int incX = -1;
float expected[] = { 0.651f };
cblas_sscal(N, alpha, X, incX);
{
 int i;
 for (i = 0; i < 1; i++) {
   gsl_test_rel(X[i], expected[i], flteps, "sscal(case 112)");
 }
};

GSL's sscal function does this as its first operation:

if (incX <= 0) { return; }

I might be mistaken, but I don't think other BLAS implementations have this 
restriction. Instead incX=-1 just means step backward in memory from &X[0], 
instead of forward.

Any thoughts?

Regards,

Albert

P. S. The complete list of failures when running GSL 1.9 test_cblas with 
Intel MKL:

FAIL: snrm2(case 28) (0.31700000166893005 observed vs 0 expected) [40]
FAIL: dnrm2(case 29) (0.070999999999999994 observed vs 0 expected) [41]
FAIL: scnrm2(case 30) (1.2523837089538574 observed vs 0 expected) [42]
FAIL: dznrm2(case 31) (0.65323579203837268 observed vs 0 expected) [43]
FAIL: snrm2(case 36) (0.92660671472549438 observed vs 0 expected) [48]
FAIL: dnrm2(case 37) (0.62676391089468442 observed vs 0 expected) [49]
FAIL: scnrm2(case 38) (1.1265496015548706 observed vs 0 expected) [50]
FAIL: dznrm2(case 39) (1.2982311042337569 observed vs 0 expected) [51]
FAIL: sasum(case 40) (0.23899999260902405 observed vs 0 expected) [52]
FAIL: dasum(case 41) (0.41299999999999998 observed vs 0 expected) [53]
FAIL: scasum(case 42) (0.11699999868869781 observed vs 0 expected) [54]
FAIL: dzasum(case 43) (0.72999999999999998 observed vs 0 expected) [55]
FAIL: sasum(case 48) (1.0019999742507935 observed vs 0 expected) [60]
FAIL: dasum(case 49) (1.0720000000000001 observed vs 0 expected) [61]
FAIL: scasum(case 50) (2.2739999294281006 observed vs 0 expected) [62]
FAIL: dzasum(case 51) (2.4770000000000003 observed vs 0 expected) [63]
FAIL: sscal(case 112) (0 observed vs 0.65100002288818359 expected) [148]
FAIL: sscal(case 113) (0.065100006759166718 observed vs 0.65100002288818359 
expected) [149]
FAIL: dscal(case 115) (0 observed vs 0.68600000000000005 expected) [151]
FAIL: dscal(case 116) (0.068600000000000008 observed vs 0.68600000000000005 
expected) [152]
FAIL: cscal(case 118) real (0 observed vs 0.98600000143051147 expected) 
[154]
FAIL: cscal(case 118) imag (0 observed vs -0.77499997615814209 expected) 
[155]
FAIL: cscal(case 119) real (0.098600000143051147 observed vs 
0.98600000143051147 expected) [156]
FAIL: cscal(case 119) imag (-0.077500000596046448 observed 
vs -0.77499997615814209 expected) [157]
FAIL: cscal(case 121) real (0.077500000596046448 observed vs 
0.98600000143051147 expected) [160]
FAIL: cscal(case 121) imag (0.098600000143051147 observed 
vs -0.77499997615814209 expected) [161]
FAIL: cscal(case 122) real (0.25360000133514404 observed vs 
0.98600000143051147 expected) [162]
FAIL: cscal(case 122) imag (0.11969999969005585 observed 
vs -0.77499997615814209 expected) [163]
FAIL: cscal(case 123) real (1.218500018119812 observed vs 
0.98600000143051147 expected) [164]
FAIL: cscal(case 123) imag (-0.47919997572898865 observed 
vs -0.77499997615814209 expected) [165]
FAIL: zscal(case 124) real (0 observed vs 0.45400000000000001 expected) 
[166]
FAIL: zscal(case 124) imag (0 observed vs -0.47799999999999998 expected) 
[167]
FAIL: zscal(case 125) real (0.045400000000000003 observed vs 
0.45400000000000001 expected) [168]
FAIL: zscal(case 125) imag (-0.047800000000000002 observed 
vs -0.47799999999999998 expected) [169]
FAIL: zscal(case 127) real (0.047800000000000002 observed vs 
0.45400000000000001 expected) [172]
FAIL: zscal(case 127) imag (0.045400000000000003 observed 
vs -0.47799999999999998 expected) [173]
FAIL: zscal(case 128) real (0.14100000000000001 observed vs 
0.45400000000000001 expected) [174]
FAIL: zscal(case 128) imag (0.043000000000000003 observed 
vs -0.47799999999999998 expected) [175]
FAIL: zscal(case 129) real (0.59740000000000004 observed vs 
0.45400000000000001 expected) [176]
FAIL: zscal(case 129) imag (-0.34179999999999999 observed 
vs -0.47799999999999998 expected) [177]
FAIL: sscal(case 148) (0 observed vs 0.62900000810623169 expected) [238]
FAIL: sscal(case 148) (-0 observed vs -0.41899999976158142 expected) [239]
FAIL: sscal(case 149) (0.06289999932050705 observed vs 0.62900000810623169 
expected) [240]
FAIL: sscal(case 149) (-0.041900001466274261 observed 
vs -0.41899999976158142 expected) [241]
FAIL: dscal(case 151) (0 observed vs 0.39800000000000002 expected) [244]
FAIL: dscal(case 151) (-0 observed vs -0.65600000000000003 expected) [245]
FAIL: dscal(case 152) (0.039800000000000002 observed vs 0.39800000000000002 
expected) [246]
FAIL: dscal(case 152) (-0.065600000000000006 observed 
vs -0.65600000000000003 expected) [247]
FAIL: cscal(case 154) real (0 observed vs 0.73600000143051147 expected) 
[250]
FAIL: cscal(case 154) imag (0 observed vs 0.33100000023841858 expected) 
[251]
FAIL: cscal(case 154) real (-0 observed vs -0.31799998879432678 expected) 
[252]
FAIL: cscal(case 154) imag (0 observed vs 0.62199997901916504 expected) 
[253]
FAIL: cscal(case 155) real (0.073600001633167267 observed vs 
0.73600000143051147 expected) [254]
FAIL: cscal(case 155) imag (0.033100001513957977 observed vs 
0.33100000023841858 expected) [255]
FAIL: cscal(case 155) real (-0.031799998134374619 observed 
vs -0.31799998879432678 expected) [256]
FAIL: cscal(case 155) imag (0.062199998646974564 observed vs 
0.62199997901916504 expected) [257]
FAIL: cscal(case 157) real (-0.033100001513957977 observed vs 
0.73600000143051147 expected) [262]
FAIL: cscal(case 157) imag (0.073600001633167267 observed vs 
0.33100000023841858 expected) [263]
FAIL: cscal(case 157) real (-0.062199998646974564 observed 
vs -0.31799998879432678 expected) [264]
FAIL: cscal(case 157) imag (-0.031799998134374619 observed vs 
0.62199997901916504 expected) [265]
FAIL: cscal(case 158) real (0.0073999986052513123 observed vs 
0.73600000143051147 expected) [266]
FAIL: cscal(case 158) imag (0.18029999732971191 observed vs 
0.33100000023841858 expected) [267]
FAIL: cscal(case 158) real (-0.15619999170303345 observed 
vs -0.31799998879432678 expected) [268]
FAIL: cscal(case 158) imag (-0.0013999976217746735 observed vs 
0.62199997901916504 expected) [269]
FAIL: cscal(case 159) real (0.63669997453689575 observed vs 
0.73600000143051147 expected) [270]
FAIL: cscal(case 159) imag (0.55180001258850098 observed vs 
0.33100000023841858 expected) [271]
FAIL: cscal(case 159) real (-0.50459998846054077 observed 
vs -0.31799998879432678 expected) [272]
FAIL: cscal(case 159) imag (0.52660000324249268 observed vs 
0.62199997901916504 expected) [273]
FAIL: zscal(case 160) real (0 observed vs 0.52100000000000002 expected) 
[274]
FAIL: zscal(case 160) imag (0 observed vs -0.81100000000000005 expected) 
[275]
FAIL: zscal(case 160) real (0 observed vs 0.55600000000000005 expected) 
[276]
FAIL: zscal(case 160) imag (0 observed vs -0.14699999999999999 expected) 
[277]
FAIL: zscal(case 161) real (0.052100000000000007 observed vs 
0.52100000000000002 expected) [278]
FAIL: zscal(case 161) imag (-0.081100000000000005 observed 
vs -0.81100000000000005 expected) [279]
FAIL: zscal(case 161) real (0.055600000000000011 observed vs 
0.55600000000000005 expected) [280]
FAIL: zscal(case 161) imag (-0.0147 observed vs -0.14699999999999999 
expected) [281]
FAIL: zscal(case 163) real (0.081100000000000005 observed vs 
0.52100000000000002 expected) [286]
FAIL: zscal(case 163) imag (0.052100000000000007 observed 
vs -0.81100000000000005 expected) [287]
FAIL: zscal(case 163) real (0.0147 observed vs 0.55600000000000005 expected) 
[288]
FAIL: zscal(case 163) imag (0.055600000000000011 observed 
vs -0.14699999999999999 expected) [289]
FAIL: zscal(case 164) real (0.21430000000000002 observed vs 
0.52100000000000002 expected) [290]
FAIL: zscal(case 164) imag (0.023100000000000009 observed 
vs -0.81100000000000005 expected) [291]
FAIL: zscal(case 164) real (0.085000000000000006 observed vs 
0.55600000000000005 expected) [292]
FAIL: zscal(case 164) imag (0.096500000000000016 observed 
vs -0.14699999999999999 expected) [293]
FAIL: zscal(case 165) real (0.76429999999999998 observed vs 
0.52100000000000002 expected) [294]
FAIL: zscal(case 165) imag (-0.65470000000000006 observed 
vs -0.81100000000000005 expected) [295]
FAIL: zscal(case 165) real (0.60010000000000008 observed vs 
0.55600000000000005 expected) [296]
FAIL: zscal(case 165) imag (0.019800000000000012 observed 
vs -0.14699999999999999 expected) [297]