| b_log.c (93211) | b_log.c (108533) |
|---|---|
| 1/* 2 * Copyright (c) 1992, 1993 3 * The Regents of the University of California. All rights reserved. 4 * 5 * Redistribution and use in source and binary forms, with or without 6 * modification, are permitted provided that the following conditions 7 * are met: 8 * 1. Redistributions of source code must retain the above copyright --- 21 unchanged lines hidden (view full) --- 30 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF 31 * SUCH DAMAGE. 32 */ 33 34#ifndef lint 35static char sccsid[] = "@(#)log.c 8.2 (Berkeley) 11/30/93"; 36#endif /* not lint */ 37#include <sys/cdefs.h> | 1/* 2 * Copyright (c) 1992, 1993 3 * The Regents of the University of California. All rights reserved. 4 * 5 * Redistribution and use in source and binary forms, with or without 6 * modification, are permitted provided that the following conditions 7 * are met: 8 * 1. Redistributions of source code must retain the above copyright --- 21 unchanged lines hidden (view full) --- 30 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF 31 * SUCH DAMAGE. 32 */ 33 34#ifndef lint 35static char sccsid[] = "@(#)log.c 8.2 (Berkeley) 11/30/93"; 36#endif /* not lint */ 37#include <sys/cdefs.h> |
| 38__FBSDID("$FreeBSD: head/lib/msun/bsdsrc/b_log.c 93211 2002-03-26 11:59:29Z bde $"); | 38__FBSDID("$FreeBSD: head/lib/msun/bsdsrc/b_log.c 108533 2003-01-01 18:49:04Z schweikh $"); |
| 39 40#include <math.h> 41#include <errno.h> 42 43#include "mathimpl.h" 44 45/* Table-driven natural logarithm. 46 * --- 390 unchanged lines hidden (view full) --- 437 u2 = (u2 + logF_tail[j]) + q; /* tiny */ 438 u2 += logF_tail[N]*m; 439 return (u1 + u2); 440} 441#endif 442 443/* 444 * Extra precision variant, returning struct {double a, b;}; | 39 40#include <math.h> 41#include <errno.h> 42 43#include "mathimpl.h" 44 45/* Table-driven natural logarithm. 46 * --- 390 unchanged lines hidden (view full) --- 437 u2 = (u2 + logF_tail[j]) + q; /* tiny */ 438 u2 += logF_tail[N]*m; 439 return (u1 + u2); 440} 441#endif 442 443/* 444 * Extra precision variant, returning struct {double a, b;}; |
| 445 * log(x) = a+b to 63 bits, with a is rounded to 26 bits. | 445 * log(x) = a+b to 63 bits, with a rounded to 26 bits. |
| 446 */ 447struct Double 448#ifdef _ANSI_SOURCE 449__log__D(double x) 450#else 451__log__D(x) double x; 452#endif 453{ --- 37 unchanged lines hidden --- | 446 */ 447struct Double 448#ifdef _ANSI_SOURCE 449__log__D(double x) 450#else 451__log__D(x) double x; 452#endif 453{ --- 37 unchanged lines hidden --- |