1 files changed, 56 insertions, 0 deletions
diff --git a/usr/src/lib/libbc/libc/gen/common/modf.c b/usr/src/lib/libbc/libc/gen/common/modf.c
new file mode 100644
index 0000000000..3ce1bb8256
--- /dev/null
+++ b/usr/src/lib/libbc/libc/gen/common/modf.c
@@ -0,0 +1,56 @@
+/*
+ * CDDL HEADER START
+ *
+ * The contents of this file are subject to the terms of the
+ * Common Development and Distribution License, Version 1.0 only
+ * (the "License").  You may not use this file except in compliance
+ * with the License.
+ *
+ * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
+ * or http://www.opensolaris.org/os/licensing.
+ * See the License for the specific language governing permissions
+ * and limitations under the License.
+ *
+ * When distributing Covered Code, include this CDDL HEADER in each
+ * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
+ * If applicable, add the following below this CDDL HEADER, with the
+ * fields enclosed by brackets "[]" replaced with your own identifying
+ * information: Portions Copyright [yyyy] [name of copyright owner]
+ *
+ * CDDL HEADER END
+ */
+#pragma ident	"%Z%%M%	%I%	%E% SMI"  /* from ATT S5R3 */
+
+/*	The following is extracted from...	*/
+
+/*	Copyright (c) 1984 AT&T	*/
+/*	  All Rights Reserved  	*/
+
+
+/*
+ * modf(value, iptr) returns the signed fractional part of value
+ * and stores the integer part indirectly through iptr.
+ *
+ */
+
+#define MAXPOWTWO	4.503599627370496000E+15
+			/* doubles >= MAXPOWTWO are already integers */
+double
+modf(value, iptr)
+double value;
+register double *iptr;
+{
+	register double absvalue;
+
+	if ((absvalue = (value >= 0.0) ? value : -value) >= MAXPOWTWO)
+		*iptr = value; /* it must be an integer */
+	else {
+		*iptr = absvalue + MAXPOWTWO; /* shift fraction off right */
+		*iptr -= MAXPOWTWO; /* shift back without fraction */
+		while (*iptr > absvalue) /* above arithmetic might round */
+			*iptr -= 1.0; /* test again just to be sure */
+		if (value < 0.0)
+			*iptr = -*iptr;
+	}
+	return (value - *iptr); /* signed fractional part */
+}