c98bd710dd6de848d5850a3de9c8d14a11d66b68
[scilab.git] / scilab / modules / linear_algebra / tests / unit_tests / det.dia.ref
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 //==============================   det        ==============================
11 //==========================================================================
12 if execstr('det([1 2;3 4;5 6])','errcatch')==0 then bugmes();quit;end
13 //Small dimension
14 //Real
15 A=[1 1; 1 2];
16 if Err(det(A)-1)>10*%eps then bugmes();quit;end
17 [e,m]=det(A);
18 if e<>0 |Err(m-1)>10*%eps then bugmes();quit;end
19 //Complex
20 A=A+%i;
21 if Err(det(A)-1-%i)>10*%eps then bugmes();quit;end
22 [e,m]=det(A);
23 if e<>0 |Err(m-1-%i)>10*%eps then bugmes();quit;end
24 //Large dimension
25 //Real
26 v=rand(1,21);
27 A=rand(21,21); A=(triu(A,1)+diag(v))*(tril(A,-1)+diag(ones(1,21)));
28 if Err(det(A)-prod(v))>400000*%eps then bugmes();quit;end
29 [e,m]=det(A);
30 if Err(m*(10^e)-prod(v))>400000*%eps then bugmes();quit;end
31 //Complex
32 v=(v+rand(v)*%i)/2;
33 A=rand(21,21); A=(triu(A,1)+diag(v))*(tril(A,-1)+diag(ones(1,21)));
34 if Err(det(A)-prod(v))>10000*%eps then bugmes();quit;end
35 [e,m]=det(A);
36 if Err(m*(10^e)-prod(v))>10000*%eps then bugmes();quit;end