aboutsummaryrefslogtreecommitdiff
path: root/externals/grill/vasp/source/ops_dft.cpp
diff options
context:
space:
mode:
authorThomas Grill <xovo@users.sourceforge.net>2003-06-13 02:45:16 +0000
committerThomas Grill <xovo@users.sourceforge.net>2003-06-13 02:45:16 +0000
commitc0c7359304c7bb54d794c74b260d9b0bd18f02ae (patch)
treec16ff1091eb663e4483e6fffe8f2db1abc152d28 /externals/grill/vasp/source/ops_dft.cpp
parentd1f912c93aef58acb6b69b16576792bbcd608399 (diff)
""
svn path=/trunk/; revision=701
Diffstat (limited to 'externals/grill/vasp/source/ops_dft.cpp')
-rw-r--r--externals/grill/vasp/source/ops_dft.cpp12
1 files changed, 6 insertions, 6 deletions
diff --git a/externals/grill/vasp/source/ops_dft.cpp b/externals/grill/vasp/source/ops_dft.cpp
index 354f126b..fd00d921 100644
--- a/externals/grill/vasp/source/ops_dft.cpp
+++ b/externals/grill/vasp/source/ops_dft.cpp
@@ -79,7 +79,7 @@ static BL fft_fwd_real_any(I cnt,F *rsdt,I _rss,F *rddt,I _rds)
#endif
}
if(ret) {
- const F nrm = 1./sqrt(cnt);
+ const F nrm = 1./sqrt((F)cnt);
const I n2 = cnt/2;
#ifndef VASP_COMPACT
@@ -158,7 +158,7 @@ static BL fft_inv_real_any(I cnt,F *rsdt,I _rss,F *rddt,I _rds)
#endif
}
if(ret) {
- const F nrm = 1./sqrt(cnt);
+ const F nrm = 1./sqrt((F)cnt);
#ifndef VASP_COMPACT
if(rds == 1)
for(i = 0; i < cnt; ++i)
@@ -230,7 +230,7 @@ static BL fft_fwd_complex_any(I cnt,F *rsdt,I _rss,F *isdt,I _iss,F *rddt,I _rds
#endif
}
if(ret) {
- const F nrm = 1./sqrt(cnt);
+ const F nrm = 1./sqrt((F)cnt);
#ifdef VASP_COMPACT
for(i = 0; i < cnt; ++i) {
@@ -359,7 +359,7 @@ static BL fft_complex_radix2(I cnt,F *rsdt,I _rss,F *isdt,I _iss,F *rddt,I _rds,
BL ret = fft_bidir_complex_radix2(cnt,rtmp,itmp,dir);
if(ret) {
- const F nrm = 1./sqrt(cnt);
+ const F nrm = 1./sqrt((F)cnt);
#ifndef VASP_COMPACT
if(rtmp == rddt)
@@ -426,7 +426,7 @@ BL fft_fwd_real_radix2(I cnt,F *src,I _sstr,F *dst,I _dstr)
register I i;
const I n2 = cnt/2;
- const F fn = (F)(1./sqrt(cnt));
+ const F fn = (F)(1./sqrt((F)cnt));
F *stmp;
if(!dst || src == dst) {
// in-place
@@ -502,7 +502,7 @@ BL fft_inv_real_radix2(I cnt,F *src,I _sstr,F *dst,I _dstr)
register I i;
const I n2 = cnt/2;
- const F fn = (F)(1./sqrt(cnt));
+ const F fn = (F)(1./sqrt((F)cnt));
F *stmp;
if(!dst || src == dst) {
// in-place