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