Hi all,
I'm not sure, but from my point of view the following abscissae values for n=6 in the quadrature.m file should be
r1pt(5) = 0.238619186083152; r1pt(6) =-0.238619186083152;
instead of
r1pt(5) = 0.238619186003152; r1pt(6) =-0.238619186003152;
?
Best regards! Frederik
I think you are correct. We will check the integration subroutines and fix this bug in the next release. Thanks for pointing it out!
Log in to post a comment.
Hi all,
I'm not sure, but from my point of view the following abscissae values for n=6 in the quadrature.m file should be
r1pt(5) = 0.238619186083152;
r1pt(6) =-0.238619186083152;
instead of
r1pt(5) = 0.238619186003152;
r1pt(6) =-0.238619186003152;
?
Best regards!
Frederik
I think you are correct. We will check the integration subroutines and fix this bug in the next release. Thanks for pointing it out!