Commit d14ce744 by Francois Gygi

test only if deltar positive

git-svn-id: http://qboxcode.org/svn/qb/trunk@1706 cba15fb0-1239-40c8-b417-11db7ca47a34
parent 793596f0
......@@ -112,8 +112,7 @@ bool Species::initialize_ncpp()
// ndft_ is a power of 2 larger than ( rdftmin / deltar_ )
// minimum outer bound in (a.u.)
const double rdftmin = 40.0;
// next line: limit small mesh sizes
assert(deltar_ > 0.0001);
assert(deltar_ > 0.0);
ndft_ = 1;
while ( ndft_ * deltar_ < rdftmin )
ndft_ *= 2;
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment