example matrices that were chosen for Smith form testing.
example matrices that were chosen for Smith form testing.
- Author
- bds & zw
Various Smith form algorithms may be used for matrices over the integers or over Z_m. Moduli greater than 2^32 are not supported here. Several types of example matrices may be constructed or the matrix be read from a file. Run the program with no arguments for a synopsis of the command line parameters.
For the "adaptive" method, the matrix must be over the integers. This is expected to work best for large matrices.
For the "2local" method, the computation is done mod 2^32.
For the "local" method, the modulus must be a prime power.
For the "ilio" method, the modulus may be arbitrary composite. If the modulus is a multiple of the integer determinant, the integer Smith form is obtained. Determinant plus ilio may be best for smaller matrices.
This example was used during the design process of the adaptive algorithm.
#include <iostream>
#include <string>
#include "matrices.h"
template <class PIR>
int main(int argc, char* argv[])
{
if (argc < 3 or argc > 4) {
cout << "usage: " << argv[0] << " type n [filename]" << endl;
cout << " type = `random', `random-rough', `tref',"
<< " 'moler', 'redheffer', "
<< " or `fib',"
<< " and n is the dimension" << endl;
cout << " If filename is present, matrix is written there, else to cout." << endl;
return 0;
}
string type = argv[1];
int n = atoi(argv[2]);
typedef Givaro::ZRing<Givaro::Integer> PIR;
PIR R;
if (M.rowdim() <= 20 && M.coldim() <= 20) {
M.write(std::clog, LinBox::Tag::FileFormat::Maple) << std::endl;
}
if (argc == 4) {
ofstream out(argv[3]);
M.write(out) << endl;
} else {
M.write(cout) << endl;
}
}
template <class PIR>
string src) {
if (src == "random-rough") RandomRoughMat(M, R, n);
else if (src == "random") RandomFromDiagMat(M, R, n);
else if (src == "fib") RandomFibMat(M, R, n);
else if (src == "tref") TrefMat(M, R, n);
else if (src == "krat") KratMat(M, R, n);
else if (src == "moler") MolerMat(M, R, n);
else if (src == "redheffer") RedhefferMat(M, R, n);
else {
}
}
void Mat(DenseMatrix< PIR > &M, PIR &R, int n, string src, string file, string format)
Output matrix is determined by src which may be: "random-rough" This mat will have s,...
Definition: blackbox/smith.C:230
Dense matrix representation.
Definition: blas-matrix.h:62
size_t rowdim() const
Get the number of rows in the matrix.
Definition: blas-matrix.h:233
std::istream & read(std::istream &file)
Read the matrix from an input stream.
Definition: blas-matrix.inl:290
linbox base configuration file
LinBox timer is Givaro's.