presentation04
fc_amat.init()
A=fc_amat.random.randnsympd(10^5,4);
D=det(A);
[L,U,P]=lu(A);
E1=abs(L*U-P*A);
fprintf('max of E1 elements: %.6e\n',max(E1(:)))
R=chol(A);
E2=abs(R'*R-A);
fprintf('max of E2 elements: %.6e\n',max(E2(:)))
b=ones(4,1);
X=A\b;
E3=abs(A*X-b);
fprintf('max of E3 elements: %.6e\n',max(E3(:)))
B=fc_amat.random.randn(10^5,4,1);
Y=A\B;
E4=abs(A*Y-B);
fprintf('max of E4 elements: %.6e\n',max(E4(:)))
whos
Using fc_amat[0.0.3] with fc_tools[0.0.25], fc_bench[0.0.6].
max of E1 elements: 1.421085e-14
max of E2 elements: 1.421085e-14
max of E3 elements: 4.547474e-13
max of E4 elements: 9.525714e-14
Variables in the current scope:
Attr Name Size Bytes Class
==== ==== ==== ===== =====
A 1x1 0 amat
B 1x1 0 amat
D 1x1 0 amat
E1 1x1 0 amat
E2 1x1 0 amat
E3 1x1 0 amat
E4 1x1 0 amat
L 1x1 0 amat
P 1x1 0 amat
R 1x1 0 amat
U 1x1 0 amat
X 1x1 0 amat
Y 1x1 0 amat
f __code__ 1x864 864 char
a argn 1x8 8 char
b 4x1 32 double
Total is 889 elements using 904 bytes