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
e4f123f9
Commit
e4f123f9
authored
Mar 25, 2015
by
Jean-Eric Campagne
Browse files
(JEC) 25/3/15 protection si NTL pas installee
parent
d7db54b0
Changes
1
Hide whitespace changes
Inline
Side-by-side
Showing
1 changed file
with
10 additions
and
6 deletions
+10
-6
tstlag.cc
tstlag.cc
+10
-6
No files found.
tstlag.cc
View file @
e4f123f9
...
...
@@ -18,7 +18,6 @@
using
namespace
NTL
;
//#include "laguerreBuilder.h"
#include "lagSphericTransform.h"
int
main
(
int
narg
,
char
*
arg
[])
{
...
...
@@ -71,6 +70,7 @@ int main(int narg, char *arg[]) {
#if TESTLEVEL == 1
#ifdef SO_NTLXDBLE
{
//Test Laguerre quadrature nodes & weights Xdouble
Timer
tm
(
"Test 1 (xdble)"
);
...
...
@@ -122,7 +122,7 @@ int main(int narg, char *arg[]) {
for
(
int
k
=
2
;
k
<=
alpha
;
k
++
)
facta
*=
(
xdouble
)(
k
);
cout
<<
"Xdble Sum weights = "
<<
sumW
<<
" ->"
<<
sumW
-
facta
<<
endl
;
}
#endif
{
//idem en r_16
Timer
tm
(
"Test 1 (r_16)"
);
...
...
@@ -182,6 +182,7 @@ int main(int narg, char *arg[]) {
#if TESTLEVEL == 2
//Single (Inverse) Laguerre Transform
#ifdef SO_NTLXDBLE
{
//xdouble case
...
...
@@ -227,6 +228,7 @@ int main(int narg, char *arg[]) {
cout
<<
"Xdble Err. Max. "
<<
err_abs
<<
" ["
<<
imax
<<
"], Err. Rel. "
<<
err_rel
<<
endl
;
cout
<<
" >>>>>>>>>>>>>>>>>>>>> <<<<<<<<<<<<<<<<<<<<<<<<<<"
<<
endl
;
}
//end xdouble
#endif
{
//r_16 case
...
...
@@ -278,7 +280,7 @@ int main(int narg, char *arg[]) {
#if TESTLEVEL == 3
//Multiple (Inverse) Laguerre Transform
#ifdef SO_NTLXDBLE
{
//xdouble
Timer
tm
(
"Test 3 (xdble)"
);
...
...
@@ -324,7 +326,7 @@ int main(int narg, char *arg[]) {
cout
<<
"Xdble Err. Max. "
<<
err_abs
<<
" ["
<<
imax
<<
"], Err. Rel. "
<<
err_rel
<<
endl
;
cout
<<
" >>>>>>>>>>>>>>>>>>>>> <<<<<<<<<<<<<<<<<<<<<<<<<<"
<<
endl
;
}
#endif
{
//r_16
...
...
@@ -379,6 +381,9 @@ int main(int narg, char *arg[]) {
#if TESTLEVEL == 4
// Spherical Harmonics & Laguerre Transform; input/output as internal computations same type 'T'
#ifdef SO_NTLXDBLE
{
// T = xdouble
...
...
@@ -462,6 +467,7 @@ int main(int narg, char *arg[]) {
cout
<<
" >>>>>>>>>>>>>>>>>>>>> <<<<<<<<<<<<<<<<<<<<<<<<<<"
<<
endl
;
}
#endif
{
// T = r_16
ThSDR48RandGen
rg
;
...
...
@@ -545,8 +551,6 @@ int main(int narg, char *arg[]) {
cout
<<
" >>>>>>>>>>>>>>>>>>>>> <<<<<<<<<<<<<<<<<<<<<<<<<<"
<<
endl
;
}
#endif
...
...
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