htswanalysis/MACS/lib/gsl/gsl-1.11/specfunc/gsl_sf_coupling.h

   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  */
  45 int 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                             );
  49 double 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  */
  59 int 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                          );
  63 double 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  */
  74 int 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                              );
  78 double 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  */
  89 int 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                          );
  94 double 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
 113 int 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                                    );
 117 double 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__ */