dp_simple.c 1.3 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849
  1. /* IEEE754 floating point arithmetic
  2. * double precision: common utilities
  3. */
  4. /*
  5. * MIPS floating point support
  6. * Copyright (C) 1994-2000 Algorithmics Ltd.
  7. *
  8. * This program is free software; you can distribute it and/or modify it
  9. * under the terms of the GNU General Public License (Version 2) as
  10. * published by the Free Software Foundation.
  11. *
  12. * This program is distributed in the hope it will be useful, but WITHOUT
  13. * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  14. * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
  15. * for more details.
  16. *
  17. * You should have received a copy of the GNU General Public License along
  18. * with this program; if not, write to the Free Software Foundation, Inc.,
  19. * 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA.
  20. */
  21. #include "ieee754dp.h"
  22. union ieee754dp ieee754dp_neg(union ieee754dp x)
  23. {
  24. unsigned int oldrm;
  25. union ieee754dp y;
  26. oldrm = ieee754_csr.rm;
  27. ieee754_csr.rm = FPU_CSR_RD;
  28. y = ieee754dp_sub(ieee754dp_zero(0), x);
  29. ieee754_csr.rm = oldrm;
  30. return y;
  31. }
  32. union ieee754dp ieee754dp_abs(union ieee754dp x)
  33. {
  34. unsigned int oldrm;
  35. union ieee754dp y;
  36. oldrm = ieee754_csr.rm;
  37. ieee754_csr.rm = FPU_CSR_RD;
  38. if (DPSIGN(x))
  39. y = ieee754dp_sub(ieee754dp_zero(0), x);
  40. else
  41. y = ieee754dp_add(ieee754dp_zero(0), x);
  42. ieee754_csr.rm = oldrm;
  43. return y;
  44. }