// -*- mode:C++ ; compile-command: "g++ -I.. -g gausspol.o sym2poly.o index.o input_lexer.o input_parser.o gaussint.o identificateur.o symbolic.o unary.o modpoly.o modfactor.o usual.o derive.o intg.o series.o subst.o vecteur.o moyal.o ifactor.o alg_EXT.o integrate.cc -lgmp" -*-

/*
 *  Copyright (C) 2000 B. Parisse, Institut Fourier, 38402 St Martin d'Heres
 *
 *  This program is free software; you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation; either version 3 of the License, or
 *  (at your option) any later version.
 *
 *  This program 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 General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with this program. If not, see <http://www.gnu.org/licenses/>.
 */

/* for profiler use -pg compile option then 
   a.out 
   then 
   gprof > test.gp 
   and 
   emacs test.gp &
*/

#include "sym2poly.h"
#include "intg.h"
#include <iostream>
#include <time.h>

using namespace giac;

int main(int ARGC, char *ARGV[]){
  vecteur v,l;
  readargs(ARGC,ARGV,l,v);
  identificateur x;
  if (ARGC>2){
    assert(v[1].type==_IDNT);
    x=*v[1]._IDNTptr;
  }
  else
    x=*(l.compptr->front()._IDNTptr);
  cout << x << endl;
  clock_t start, end;  
  start = clock();
  gen r;
  gen primitive(integrate(v.front(),x,r));
  if (r==0)
    cout << primitive << endl;
  else
    cout << primitive << "+int(" << r << "," << x << ")" << endl;
  end = clock();
  cout << "# Time for (integrate)" << double(end-start)/CLOCKS_PER_SEC << endl;
  // cout << "Verifying: " << simplify(r+primitive.diff(x)-v.front()) << endl ;
}