1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
import TestLib;
import math;
StartTest("determinant");
assert(determinant(new real[][] {{0}}) == 0);
assert(determinant(new real[][] {{1}}) == 1);
assert(determinant(new real[][] {{1,2},{3,4}}) == -2);
real e=1e-20;
assert(close(determinant(new real[][] {{1e,2e},{3e,4e}}),-2e-40));
assert(close(determinant(new real[][] {{1,2,3},{4,5,6},{7,8,9}}),0));
assert(close(determinant(new real[][] {{1,2},{1,2}}),0));
assert(close(determinant(new real[][] {{1,2,3,4},
{5,6,7,8},{9,10,11,12},{13,14,15,16}}),0));
assert(close(determinant(new real[][] {{1,2,3,4},
{5,0,7,8},{9,10,0,12},{13,14,15,16}}),-2376));
assert(close(determinant(new real[][]{{1,-2,3,0},{4,-5,6,2},{-7,-8,10,5},
{1,50,1,-2}}),-4588));
EndTest();
|
