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
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
|
// =============================================================================
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) ????-2008 - INRIA Michael Baudin
//
// This file is distributed under the same license as the Scilab package.
// =============================================================================
//define tools
function A=testmat1(a,n)
//eigen values are given by a dilation of nth roots of 1
A=diag(a*ones(1,n-1),1)+diag((1/a)*ones(1,n-1),-1)
A(1,n)=1/a;A(n,1)=a
endfunction
function r=Err(x)
r=norm(x,1)
endfunction
rand('normal')
//==========================================================================
//============================== inv ==============================
//==========================================================================
//Empty matrix
A=[];
if inv(A)<>[] then bugmes();quit;end
//Singular matrix
if execstr('inv([0 0;2 3])','errcatch')==0 then bugmes();quit;end
if execstr('inv([0 0;%i 3])','errcatch')==0 then bugmes();quit;end
//Rectangular matrix
if execstr('inv(rand(2,3))','errcatch')==0 then bugmes();quit;end
if execstr('inv(rand(2,3)+%i*eye())','errcatch')==0 then bugmes();quit;end
//Small dimension
//---------------
//Unsymetric
A=testmat1(3,5);Ac=testmat1(3+%i,5);
//Real Case
if Err(A*inv(A)-eye(A)) >200*%eps then bugmes();quit;end
//Complex Case
if Err(Ac*inv(Ac)-eye(A)) >200*%eps then bugmes();quit;end
//Symetric
A=A*A';Ac=Ac*Ac';
//Real Case
if Err(A*inv(A)-eye(A)) >1000*%eps then bugmes();quit;end
//Complex Case
if Err(Ac*inv(Ac)-eye(A)) >1000*%eps then bugmes();quit;end
//Large dimension
//---------------
//Unsymetric
A=testmat1(3,50);Ac=testmat1(3+%i,50);
//Real Case
if Err(A*inv(A)-eye(A)) >1000*%eps then bugmes();quit;end
//Complex Case
if Err(Ac*inv(Ac)-eye(A)) >2000*%eps then bugmes();quit;end
//Symetric
A=A*A';Ac=Ac*Ac';
//Real Case
if Err(A*inv(A)-eye(A)) >1.d-10 then bugmes();quit;end
//Complex Case
if Err(Ac*inv(Ac)-eye(A)) >4.d-10 then bugmes();quit;end
|
