c8f9c397bee8c2c9404af1037063899c7763495f
[scilab.git] / scilab / modules / linear_algebra / tests / unit_tests / det.tst
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) ????-2008 - INRIA Michael Baudin
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7 function r=Err(x),r=norm(x,1),endfunction
8 rand('normal')
9
10 //==========================================================================
11 //==============================   det        ============================== 
12 //==========================================================================
13 if execstr('det([1 2;3 4;5 6])','errcatch')==0 then pause,end
14 //Small dimension 
15 //Real
16 A=[1 1; 1 2];
17 if Err(det(A)-1)>10*%eps then pause,end
18 [e,m]=det(A);
19 if e<>0 |Err(m-1)>10*%eps then pause,end
20 //Complex
21 A=A+%i;
22 if Err(det(A)-1-%i)>10*%eps then pause,end
23 [e,m]=det(A);
24 if e<>0 |Err(m-1-%i)>10*%eps then pause,end
25 //Large dimension
26 //Real
27 v=rand(1,21);
28 A=rand(21,21); A=(triu(A,1)+diag(v))*(tril(A,-1)+diag(ones(1,21)));
29 if Err(det(A)-prod(v))>400000*%eps then pause,end
30 [e,m]=det(A);
31 if Err(m*(10^e)-prod(v))>400000*%eps then pause,end
32 //Complex
33 v=(v+rand(v)*%i)/2;
34 A=rand(21,21); A=(triu(A,1)+diag(v))*(tril(A,-1)+diag(ones(1,21)));
35 if Err(det(A)-prod(v))>10000*%eps then pause,end
36 [e,m]=det(A);
37 if Err(m*(10^e)-prod(v))>10000*%eps then pause,end
38