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 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129
|
/* Copyright (C) 2005-2008 Damien Stehle.
Copyright (C) 2007 David Cade.
This file is part of the fplll Library.
The fplll Library is free software; you can redistribute it and/or modify
it under the terms of the GNU Lesser General Public License as published by
the Free Software Foundation; either version 2.1 of the License, or (at your
option) any later version.
The fplll Library is distributed in the hope that it will be useful, but
WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
License for more details.
You should have received a copy of the GNU Lesser General Public License
along with the fplll Library; see the file COPYING. If not, write to
the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
MA 02111-1307, USA. */
#include "util.h"
using namespace fplll;
template <class ZT> int llldiff(ZZ_mat<ZT> &B1, ZZ_mat<ZT> &B2, int c, int r)
{
int test = 1, i, j, sg;
Z_NR<ZT> tmp1;
Z_NR<ZT> tmp2;
for (i = 0; i < r; i++)
{
sg = 1;
tmp1.abs(B1(i, 0));
tmp2.abs(B2(i, 0));
if (tmp1.cmp(tmp2) != 0)
{
// cerr << r << ", 0\n";
test = 0;
}
if (tmp1.cmp(B1(i, 0)) != 0)
sg *= -1;
if (tmp1.cmp(B2(i, 0)) != 0)
sg *= -1;
if (sg == 1)
{
for (j = 1; j < c; j++)
{
if (B1(i, j).cmp(B2(i, j)) != 0)
{
// cerr << i << " " << j << "\n";
test = 0;
}
}
}
else
{
for (j = 1; j < c; j++)
{
tmp1.mul_si(B1(i, j), -1);
if (tmp1.cmp(B2(i, j)) != 0)
{
// cerr << i << " " << j << "\n";
test = 0;
}
}
}
}
return (test);
}
/* ********************** */
/* MAIN **************** */
/* ********************** */
int main(int argc, char **argv)
{
int c, r, ac;
ZZ_mat<mpz_t> mat1, mat2;
for (ac = 1; ac < argc; ac++)
{
if (strcmp(argv[ac], "-c") == 0 || strcmp(argv[ac], "-r") == 0)
{
if (ac < argc - 1)
++ac;
}
else if (strcmp(argv[ac], "--help") == 0)
{
cout << "Usage: cat matrix1 matrix2 | " << argv[0] << endl;
return 0;
}
else if (argv[ac][0] == '-')
{
cerr << "llldiff: invalid option '" << argv[ac] << "'" << endl;
return 1;
}
else
{
break;
}
}
istream *inputStream;
if (argv[ac])
inputStream = new ifstream(argv[ac]);
else
inputStream = &cin;
*inputStream >> mat1;
*inputStream >> mat2;
if (argv[ac])
delete inputStream;
r = mat1.get_rows();
c = mat1.get_cols();
int difference = !llldiff<mpz_t>(mat1, mat2, c, r);
if (difference)
{
cerr << "===INVALID RESULT===" << endl;
}
return difference;
}
|