Back to index

glibc  2.9
s_remquo.c
Go to the documentation of this file.
00001 /* Compute remainder and a congruent to the quotient.  m68k fpu version
00002    Copyright (C) 1997 Free Software Foundation, Inc.
00003    This file is part of the GNU C Library.
00004    Contributed by Andreas Schwab <schwab@issan.informatik.uni-dortmund.de>
00005 
00006    The GNU C Library is free software; you can redistribute it and/or
00007    modify it under the terms of the GNU Lesser General Public
00008    License as published by the Free Software Foundation; either
00009    version 2.1 of the License, or (at your option) any later version.
00010 
00011    The GNU C Library is distributed in the hope that it will be useful,
00012    but WITHOUT ANY WARRANTY; without even the implied warranty of
00013    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
00014    Lesser General Public License for more details.
00015 
00016    You should have received a copy of the GNU Lesser General Public
00017    License along with the GNU C Library; if not, write to the Free
00018    Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
00019    02111-1307 USA.  */
00020 
00021 #include <math.h>
00022 
00023 #ifndef SUFF
00024 #define SUFF
00025 #endif
00026 #ifndef float_type
00027 #define float_type double
00028 #endif
00029 
00030 #define CONCATX(a,b) __CONCAT(a,b)
00031 #define s(name) CONCATX(name,SUFF)
00032 
00033 float_type
00034 s(__remquo) (float_type x, float_type y, int *quo)
00035 {
00036   float_type result;
00037   int cquo, fpsr;
00038 
00039   __asm ("frem%.x %2,%0\n\tfmove%.l %/fpsr,%1"
00040         : "=f" (result), "=dm" (fpsr) : "f" (y), "0" (x));
00041   cquo = (fpsr >> 16) & 0x7f;
00042   if (fpsr & (1 << 23))
00043     cquo = -cquo;
00044   *quo = cquo;
00045   return result;
00046 }
00047 #define weak_aliasx(a,b) weak_alias(a,b)
00048 weak_aliasx (s(__remquo), s(remquo))