as it should be.
SCCS-vsn: lib/libm/common_source/cosh.c 1.3
SCCS-vsn: lib/libm/common_source/exp.c 1.3
SCCS-vsn: lib/libm/common_source/exp__E.c 1.3
SCCS-vsn: lib/libm/common_source/expm1.c 1.3
SCCS-vsn: lib/libm/common_source/log.c 1.3
SCCS-vsn: lib/libm/common_source/log1p.c 1.3
SCCS-vsn: lib/libm/common_source/pow.c 1.3
SCCS-vsn: lib/libm/common_source/sinh.c 1.3
#ifndef lint
static char sccsid[] =
#ifndef lint
static char sccsid[] =
-"@(#)cosh.c 1.2 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%";
+"@(#)cosh.c 1.2 (Berkeley) 8/21/85; 1.3 (ucb.elefunt) %G%";
#endif not lint
/* COSH(X)
#endif not lint
/* COSH(X)
/* mln2hi = 8.8029691931113054792E1 , Hex 2^ 7 * .B00F33C7E22BDB */
/* mln2lo = -4.9650192275318476525E-16 , Hex 2^-50 * -.8F1B60279E582A */
/* lnovfl = 8.8029691931113053016E1 ; Hex 2^ 7 * .B00F33C7E22BDA */
/* mln2hi = 8.8029691931113054792E1 , Hex 2^ 7 * .B00F33C7E22BDB */
/* mln2lo = -4.9650192275318476525E-16 , Hex 2^-50 * -.8F1B60279E582A */
/* lnovfl = 8.8029691931113053016E1 ; Hex 2^ 7 * .B00F33C7E22BDA */
#define mln2lo (*(double*)mln2lox)
#define lnovfl (*(double*)lnovflx)
#else /* IEEE double */
#define mln2lo (*(double*)mln2lox)
#define lnovfl (*(double*)lnovflx)
#else /* IEEE double */
mln2hi = 7.0978271289338397310E2 , /*Hex 2^ 10 * 1.62E42FEFA39EF */
mln2lo = 2.3747039373786107478E-14 , /*Hex 2^-45 * 1.ABC9E3B39803F */
lnovfl = 7.0978271289338397310E2 ; /*Hex 2^ 9 * 1.62E42FEFA39EF */
mln2hi = 7.0978271289338397310E2 , /*Hex 2^ 10 * 1.62E42FEFA39EF */
mln2lo = 2.3747039373786107478E-14 , /*Hex 2^-45 * 1.ABC9E3B39803F */
lnovfl = 7.0978271289338397310E2 ; /*Hex 2^ 9 * 1.62E42FEFA39EF */
#ifndef lint
static char sccsid[] =
#ifndef lint
static char sccsid[] =
-"@(#)exp.c 4.3 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%";
+"@(#)exp.c 4.3 (Berkeley) 8/21/85; 1.3 (ucb.elefunt) %G%";
#endif not lint
/* EXP(X)
#endif not lint
/* EXP(X)
*/
#ifdef VAX /* VAX D format */
*/
#ifdef VAX /* VAX D format */
/* ln2hi = 6.9314718055829871446E-1 , Hex 2^ 0 * .B17217F7D00000 */
/* ln2lo = 1.6465949582897081279E-12 , Hex 2^-39 * .E7BCD5E4F1D9CC */
/* lnhuge = 9.4961163736712506989E1 , Hex 2^ 7 * .BDEC1DA73E9010 */
/* ln2hi = 6.9314718055829871446E-1 , Hex 2^ 0 * .B17217F7D00000 */
/* ln2lo = 1.6465949582897081279E-12 , Hex 2^-39 * .E7BCD5E4F1D9CC */
/* lnhuge = 9.4961163736712506989E1 , Hex 2^ 7 * .BDEC1DA73E9010 */
#define lntiny (*(double*)lntinyx)
#define invln2 (*(double*)invln2x)
#else /* IEEE double */
#define lntiny (*(double*)lntinyx)
#define invln2 (*(double*)invln2x)
#else /* IEEE double */
ln2hi = 6.9314718036912381649E-1 , /*Hex 2^ -1 * 1.62E42FEE00000 */
ln2lo = 1.9082149292705877000E-10 , /*Hex 2^-33 * 1.A39EF35793C76 */
lnhuge = 7.1602103751842355450E2 , /*Hex 2^ 9 * 1.6602B15B7ECF2 */
ln2hi = 6.9314718036912381649E-1 , /*Hex 2^ -1 * 1.62E42FEE00000 */
ln2lo = 1.9082149292705877000E-10 , /*Hex 2^-33 * 1.A39EF35793C76 */
lnhuge = 7.1602103751842355450E2 , /*Hex 2^ 9 * 1.6602B15B7ECF2 */
#ifndef lint
static char sccsid[] =
#ifndef lint
static char sccsid[] =
-"@(#)exp__E.c 1.2 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%";
+"@(#)exp__E.c 1.2 (Berkeley) 8/21/85; 1.3 (ucb.elefunt) %G%";
#endif not lint
/* exp__E(x,c)
#endif not lint
/* exp__E(x,c)
double exp__E(x,c)
double x,c;
{
double exp__E(x,c)
double x,c;
{
- double static zero=0.0, one=1.0, half=1.0/2.0, small=1.0E-19;
+ static double zero=0.0, one=1.0, half=1.0/2.0, small=1.0E-19;
double copysign(),z,p,q,xp,xh,w;
if(copysign(x,one)>small) {
z = x*x ;
double copysign(),z,p,q,xp,xh,w;
if(copysign(x,one)>small) {
z = x*x ;
#ifndef lint
static char sccsid[] =
#ifndef lint
static char sccsid[] =
-"@(#)expm1.c 1.2 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%";
+"@(#)expm1.c 1.2 (Berkeley) 8/21/85; 1.3 (ucb.elefunt) %G%";
#endif not lint
/* EXPM1(X)
#endif not lint
/* EXPM1(X)
*/
#ifdef VAX /* VAX D format */
*/
#ifdef VAX /* VAX D format */
/* ln2hi = 6.9314718055829871446E-1 , Hex 2^ 0 * .B17217F7D00000 */
/* ln2lo = 1.6465949582897081279E-12 , Hex 2^-39 * .E7BCD5E4F1D9CC */
/* lnhuge = 9.4961163736712506989E1 , Hex 2^ 7 * .BDEC1DA73E9010 */
/* ln2hi = 6.9314718055829871446E-1 , Hex 2^ 0 * .B17217F7D00000 */
/* ln2lo = 1.6465949582897081279E-12 , Hex 2^-39 * .E7BCD5E4F1D9CC */
/* lnhuge = 9.4961163736712506989E1 , Hex 2^ 7 * .BDEC1DA73E9010 */
#define lnhuge (*(double*)lnhugex)
#define invln2 (*(double*)invln2x)
#else /* IEEE double */
#define lnhuge (*(double*)lnhugex)
#define invln2 (*(double*)invln2x)
#else /* IEEE double */
ln2hi = 6.9314718036912381649E-1 , /*Hex 2^ -1 * 1.62E42FEE00000 */
ln2lo = 1.9082149292705877000E-10 , /*Hex 2^-33 * 1.A39EF35793C76 */
lnhuge = 7.1602103751842355450E2 , /*Hex 2^ 9 * 1.6602B15B7ECF2 */
ln2hi = 6.9314718036912381649E-1 , /*Hex 2^ -1 * 1.62E42FEE00000 */
ln2lo = 1.9082149292705877000E-10 , /*Hex 2^-33 * 1.A39EF35793C76 */
lnhuge = 7.1602103751842355450E2 , /*Hex 2^ 9 * 1.6602B15B7ECF2 */
double expm1(x)
double x;
{
double expm1(x)
double x;
{
- double static one=1.0, half=1.0/2.0;
+ static double one=1.0, half=1.0/2.0;
double scalb(), copysign(), exp__E(), z,hi,lo,c;
int k,finite();
#ifdef VAX
double scalb(), copysign(), exp__E(), z,hi,lo,c;
int k,finite();
#ifdef VAX
#ifndef lint
static char sccsid[] =
#ifndef lint
static char sccsid[] =
-"@(#)log.c 4.5 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%";
+"@(#)log.c 4.5 (Berkeley) 8/21/85; 1.3 (ucb.elefunt) %G%";
#endif not lint
/* LOG(X)
#endif not lint
/* LOG(X)
#ifdef VAX /* VAX D format */
#include <errno.h>
#ifdef VAX /* VAX D format */
#include <errno.h>
/* ln2hi = 6.9314718055829871446E-1 , Hex 2^ 0 * .B17217F7D00000 */
/* ln2lo = 1.6465949582897081279E-12 , Hex 2^-39 * .E7BCD5E4F1D9CC */
/* sqrt2 = 1.4142135623730950622E0 ; Hex 2^ 1 * .B504F333F9DE65 */
/* ln2hi = 6.9314718055829871446E-1 , Hex 2^ 0 * .B17217F7D00000 */
/* ln2lo = 1.6465949582897081279E-12 , Hex 2^-39 * .E7BCD5E4F1D9CC */
/* sqrt2 = 1.4142135623730950622E0 ; Hex 2^ 1 * .B504F333F9DE65 */
#define ln2lo (*(double*)ln2lox)
#define sqrt2 (*(double*)sqrt2x)
#else /* IEEE double */
#define ln2lo (*(double*)ln2lox)
#define sqrt2 (*(double*)sqrt2x)
#else /* IEEE double */
ln2hi = 6.9314718036912381649E-1 , /*Hex 2^ -1 * 1.62E42FEE00000 */
ln2lo = 1.9082149292705877000E-10 , /*Hex 2^-33 * 1.A39EF35793C76 */
sqrt2 = 1.4142135623730951455E0 ; /*Hex 2^ 0 * 1.6A09E667F3BCD */
ln2hi = 6.9314718036912381649E-1 , /*Hex 2^ -1 * 1.62E42FEE00000 */
ln2lo = 1.9082149292705877000E-10 , /*Hex 2^-33 * 1.A39EF35793C76 */
sqrt2 = 1.4142135623730951455E0 ; /*Hex 2^ 0 * 1.6A09E667F3BCD */
#ifndef lint
static char sccsid[] =
#ifndef lint
static char sccsid[] =
-"@(#)log1p.c 1.3 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%";
+"@(#)log1p.c 1.3 (Berkeley) 8/21/85; 1.3 (ucb.elefunt) %G%";
#endif not lint
/* LOG1P(x)
#endif not lint
/* LOG1P(x)
#ifdef VAX /* VAX D format */
#include <errno.h>
#ifdef VAX /* VAX D format */
#include <errno.h>
/* ln2hi = 6.9314718055829871446E-1 , Hex 2^ 0 * .B17217F7D00000 */
/* ln2lo = 1.6465949582897081279E-12 , Hex 2^-39 * .E7BCD5E4F1D9CC */
/* sqrt2 = 1.4142135623730950622E0 ; Hex 2^ 1 * .B504F333F9DE65 */
/* ln2hi = 6.9314718055829871446E-1 , Hex 2^ 0 * .B17217F7D00000 */
/* ln2lo = 1.6465949582897081279E-12 , Hex 2^-39 * .E7BCD5E4F1D9CC */
/* sqrt2 = 1.4142135623730950622E0 ; Hex 2^ 1 * .B504F333F9DE65 */
#define ln2lo (*(double*)ln2lox)
#define sqrt2 (*(double*)sqrt2x)
#else /* IEEE double */
#define ln2lo (*(double*)ln2lox)
#define sqrt2 (*(double*)sqrt2x)
#else /* IEEE double */
ln2hi = 6.9314718036912381649E-1 , /*Hex 2^ -1 * 1.62E42FEE00000 */
ln2lo = 1.9082149292705877000E-10 , /*Hex 2^-33 * 1.A39EF35793C76 */
sqrt2 = 1.4142135623730951455E0 ; /*Hex 2^ 0 * 1.6A09E667F3BCD */
ln2hi = 6.9314718036912381649E-1 , /*Hex 2^ -1 * 1.62E42FEE00000 */
ln2lo = 1.9082149292705877000E-10 , /*Hex 2^-33 * 1.A39EF35793C76 */
sqrt2 = 1.4142135623730951455E0 ; /*Hex 2^ 0 * 1.6A09E667F3BCD */
#ifndef lint
static char sccsid[] =
#ifndef lint
static char sccsid[] =
-"@(#)pow.c 4.5 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%";
+"@(#)pow.c 4.5 (Berkeley) 8/21/85; 1.3 (ucb.elefunt) %G%";
#endif not lint
/* POW(X,Y)
#endif not lint
/* POW(X,Y)
#include <errno.h>
extern double infnan();
#include <errno.h>
extern double infnan();
/* ln2hi = 6.9314718055829871446E-1 , Hex 2^ 0 * .B17217F7D00000 */
/* ln2lo = 1.6465949582897081279E-12 , Hex 2^-39 * .E7BCD5E4F1D9CC */
/* invln2 = 1.4426950408889634148E0 , Hex 2^ 1 * .B8AA3B295C17F1 */
/* ln2hi = 6.9314718055829871446E-1 , Hex 2^ 0 * .B17217F7D00000 */
/* ln2lo = 1.6465949582897081279E-12 , Hex 2^-39 * .E7BCD5E4F1D9CC */
/* invln2 = 1.4426950408889634148E0 , Hex 2^ 1 * .B8AA3B295C17F1 */
#define invln2 (*(double*)invln2x)
#define sqrt2 (*(double*)sqrt2x)
#else /* IEEE double */
#define invln2 (*(double*)invln2x)
#define sqrt2 (*(double*)sqrt2x)
#else /* IEEE double */
ln2hi = 6.9314718036912381649E-1 , /*Hex 2^ -1 * 1.62E42FEE00000 */
ln2lo = 1.9082149292705877000E-10 , /*Hex 2^-33 * 1.A39EF35793C76 */
invln2 = 1.4426950408889633870E0 , /*Hex 2^ 0 * 1.71547652B82FE */
sqrt2 = 1.4142135623730951455E0 ; /*Hex 2^ 0 * 1.6A09E667F3BCD */
#endif
ln2hi = 6.9314718036912381649E-1 , /*Hex 2^ -1 * 1.62E42FEE00000 */
ln2lo = 1.9082149292705877000E-10 , /*Hex 2^-33 * 1.A39EF35793C76 */
invln2 = 1.4426950408889633870E0 , /*Hex 2^ 0 * 1.71547652B82FE */
sqrt2 = 1.4142135623730951455E0 ; /*Hex 2^ 0 * 1.6A09E667F3BCD */
#endif
-double static zero=0.0, half=1.0/2.0, one=1.0, two=2.0, negone= -1.0;
+static double zero=0.0, half=1.0/2.0, one=1.0, two=2.0, negone= -1.0;
double pow(x,y)
double x,y;
double pow(x,y)
double x,y;
#ifndef lint
static char sccsid[] =
#ifndef lint
static char sccsid[] =
-"@(#)sinh.c 4.3 (Berkeley) 8/21/85; 1.2 (ucb.elefunt) %G%";
+"@(#)sinh.c 4.3 (Berkeley) 8/21/85; 1.3 (ucb.elefunt) %G%";
#endif not lint
/* SINH(X)
#endif not lint
/* SINH(X)
/* mln2hi = 8.8029691931113054792E1 , Hex 2^ 7 * .B00F33C7E22BDB */
/* mln2lo = -4.9650192275318476525E-16 , Hex 2^-50 * -.8F1B60279E582A */
/* lnovfl = 8.8029691931113053016E1 ; Hex 2^ 7 * .B00F33C7E22BDA */
/* mln2hi = 8.8029691931113054792E1 , Hex 2^ 7 * .B00F33C7E22BDB */
/* mln2lo = -4.9650192275318476525E-16 , Hex 2^-50 * -.8F1B60279E582A */
/* lnovfl = 8.8029691931113053016E1 ; Hex 2^ 7 * .B00F33C7E22BDA */
#define mln2lo (*(double*)mln2lox)
#define lnovfl (*(double*)lnovflx)
#else /* IEEE double */
#define mln2lo (*(double*)mln2lox)
#define lnovfl (*(double*)lnovflx)
#else /* IEEE double */
mln2hi = 7.0978271289338397310E2 , /*Hex 2^ 10 * 1.62E42FEFA39EF */
mln2lo = 2.3747039373786107478E-14 , /*Hex 2^-45 * 1.ABC9E3B39803F */
lnovfl = 7.0978271289338397310E2 ; /*Hex 2^ 9 * 1.62E42FEFA39EF */
mln2hi = 7.0978271289338397310E2 , /*Hex 2^ 10 * 1.62E42FEFA39EF */
mln2lo = 2.3747039373786107478E-14 , /*Hex 2^-45 * 1.ABC9E3B39803F */
lnovfl = 7.0978271289338397310E2 ; /*Hex 2^ 9 * 1.62E42FEFA39EF */