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
50
51
52
53
54
55
56
57
|
// -*- 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 ;
}
|