document and mlink long double functions. ok jmc@

15 files changed, 196 insertions, 61 deletions

diff --git a/lib/libm/Makefile b/lib/libm/Makefile
index 5cbdedb47fc..b13ae38a514 100644
--- a/lib/libm/Makefile
+++ b/lib/libm/Makefile
@@ -1,4 +1,4 @@
-#  $OpenBSD: Makefile,v 1.60 2008/12/09 23:13:31 martynas Exp $
+#  $OpenBSD: Makefile,v 1.61 2008/12/12 00:10:26 martynas Exp $
 #  $NetBSD: Makefile,v 1.28 1995/11/20 22:06:19 jtc Exp $
 #
 #  @(#)Makefile 5.1beta 93/09/24 
@@ -236,6 +236,24 @@ MLINKS+=tan.3 tanf.3
 MLINKS+=tanh.3 tanhf.3
 MLINKS+=trunc.3 truncf.3
 
+# long double versions
+MLINKS+=acos.3 acosl.3
+MLINKS+=asin.3 asinl.3
+MLINKS+=atan.3 atanl.3
+MLINKS+=atan2.3 atan2l.3
+MLINKS+=cos.3 cosl.3
+MLINKS+=exp.3 exp2l.3
+MLINKS+=fabs.3 fabsl.3
+MLINKS+=ieee.3 copysignl.3
+MLINKS+=ieee.3 ilogbl.3
+MLINKS+=ieee.3 scalbnl.3
+MLINKS+=logb.3 logbl.3
+MLINKS+=nan.3 nanl.3
+MLINKS+=rint.3 rintl.3
+MLINKS+=sin.3 sinl.3
+MLINKS+=sqrt.3 sqrtl.3
+MLINKS+=tan.3 tanl.3
+
 #.if (${MACHINE_ARCH} == "i386")
 #LINKS=	${LIBDIR}/libm.a ${LIBDIR}/libm_i387.a
 #.if !defined(NOPROFILE)
diff --git a/lib/libm/man/acos.3 b/lib/libm/man/acos.3
index 67f895dc811..678becf86d8 100644
--- a/lib/libm/man/acos.3
+++ b/lib/libm/man/acos.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: acos.3,v 1.11 2007/05/31 19:19:35 jmc Exp $
+.\"	$OpenBSD: acos.3,v 1.12 2008/12/12 00:10:26 martynas Exp $
 .\" Copyright (c) 1991 The Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)acos.3	5.1 (Berkeley) 5/2/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: December 12 2008 $
 .Dt ACOS 3
 .Os
 .Sh NAME
 .Nm acos ,
-.Nm acosf
+.Nm acosf ,
+.Nm acosl
 .Nd arc cosine functions
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -41,6 +42,8 @@
 .Fn acos "double x"
 .Ft float
 .Fn acosf "float x"
+.Ft long double
+.Fn acosl "long double x"
 .Sh DESCRIPTION
 The
 .Fn acos
@@ -52,11 +55,16 @@ The
 .Fn acosf
 is a single precision version of
 .Fn acos .
+The
+.Fn acosl
+is an extended precision version of
+.Fn acos .
 .Sh RETURN VALUES
 If |x|>1,
-.Fn acos "x"
-and
+.Fn acos "x" ,
 .Fn acosf "x"
+and
+.Fn acosl "x"
 set the global variable
 .Va errno
 to EDOM.
diff --git a/lib/libm/man/asin.3 b/lib/libm/man/asin.3
index b661246123e..5c247e318e3 100644
--- a/lib/libm/man/asin.3
+++ b/lib/libm/man/asin.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: asin.3,v 1.11 2007/05/31 19:19:35 jmc Exp $
+.\"	$OpenBSD: asin.3,v 1.12 2008/12/12 00:10:26 martynas Exp $
 .\" Copyright (c) 1991 The Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)asin.3	5.1 (Berkeley) 5/2/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: December 12 2008 $
 .Dt ASIN 3
 .Os
 .Sh NAME
 .Nm asin ,
-.Nm asinf
+.Nm asinf ,
+.Nm asinl
 .Nd arc sine functions
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -41,6 +42,8 @@
 .Fn asin "double x"
 .Ft float
 .Fn asinf "float x"
+.Ft long double
+.Fn asinl "long double x"
 .Sh DESCRIPTION
 The
 .Fn asin
@@ -54,11 +57,16 @@ The
 .Fn asinf
 function is a single precision version of
 .Fn asin .
+The
+.Fn asinl
+function is an extended precision version of
+.Fn asin .
 .Sh RETURN VALUES
 If |x|>1,
-.Fn asin "x"
-and
+.Fn asin "x" ,
 .Fn asinf "x"
+and
+.Fn asinl "x"
 return NaN and set the global variable
 .Va errno
 to EDOM.
diff --git a/lib/libm/man/atan.3 b/lib/libm/man/atan.3
index 7be5974fe98..1976b7c463c 100644
--- a/lib/libm/man/atan.3
+++ b/lib/libm/man/atan.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: atan.3,v 1.10 2007/05/31 19:19:35 jmc Exp $
+.\"	$OpenBSD: atan.3,v 1.11 2008/12/12 00:10:26 martynas Exp $
 .\" Copyright (c) 1991 The Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)atan.3	5.1 (Berkeley) 5/2/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: December 12 2008 $
 .Dt ATAN 3
 .Os
 .Sh NAME
 .Nm atan ,
-.Nm atanf
+.Nm atanf ,
+.Nm atanl
 .Nd arc tangent functions of one variable
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -41,6 +42,8 @@
 .Fn atan "double x"
 .Ft float
 .Fn atanf "float x"
+.Ft long double
+.Fn atanl "long double x"
 .Sh DESCRIPTION
 The
 .Fn atan
@@ -54,6 +57,10 @@ The
 .Fn atanf
 function is a single precision version of
 .Fn atan .
+The
+.Fn atanl
+function is an extended precision version of
+.Fn atan .
 .Sh SEE ALSO
 .Xr acos 3 ,
 .Xr asin 3 ,
diff --git a/lib/libm/man/atan2.3 b/lib/libm/man/atan2.3
index 795db059a1d..da783d74613 100644
--- a/lib/libm/man/atan2.3
+++ b/lib/libm/man/atan2.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: atan2.3,v 1.11 2007/05/31 19:19:35 jmc Exp $
+.\"	$OpenBSD: atan2.3,v 1.12 2008/12/12 00:10:26 martynas Exp $
 .\" Copyright (c) 1991 The Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)atan2.3	5.1 (Berkeley) 5/2/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: December 12 2008 $
 .Dt ATAN2 3
 .Os
 .Sh NAME
 .Nm atan2 ,
-.Nm atan2f
+.Nm atan2f ,
+.Nm atan2l
 .Nd arc tangent functions of two variables
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -41,6 +42,8 @@
 .Fn atan2 "double y" "double x"
 .Ft float
 .Fn atan2f "float y" "float x"
+.Ft long double
+.Fn atan2l "long double y" "long double x"
 .Sh DESCRIPTION
 The
 .Fn atan2
@@ -52,11 +55,16 @@ The
 .Fn atan2f
 function is a single precision version of
 .Fn atan2 .
+The
+.Fn atan2l
+function is an extended precision version of
+.Fn atan2 .
 .Sh RETURN VALUES
 The
-.Fn atan2
-and
+.Fn atan2 ,
 .Fn atan2f
+and
+.Fn atan2l
 functions, if successful,
 return the arc tangent of
 .Ar y/ Ns Ar x
diff --git a/lib/libm/man/cos.3 b/lib/libm/man/cos.3
index 83f9fe22694..74431c4a700 100644
--- a/lib/libm/man/cos.3
+++ b/lib/libm/man/cos.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: cos.3,v 1.11 2007/05/31 19:19:35 jmc Exp $
+.\"	$OpenBSD: cos.3,v 1.12 2008/12/12 00:10:26 martynas Exp $
 .\" Copyright (c) 1991 The Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)cos.3	5.1 (Berkeley) 5/2/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: December 12 2008 $
 .Dt COS 3
 .Os
 .Sh NAME
 .Nm cos ,
-.Nm cosf
+.Nm cosf ,
+.Nm cosl
 .Nd cosine functions
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -41,6 +42,8 @@
 .Fn cos "double x"
 .Ft float
 .Fn cosf "float x"
+.Ft long double
+.Fn cosl "long double x"
 .Sh DESCRIPTION
 The
 .Fn cos
@@ -51,15 +54,20 @@ The
 .Fn cosf
 function is a single precision version of
 .Fn cos .
+The
+.Fn cosl
+function is an extended precision version of
+.Fn cos .
 A large magnitude argument may yield a result with little or no
 significance.
 For a discussion of error due to roundoff, see
 .Xr math 3 .
 .Sh RETURN VALUES
 The
-.Fn cos
-and
+.Fn cos ,
 .Fn cosf
+and
+.Fn cosl
 functions return the cosine value.
 .Sh SEE ALSO
 .Xr acos 3 ,
diff --git a/lib/libm/man/exp.3 b/lib/libm/man/exp.3
index ac411e2f24c..37ce6a054c3 100644
--- a/lib/libm/man/exp.3
+++ b/lib/libm/man/exp.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: exp.3,v 1.19 2008/07/30 08:02:35 jmc Exp $
+.\"	$OpenBSD: exp.3,v 1.20 2008/12/12 00:10:26 martynas Exp $
 .\"
 .\" Copyright (c) 1985, 1991 Regents of the University of California.
 .\" All rights reserved.
@@ -29,7 +29,7 @@
 .\"
 .\"     from: @(#)exp.3	6.12 (Berkeley) 7/31/91
 .\"
-.Dd $Mdocdate: July 30 2008 $
+.Dd $Mdocdate: December 12 2008 $
 .Dt EXP 3
 .Os
 .Sh NAME
@@ -37,6 +37,7 @@
 .Nm expf ,
 .Nm exp2 ,
 .Nm exp2f ,
+.Nm exp2l ,
 .Nm expm1 ,
 .Nm expm1f ,
 .Nm log ,
@@ -58,6 +59,8 @@
 .Fn exp2 "double x"
 .Ft float
 .Fn exp2f "float x"
+.Ft long double
+.Fn exp2l "long double x"
 .Ft double
 .Fn expm1 "double x"
 .Ft float
@@ -98,6 +101,10 @@ The
 .Fn exp2f
 function is a single precision version of
 .Fn exp2 .
+The
+.Fn exp2l
+function is an extended precision version of
+.Fn exp2 .
 .Pp
 The
 .Fn expm1
diff --git a/lib/libm/man/fabs.3 b/lib/libm/man/fabs.3
index b7428a04140..7232a07ef2a 100644
--- a/lib/libm/man/fabs.3
+++ b/lib/libm/man/fabs.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: fabs.3,v 1.9 2007/05/31 19:19:35 jmc Exp $
+.\"	$OpenBSD: fabs.3,v 1.10 2008/12/12 00:10:26 martynas Exp $
 .\" Copyright (c) 1991 The Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -29,12 +29,13 @@
 .\"
 .\"     from: @(#)fabs.3	5.1 (Berkeley) 5/2/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: December 12 2008 $
 .Dt FABS 3
 .Os
 .Sh NAME
 .Nm fabs ,
-.Nm fabsf
+.Nm fabsf ,
+.Nm fabsl
 .Nd floating-point absolute value functions
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -42,6 +43,8 @@
 .Fn fabs "double x"
 .Ft float
 .Fn fabsf "float x"
+.Ft long double
+.Fn fabsl "long double x"
 .Sh DESCRIPTION
 The
 .Fn fabs
@@ -51,11 +54,16 @@ The
 .Fn fabsf
 function is a single precision version of
 .Fn fabs .
+The
+.Fn fabsl
+function is an extended precision version of
+.Fn fabs .
 .Sh RETURN VALUES
 The
-.Fn fabs
-and
+.Fn fabs ,
 .Fn fabsf
+and
+.Fn fabsl
 functions return the absolute value of
 .Fa x .
 .Sh SEE ALSO
diff --git a/lib/libm/man/ieee.3 b/lib/libm/man/ieee.3
index 3b8219dd38d..27cd81767a9 100644
--- a/lib/libm/man/ieee.3
+++ b/lib/libm/man/ieee.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: ieee.3,v 1.16 2008/07/29 22:08:43 martynas Exp $
+.\"	$OpenBSD: ieee.3,v 1.17 2008/12/12 00:10:26 martynas Exp $
 .\" Copyright (c) 1985, 1991 Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,18 +28,21 @@
 .\"
 .\"     from: @(#)ieee.3	6.4 (Berkeley) 5/6/91
 .\"
-.Dd $Mdocdate: July 29 2008 $
+.Dd $Mdocdate: December 12 2008 $
 .Dt IEEE 3
 .Os
 .Sh NAME
 .Nm copysign ,
 .Nm copysignf ,
+.Nm copysignl ,
 .Nm ilogb ,
 .Nm ilogbf ,
+.Nm ilogbl ,
 .Nm nextafter ,
 .Nm nextafterf ,
 .Nm scalbn ,
-.Nm scalbnf
+.Nm scalbnf ,
+.Nm scalbnl
 .Nd functions for IEEE arithmetic
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -47,10 +50,14 @@
 .Fn copysign "double x" "double y"
 .Ft float
 .Fn copysignf "float x" "float y"
+.Ft long double
+.Fn copysignl "long double x" "long double y"
 .Ft int
 .Fn ilogb "double x"
 .Ft int
 .Fn ilogbf "float x"
+.Ft int
+.Fn ilogbl "long double x"
 .Ft double
 .Fn nextafter "double x" "double y"
 .Ft float
@@ -59,6 +66,8 @@
 .Fn scalbn "double x" "int n"
 .Ft float
 .Fn scalbnf "float x" "int n"
+.Ft long double
+.Fn scalbnl "long double x" "int n"
 .Sh DESCRIPTION
 These functions are required or recommended by
 .St -ieee754 .
@@ -72,6 +81,10 @@ The
 .Fn copysignf
 function is a single precision version of
 .Fn copysign .
+The
+.Fn copysignl
+function is an extended precision version of
+.Fn copysign .
 .Pp
 .Fn ilogb
 returns
@@ -89,6 +102,10 @@ The
 .Fn ilogbf
 function is a single precision version of
 .Fn ilogb .
+The
+.Fn ilogbl
+function is an extended precision version of
+.Fn ilogb .
 .Pp
 .Fn nextafter
 returns the next machine representable number from
@@ -108,6 +125,10 @@ The
 .Fn scalbnf
 function is a single precision version of
 .Fn scalbn .
+The
+.Fn scalbnl
+function is an extended precision version of
+.Fn scalbn .
 .Sh SEE ALSO
 .Xr math 3
 .Sh STANDARDS
diff --git a/lib/libm/man/logb.3 b/lib/libm/man/logb.3
index 430c54374bc..6cdd11514e0 100644
--- a/lib/libm/man/logb.3
+++ b/lib/libm/man/logb.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: logb.3,v 1.3 2007/05/31 19:19:35 jmc Exp $
+.\"	$OpenBSD: logb.3,v 1.4 2008/12/12 00:10:26 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: May 31 2007 $
+.Dd $Mdocdate: December 12 2008 $
 .Dt IEEE_TEST 3
 .Os
 .Sh NAME
 .Nm logb ,
 .Nm logbf ,
+.Nm logbl ,
 .Nm scalb ,
 .Nm scalbf ,
 .Nm significand ,
@@ -45,6 +46,8 @@
 .Fn logb "double x"
 .Ft float
 .Fn logbf "float x"
+.Ft long double
+.Fn logbl "long double x"
 .Ft double
 .Fn scalb "double x" "double n"
 .Ft float
@@ -70,6 +73,9 @@ a signed integer converted to double\-precision floating\-point.
 .Fn logbf
 is the single precision form of
 .Fn logb .
+.Fn logbl
+is the extended precision form of
+.Fn logb .
 .Pp
 .Fn scalb x n
 returns
diff --git a/lib/libm/man/nan.3 b/lib/libm/man/nan.3
index 0273361e902..48a10c97741 100644
--- a/lib/libm/man/nan.3
+++ b/lib/libm/man/nan.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: nan.3,v 1.2 2008/07/30 08:02:35 jmc Exp $
+.\"	$OpenBSD: nan.3,v 1.3 2008/12/12 00:10:26 martynas Exp $
 .\"
 .\" Copyright (c) 2007 David Schultz <das@FreeBSD.org>
 .\" All rights reserved.
@@ -26,12 +26,13 @@
 .\"
 .\" $FreeBSD: src/lib/msun/man/nan.3,v 1.1 2007/12/16 21:19:28 das Exp $
 .\"
-.Dd $Mdocdate: July 30 2008 $
+.Dd $Mdocdate: December 12 2008 $
 .Dt NAN 3
 .Os
 .Sh NAME
 .Nm nan ,
-.Nm nanf
+.Nm nanf ,
+.Nm nanl
 .Nd quiet \*(Nas
 .Sh SYNOPSIS
 .In math.h
@@ -39,14 +40,17 @@
 .Fn nan "const char *s"
 .Ft float
 .Fn nanf "const char *s"
+.Ft long double
+.Fn nanl "const char *s"
 .Sh DESCRIPTION
 The
 .Dv NAN
 macro expands to a quiet \*(Na (Not A Number).
 Similarly both the
-.Fn nan
-and
+.Fn nan ,
 .Fn nanf
+and
+.Fn nanl
 functions generate a quiet \*(Na value without raising an invalid exception.
 The argument
 .Fa s
@@ -84,9 +88,10 @@ macro instead.
 .Xr strtod 3
 .Sh STANDARDS
 The
-.Fn nan
-and
+.Fn nan ,
 .Fn nanf
+and
+.Fn nanl
 functions and the
 .Dv NAN
 macro conform to
diff --git a/lib/libm/man/rint.3 b/lib/libm/man/rint.3
index b61023bf77e..0c12cebc0be 100644
--- a/lib/libm/man/rint.3
+++ b/lib/libm/man/rint.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: rint.3,v 1.8 2007/05/31 19:19:36 jmc Exp $
+.\"	$OpenBSD: rint.3,v 1.9 2008/12/12 00:10:26 martynas Exp $
 .\" Copyright (c) 1985, 1991 Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)rint.3	5.1 (Berkeley) 5/2/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: December 12 2008 $
 .Dt RINT 3
 .Os
 .Sh NAME
 .Nm rint ,
-.Nm rintf
+.Nm rintf ,
+.Nm rintl
 .Nd round to integral value in floating-point format
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -41,6 +42,8 @@
 .Fn rint "double x"
 .Ft float
 .Fn rintf "float x"
+.Ft long double
+.Fn rintl "long double x"
 .Sh DESCRIPTION
 The
 .Fn rint
@@ -52,6 +55,10 @@ The
 .Fn rintf
 function is a single precision version of
 .Fn rint .
+The
+.Fn rintl
+function is an extended precision version of
+.Fn rint .
 .Sh SEE ALSO
 .Xr abs 3 ,
 .Xr ceil 3 ,
diff --git a/lib/libm/man/sin.3 b/lib/libm/man/sin.3
index 7837fa26809..5dd252d9b9f 100644
--- a/lib/libm/man/sin.3
+++ b/lib/libm/man/sin.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: sin.3,v 1.11 2007/05/31 19:19:36 jmc Exp $
+.\"	$OpenBSD: sin.3,v 1.12 2008/12/12 00:10:26 martynas Exp $
 .\" Copyright (c) 1991 The Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)cos.3	5.1 (Berkeley) 5/2/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: December 12 2008 $
 .Dt COS 3
 .Os
 .Sh NAME
 .Nm sin ,
-.Nm sinf
+.Nm sinf ,
+.Nm sinl
 .Nd sine functions
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -41,6 +42,8 @@
 .Fn sin "double x"
 .Ft float
 .Fn sinf "float x"
+.Ft long double
+.Fn sinl "long double x"
 .Sh DESCRIPTION
 The
 .Fn sin
@@ -51,15 +54,20 @@ The
 .Fn sinf
 function is a single precision version of
 .Fn sin .
+The
+.Fn sinl
+function is an extended precision version of
+.Fn sin .
 A large magnitude argument may yield a result with little or no
 significance.
 For a discussion of error due to roundoff, see
 .Xr math 3 .
 .Sh RETURN VALUES
 The
-.Fn sin
-and
+.Fn sin ,
 .Fn sinf
+and
+.Fn sinl
 functions return the sine value.
 .Sh SEE ALSO
 .Xr acos 3 ,
diff --git a/lib/libm/man/sqrt.3 b/lib/libm/man/sqrt.3
index 8b6c63fe3eb..7b6ae4d8a69 100644
--- a/lib/libm/man/sqrt.3
+++ b/lib/libm/man/sqrt.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: sqrt.3,v 1.10 2007/05/31 19:19:36 jmc Exp $
+.\"	$OpenBSD: sqrt.3,v 1.11 2008/12/12 00:10:26 martynas Exp $
 .\" Copyright (c) 1985, 1991 Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,14 +28,15 @@
 .\"
 .\"     from: @(#)sqrt.3	6.4 (Berkeley) 5/6/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: December 12 2008 $
 .Dt SQRT 3
 .Os
 .Sh NAME
 .Nm cbrt ,
 .Nm cbrtf ,
 .Nm sqrt ,
-.Nm sqrtf
+.Nm sqrtf ,
+.Nm sqrtl
 .Nd cube root and square root functions
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -47,6 +48,8 @@
 .Fn sqrt "double x"
 .Ft float
 .Fn sqrtf "float x"
+.Ft long double
+.Fn sqrtl "long double x"
 .Sh DESCRIPTION
 The
 .Fn cbrt
@@ -65,11 +68,16 @@ The
 .Fn sqrtf
 function is a single precision version of
 .Fn sqrt .
+The
+.Fn sqrtl
+function is an extended precision version of
+.Fn sqrt .
 .Sh RETURN VALUES
 If x is negative,
-.Fn sqrt "x"
-and
+.Fn sqrt "x" ,
 .Fn sqrtf "x"
+and
+.Fn sqrtl "x"
 set the global variable
 .Va errno
 to EDOM.
diff --git a/lib/libm/man/tan.3 b/lib/libm/man/tan.3
index 585cc90cc61..11b097fc5bb 100644
--- a/lib/libm/man/tan.3
+++ b/lib/libm/man/tan.3
@@ -1,4 +1,4 @@
-.\"	$OpenBSD: tan.3,v 1.9 2007/05/31 19:19:36 jmc Exp $
+.\"	$OpenBSD: tan.3,v 1.10 2008/12/12 00:10:26 martynas Exp $
 .\" Copyright (c) 1991 The Regents of the University of California.
 .\" All rights reserved.
 .\"
@@ -28,12 +28,13 @@
 .\"
 .\"     from: @(#)tan.3	5.1 (Berkeley) 5/2/91
 .\"
-.Dd $Mdocdate: May 31 2007 $
+.Dd $Mdocdate: December 12 2008 $
 .Dt TAN 3
 .Os
 .Sh NAME
 .Nm tan ,
-.Nm tanf
+.Nm tanf ,
+.Nm tanl
 .Nd tangent functions
 .Sh SYNOPSIS
 .Fd #include <math.h>
@@ -41,6 +42,8 @@
 .Fn tan "double x"
 .Ft float
 .Fn tanf "float x"
+.Ft long double
+.Fn tanl "long double x"
 .Sh DESCRIPTION
 The
 .Fn tan
@@ -51,15 +54,20 @@ The
 .Fn tanf
 function is a single precision version of
 .Fn tan .
+The
+.Fn tanl
+function is an extended precision version of
+.Fn tan .
 A large magnitude argument may yield a result
 with little or no significance.
 For a discussion of error due to roundoff, see
 .Xr math 3 .
 .Sh RETURN VALUES
 The
-.Fn tan
-and
+.Fn tan ,
 .Fn tanf
+and
+.Fn tanl
 functions return the tangent value.
 .Sh SEE ALSO
 .Xr acos 3 ,