home *** CD-ROM | disk | FTP | other *** search
/ Linux Cubed Series 3: Developer Tools / Linux Cubed Series 3 - Developer Tools.iso / devel / lang / lisp / gcl-1.000 / gcl-1 / gcl-1.0 / mp / mp_dblrsl3.c < prev    next >
Encoding:
C/C++ Source or Header  |  1994-05-07  |  1.2 KB  |  46 lines
  1.  
  2. /*          Copyright (C) 1994 W. Schelter
  3.  
  4. This file is part of GNU Common Lisp, herein referred to as GCL  
  5.  
  6. GCL is free software; you can redistribute it and/or modify it under
  7. the terms of the GNU LIBRARY GENERAL PUBLIC LICENSE as published by
  8. the Free Software Foundation; either version 2, or (at your option)
  9. any later version.
  10.  
  11. GCL is distributed in the hope that it will be useful, but WITHOUT
  12. ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  13. FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  14. for more details.
  15.  
  16. You should have received a copy of the GNU library general public
  17. license along with GCL; see the file COPYING.  If not, write to the
  18. Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.
  19. */
  20.  
  21. #include "include.h"
  22. #include "genpari.h"
  23. #include "arith.h"
  24.  
  25. int dblremsl3(x,y,z)
  26.      int x,y,z;
  27. { unsigned long h;
  28.   unsigned long w;
  29.   if (x>= 0 && y>= 0 && z>0)
  30.      {w = mulul(x,y,h);
  31.       divul(x,z,h);
  32.       return h;}
  33.   else
  34.     { long save = avma;
  35.       GEN yy = stoi(y);
  36.       GEN xx = stoi(x);
  37.       GEN ans = mulss(xx,yy);
  38.       ans = dvmdii(ans,stoi(z),-1);
  39.       avma = save;
  40.       if (signe(ans) > 0) return ans[2];
  41.       if (signe(ans) < 0) return -ans[2];
  42.       return 0;}
  43. }
  44.   
  45.  
  46.