Home | History | Annotate | Line # | Download | only in complex
ctanl.c revision 1.1.2.2 
      1  1.1.2.2  martin /* $NetBSD: ctanl.c,v 1.1.2.2 2014/10/13 19:34:58 martin Exp $ */
      2  1.1.2.2  martin 
      3  1.1.2.2  martin /*-
      4  1.1.2.2  martin  * Copyright (c) 2007 The NetBSD Foundation, Inc.
      5  1.1.2.2  martin  * All rights reserved.
      6  1.1.2.2  martin  *
      7  1.1.2.2  martin  * This code is derived from software written by Stephen L. Moshier.
      8  1.1.2.2  martin  * It is redistributed by the NetBSD Foundation by permission of the author.
      9  1.1.2.2  martin  *
     10  1.1.2.2  martin  * Redistribution and use in source and binary forms, with or without
     11  1.1.2.2  martin  * modification, are permitted provided that the following conditions
     12  1.1.2.2  martin  * are met:
     13  1.1.2.2  martin  * 1. Redistributions of source code must retain the above copyright
     14  1.1.2.2  martin  *    notice, this list of conditions and the following disclaimer.
     15  1.1.2.2  martin  * 2. Redistributions in binary form must reproduce the above copyright
     16  1.1.2.2  martin  *    notice, this list of conditions and the following disclaimer in the
     17  1.1.2.2  martin  *    documentation and/or other materials provided with the distribution.
     18  1.1.2.2  martin  *
     19  1.1.2.2  martin  * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS
     20  1.1.2.2  martin  * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
     21  1.1.2.2  martin  * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
     22  1.1.2.2  martin  * PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS
     23  1.1.2.2  martin  * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
     24  1.1.2.2  martin  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
     25  1.1.2.2  martin  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
     26  1.1.2.2  martin  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
     27  1.1.2.2  martin  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
     28  1.1.2.2  martin  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
     29  1.1.2.2  martin  * POSSIBILITY OF SUCH DAMAGE.
     30  1.1.2.2  martin  */
     31  1.1.2.2  martin 
     32  1.1.2.2  martin #include "../src/namespace.h"
     33  1.1.2.2  martin #include <complex.h>
     34  1.1.2.2  martin #include <math.h>
     35  1.1.2.2  martin #include <float.h>
     36  1.1.2.2  martin #include "cephes_subrl.h"
     37  1.1.2.2  martin 
     38  1.1.2.2  martin #define MAXNUM LDBL_MAX
     39  1.1.2.2  martin 
     40  1.1.2.2  martin long double complex
     41  1.1.2.2  martin ctanl(long double complex z)
     42  1.1.2.2  martin {
     43  1.1.2.2  martin 	long double complex w;
     44  1.1.2.2  martin 	long double d;
     45  1.1.2.2  martin 
     46  1.1.2.2  martin 	d = cosl(2.0L * creall(z)) + coshl(2.0L * cimagl(z));
     47  1.1.2.2  martin 
     48  1.1.2.2  martin 	if (fabsl(d) < 0.25L)
     49  1.1.2.2  martin 		d = _ctansl(z);
     50  1.1.2.2  martin 
     51  1.1.2.2  martin 	if (d == 0.0L) {
     52  1.1.2.2  martin 		/* mtherr ("ctan", OVERFLOW); */
     53  1.1.2.2  martin 		w = MAXNUM + MAXNUM * I;
     54  1.1.2.2  martin 		return w;
     55  1.1.2.2  martin 	}
     56  1.1.2.2  martin 
     57  1.1.2.2  martin 	w = sinl(2.0L * creall(z)) / d + (sinhl(2.0L * cimagl(z)) / d) * I;
     58  1.1.2.2  martin 	return w;
     59  1.1.2.2  martin }
     60