include/gsl/gsl_sf_coupling.h

0001 /* specfunc/gsl_sf_coupling.h
0002  *
0003  * Copyright (C) 1996,1997,1998,1999,2000,2001,2002 Gerard Jungman
0004  *
0005  * This program is free software; you can redistribute it and/or modify
0006  * it under the terms of the GNU General Public License as published by
0007  * the Free Software Foundation; either version 3 of the License, or (at
0008  * your option) any later version.
0009  *
0010  * This program is distributed in the hope that it will be useful, but
0011  * WITHOUT ANY WARRANTY; without even the implied warranty of
0012  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
0013  * General Public License for more details.
0014  *
0015  * You should have received a copy of the GNU General Public License
0016  * along with this program; if not, write to the Free Software
0017  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
0018  */
0019
0020 /* Author:  G. Jungman */
0021
0022 #ifndef __GSL_SF_COUPLING_H__
0023 #define __GSL_SF_COUPLING_H__
0024
0025 #include <gsl/gsl_sf_result.h>
0026
0027 #undef __BEGIN_DECLS
0028 #undef __END_DECLS
0029 #ifdef __cplusplus
0030 # define __BEGIN_DECLS extern "C" {
0031 # define __END_DECLS }
0032 #else
0033 # define __BEGIN_DECLS /* empty */
0034 # define __END_DECLS /* empty */
0035 #endif
0036
0037 __BEGIN_DECLS
0038
0039
0040 /* 3j Symbols:  / ja jb jc \
0041  *              \ ma mb mc /
0042  *
0043  * exceptions: GSL_EDOM, GSL_EOVRFLW
0044  */
0045 int gsl_sf_coupling_3j_e(int two_ja, int two_jb, int two_jc,
0046                             int two_ma, int two_mb, int two_mc,
0047                             gsl_sf_result * result
0048                             );
0049 double gsl_sf_coupling_3j(int two_ja, int two_jb, int two_jc,
0050                           int two_ma, int two_mb, int two_mc
0051                           );
0052
0053
0054 /* 6j Symbols:  / ja jb jc \
0055  *              \ jd je jf /
0056  *
0057  * exceptions: GSL_EDOM, GSL_EOVRFLW
0058  */
0059 int gsl_sf_coupling_6j_e(int two_ja, int two_jb, int two_jc,
0060                          int two_jd, int two_je, int two_jf,
0061                          gsl_sf_result * result
0062                          );
0063 double gsl_sf_coupling_6j(int two_ja, int two_jb, int two_jc,
0064                           int two_jd, int two_je, int two_jf
0065                           );
0066
0067 /* Racah W coefficients:
0068  *
0069  *   W(a b c d; e f) = (-1)^{a+b+c+d} / a b e \
0070  *                                    \ d c f /
0071  *
0072  * exceptions: GSL_EDOM, GSL_EOVRFLW
0073  */
0074 int gsl_sf_coupling_RacahW_e(int two_ja, int two_jb, int two_jc,
0075                              int two_jd, int two_je, int two_jf,
0076                              gsl_sf_result * result
0077                              );
0078 double gsl_sf_coupling_RacahW(int two_ja, int two_jb, int two_jc,
0079                               int two_jd, int two_je, int two_jf
0080                               );
0081
0082
0083 /* 9j Symbols:  / ja jb jc \
0084  *              | jd je jf |
0085  *              \ jg jh ji /
0086  *
0087  * exceptions: GSL_EDOM, GSL_EOVRFLW
0088  */
0089 int gsl_sf_coupling_9j_e(int two_ja, int two_jb, int two_jc,
0090                          int two_jd, int two_je, int two_jf,
0091                          int two_jg, int two_jh, int two_ji,
0092                          gsl_sf_result * result
0093                          );
0094 double gsl_sf_coupling_9j(int two_ja, int two_jb, int two_jc,
0095                           int two_jd, int two_je, int two_jf,
0096                           int two_jg, int two_jh, int two_ji
0097                           );
0098
0099
0100 /* INCORRECT version of 6j Symbols:
0101  * This function actually calculates
0102  *              / ja jb je \
0103  *              \ jd jc jf /
0104  * It represents the original implementation,
0105  * which had the above permutation of the
0106  * arguments. This was wrong and confusing,
0107  * and I had to fix it. Sorry for the trouble.
0108  * [GJ] Tue Nov 26 12:53:39 MST 2002
0109  *
0110  * exceptions: GSL_EDOM, GSL_EOVRFLW
0111  */
0112 #ifndef GSL_DISABLE_DEPRECATED
0113 int gsl_sf_coupling_6j_INCORRECT_e(int two_ja, int two_jb, int two_jc,
0114                                    int two_jd, int two_je, int two_jf,
0115                                    gsl_sf_result * result
0116                                    );
0117 double gsl_sf_coupling_6j_INCORRECT(int two_ja, int two_jb, int two_jc,
0118                                     int two_jd, int two_je, int two_jf
0119                                     );
0120 #endif /* !GSL_DISABLE_DEPRECATED */
0121
0122
0123 __END_DECLS
0124
0125 #endif /* __GSL_SF_COUPLING_H__ */