1function test19
2%TEST19 look for NaN's from lchol (caused by Intel MKL 7.x bug)
3% Example:
4%   test19
5% See also cholmod_test
6
7% Copyright 2007, Timothy A. Davis, http://www.suitesparse.com
8
9fprintf ('=================================================================\n');
10fprintf ('test19: look for NaN''s from lchol (caused by Intel MKL 7.x bug)\n') ;
11
12Prob = ssget (936)							    %#ok
13A = Prob.A ;
14[p count] = analyze (A) ;
15A = A (p,p) ;
16tic
17L = lchol (A) ;
18t = toc ;
19fl = sum (count.^2) ;
20fprintf ('mflop rate: %8.2f\n', 1e-6*fl/t) ;
21n = size (L,1) ;
22for k = 1:n
23    if (any (isnan (L (:,k))))
24	k								    %#ok
25	error ('!') ;
26    end
27end
28
29fprintf ('test19 passed; you have a NaN-free BLAS (must not be MKL 7.x...)\n') ;
30