1/* specfunc/gsl_sf_coupling.h
2 *
3 * Copyright (C) 1996,1997,1998,1999,2000,2001,2002 Gerard Jungman
4 *
5 * This program is free software; you can redistribute it and/or modify
6 * it under the terms of the GNU General Public License as published by
7 * the Free Software Foundation; either version 3 of the License, or (at
8 * your option) any later version.
9 *
10 * This program is distributed in the hope that it will be useful, but
11 * WITHOUT ANY WARRANTY; without even the implied warranty of
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 * General Public License for more details.
14 *
15 * You should have received a copy of the GNU General Public License
16 * along with this program; if not, write to the Free Software
17 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
18 */
19
20/* Author: G. Jungman */
21
22#ifndef __GSL_SF_COUPLING_H__
23#define __GSL_SF_COUPLING_H__
24
25#include <gsl/gsl_sf_result.h>
26
27#undef __BEGIN_DECLS
28#undef __END_DECLS
29#ifdef __cplusplus
30# define __BEGIN_DECLS extern "C" {
31# define __END_DECLS }
32#else
33# define __BEGIN_DECLS /* empty */
34# define __END_DECLS /* empty */
35#endif
36
37__BEGIN_DECLS
38
39
40/* 3j Symbols: / ja jb jc \
41 * \ ma mb mc /
42 *
43 * exceptions: GSL_EDOM, GSL_EOVRFLW
44 */
45int gsl_sf_coupling_3j_e(int two_ja, int two_jb, int two_jc,
46 int two_ma, int two_mb, int two_mc,
47 gsl_sf_result * result
48 );
49double gsl_sf_coupling_3j(int two_ja, int two_jb, int two_jc,
50 int two_ma, int two_mb, int two_mc
51 );
52
53
54/* 6j Symbols: / ja jb jc \
55 * \ jd je jf /
56 *
57 * exceptions: GSL_EDOM, GSL_EOVRFLW
58 */
59int gsl_sf_coupling_6j_e(int two_ja, int two_jb, int two_jc,
60 int two_jd, int two_je, int two_jf,
61 gsl_sf_result * result
62 );
63double gsl_sf_coupling_6j(int two_ja, int two_jb, int two_jc,
64 int two_jd, int two_je, int two_jf
65 );
66
67/* Racah W coefficients:
68 *
69 * W(a b c d; e f) = (-1)^{a+b+c+d} / a b e \
70 * \ d c f /
71 *
72 * exceptions: GSL_EDOM, GSL_EOVRFLW
73 */
74int gsl_sf_coupling_RacahW_e(int two_ja, int two_jb, int two_jc,
75 int two_jd, int two_je, int two_jf,
76 gsl_sf_result * result
77 );
78double gsl_sf_coupling_RacahW(int two_ja, int two_jb, int two_jc,
79 int two_jd, int two_je, int two_jf
80 );
81
82
83/* 9j Symbols: / ja jb jc \
84 * | jd je jf |
85 * \ jg jh ji /
86 *
87 * exceptions: GSL_EDOM, GSL_EOVRFLW
88 */
89int gsl_sf_coupling_9j_e(int two_ja, int two_jb, int two_jc,
90 int two_jd, int two_je, int two_jf,
91 int two_jg, int two_jh, int two_ji,
92 gsl_sf_result * result
93 );
94double gsl_sf_coupling_9j(int two_ja, int two_jb, int two_jc,
95 int two_jd, int two_je, int two_jf,
96 int two_jg, int two_jh, int two_ji
97 );
98
99
100/* INCORRECT version of 6j Symbols:
101 * This function actually calculates
102 * / ja jb je \
103 * \ jd jc jf /
104 * It represents the original implementation,
105 * which had the above permutation of the
106 * arguments. This was wrong and confusing,
107 * and I had to fix it. Sorry for the trouble.
108 * [GJ] Tue Nov 26 12:53:39 MST 2002
109 *
110 * exceptions: GSL_EDOM, GSL_EOVRFLW
111 */
112#ifndef GSL_DISABLE_DEPRECATED
113int gsl_sf_coupling_6j_INCORRECT_e(int two_ja, int two_jb, int two_jc,
114 int two_jd, int two_je, int two_jf,
115 gsl_sf_result * result
116 );
117double gsl_sf_coupling_6j_INCORRECT(int two_ja, int two_jb, int two_jc,
118 int two_jd, int two_je, int two_jf
119 );
120#endif /* !GSL_DISABLE_DEPRECATED */
121
122
123__END_DECLS
124
125#endif /* __GSL_SF_COUPLING_H__ */
126