src/math/Wigner3jm.hpp

/*
 * Matpack Wigner3jm special function imported and modified for use in
 * OpenMD
 *                                                                              
 * Matpack Library Release 1.9.0                                                
 * Copyright (C) 1991-2003 by Berndt M. Gammel. All rights reserved.            
 *
 * Permission to use, copy, and distribute Matpack in its entirety
 * and its documentation for non-commercial purpose and without fee
 * is hereby granted, provided that this license information and
 * copyright notice appear unmodified in all copies.  This software
 * is provided 'as is' without express or implied warranty.  In no
 * event will the author be held liable for any damages arising from
 * the use of this software.
 *
 * Note that distributing Matpack 'bundled' in with any product is
 * considered to be a 'commercial purpose'.
 *
 * The software may be modified for your own purposes, but modified
 * versions may not be distributed without prior consent of the
 * author.
 *                                                                  
 * Read the COPYRIGHT and README files in this distribution about
 * registration and installation of Matpack.
 */

#ifndef MATH_WIGNER3JM_HPP
#define MATH_WIGNER3JM_HPP

#include "config.h"

namespace MATPACK {
    
/// return sign of number
template <class T> inline int sign (T x)
{ return (x > 0) ? 1 : (x < 0) ? -1 : 0; }
  
  /*
   * MpMin(), MpMin(): min and max templates for 2 arguments (all that
   * is required for Wigner3jm )
   */
template <class T> inline T MpMin (T x, T y) 
{ return (x<y)?x:y; }
template <class T> inline T MpMax (T x, T y) 
{ return (x>y)?x:y; }

// even and odd

inline int even(int x){return!(x&1);}
inline int odd(int x){return(x&1);}

void Wigner3jm(RealType l1, RealType l2, RealType l3, RealType m1, 
               RealType &m2min, RealType &m2max, RealType *thrcof, int ndim, 
               int &errflag);
}

#endif
Revision:	1600
Committed:	Wed Aug 3 20:20:37 2011 UTC (13 years, 9 months ago) by gezelter
File size:	1952 byte(s)
Log Message:	Completing the Fortran removal project. Fixes for compilation with clang / llvm, debugging, removing code that we'd never used.
#	User	Rev	Content
1	gezelter	1600	/*
2			* Matpack Wigner3jm special function imported and modified for use in
3			* OpenMD
4			*
5			* Matpack Library Release 1.9.0
6			* Copyright (C) 1991-2003 by Berndt M. Gammel. All rights reserved.
7			*
8			* Permission to use, copy, and distribute Matpack in its entirety
9			* and its documentation for non-commercial purpose and without fee
10			* is hereby granted, provided that this license information and
11			* copyright notice appear unmodified in all copies. This software
12			* is provided 'as is' without express or implied warranty. In no
13			* event will the author be held liable for any damages arising from
14			* the use of this software.
15			*
16			* Note that distributing Matpack 'bundled' in with any product is
17			* considered to be a 'commercial purpose'.
18			*
19			* The software may be modified for your own purposes, but modified
20			* versions may not be distributed without prior consent of the
21			* author.
22			*
23			* Read the COPYRIGHT and README files in this distribution about
24			* registration and installation of Matpack.
25			*/
26
27			#ifndef MATH_WIGNER3JM_HPP
28			#define MATH_WIGNER3JM_HPP
29
30			#include "config.h"
31
32			namespace MATPACK {
33
34			/// return sign of number
35			template <class T> inline int sign (T x)
36			{ return (x > 0) ? 1 : (x < 0) ? -1 : 0; }
37
38			/*
39			* MpMin(), MpMin(): min and max templates for 2 arguments (all that
40			* is required for Wigner3jm )
41			*/
42			template <class T> inline T MpMin (T x, T y)
43			{ return (x<y)?x:y; }
44			template <class T> inline T MpMax (T x, T y)
45			{ return (x>y)?x:y; }
46
47			// even and odd
48
49			inline int even(int x){return!(x&1);}
50			inline int odd(int x){return(x&1);}
51
52			void Wigner3jm(RealType l1, RealType l2, RealType l3, RealType m1,
53			RealType &m2min, RealType &m2max, RealType *thrcof, int ndim,
54			int &errflag);
55			}
56
57			#endif