Document extended-precision routines.

22 files changed, 341 insertions, 117 deletions

diff --git a/lib/libm/Makefile b/lib/libm/Makefile
index 7023ce8dce1..783a2ea308c 100644
--- a/lib/libm/Makefile
+++ b/lib/libm/Makefile
@@ -1,4 +1,4 @@
-#  $OpenBSD: Makefile,v 1.84 2011/07/06 03:50:23 martynas Exp $
+#  $OpenBSD: Makefile,v 1.85 2011/07/07 00:54:15 martynas Exp $
 #  $NetBSD: Makefile,v 1.28 1995/11/20 22:06:19 jtc Exp $
 #
 #  @(#)Makefile 5.1beta 93/09/24 
@@ -299,25 +299,57 @@ MLINKS+=trunc.3 truncf.3
 
 # long double versions
 MLINKS+=acos.3 acosl.3
+MLINKS+=acosh.3 acoshl.3
 MLINKS+=asin.3 asinl.3
+MLINKS+=asinh.3 asinhl.3
 MLINKS+=atan.3 atanl.3
 MLINKS+=atan2.3 atan2l.3
+MLINKS+=atanh.3 atanhl.3
+MLINKS+=ceil.3 ceill.3
 MLINKS+=copysign.3 copysignl.3
 MLINKS+=cos.3 cosl.3
+MLINKS+=cosh.3 coshl.3
+MLINKS+=erf.3 erfcl.3
+MLINKS+=erf.3 erfl.3
 MLINKS+=exp.3 exp2l.3
+MLINKS+=exp.3 expl.3
+MLINKS+=exp.3 expm1l.3
+MLINKS+=exp.3 log10l.3
+MLINKS+=exp.3 log1pl.3
+MLINKS+=exp.3 log2l.3
+MLINKS+=exp.3 logl.3
+MLINKS+=exp.3 powl.3
 MLINKS+=fabs.3 fabsl.3
 MLINKS+=fdim.3 fdiml.3
+MLINKS+=floor.3 floorl.3
 MLINKS+=fmax.3 fmaxl.3
 MLINKS+=fmax.3 fminl.3
+MLINKS+=fmod.3 fmodl.3
+MLINKS+=hypot.3 hypotl.3
 MLINKS+=ilogb.3 ilogbl.3
+MLINKS+=lgamma.3 lgammal.3
+MLINKS+=lgamma.3 tgammal.3
 MLINKS+=logb.3 logbl.3
+MLINKS+=logb.3 scalbl.3
+MLINKS+=lrint.3 llrintl.3
+MLINKS+=lrint.3 lrintl.3
+MLINKS+=lround.3 llroundl.3
+MLINKS+=lround.3 lroundl.3
 MLINKS+=nan.3 nanl.3
 MLINKS+=nearbyint.3 nearbyintl.3
+MLINKS+=nextafter.3 nextafterl.3
+MLINKS+=remainder.3 remainderl.3
+MLINKS+=remainder.3 remquol.3
 MLINKS+=rint.3 rintl.3
+MLINKS+=round.3 roundl.3
 MLINKS+=scalbln.3 scalblnl.3
 MLINKS+=scalbn.3 scalbnl.3
 MLINKS+=sin.3 sinl.3
+MLINKS+=sinh.3 sinhl.3
+MLINKS+=sqrt.3 cbrtl.3
 MLINKS+=sqrt.3 sqrtl.3
 MLINKS+=tan.3 tanl.3
+MLINKS+=tanh.3 tanhl.3
+MLINKS+=trunc.3 truncl.3
 
 .include <bsd.lib.mk>
diff --git a/lib/libm/man/acosh.3 b/lib/libm/man/acosh.3
index d2c0fe775e5..0e87b13772c 100644
--- a/lib/libm/man/acosh.3
+++ b/lib/libm/man/acosh.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: acosh.3,v 1.11 2007/05/31 19:19:35 jmc Exp $
+.\"	$OpenBSD: acosh.3,v 1.12 2011/07/07 00:54:16 martynas Exp $
 .\" Copyright (c) 1991 Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)acosh.3	5.2 (Berkeley) 5/6/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt ACOSH 3
 .Os
 .Sh NAME
 .Nm acosh ,
-.Nm acoshf
+.Nm acoshf ,
+.Nm acoshl
 .Nd inverse hyperbolic cosine functions
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -41,6 +42,8 @@
 .Fn acosh "double x"
 .Ft float
 .Fn acoshf "float x"
+.Ft long double
+.Fn acoshl "long double x"
 .Sh DESCRIPTION
 The
 .Fn acosh
@@ -52,11 +55,16 @@ The
 .Fn acoshf
 function is a single precision version of
 .Fn acosh .
+The
+.Fn acoshl
+function is an extended precision version of
+.Fn acosh .
 .Sh RETURN VALUES
 If x is less than one,
-.Fn acosh "x"
-and
+.Fn acosh "x" ,
 .Fn acoshf "x"
+and
+.Fn acoshl "x"
 return NaN and set the global variable
 .Va errno
 to EDOM.
diff --git a/lib/libm/man/asinh.3 b/lib/libm/man/asinh.3
index ce7223fcc1b..17f72957391 100644
--- a/lib/libm/man/asinh.3
+++ b/lib/libm/man/asinh.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: asinh.3,v 1.11 2007/05/31 19:19:35 jmc Exp $
+.\"	$OpenBSD: asinh.3,v 1.12 2011/07/07 00:54:16 martynas Exp $
 .\" Copyright (c) 1985, 1991 Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)asinh.3	6.4 (Berkeley) 5/6/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt ASINH 3
 .Os
 .Sh NAME
 .Nm asinh ,
-.Nm asinhf
+.Nm asinhf ,
+.Nm asinhl
 .Nd inverse hyperbolic sine functions
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -41,6 +42,8 @@
 .Fn asinh "double x"
 .Ft float
 .Fn asinhf "float x"
+.Ft long double
+.Fn asinhl "long double x"
 .Sh DESCRIPTION
 The
 .Fn asinh
@@ -50,11 +53,16 @@ The
 .Fn asinhf
 function is a single precision version of
 .Fn asinh .
+The
+.Fn asinhl
+function is an extended precision version of
+.Fn asinh .
 .Sh RETURN VALUES
 The
-.Fn asinh
-and
+.Fn asinh ,
 .Fn asinhf
+and
+.Fn asinhl
 functions return the inverse hyperbolic sine of
 .Ar x .
 .Sh SEE ALSO
diff --git a/lib/libm/man/atanh.3 b/lib/libm/man/atanh.3
index 03ac2e4bd99..e1c8a770335 100644
--- a/lib/libm/man/atanh.3
+++ b/lib/libm/man/atanh.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: atanh.3,v 1.11 2007/05/31 19:19:35 jmc Exp $
+.\"	$OpenBSD: atanh.3,v 1.12 2011/07/07 00:54:16 martynas Exp $
 .\" Copyright (c) 1985, 1991 Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)atanh.3	5.2 (Berkeley) 5/6/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt ATANH 3
 .Os
 .Sh NAME
 .Nm atanh ,
-.Nm atanhf
+.Nm atanhf ,
+.Nm atanhl
 .Nd inverse hyperbolic tangent functions
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -41,6 +42,8 @@
 .Fn atanh "double x"
 .Ft float
 .Fn atanhf "float x"
+.Ft long double
+.Fn atanhl "long double x"
 .Sh DESCRIPTION
 The
 .Fn atanh
@@ -52,11 +55,16 @@ The
 .Fn atanhf
 function is a single precision version of
 .Fn atanh .
+The
+.Fn atanhl
+function is an extended precision version of
+.Fn atanh .
 .Sh RETURN VALUES
 If |x|>=1,
-.Fn atanh "x"
-and
+.Fn atanh "x" ,
 .Fn atanhf "x"
+and
+.Fn atanhl "x"
 return +inf, -inf or NaN, and set the global variable
 .Va errno
 to EDOM.
diff --git a/lib/libm/man/ceil.3 b/lib/libm/man/ceil.3
index ff0f95905c0..1ef2c0a11ec 100644
--- a/lib/libm/man/ceil.3
+++ b/lib/libm/man/ceil.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: ceil.3,v 1.11 2009/08/03 21:26:07 martynas Exp $
+.\"	$OpenBSD: ceil.3,v 1.12 2011/07/07 00:54:16 martynas Exp $
 .\" Copyright (c) 1991 The Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)ceil.3	5.1 (Berkeley) 5/2/91
 .\"
-.Dd $Mdocdate: August 3 2009 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt CEIL 3
 .Os
 .Sh NAME
 .Nm ceil ,
-.Nm ceilf
+.Nm ceilf ,
+.Nm ceill
 .Nd "round to smallest integral value greater than or equal to x"
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -41,6 +42,8 @@
 .Fn ceil "double x"
 .Ft float
 .Fn ceilf "float x"
+.Ft long double
+.Fn ceill "long double x"
 .Sh DESCRIPTION
 The
 .Fn ceil
@@ -51,6 +54,10 @@ The
 .Fn ceilf
 function is a single precision version of
 .Fn ceil .
+The
+.Fn ceill
+function is an extended precision version of
+.Fn ceil .
 .Sh SEE ALSO
 .Xr abs 3 ,
 .Xr fabs 3 ,
diff --git a/lib/libm/man/cosh.3 b/lib/libm/man/cosh.3
index d3ed3e15e84..59b09faf09e 100644
--- a/lib/libm/man/cosh.3
+++ b/lib/libm/man/cosh.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: cosh.3,v 1.10 2007/05/31 19:19:35 jmc Exp $
+.\"	$OpenBSD: cosh.3,v 1.11 2011/07/07 00:54:16 martynas Exp $
 .\" Copyright (c) 1989, 1991 The Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)cosh.3	5.1 (Berkeley) 5/2/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt COSH 3
 .Os
 .Sh NAME
 .Nm cosh ,
-.Nm coshf
+.Nm coshf ,
+.Nm coshl
 .Nd hyperbolic cosine functions
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -41,6 +42,8 @@
 .Fn cosh "double x"
 .Ft float
 .Fn coshf "float x"
+.Ft long double
+.Fn coshl "long double x"
 .Sh DESCRIPTION
 The
 .Fn cosh
@@ -50,11 +53,16 @@ The
 .Fn coshf
 function is a single precision version of
 .Fn cosh .
+The
+.Fn coshl
+function is an extended precision version of
+.Fn cosh .
 .Sh RETURN VALUES
 If the magnitude of x is too large,
-.Fn cosh "x"
-and
+.Fn cosh "x" ,
 .Fn coshf "x"
+and
+.Fn coshl "x"
 return Inf and set the global variable
 .Va errno
 to ERANGE.
diff --git a/lib/libm/man/erf.3 b/lib/libm/man/erf.3
index 2414dc90cb0..732307f38da 100644
--- a/lib/libm/man/erf.3
+++ b/lib/libm/man/erf.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: erf.3,v 1.11 2007/05/31 19:19:35 jmc Exp $
+.\"	$OpenBSD: erf.3,v 1.12 2011/07/07 00:54:16 martynas Exp $
 .\" Copyright (c) 1985, 1991 Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,14 +28,16 @@
 .\"
 .\"     from: @(#)erf.3	6.4 (Berkeley) 4/20/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt ERFC 3
 .Os
 .Sh NAME
 .Nm erf ,
 .Nm erff ,
+.Nm erfl ,
 .Nm erfc ,
-.Nm erfcf
+.Nm erfcf ,
+.Nm erfcl
 .Nd error function operators
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -43,10 +45,14 @@
 .Fn erf "double x"
 .Ft float
 .Fn erff "float x"
+.Ft long double
+.Fn erfl "long double x"
 .Ft double
 .Fn erfc "double x"
 .Ft float
 .Fn erfcf "float x"
+.Ft long double
+.Fn erfcl "long double x"
 .Sh DESCRIPTION
 These functions calculate the error function of
 .Fa x .
@@ -56,7 +62,9 @@ The
 calculates the error function of x, and
 the
 .Fn erff
-function is a single precision version of
+and
+.Fn erfl
+functions are single and double precision versions of
 .Fn erf .
 The error function is defined as:
 .Bd -filled -offset indent
@@ -81,7 +89,9 @@ This is useful, since for large
 places disappear.
 The
 .Fn erfcf
-function is a single precision version of
+and
+.Fn erfcl
+functions are single and double precision version of
 .Fn erfc .
 .Sh SEE ALSO
 .Xr math 3
diff --git a/lib/libm/man/exp.3 b/lib/libm/man/exp.3
index 39037e299d5..c8a625a9a88 100644
--- a/lib/libm/man/exp.3
+++ b/lib/libm/man/exp.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: exp.3,v 1.22 2009/10/27 19:26:14 jmc Exp $
+.\"	$OpenBSD: exp.3,v 1.23 2011/07/07 00:54:16 martynas Exp $
 .\"
 .\" Copyright (c) 1985, 1991 Regents of the University of California.
 .\" All rights reserved.
@@ -29,27 +29,34 @@
 .\"
 .\"     from: @(#)exp.3	6.12 (Berkeley) 7/31/91
 .\"
-.Dd $Mdocdate: October 27 2009 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt EXP 3
 .Os
 .Sh NAME
 .Nm exp ,
 .Nm expf ,
+.Nm expl ,
 .Nm exp2 ,
 .Nm exp2f ,
 .Nm exp2l ,
 .Nm expm1 ,
 .Nm expm1f ,
+.Nm expm1l ,
 .Nm log ,
 .Nm logf ,
+.Nm logl ,
 .Nm log2 ,
 .Nm log2f ,
+.Nm log2l ,
 .Nm log10 ,
 .Nm log10f ,
+.Nm log10l ,
 .Nm log1p ,
 .Nm log1pf ,
+.Nm log1pl ,
 .Nm pow ,
-.Nm powf
+.Nm powf ,
+.Nm powl
 .Nd exponential, logarithm, power functions
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -57,6 +64,8 @@
 .Fn exp "double x"
 .Ft float
 .Fn expf "float x"
+.Ft long double
+.Fn expl "long double x"
 .Ft double
 .Fn exp2 "double x"
 .Ft float
@@ -67,26 +76,38 @@
 .Fn expm1 "double x"
 .Ft float
 .Fn expm1f "float x"
+.Ft long double
+.Fn expm1l "long double x"
 .Ft double
 .Fn log "double x"
 .Ft float
 .Fn logf "float x"
+.Ft long double
+.Fn logl "long double x"
 .Ft double
 .Fn log2 "double x"
 .Ft float
 .Fn log2f "float x"
+.Ft long double
+.Fn log2l "long double x"
 .Ft double
 .Fn log10 "double x"
 .Ft float
 .Fn log10f "float x"
+.Ft long double
+.Fn log10l "long double x"
 .Ft double
 .Fn log1p "double x"
 .Ft float
 .Fn log1pf "float x"
+.Ft long double
+.Fn log1pl "long double x"
 .Ft double
 .Fn pow "double x" "double y"
 .Ft float
 .Fn powf "float x" "float y"
+.Ft long double
+.Fn powl "long double x" "long double y"
 .Sh DESCRIPTION
 The
 .Fn exp
@@ -98,6 +119,10 @@ The
 .Fn expf
 function is a single precision version of
 .Fn exp .
+The
+.Fn expl
+function is an extended precision version of
+.Fn exp .
 .Pp
 The
 .Fn exp2
@@ -120,6 +145,10 @@ The
 .Fn expm1f
 function is a single precision version of
 .Fn expm1 .
+The
+.Fn expm1l
+function is an extended precision version of
+.Fn expm1 .
 .Pp
 The
 .Fn log
@@ -129,6 +158,10 @@ The
 .Fn logf
 function is a single precision version of
 .Fn log .
+The
+.Fn logl
+function is an extended precision version of
+.Fn log .
 .Pp
 The
 .Fn log2
@@ -139,6 +172,10 @@ The
 .Fn log2f
 function is a single precision version of
 .Fn log2 .
+The
+.Fn log2l
+function is an extended precision version of
+.Fn log2 .
 .Pp
 The
 .Fn log10
@@ -149,6 +186,10 @@ The
 .Fn log10f
 function is a single precision version of
 .Fn log10 .
+The
+.Fn log10l
+function is an extended precision version of
+.Fn log10 .
 .Pp
 The
 .Fn log1p
@@ -159,6 +200,10 @@ The
 .Fn log1pf
 function is a single precision version of
 .Fn log1p .
+The
+.Fn log1pl
+function is an extended precision version of
+.Fn log1p .
 .Pp
 The
 .Fn pow
@@ -170,6 +215,10 @@ The
 .Fn powf
 function is a single precision version of
 .Fn pow .
+The
+.Fn powl
+function is an extended precision version of
+.Fn pow .
 .Sh RETURN VALUES
 These functions will return the appropriate computation unless an error
 occurs or an argument is out of range.
diff --git a/lib/libm/man/floor.3 b/lib/libm/man/floor.3
index 6f6972d2783..a62dac42f05 100644
--- a/lib/libm/man/floor.3
+++ b/lib/libm/man/floor.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: floor.3,v 1.10 2009/08/03 21:26:07 martynas Exp $
+.\"	$OpenBSD: floor.3,v 1.11 2011/07/07 00:54:16 martynas Exp $
 .\" Copyright (c) 1985, 1991 The Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)floor.3	6.5 (Berkeley) 4/19/91
 .\"
-.Dd $Mdocdate: August 3 2009 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt FLOOR 3
 .Os
 .Sh NAME
 .Nm floor ,
-.Nm floorf
+.Nm floorf ,
+.Nm floorl
 .Nd round to largest integral value not greater than x
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -41,6 +42,8 @@
 .Fn floor "double x"
 .Ft float
 .Fn floorf "float x"
+.Ft long double
+.Fn floorl "long double x"
 .Sh DESCRIPTION
 The
 .Fn floor
@@ -51,6 +54,10 @@ The
 .Fn floorf
 function is a single precision version of
 .Fn floor .
+The
+.Fn floorl
+function is an extended precision version of
+.Fn floor .
 .Sh SEE ALSO
 .Xr abs 3 ,
 .Xr ceil 3 ,
diff --git a/lib/libm/man/fmod.3 b/lib/libm/man/fmod.3
index 0f1232bba5e..ff8ec95d286 100644
--- a/lib/libm/man/fmod.3
+++ b/lib/libm/man/fmod.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: fmod.3,v 1.8 2007/05/31 19:19:35 jmc Exp $
+.\"	$OpenBSD: fmod.3,v 1.9 2011/07/07 00:54:16 martynas Exp $
 .\" Copyright (c) 1991 The Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)fmod.3	5.1 (Berkeley) 5/2/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt FMOD 3
 .Os
 .Sh NAME
 .Nm fmod ,
-.Nm fmodf
+.Nm fmodf ,
+.Nm fmodl
 .Nd floating-point remainder functions
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -41,6 +42,8 @@
 .Fn fmod "double x" "double y"
 .Ft float
 .Fn fmodf "float x" "float y"
+.Ft long double
+.Fn fmodl "long double x" "long double y"
 .Sh DESCRIPTION
 The
 .Fn fmod
@@ -50,11 +53,16 @@ The
 .Fn fmodf
 function is a single precision version of
 .Fn fmod .
+The
+.Fn fmodl
+function is an extended precision version of
+.Fn fmod .
 .Sh RETURN VALUES
 The
-.Fn fmod
-and
+.Fn fmod ,
 .Fn fmodf
+and
+.Fn fmodl
 functions return the value
 .Sm off
 .Fa x - Em i * Fa y ,
diff --git a/lib/libm/man/hypot.3 b/lib/libm/man/hypot.3
index 0c0de8a39c0..c853e6bb9a2 100644
--- a/lib/libm/man/hypot.3
+++ b/lib/libm/man/hypot.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: hypot.3,v 1.15 2007/05/31 19:19:35 jmc Exp $
+.\"	$OpenBSD: hypot.3,v 1.16 2011/07/07 00:54:16 martynas Exp $
 .\" Copyright (c) 1985, 1991 Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)hypot.3	6.7 (Berkeley) 5/6/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt HYPOT 3
 .Os
 .Sh NAME
 .Nm hypot ,
 .Nm hypotf ,
+.Nm hypotl ,
 .Nm cabs ,
 .Nm cabsf
 .Nd Euclidean distance and complex absolute value functions
@@ -43,6 +44,8 @@
 .Fn hypot "double x" "double y"
 .Ft float
 .Fn hypotf "float x" "float y"
+.Ft long double
+.Fn hypotl "long double x" "long double y"
 .Ft double
 .Fn cabs "struct complex { double x; double y; } z"
 .Ft float
@@ -66,6 +69,10 @@ functions are single precision versions of
 and
 .Fn cabs ,
 respectively.
+The
+.Fn hypotl
+function is an extended precision version of
+.Fn hypot .
 .Pp
 .Fn hypot "\*(If" "v"
 =
diff --git a/lib/libm/man/lgamma.3 b/lib/libm/man/lgamma.3
index ae00e804eaf..9fb55417ee5 100644
--- a/lib/libm/man/lgamma.3
+++ b/lib/libm/man/lgamma.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: lgamma.3,v 1.16 2008/06/12 15:51:42 jmc Exp $
+.\"	$OpenBSD: lgamma.3,v 1.17 2011/07/07 00:54:16 martynas Exp $
 .\" Copyright (c) 1985, 1991 Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,14 +28,16 @@
 .\"
 .\"     from: @(#)lgamma.3	6.6 (Berkeley) 12/3/92
 .\"
-.Dd $Mdocdate: June 12 2008 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt LGAMMA 3
 .Os
 .Sh NAME
 .Nm lgamma ,
 .Nm lgammaf ,
+.Nm lgammal ,
 .Nm tgamma ,
-.Nm tgammaf
+.Nm tgammaf ,
+.Nm tgammal
 .Nd log gamma functions
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -46,10 +48,14 @@
 .Fn lgamma "double x"
 .Ft float
 .Fn lgammaf "float x"
+.Ft long double
+.Fn lgammal "long double x"
 .Ft double
 .Fn tgamma "double x"
 .Ft float
 .Fn tgammaf "float x"
+.Ft long double
+.Fn tgammal "long double x"
 .Sh DESCRIPTION
 .Fn lgamma x
 .if t \{\
@@ -70,11 +76,16 @@ The
 .Fn lgammaf
 function is a single precision version of
 .Fn lgamma .
+The
+.Fn lgammal
+function is an extended precision version of
+.Fn lgamma .
 .Pp
 The
-.Fn tgamma x
-and
+.Fn tgamma x ,
 .Fn tgammaf x
+and
+.Fn tgammal x
 functions return \(*G(x), with no effect on
 .Fa signgam .
 .Sh IDIOSYNCRASIES
@@ -116,9 +127,11 @@ is set to
 The
 .Fn lgamma ,
 .Fn lgammaf ,
+.Fn lgammal ,
 .Fn tgamma ,
+.Fn tgammaf ,
 and
-.Fn tgammaf
+.Fn tgammal
 functions are expected to conform to
 .St -isoC-99 .
 .Pp
diff --git a/lib/libm/man/logb.3 b/lib/libm/man/logb.3
index 950e77f1a2b..c4e8db32916 100644
--- a/lib/libm/man/logb.3
+++ b/lib/libm/man/logb.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: logb.3,v 1.6 2009/08/03 21:26:07 martynas Exp $
+.\"	$OpenBSD: logb.3,v 1.7 2011/07/07 00:54:16 martynas Exp $
 .\" Copyright (c) 1985, 1991 Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,7 +28,7 @@
 .\"
 .\"     from: @(#)ieee.3	6.4 (Berkeley) 5/6/91
 .\"
-.Dd $Mdocdate: August 3 2009 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt LOGB 3
 .Os
 .Sh NAME
@@ -37,6 +37,7 @@
 .Nm logbl ,
 .Nm scalb ,
 .Nm scalbf ,
+.Nm scalbl ,
 .Nm significand ,
 .Nm significandf
 .Nd IEEE test functions
@@ -52,6 +53,8 @@
 .Fn scalb "double x" "double n"
 .Ft float
 .Fn scalbf "float x" "float n"
+.Ft long double
+.Fn scalbl "long double x" "long double n"
 .Ft double
 .Fn significand "double x"
 .Ft float
@@ -84,6 +87,9 @@ computed by exponent manipulation.
 .Fn scalbf
 is the single precision form of
 .Fn scalb .
+.Fn scalbl
+is the extended precision form of
+.Fn scalb .
 .Pp
 .Fn significand x
 returns
diff --git a/lib/libm/man/lrint.3 b/lib/libm/man/lrint.3
index 0e34d73b1aa..a8632138f1e 100644
--- a/lib/libm/man/lrint.3
+++ b/lib/libm/man/lrint.3
@@ -1,4 +1,4 @@
-.\" 	$OpenBSD: lrint.3,v 1.2 2007/05/31 19:19:35 jmc Exp $
+.\" 	$OpenBSD: lrint.3,v 1.3 2011/07/07 00:54:16 martynas Exp $
 .\" $NetBSD: lrint.3,v 1.1 2005/09/16 15:26:47 wiz Exp $
 .\"
 .\" Copyright (c) 2005 David Schultz <das@FreeBSD.org>
@@ -27,14 +27,16 @@
 .\"
 .\" $FreeBSD: /repoman/r/ncvs/src/lib/msun/man/lrint.3,v 1.2.2.2 2005/03/01 16:18:39 brueffer Exp $
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt LRINT 3
 .Os
 .Sh NAME
 .Nm llrint ,
 .Nm llrintf ,
+.Nm llrintl ,
 .Nm lrint ,
-.Nm lrintf
+.Nm lrintf ,
+.Nm lrintl
 .Nd convert to integer
 .Sh SYNOPSIS
 .In math.h
@@ -42,10 +44,14 @@
 .Fn llrint "double x"
 .Ft long long
 .Fn llrintf "float x"
+.Ft long long
+.Fn llrintl "long double x"
 .Ft long
 .Fn lrint "double x"
 .Ft long
 .Fn lrintf "float x"
+.Ft long
+.Fn lrintl "long double x"
 .Sh DESCRIPTION
 The
 .Fn lrint
@@ -56,8 +62,10 @@ according to the current rounding mode.
 The
 .Fn llrint ,
 .Fn llrintf ,
+.Fn llrintl ,
+.Fn lrintf ,
 and
-.Fn lrintf
+.Fn lrintl
 functions differ from
 .Fn lrint
 only in their input and output types.
@@ -65,9 +73,11 @@ only in their input and output types.
 The
 .Nm llrint ,
 .Nm llrintf ,
+.Nm llrintl ,
 .Nm lrint ,
+.Nm lrintf ,
 and
-.Nm lrintf
+.Nm lrintl
 functions return the integer nearest to their argument
 .Fa x
 according to the current rounding mode.
@@ -94,8 +104,10 @@ the return value is undefined.
 The
 .Fn llrint ,
 .Fn llrintf ,
+.Fn llrintl ,
 .Fn lrint ,
+.Fn lrintf ,
 and
-.Fn lrintf
+.Fn lrintl
 functions conform to
 .St -isoC-99 .
diff --git a/lib/libm/man/lround.3 b/lib/libm/man/lround.3
index 0a2ce7f6cd3..e3feaf16a8f 100644
--- a/lib/libm/man/lround.3
+++ b/lib/libm/man/lround.3
@@ -1,4 +1,4 @@
-.\" 	$OpenBSD: lround.3,v 1.2 2007/05/31 19:19:36 jmc Exp $
+.\" 	$OpenBSD: lround.3,v 1.3 2011/07/07 00:54:16 martynas Exp $
 .\"
 .\" Copyright (c) 2005 David Schultz <das@FreeBSD.org>
 .\" All rights reserved.
@@ -26,16 +26,16 @@
 .\"
 .\" $FreeBSD: /repoman/r/ncvs/src/lib/msun/man/lround.3,v 1.4 2005/06/15 19:04:04 ru Exp $
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt LROUND 3
 .Os
 .Sh NAME
 .Nm llround ,
 .Nm llroundf ,
-.\" .Nm llroundl ,
+.Nm llroundl ,
 .Nm lround ,
-.Nm lroundf
-.\" .Nm lroundl
+.Nm lroundf ,
+.Nm lroundl
 .Nd "convert to nearest integral value"
 .Sh SYNOPSIS
 .In math.h
@@ -43,14 +43,14 @@
 .Fn llround "double x"
 .Ft "long long"
 .Fn llroundf "float x"
-.\" .Ft "long long"
-.\" .Fn llroundl "long double x"
+.Ft "long long"
+.Fn llroundl "long double x"
 .Ft long
 .Fn lround "double x"
 .Ft long
 .Fn lroundf "float x"
-.\" .Ft long
-.\" .Fn lroundl "long double x"
+.Ft long
+.Fn lroundl "long double x"
 .Sh DESCRIPTION
 The
 .Fn lround
@@ -76,10 +76,10 @@ is equivalent to
 The
 .Fn llround ,
 .Fn llroundf ,
-.\" .Fn llroundl ,
-.\" .Fn lroundl
-and
+.Fn llroundl ,
 .Fn lroundf
+and
+.Fn lroundl
 functions differ from
 .Fn lround
 only in their input and output types.
@@ -91,10 +91,10 @@ only in their input and output types.
 The
 .Fn llround ,
 .Fn llroundf ,
-.\" .Fn lroundl
-.\" .Fn llroundl ,
+.Fn llroundl ,
 .Fn lround ,
+.Fn lroundf ,
 and
-.Fn lroundf
+.Fn lroundl
 functions conform to
 .St -isoC-99 .
diff --git a/lib/libm/man/nextafter.3 b/lib/libm/man/nextafter.3
index 83352cc9baf..5406bce0556 100644
--- a/lib/libm/man/nextafter.3
+++ b/lib/libm/man/nextafter.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: nextafter.3,v 1.3 2011/04/16 22:01:00 martynas Exp $
+.\"	$OpenBSD: nextafter.3,v 1.4 2011/07/07 00:54:16 martynas Exp $
 .\" Copyright (c) 1985, 1991 Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)ieee.3	6.4 (Berkeley) 5/6/91
 .\"
-.Dd $Mdocdate: April 16 2011 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt NEXTAFTER 3
 .Os
 .Sh NAME
 .Nm nextafter ,
-.Nm nextafterf
+.Nm nextafterf ,
+.Nm nextafterl
 .Nd next representable value
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -41,6 +42,8 @@
 .Fn nextafter "double x" "double y"
 .Ft float
 .Fn nextafterf "float x" "float y"
+.Ft long double
+.Fn nextafterl "long double x" "long double y"
 .Sh DESCRIPTION
 .Fn nextafter
 returns the next machine representable number from
@@ -58,6 +61,10 @@ The
 .Fn nextafterf
 function is a single precision version of
 .Fn nextafter .
+The
+.Fn nextafterl
+function is an extended precision version of
+.Fn nextafter .
 .Sh SEE ALSO
 .Xr math 3
 .Sh STANDARDS
diff --git a/lib/libm/man/remainder.3 b/lib/libm/man/remainder.3
index 3a84e5c4f80..8fdf7d8194e 100644
--- a/lib/libm/man/remainder.3
+++ b/lib/libm/man/remainder.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: remainder.3,v 1.4 2010/03/26 19:30:40 jmc Exp $
+.\"	$OpenBSD: remainder.3,v 1.5 2011/07/07 00:54:16 martynas Exp $
 .\"
 .\" Copyright (c) 1985, 1991 Regents of the University of California.
 .\" All rights reserved.
@@ -30,14 +30,16 @@
 .\"     from: @(#)ieee.3	6.4 (Berkeley) 5/6/91
 .\" $FreeBSD: src/lib/msun/man/remainder.3,v 1.6 2008/03/30 20:48:02 das Exp $
 .\"
-.Dd $Mdocdate: March 26 2010 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt REMAINDER 3
 .Os
 .Sh NAME
 .Nm remainder ,
 .Nm remainderf ,
+.Nm remainderl ,
 .Nm remquo ,
-.Nm remquof
+.Nm remquof ,
+.Nm remquol
 .Nd minimal residue functions
 .Sh SYNOPSIS
 .In math.h
@@ -45,16 +47,22 @@
 .Fn remainder "double x" "double y"
 .Ft float
 .Fn remainderf "float x" "float y"
+.Ft long double
+.Fn remainderl "long double x" "long double y"
 .Ft double
 .Fn remquo "double x" "double y" "int *quo"
 .Ft float
 .Fn remquof "float x" "float y" "int *quo"
+.Ft long double
+.Fn remquol "long double x" "long double y" "int *quo"
 .Sh DESCRIPTION
 .Fn remainder ,
 .Fn remainderf ,
+.Fn remainderl ,
 .Fn remquo ,
+.Fn remquof ,
 and
-.Fn remquof
+.Fn remquol
 return the remainder
 .Fa r
 :=
@@ -94,9 +102,10 @@ is 0 or
 is \*(Pm\*(If is an invalid operation that produces a \*(Na.
 .Pp
 The
-.Fn remquo
-and
+.Fn remquo ,
 .Fn remquof
+and
+.Fn remquol
 functions also store the last
 .Va k
 bits of
@@ -117,9 +126,11 @@ is platform-specific, but is guaranteed to be at least 3.
 The
 .Fn remainder ,
 .Fn remainderf ,
+.Fn remainderl ,
 .Fn remquo ,
+.Fn remquof ,
 and
-.Fn remquof
+.Fn remquol
 routines conform to
 .St -isoC-99 .
 The remainder is as defined in
diff --git a/lib/libm/man/round.3 b/lib/libm/man/round.3
index a9023a5cf71..afe352bce0b 100644
--- a/lib/libm/man/round.3
+++ b/lib/libm/man/round.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: round.3,v 1.3 2009/08/03 21:26:07 martynas Exp $
+.\"	$OpenBSD: round.3,v 1.4 2011/07/07 00:54:16 martynas Exp $
 .\" Copyright (c) 2003, Steven G. Kargl
 .\" All rights reserved.
 .\"
@@ -25,13 +25,13 @@
 .\"
 .\" $FreeBSD: src/lib/msun/man/round.3,v 1.6 2005/06/15 19:04:04 ru Exp $
 .\"
-.Dd $Mdocdate: August 3 2009 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt ROUND 3
 .Os
 .Sh NAME
 .Nm round ,
-.Nm roundf
-.\" .Nm roundl
+.Nm roundf ,
+.Nm roundl
 .Nd round to nearest integral value
 .Sh SYNOPSIS
 .In math.h
@@ -39,14 +39,14 @@
 .Fn round "double x"
 .Ft float
 .Fn roundf "float x"
-.\" .Ft "long double"
-.\" .Fn roundl "long double x"
+.Ft long double
+.Fn roundl "long double x"
 .Sh DESCRIPTION
 The
-.Fn round
-and
+.Fn round ,
 .Fn roundf
-.\" .Fn roundl
+and
+.Fn roundl
 functions return the nearest integral value to
 .Fa x ;
 if
diff --git a/lib/libm/man/sinh.3 b/lib/libm/man/sinh.3
index 1aee5d86ae5..977fcef9946 100644
--- a/lib/libm/man/sinh.3
+++ b/lib/libm/man/sinh.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: sinh.3,v 1.11 2007/05/31 19:19:36 jmc Exp $
+.\"	$OpenBSD: sinh.3,v 1.12 2011/07/07 00:54:16 martynas Exp $
 .\" Copyright (c) 1991 The Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -27,12 +27,13 @@
 .\" SUCH DAMAGE.
 .\"
 .\"	from: @(#)sinh.3	6.6 (Berkeley) 4/19/91
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt SINH 3
 .Os
 .Sh NAME
 .Nm sinh ,
-.Nm sinhf
+.Nm sinhf ,
+.Nm sinhl
 .Nd hyperbolic sine functions
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -40,6 +41,8 @@
 .Fn sinh "double x"
 .Ft float
 .Fn sinhf "float x"
+.Ft long double
+.Fn sinhl "long double x"
 .Sh DESCRIPTION
 The
 .Fn sinh
@@ -49,11 +52,16 @@ The
 .Fn sinhf
 function is a single precision version of
 .Fn sinh .
+The
+.Fn sinhl
+function is an extended precision version of
+.Fn sinh .
 .Sh RETURN VALUES
 The
-.Fn sinh
-and
+.Fn sinh ,
 .Fn sinhf
+and
+.Fn sinhl
 functions return the hyperbolic sine value unless
 the magnitude
 of
diff --git a/lib/libm/man/sqrt.3 b/lib/libm/man/sqrt.3
index 7b6ae4d8a69..84379faf1e1 100644
--- a/lib/libm/man/sqrt.3
+++ b/lib/libm/man/sqrt.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: sqrt.3,v 1.11 2008/12/12 00:10:26 martynas Exp $
+.\"	$OpenBSD: sqrt.3,v 1.12 2011/07/07 00:54:16 martynas Exp $
 .\" Copyright (c) 1985, 1991 Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)sqrt.3	6.4 (Berkeley) 5/6/91
 .\"
-.Dd $Mdocdate: December 12 2008 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt SQRT 3
 .Os
 .Sh NAME
 .Nm cbrt ,
 .Nm cbrtf ,
+.Nm cbrtl ,
 .Nm sqrt ,
 .Nm sqrtf ,
 .Nm sqrtl
@@ -44,6 +45,8 @@
 .Fn cbrt "double x"
 .Ft float
 .Fn cbrtf "float x"
+.Ft long double
+.Fn cbrtl "long double x"
 .Ft double
 .Fn sqrt "double x"
 .Ft float
@@ -59,6 +62,10 @@ The
 .Fn cbrtf
 function is a single precision version of
 .Fn cbrt .
+The
+.Fn cbrtl
+function is an extended precision version of
+.Fn cbrt .
 .Pp
 The
 .Fn sqrt
diff --git a/lib/libm/man/tanh.3 b/lib/libm/man/tanh.3
index 97ced3cd016..73d4e89e793 100644
--- a/lib/libm/man/tanh.3
+++ b/lib/libm/man/tanh.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: tanh.3,v 1.8 2007/05/31 19:19:36 jmc Exp $
+.\"	$OpenBSD: tanh.3,v 1.9 2011/07/07 00:54:16 martynas Exp $
 .\" Copyright (c) 1991 The Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)tanh.3	5.1 (Berkeley) 5/2/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt TANH 3
 .Os
 .Sh NAME
 .Nm tanh ,
-.Nm tanhf
+.Nm tanhf ,
+.Nm tanhl
 .Nd hyperbolic tangent functions
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -41,6 +42,8 @@
 .Fn tanh "double x"
 .Ft float
 .Fn tanhf "float x"
+.Ft long double
+.Fn tanhl "long double x"
 .Sh DESCRIPTION
 The
 .Fn tanh
@@ -50,13 +53,18 @@ The
 .Fn tanhf
 function is a single precision version of
 .Fn tanh .
+The
+.Fn tanhl
+function is an extended precision version of
+.Fn tanh .
 For a discussion of error due to roundoff, see
 .Xr math 3 .
 .Sh RETURN VALUES
 The
 .Fn tanh
+.Fn tanhf ,
 and
-.Fn tanhf
+.Fn tanhl
 functions return the hyperbolic tangent value.
 .Sh SEE ALSO
 .Xr acos 3 ,
diff --git a/lib/libm/man/trunc.3 b/lib/libm/man/trunc.3
index 84fdc004c98..142771c5b0b 100644
--- a/lib/libm/man/trunc.3
+++ b/lib/libm/man/trunc.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: trunc.3,v 1.4 2009/07/06 00:09:48 martynas Exp $
+.\"	$OpenBSD: trunc.3,v 1.5 2011/07/07 00:54:16 martynas Exp $
 .\"
 .\" Copyright (c) 2004, 2005 David Schultz <das@FreeBSD.org>
 .\" All rights reserved.
@@ -26,13 +26,13 @@
 .\"
 .\" $FreeBSD: src/lib/msun/man/trunc.3,v 1.3 2005/06/15 19:04:04 ru Exp $
 .\"
-.Dd $Mdocdate: July 6 2009 $
+.Dd $Mdocdate: July 7 2011 $
 .Dt TRUNC 3
 .Os
 .Sh NAME
 .Nm trunc ,
-.Nm truncf
-.\" .Nm truncl
+.Nm truncf ,
+.Nm truncl
 .Nd "nearest integral value with magnitude less than or equal to |x|"
 .Sh SYNOPSIS
 .In math.h
@@ -40,26 +40,26 @@
 .Fn trunc "double x"
 .Ft float
 .Fn truncf "float x"
-.\" .Ft "long double"
-.\" .Fn truncl "long double x"
+.Ft long double
+.Fn truncl "long double x"
 .Sh DESCRIPTION
 The
-.Fn trunc
-and
+.Fn trunc ,
 .Fn truncf
-.\" .Fn truncl
+and
+.Fn truncl
 functions return the nearest integral value with magnitude less than
 or equal to
 .Pf | Fa x Ns | .
 They are equivalent to
-.Fn rint
-and
+.Fn rint ,
 .Fn rintf
-.\" .Fn rintl
+and
+.Fn rintl
 respectively, in the rounding towards zero mode.
 .Sh SEE ALSO
 .Xr ceil 3 ,
-.\" .Xr fesetround 3
+.Xr fesetround 3
 .Xr floor 3 ,
 .Xr math 3 ,
 .Xr nextafter 3 ,
@@ -67,9 +67,9 @@ respectively, in the rounding towards zero mode.
 .Xr round 3
 .Sh STANDARDS
 The
-.Fn trunc
-and
+.Fn trunc ,
 .Fn truncf
-.\" .Fn truncl
+and
+.Fn truncl
 functions conform to
 .St -isoC-99 .