Guido van Rossum | f70e43a | 1991-02-19 12:39:46 +0000 | [diff] [blame] | 1 | /*********************************************************** |
Guido van Rossum | 6d023c9 | 1995-01-04 19:12:13 +0000 | [diff] [blame] | 2 | Copyright 1991-1995 by Stichting Mathematisch Centrum, Amsterdam, |
| 3 | The Netherlands. |
Guido van Rossum | f70e43a | 1991-02-19 12:39:46 +0000 | [diff] [blame] | 4 | |
| 5 | All Rights Reserved |
| 6 | |
| 7 | Permission to use, copy, modify, and distribute this software and its |
| 8 | documentation for any purpose and without fee is hereby granted, |
| 9 | provided that the above copyright notice appear in all copies and that |
| 10 | both that copyright notice and this permission notice appear in |
| 11 | supporting documentation, and that the names of Stichting Mathematisch |
| 12 | Centrum or CWI not be used in advertising or publicity pertaining to |
| 13 | distribution of the software without specific, written prior permission. |
| 14 | |
| 15 | STICHTING MATHEMATISCH CENTRUM DISCLAIMS ALL WARRANTIES WITH REGARD TO |
| 16 | THIS SOFTWARE, INCLUDING ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND |
| 17 | FITNESS, IN NO EVENT SHALL STICHTING MATHEMATISCH CENTRUM BE LIABLE |
| 18 | FOR ANY SPECIAL, INDIRECT OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES |
| 19 | WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN |
| 20 | ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT |
| 21 | OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE. |
| 22 | |
| 23 | ******************************************************************/ |
| 24 | |
Guido van Rossum | 85a5fbb | 1990-10-14 12:07:46 +0000 | [diff] [blame] | 25 | /* Portable fmod(x, y) implementation for systems that don't have it */ |
| 26 | |
Guido van Rossum | 464a0a1 | 1995-02-27 10:15:10 +0000 | [diff] [blame] | 27 | #include "config.h" |
Guido van Rossum | 464a0a1 | 1995-02-27 10:15:10 +0000 | [diff] [blame] | 28 | |
| 29 | #include "mymath.h" |
Guido van Rossum | 85a5fbb | 1990-10-14 12:07:46 +0000 | [diff] [blame] | 30 | #include <errno.h> |
| 31 | |
Guido van Rossum | 85a5fbb | 1990-10-14 12:07:46 +0000 | [diff] [blame] | 32 | double |
Guido van Rossum | 464a0a1 | 1995-02-27 10:15:10 +0000 | [diff] [blame] | 33 | fmod(double x, double y) |
Guido van Rossum | 85a5fbb | 1990-10-14 12:07:46 +0000 | [diff] [blame] | 34 | { |
| 35 | double i, f; |
| 36 | |
| 37 | if (y == 0.0) { |
| 38 | errno = EDOM; |
| 39 | return 0.0; |
| 40 | } |
| 41 | |
| 42 | /* return f such that x = i*y + f for some integer i |
| 43 | such that |f| < |y| and f has the same sign as x */ |
| 44 | |
| 45 | i = floor(x/y); |
| 46 | f = x - i*y; |
| 47 | if ((x < 0.0) != (y < 0.0)) |
| 48 | f = f-y; |
| 49 | return f; |
| 50 | } |