We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 56822b5 commit 1a914ddCopy full SHA for 1a914dd
src/stdlib_quadrature_gauss.f90
@@ -26,7 +26,7 @@ pure module subroutine gauss_legendre_fp64 (x, w, interval)
26
integer :: i,j
27
real(dp) :: leg, dleg, delta
28
29
- do i = 0, int(floor((N+1)/2._dp)-1)
+ do i = 0, (N+1)/2-1
30
x(i+1) = -cos((2*i+1)/(2._dp*N+2._dp) * PI)
31
do j = 0, newton_iters-1
32
leg = legendre(N+1,x(i+1))
0 commit comments