Blame view

Giac_maj/giac-1.4.9/src/factor.cc 1.55 KB
6663b6c9   adorian   projet complet av...
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
  // without NTL -*- mode:C++ ; compile-command: "g++ -I.. -g gausspol.o sym2poly.o index.o input_lexer.o input_parser.o gaussint.o identificateur.o unary.o modpoly.o modfactor.o symbolic.o derive.o intg.o series.o subst.o usual.o vecteur.o moyal.o ifactor.o alg_ext.o factor.cc -lgmp" -*-
  
  
  /* for profiler use -pg compile option then 
     a.out 
     then 
     gprof > test.gp 
     and 
     emacs test.gp &
  */
  
  /*
   *  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/>.
   */
  
  
  #include "sym2poly.h"
  #include <iostream>
  //#include <time.h>
  
  using namespace giac;
  
  // int tensor_dim;
  
  int main(int ARGC, char *ARGV[]){
    vecteur v,l;
    readargs(ARGC,ARGV,l,v);
    clock_t start, end;  
    vecteur::const_iterator it=v.begin(),itend=v.end();
    for (;it!=itend;++it){
      start = clock();
      cout << factor(*it) << " " <<endl ;
      end = clock();
      cout << "# Time" << double(end-start)/CLOCKS_PER_SEC << endl;
    }
  }