correct error in legendre polinomial evaluation

nomad-coe · Nov 6, 2024 · b99e7b5 · b99e7b5
1 parent 87b81cd
commit b99e7b5
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/GX-NAOIntegrals/src/legendre_polynomial.f90 b/GX-NAOIntegrals/src/legendre_polynomial.f90
@@ -35,7 +35,7 @@ subroutine evaluate_legendre_polinomial(n, x, p)
         p(1) = 1.0_dp 
         p(2) = x
         do i = 2, n 
-            p(i+1) = (2.0_dp*n-1.0_dp)/dble(n) * x * p(i) - (n-1.0_dp)/dble(n) * p(i-1)
+            p(i+1) = (2.0_dp*i-1.0_dp)/dble(i) * x * p(i) - (i-1.0_dp)/dble(i) * p(i-1)
         end do 
 
     end subroutine evaluate_legendre_polinomial
@@ -67,7 +67,7 @@ subroutine evaluate_legendre_polinomial_batch(n, n_points, x, p)
         p(:, 1) = 1.0_dp 
         p(:, 2) = x(:)
         do i = 2, n 
-            p(:, i+1) = (2.0_dp*n-1.0_dp)/dble(n) * x(:) * p(:, i) - (n-1.0_dp)/dble(n) * p(:, i-1)
+            p(:, i+1) = (2.0_dp*i-1.0_dp)/dble(i) * x(:) * p(:, i) - (i-1.0_dp)/dble(i) * p(:, i-1)
         end do 
 
     end subroutine evaluate_legendre_polinomial_batch