00001 /* 00002 * Copyright 2004 PathScale, Inc. All Rights Reserved. 00003 */ 00004 00005 /* 00006 00007 Copyright (C) 2000, 2001 Silicon Graphics, Inc. All Rights Reserved. 00008 00009 This program is free software; you can redistribute it and/or modify it 00010 under the terms of version 2.1 of the GNU Lesser General Public License 00011 as published by the Free Software Foundation. 00012 00013 This program is distributed in the hope that it would be useful, but 00014 WITHOUT ANY WARRANTY; without even the implied warranty of 00015 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. 00016 00017 Further, this software is distributed without any warranty that it is 00018 free of the rightful claim of any third person regarding infringement 00019 or the like. Any license provided herein, whether implied or 00020 otherwise, applies only to this software file. Patent licenses, if 00021 any, provided herein do not apply to combinations of this program with 00022 other software, or any other product whatsoever. 00023 00024 You should have received a copy of the GNU Lesser General Public 00025 License along with this program; if not, write the Free Software 00026 Foundation, Inc., 59 Temple Place - Suite 330, Boston MA 02111-1307, 00027 USA. 00028 00029 Contact information: Silicon Graphics, Inc., 1600 Amphitheatre Pky, 00030 Mountain View, CA 94043, or: 00031 00032 http://www.sgi.com 00033 00034 For further information regarding this notice, see: 00035 00036 http://oss.sgi.com/projects/GenInfo/NoticeExplan 00037 00038 */ 00039 00040 00041 /* --------------------------------------------------- */ 00042 /* | All Rights Reserved. | */ 00043 /* --------------------------------------------------- */ 00044 /* $Header: /proj/osprey/CVS/open64/osprey1.0/libF77/r_mod.c,v 1.1.1.1 2005/10/21 19:00:00 marcel Exp $ */ 00045 #include "cmplrs/host.h" 00046 #include <math.h> 00047 #include "moremath.h" 00048 00049 #ifdef KEY 00050 float __rmod(float x, float y) 00051 #else 00052 float __rmod(float x, float_t y) 00053 #endif // KEY 00054 { 00055 float q1, q2; 00056 00057 if ((q1 = x/y) >= 0.0F) 00058 q2 = floorf(q1); 00059 else 00060 q2 = -floorf(-q1); 00061 return (x - y * q2); 00062 } 00063 00064 #ifdef KEY 00065 float r_mod(float *x, float *y) 00066 #else 00067 float r_mod(float *x, float_t *y) 00068 #endif // KEY 00069 { 00070 return (__rmod(*x, *y)); 00071 }
1.5.6