diff options
Diffstat (limited to 'gsl-1.9/specfunc/test_dilog.c')
| -rw-r--r-- | gsl-1.9/specfunc/test_dilog.c | 310 |
1 files changed, 310 insertions, 0 deletions
diff --git a/gsl-1.9/specfunc/test_dilog.c b/gsl-1.9/specfunc/test_dilog.c new file mode 100644 index 0000000..7e73121 --- /dev/null +++ b/gsl-1.9/specfunc/test_dilog.c @@ -0,0 +1,310 @@ +/* specfunc/test_dilog.c + * + * Copyright (C) 1996, 1997, 1998, 1999, 2000, 2004 Gerard Jungman + * + * 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 2 of the License, or (at + * your option) any later version. + * + * 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 General Public License + * along with this program; if not, write to the Free Software + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. + */ + +/* Author: G. Jungman */ + +#include <config.h> +#include <gsl/gsl_test.h> +#include <gsl/gsl_sf.h> +#include "test_sf.h" + + +int test_dilog(void) +{ + gsl_sf_result r; + gsl_sf_result r1, r2; + int s = 0; + + /* real dilog */ + + TEST_SF(s, gsl_sf_dilog_e, (-3.0, &r), -1.9393754207667089531, TEST_TOL0, GSL_SUCCESS); + TEST_SF(s, gsl_sf_dilog_e, (-0.5, &r), -0.4484142069236462024, TEST_TOL0, GSL_SUCCESS); + TEST_SF(s, gsl_sf_dilog_e, (-0.001, &r), -0.0009997501110486510834, TEST_TOL0, GSL_SUCCESS); + TEST_SF(s, gsl_sf_dilog_e, (0.1, &r), 0.1026177910993911, TEST_TOL0, GSL_SUCCESS); + TEST_SF(s, gsl_sf_dilog_e, (0.7, &r), 0.8893776242860387386, TEST_TOL0, GSL_SUCCESS); + TEST_SF(s, gsl_sf_dilog_e, (1.0, &r), 1.6449340668482260, TEST_TOL0, GSL_SUCCESS); + TEST_SF(s, gsl_sf_dilog_e, (1.5, &r), 2.3743952702724802007, TEST_TOL0, GSL_SUCCESS); + TEST_SF(s, gsl_sf_dilog_e, (2.0, &r), 2.4674011002723397, TEST_TOL0, GSL_SUCCESS); + TEST_SF(s, gsl_sf_dilog_e, ( 5.0, &r), 1.7837191612666306277, TEST_TOL0, GSL_SUCCESS); + TEST_SF(s, gsl_sf_dilog_e, ( 11.0, &r), 0.3218540439999117111, TEST_TOL1, GSL_SUCCESS); + TEST_SF(s, gsl_sf_dilog_e, (12.59, &r), 0.0010060918167266208634, TEST_TOL3, GSL_SUCCESS); + TEST_SF(s, gsl_sf_dilog_e, (12.595, &r), 0.00003314826006436236810, TEST_TOL5, GSL_SUCCESS); + TEST_SF(s, gsl_sf_dilog_e, (13.0, &r), -0.07806971248458575855, TEST_TOL2, GSL_SUCCESS); + TEST_SF(s, gsl_sf_dilog_e, (20.0, &r), -1.2479770861745251168, TEST_TOL2, GSL_SUCCESS); + TEST_SF(s, gsl_sf_dilog_e, (150.0, &r), -9.270042702348657270, TEST_TOL0, GSL_SUCCESS); + TEST_SF(s, gsl_sf_dilog_e, (1100.0, &r), -21.232504073931749553, TEST_TOL0, GSL_SUCCESS); + + + /* complex dilog */ + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.99999, M_PI/2.0, &r1, &r2), + -0.20561329262779687646, TEST_TOL0, + 0.91595774018131512060, TEST_TOL0, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.991, M_PI/2.0, &r1, &r2), + -0.20250384721077806127, TEST_TOL0, + 0.90888544355846447810, TEST_TOL0, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.98, M_PI/2.0, &r1, &r2), + -0.19871638377785918403, TEST_TOL2, + 0.90020045882981847610, TEST_TOL2, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.98, -M_PI/2.0, &r1, &r2), + -0.19871638377785918403, TEST_TOL2, + -0.90020045882981847610, TEST_TOL2, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.95, M_PI/2.0, &r1, &r2), + -0.18848636456893572091, TEST_TOL1, + 0.87633754133420277830, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.8, M_PI/2.0, &r1, &r2), + -0.13980800855429037810, TEST_TOL0, + 0.75310609092419884460, TEST_TOL0, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.8, -M_PI/2.0, &r1, &r2), + -0.13980800855429037810, TEST_TOL0, + -0.75310609092419884460, TEST_TOL0, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.5, M_PI/2.0, &r1, &r2), + -0.05897507442156586346, TEST_TOL1, + 0.48722235829452235710, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.5, -M_PI/2.0, &r1, &r2), + -0.05897507442156586346, TEST_TOL1, + -0.48722235829452235710, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.01, M_PI/2.0, &r1, &r2), + -0.000024999375027776215378, TEST_TOL3, + 0.009999888892888684820, TEST_TOL3, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.01, -M_PI/2.0, &r1, &r2), + -0.000024999375027776215378, TEST_TOL3, + -0.009999888892888684820, TEST_TOL3, + GSL_SUCCESS); + + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.99, M_PI/4.0, &r1, &r2), + 0.56273366219795547757, TEST_TOL3, + 0.97009284079274560384, TEST_TOL3, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.99, -M_PI/4.0, &r1, &r2), + 0.56273366219795547757, TEST_TOL3, + -0.97009284079274560384, TEST_TOL3, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.99, 3.0*M_PI/4.0, &r1, &r2), + -0.66210902664245926235, TEST_TOL1, + 0.51995305609998319025, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.99, 5.0*M_PI/4.0, &r1, &r2), + -0.66210902664245926235, TEST_TOL1, + -0.51995305609998319025, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.99, 3.0*M_PI/2.0, &r1, &r2), + -0.20215874509123277909, TEST_TOL1, + -0.90809733095648731408, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.25, 3.0*M_PI/2.0, &r1, &r2), + -0.01538741178141053563, TEST_TOL1, + -0.24830175098230686908, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.25, 15.0/8.0*M_PI, &r1, &r2), + 0.24266162342377302235, TEST_TOL1, + -0.10860883369274445067, TEST_TOL1, + GSL_SUCCESS); + + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.99, M_PI/8.0, &r1, &r2), + 1.0571539648820244720, TEST_TOL0, + 0.7469145254610851318, TEST_TOL0, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.99, M_PI/64.0, &r1, &r2), + 1.5381800285902999666, TEST_TOL0, + 0.1825271634987756651, TEST_TOL0, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (0.99, -M_PI/8.0, &r1, &r2), + 1.05715396488202447202, TEST_TOL1, + -0.74691452546108513176, TEST_TOL1, + GSL_SUCCESS); + + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (1.00001, M_PI/2.0, &r1, &r2), + -0.20562022409960237363, TEST_TOL1, + 0.91597344814458309320, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (10.0, M_PI/2.0, &r1, &r2), + -3.0596887943287347304, TEST_TOL0, + 3.7167814930680685900, TEST_TOL0, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (100.0, M_PI/2.0, &r1, &r2), + -11.015004738293824854, TEST_TOL0, + 7.2437843013083534970, TEST_TOL0, + GSL_SUCCESS); + + + /** tests brought up by Jim McElwaine bug report */ + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (1.1, -M_PI/2.0, &r1, &r2), + -0.24099184177382733037, TEST_TOL1, + -0.99309132538137822631, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (1.1, 3.0*M_PI/2.0, &r1, &r2), + -0.24099184177382733037, TEST_TOL1, + -0.99309132538137822631, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (1.1, -3.0*M_PI/2.0, &r1, &r2), + -0.24099184177382733037, TEST_TOL1, + 0.99309132538137822631, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (1.1, -M_PI - 0.25*M_PI, &r1, &r2), + -0.72908565537087935118, TEST_TOL1, + 0.56225783937234862649, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (1.1, M_PI + 0.25*M_PI, &r1, &r2), + -0.72908565537087935118, TEST_TOL1, + -0.56225783937234862649, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (1.1, -M_PI/128.0, &r1, &r2), + 1.8881719454909716580, TEST_TOL1, + -0.3556738764969238976, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (1.1, M_PI/128.0, &r1, &r2), + 1.8881719454909716580, TEST_TOL1, + 0.3556738764969238976, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (1.5, M_PI/8.0, &r1, &r2), + 1.3498525763442498343, TEST_TOL1, + 1.4976532712229749493, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (1.5, -M_PI/8.0, &r1, &r2), + 1.3498525763442498343, TEST_TOL1, + -1.4976532712229749493, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (1.5, 2.0*M_PI + M_PI/8.0, &r1, &r2), + 1.3498525763442498343, TEST_TOL1, + 1.4976532712229749493, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_e, (1.5, 2.0*M_PI - M_PI/8.0, &r1, &r2), + 1.3498525763442498343, TEST_TOL1, + -1.4976532712229749493, TEST_TOL1, + GSL_SUCCESS); + + + /* tests of the (x,y) function, which is now the underlying implementation */ + + TEST_SF_2(s, gsl_sf_complex_dilog_xy_e, (0.0, 0.5, &r1, &r2), + -0.05897507442156586346, TEST_TOL1, + 0.48722235829452235710, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_xy_e, (0.0, -0.5, &r1, &r2), + -0.05897507442156586346, TEST_TOL1, + -0.48722235829452235710, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_xy_e, (0.91464073718617389108, 0.37885659804143889673, &r1, &r2), + 1.0571539648820244720, TEST_TOL0, + 0.7469145254610851318, TEST_TOL0, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_xy_e, (0.91464073718617389108, -0.37885659804143889673, &r1, &r2), + 1.05715396488202447202, TEST_TOL1, + -0.74691452546108513176, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_xy_e, (-1.5, 0.0, &r1, &r2), + -1.1473806603755707541, TEST_TOL1, + 0.0, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_xy_e, (0.5, 0.0, &r1, &r2), + 0.58224052646501250590, TEST_TOL1, + 0.0, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_dilog_xy_e, (1.5, 0.0, &r1, &r2), + 2.3743952702724802007, TEST_TOL1, + -1.2738062049196005309, TEST_TOL1, + GSL_SUCCESS); + + + /* small set of spence tests, mostly to check the value on the cut */ + + TEST_SF_2(s, gsl_sf_complex_spence_xy_e, (1.5, 0.0, &r1, &r2), + -0.44841420692364620244, TEST_TOL1, + 0.0, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_spence_xy_e, (0.5, 0.0, &r1, &r2), + 0.58224052646501250590, TEST_TOL1, + 0.0, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_spence_xy_e, (0.0, 0.0, &r1, &r2), + 1.6449340668482264365, TEST_TOL1, + 0.0, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_spence_xy_e, (-0.5, 0.0, &r1, &r2), + 2.3743952702724802007, TEST_TOL1, + -1.2738062049196005309, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_spence_xy_e, (-0.5, 1.0/1024.0, &r1, &r2), + 2.3723507455234125018, TEST_TOL1, + -1.2742581376517839070, TEST_TOL1, + GSL_SUCCESS); + + TEST_SF_2(s, gsl_sf_complex_spence_xy_e, (-0.5, -1.0/1024.0, &r1, &r2), + 2.3723507455234125018, TEST_TOL1, + 1.2742581376517839070, TEST_TOL1, + GSL_SUCCESS); + + + return s; +} |