Blame view

Giac_maj/giac-1.4.9/examples/demo/pgcd.cpp 972 Bytes
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
  // -*- mode:C++ ; compile-command: "g++ -I.. -fPIC -DPIC -g -c pgcd.cpp -o pgcd.lo && ln -sf pgcd.lo pgcd.o && gcc -shared pgcd.lo -lc  -Wl,-soname -Wl,libpgcd.so.0 -o libpgcd.so.0.0.0 && ln -sf libpgcd.so.0.0.0 libpgcd.so.0 && ln -sf libpgcd.so.0.0.0 libpgcd.so" -*-
  using namespace std;
  #include <stdexcept>
  #include <cmath>
  #include <cstdlib>
  #include <giac/giac.h>
  //#include "pgcd.h"
  
  #ifndef NO_NAMESPACE_GIAC
  namespace giac {
  #endif // ndef NO_NAMESPACE_GIAC
  
    gen pgcd(gen a,gen b){
      gen q,r;
      for (;b!=0;){
        r=irem(a,b,q);
        a=b;
        b=r;
      }
      return a;
    }
    gen _pgcd(const gen & args){
      if ( (args.type!=_VECT) || (args._VECTptr->size()!=2))
        setsizeerr();
      vecteur &v=*args._VECTptr;
      return pgcd(v[0],v[1]);
    }
    const string _pgcd_s("pgcd");
    unary_function_unary __pgcd(&_pgcd,_pgcd_s);
    unary_function_ptr at_pgcd (&__pgcd,0,true);
    
  
  #ifndef NO_NAMESPACE_GIAC
  } // namespace giac
  #endif // ndef NO_NAMESPACE_GIAC