X-Git-Url: http://erislabs.net/gitweb/?a=blobdiff_plain;f=lib%2Ftrigl.h;h=79dc38b5f34518b5311e5ad2acc12b0696df92af;hb=3213c107e7632d3b7e462c06acaefde7d50e5067;hp=b4a7a154068c4ed06121cd5f82aeba2a9b03d98a;hpb=5657b2536d242966c41c268b9ab4fe8df727bb13;p=gnulib.git diff --git a/lib/trigl.h b/lib/trigl.h index b4a7a1540..79dc38b5f 100644 --- a/lib/trigl.h +++ b/lib/trigl.h @@ -1,27 +1,35 @@ /* Declarations for sinl, cosl, tanl internal functions Contributed by Paolo Bonzini - Copyright 2002, 2003 Free Software Foundation, Inc. + Copyright 2002, 2003, 2009, 2010 Free Software Foundation, Inc. This file is part of gnulib. - gnulib is free software; you can redistribute it and/or modify it - under the terms of the GNU Lesser General Public License as published - by the Free Software Foundation; either version 2.1, or (at your option) - any later version. + This program is free software: you can redistribute it and/or modify + it under the terms of the GNU General Public License as published by + the Free Software Foundation; either version 3 of the License, or + (at your option) any later version. - gnulib is distributed in the hope that it will be useful, but WITHOUT - ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or - FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public - License for more details. + This program is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the + GNU General Public License for more details. - You should have received a copy of the GNU Lesser General Public License - along with gnulib; see the file COPYING.LIB. If not, write to the Free - Software Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307, - USA. - */ + You should have received a copy of the GNU General Public License + along with this program. If not, see . */ +/* Decompose x into x = k * π/2 + r + where k is an integer and abs(r) <= π/4. + Store r in y[0] and y[1] (main part in y[0], small additional part in + y[1], r = y[0] + y[1]). + Return k. */ extern int ieee754_rem_pio2l (long double x, long double *y); + +/* Compute and return sinl (x + y), where x is the main part and y is the + small additional part of a floating-point number. + iy is 0 when y is known to be 0.0, otherwise iy is 1. */ extern long double kernel_sinl (long double x, long double y, int iy); -extern long double kernel_cosl (long double x, long double y); +/* Compute and return cosl (x + y), where x is the main part and y is the + small additional part of a floating-point number. */ +extern long double kernel_cosl (long double x, long double y);