Skip to content
GitLab
Projects
Groups
Snippets
Help
Loading...
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in / Register
Toggle navigation
Open sidebar
Jean-Eric Campagne
LagSHT
Commits
4fef23e3
Commit
4fef23e3
authored
Feb 03, 2016
by
Jean-Eric Campagne
Browse files
(JEC) 3/2/16 update testing suite
parent
7a98bebd
Changes
2
Hide whitespace changes
Inline
Side-by-side
Showing
2 changed files
with
19 additions
and
7 deletions
+19
-7
src/lagsht_testsuite.cc
src/lagsht_testsuite.cc
+16
-4
src/laguerre2bessel.cc
src/laguerre2bessel.cc
+3
-3
No files found.
src/lagsht_testsuite.cc
View file @
4fef23e3
...
...
@@ -834,7 +834,7 @@ int main(int narg, char *arg[]) {
char
*
pLAGSHTROOTDIR
;
pLAGSHTROOTDIR
=
getenv
(
"LAGSHTROOTDIR"
);
string
clenshawDir
=
string
(
pLAGSHTROOTDIR
)
+
"/data"
;
string
jlnpDir
=
string
(
pLAGSHTROOTDIR
)
+
"/data"
;
bool
recomputeJlnp
=
tru
e
;
bool
recomputeJlnp
=
fals
e
;
//JEC 12/1/16 alpha as input
...
...
@@ -973,7 +973,20 @@ int main(int narg, char *arg[]) {
}
break
;
case
5
:
case
50
:
{
cout
<<
"rebuild Jlnp:
\n
Clenshaw Path<"
<<
param
.
clenshawDir
<<
">
\n
and jlnp Path<"
<<
param
.
jlnpDir
<<
">"
<<
endl
;
if
(
Pmax
<
N
){
param
.
Pmax
=
N
;
cout
<<
"(Warning): Pmax<Nmax => modify the value to "
<<
param
.
Pmax
<<
endl
;}
tstack_push
(
"TestJlnpComputation"
);
TestJlnpComputation
();
tstack_pop
(
"TestLaguerre2Bessel"
);
tstack_report
(
"TestLaguerre2Bessel"
);
}
break
;
case
51
:
{
if
(
recomputeJlnp
){
cout
<<
"rebuild Jlnp:
\n
Clenshaw Path<"
<<
param
.
clenshawDir
<<
">
\n
and jlnp Path<"
...
...
@@ -984,8 +997,7 @@ int main(int narg, char *arg[]) {
}
if
(
Pmax
<
N
){
param
.
Pmax
=
16
*
N
;
cout
<<
"(Warning): Pmax<Nmax => modify the value to "
<<
param
.
Pmax
<<
endl
;}
tstack_push
(
"TestLaguerre2Bessel"
);
// TestLaguerre2Bessel();
TestJlnpComputation
();
TestLaguerre2Bessel
();
tstack_pop
(
"TestLaguerre2Bessel"
);
tstack_report
(
"TestLaguerre2Bessel"
);
}
...
...
src/laguerre2bessel.cc
View file @
4fef23e3
...
...
@@ -318,9 +318,9 @@ void Laguerre2Bessel::ComputeJInteg(string clenshawDir, string jlnpDir){
integ_val
=
theQuad
.
GlobalAdapStrat
(
tol
);
Jlnp_
[
ploff7
+
n
*
Jstride
]
=
integ_val
.
first
*
(
facts
[
n
]
*
sqrt
(
2.0
/
M_PI
));
cout
<<
l
<<
" "
<<
n
<<
" "
<<
p
<<
" "
<<
setprecision
(
30
)
<<
Jlnp_
[
ploff7
+
n
*
Jstride
]
<<
endl
;
//
cout <<l<<" "
//
<<n<<" "
//
<<p<< " " << setprecision(30) << Jlnp_[ploff7 + n*Jstride] << endl;
ofs
<<
l
<<
" "
<<
n
<<
" "
<<
p
<<
" "
<<
setprecision
(
30
)
<<
Jlnp_
[
ploff7
+
n
*
Jstride
]
<<
endl
;
...
...
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment