Commit c443d6ff authored by ulrich_y's avatar ulrich_y

Initialised allocated Bernoulli

parent 3b568954
......@@ -71,6 +71,7 @@ CONTAINS
if(verb >= 150) print*, 'initialising Bernoulli number system for m=0,', cachecontr(1)
#endif
allocate(bernoulli(0:cachecontr(1)))
bernoulli = 0._prec
bernoulli(0) = 1.
bernoulli(1) = -0.5
endif
......@@ -83,6 +84,7 @@ CONTAINS
deallocate(bernoulli)
allocate(bernoulli(0:n + cachecontr(2)))
bernoulli = 0._prec
bernoulli(0:k) = buffer(0:k)
#ifdef DEBUG
......
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