doc/exemples/exemple1/exec2/alglin.cxx

   1 #include "alglin.hxx"
   2
   3 extern "C" {
   4   double * valloc(long n);
   5   void vfree(double *x);
   6   void addvec_(double *C, double *A, double *B, long *n);
   7   double prdscl_(double *A, double *B, long *n);
   8 }
   9
  10 void   alglin::addvec(vecteur *C, vecteur *A, vecteur *B) {
  11   long n = A->n();
  12   addvec_(C->x(), A->x(), B->x(), &n);
  13 }
  14
  15 double alglin::prdscl(vecteur *A, vecteur *B) {
  16   long n = A->n();
  17   return prdscl_(A->x(), B->x(), &n);
  18 }
  19
  20 vecteur * alglin::create_vector(long n) {
  21   return new vecteur(n);
  22 }
  23
  24 void   alglin::destroy_vector(vecteur *V) {
  25   delete V;
  26 }