-
Notifications
You must be signed in to change notification settings - Fork 15
/
gh-issue-1.cpp
82 lines (71 loc) · 4.41 KB
/
gh-issue-1.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
/*******************************************************-/
* This source code is subject to the terms of the GNU -/
* Lesser Public License. If a copy of the LGPL was not -/
* distributed with this file, you can obtain one at -/
* https://www.gnu.org/licenses/lgpl.html -/
********************************************************/
/*! \file gh-issue-1.cpp
* \author Joey Dumont <[email protected]>
* \since 2014-07-11
* \brief Tests Issue 1 of WignerSymbols (see valandil.github.com/wignerSymbols.
* \copyright LGPL
* This file tests the bug described in Issue #1 of WignerSymbols (GitHub).
* It seems that for some values of l1, the algorithm returns a set
* of zeros even though it should return proper coefficients.
* Permuting the l's lead to the proper output.
*/
#include <wignerSymbols.h>
int main (int argc, char* argv[])
{
double test11 = WignerSymbols::wigner3j(325, 999, 1221, 280, 899, -1179);
double test12 = WignerSymbols::wigner3j(999, 1221, 325, 899, -1179, 280);
double test13 = WignerSymbols::wigner3j(1221, 325, 999, -1179, 280, 899);
double test14 = WignerSymbols::wigner3j_f(325, 999, 1221, 280, 899, -1179);
double test15 = WignerSymbols::wigner3j_f(999, 1221, 325, 899, -1179, 280);
double test16 = WignerSymbols::wigner3j_f(1221, 325, 999, -1179, 280, 899);
double test21 = WignerSymbols::wigner3j(693, 896, 1371, 513, -838, 325);
double test22 = WignerSymbols::wigner3j(896, 1371, 693, -838, 325, 513);
double test23 = WignerSymbols::wigner3j(1371, 693, 896, 325, 513, -838);
double test24 = WignerSymbols::wigner3j_f(693, 896, 1371, 513, -838, 325);
double test25 = WignerSymbols::wigner3j_f(896, 1371, 693, -838, 325, 513);
double test26 = WignerSymbols::wigner3j_f(1371, 693, 896, 325, 513, -838);
double test31 = WignerSymbols::wigner3j(772, 874, 1231, 442, 756, -1198);
double test32 = WignerSymbols::wigner3j(874, 1231, 772, 756, -1198, 442);
double test33 = WignerSymbols::wigner3j(1231, 772, 874, -1198, 442, 756);
double test34 = WignerSymbols::wigner3j_f(772, 874, 1231, 442, 756, -1198);
double test35 = WignerSymbols::wigner3j_f(874, 1231, 772, 756, -1198, 442);
double test36 = WignerSymbols::wigner3j_f(1231, 772, 874, -1198, 442, 756);
std::cout << test11 << std::endl;
std::cout << WignerSymbols::wigner3j(693, 896, 1371, 513, -838, 325) << std::endl;
std::cout << WignerSymbols::wigner3j(772, 874, 1231, 442, 756, -1198) << std::endl;
std::cout << WignerSymbols::wigner3j(842, 996, 1714, 376, 979, -1355) << std::endl;
std::cout << WignerSymbols::wigner3j(101, 987, 1084, -80, -935, 1015) << std::endl;
std::cout << WignerSymbols::wigner3j( 51, 1003, 978, -32, 993, -961) << std::endl;
std::cout << WignerSymbols::wigner3j(217, 1008, 1107, -76, -987, 1063)<< std::endl;
std::cout << WignerSymbols::wigner3j(408, 894, 954, -114, -840, 954) << std::endl;
std::cout << WignerSymbols::wigner3j(808, 980, 1734, 341, 954, -1295) << std::endl;
std::cout << WignerSymbols::wigner3j(827, 1020, 1570, 590, 902, -1492)<< std::endl;
std::cout << "C++ impl.: " << test11 << ", " << test12 << ", " << test13 << std::endl;
std::cout << "FOR impl.: " << test14 << ", " << test15 << ", " << test16 << std::endl;
std::cout << "C++ impl.: " << test21 << ", " << test22 << ", " << test23 << std::endl;
std::cout << "FOR impl.: " << test24 << ", " << test25 << ", " << test26 << std::endl;
std::cout << "C++ impl.: " << test31 << ", " << test32 << ", " << test33 << std::endl;
std::cout << "FOR impl.: " << test34 << ", " << test35 << ", " << test36 << std::endl;
std::vector<double> test4 = WignerSymbols::wigner3j(992, 1243, 196, -901, 705);
double test41 = WignerSymbols::wigner3j(529, 992, 1243, 196, -901, 705);
double test42 = WignerSymbols::wigner3j_f(529, 992, 1243, 196, -901, 705);
std::cout << "C++ impl.: " << test41 << std::endl;
std::cout << "FOR impl.: " << test42 << std::endl;
std::cout << std::endl;
std::vector<double> test5 = WignerSymbols::wigner3j(856, 1200, 464, -828, 364);
double test51 = WignerSymbols::wigner3j(751, 856, 1200, 464, -828, 364);
double test52 = WignerSymbols::wigner3j_f(751, 856, 1200, 464, -828, 364);
/*
[ 841 379 1011 -631 313 318] -2.44096504011e-41 -0.0
[ 570 1007 1392 327 -933 606] -1.74376347733e-98 0.0
[ 970 727 1202 533 -663 130] -6.93009562166e-12 0.0
[ 905 919 1670 869 -594 -275] -3.48516309858e-195 -0.0
[ 895 574 1392 793 -365 -428] -1.41868655509e-146 -0.0
*/
return 0;
}