rrspace.cc 2.17 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
/* Common header file for the rrspace software

This file is part of the rrspace project.

This 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.

This 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 this library; if not, write to the Free Software Foundation, Inc., 51
Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA

*/

#include <iostream>
#include <sstream>
#include "algos.h"
24
#include <unistd.h>
25 26 27 28

using namespace NTL;
using namespace std;

29
int main(int argc, char** argv) {
30 31 32 33 34 35
  long p;
  BivPol eq_curve;
  cin >> p;
  ZZ_p::init(ZZ(p));
  cin >> eq_curve;
  Curve C(&eq_curve);
SPAENLEHAUER Pierre-Jean's avatar
SPAENLEHAUER Pierre-Jean committed
36

37 38 39
  // Singular divisor for nodal curves
  EffectiveDivisor E(C);
  cin >> E;
SPAENLEHAUER Pierre-Jean's avatar
SPAENLEHAUER Pierre-Jean committed
40

41 42 43
  Divisor D(C);
  cin >> D;

44
  RRspace basisRR = RiemannRochBasis(D, E);
45 46 47 48 49 50 51 52 53 54 55 56 57

  bool MAGMACHECK_FLAG = false;

  int c;
  while ((c = getopt(argc, argv, "chvo:s:")) != -1)
    switch (c) {
      case 'c':
        MAGMACHECK_FLAG = true;
        break;
      default:
        std::cerr << "option not recognized, aborting." << std::endl;
        abort();
    }
58
  
59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80
  if (!MAGMACHECK_FLAG) {
    cout << "Dimension: " << basisRR.num_basis.size() << endl;
    cout << "Denominator: " << endl;
    PrintMagma(std::cout, basisRR.denom) << endl;
    cout << "Numerators: " << endl;
    for (size_t i = 0; i < basisRR.num_basis.size(); ++i)
      PrintMagma(std::cout, basisRR.num_basis[i]) << endl;
  } else {
    cout << "[" << endl;
    for (size_t i = 0; i < basisRR.num_basis.size(); ++i) {
      cout << "(";
      PrintMagma(std::cout, basisRR.num_basis[i]);
      cout << ")/(";
      PrintMagma(std::cout, basisRR.denom);
      cout << ")";
      if (i != basisRR.num_basis.size()-1)
        cout << ",";
      cout << endl;
    }
    cout << "]" << endl;
  }
  return EXIT_SUCCESS;
81
}