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
|
// Copyright (C) 1995 The New York Group Theory Cooperative
// See magnus/doc/COPYRIGHT for the full notice.
//
// Contents: Informal tests of class RandomMatrix.
//
// Principal Authors: Dmitry Bormotov
//
#include "RandomMatrix.h"
#include "GaussTransformation.h"
int main( )
{
/*
srand(1335);
RandomMatrix<Integer> RM(2);
for( int i = 0; i < 10; i++ ) {
Matrix<Integer> M = RM.getRandomMatrix();
cout << "M = " << M << endl;
Determinant<Integer> D(M);
D.startComputation();
while (!D.done()) D.continueComputation();
cout << "determinant = " << D.getDeterminant() << endl;
}
*/
srand(1335);
RandomMatrix<Rational> RM(2);
for( int i = 0; i < 10; i++ ) {
Matrix<Rational> M = RM.getRandomMatrix();
cout << "M = " << M << endl;
Determinant<Rational> D(M);
D.startComputation();
while (!D.done()) D.continueComputation();
cout << "determinant = " << D.getDeterminant() << endl;
}
}
|
