- #1

rabbed

- 243

- 3

Page 51 in

<Moderator's note: link to copyrighted material removed, see instead http://dx.doi.org/10.1119/1.2967703>

The image is showing functions:

y1=e^(ln(n0)-x)

y2=invDigamma(digamma(n0+1)-x)-1

where y2 is the "exact" version. How exactly is it exact?

For this example (with 7 particles):

https://bouman.chem.georgetown.edu/S02/lect21/lect21.htm

if I use

y2 = Round( invDigamma(digamma(n0+1)-0.4*x)-1 )

for x=0, 1, 2, 3

i get (n0, n1, n2, n3) = (3, 2, 1, 1) which agrees with the example.

But for 8 particles I think the correct answer should be (3, 3, 1, 1) which I cannot get with the "exact" function.

Am I thinking about this incorrectly?

<Moderator's note: link to copyrighted material removed, see instead http://dx.doi.org/10.1119/1.2967703>

The image is showing functions:

y1=e^(ln(n0)-x)

y2=invDigamma(digamma(n0+1)-x)-1

where y2 is the "exact" version. How exactly is it exact?

For this example (with 7 particles):

https://bouman.chem.georgetown.edu/S02/lect21/lect21.htm

if I use

y2 = Round( invDigamma(digamma(n0+1)-0.4*x)-1 )

for x=0, 1, 2, 3

i get (n0, n1, n2, n3) = (3, 2, 1, 1) which agrees with the example.

But for 8 particles I think the correct answer should be (3, 3, 1, 1) which I cannot get with the "exact" function.

Am I thinking about this incorrectly?

Last edited by a moderator: