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