changed stdlib_functions to stdlib_specialfunctions

Author: hsnyder

src/CMakeLists.txt

@@ -41,8 +41,8 @@ set(SRC
     stdlib_kinds.f90
     stdlib_logger.f90
     stdlib_system.F90
-    stdlib_functions.f90
-    stdlib_functions_legendre.f90
+    stdlib_specialfunctions.f90
+    stdlib_specialfunctions_legendre.f90
     stdlib_quadrature_gauss.f90
     ${outFiles}
 )

src/Makefile.manual

@@ -4,8 +4,8 @@ SRC = f18estop.f90 \
       stdlib_bitsets_64.f90 \
       stdlib_bitsets_large.f90 \
       stdlib_error.f90 \
-      stdlib_functions.f90 \
-      stdlib_functions_legendre.f90 \
+      stdlib_specialfunctions.f90 \
+      stdlib_specialfunctions_legendre.f90 \
       stdlib_io.f90 \
       stdlib_kinds.f90 \
       stdlib_linalg.f90 \
@@ -53,8 +53,8 @@ stdlib_bitsets.o: stdlib_kinds.o
 stdlib_bitsets_64.o: stdlib_bitsets.o
 stdlib_bitsets_large.o: stdlib_bitsets.o
 stdlib_error.o: stdlib_optval.o
-stdlib_functions.o: stdlib_kinds.o
-stdlib_functions_legendre.o: stdlib_kinds.o
+stdlib_specialfunctions.o: stdlib_kinds.o
+stdlib_specialfunctions_legendre.o: stdlib_kinds.o
 stdlib_io.o: \
 	stdlib_error.o \
 	stdlib_optval.o \

src/stdlib_quadrature_gauss.f90

@@ -1,6 +1,6 @@
 submodule (stdlib_quadrature) stdlib_quadrature_gauss
     use stdlib_kinds, only: sp, dp, qp
-    use stdlib_functions, only: legendre, dlegendre
+    use stdlib_specialfunctions, only: legendre, dlegendre
     implicit none
 
     real(dp), parameter :: pi = acos(-1._dp)
@@ -84,6 +84,7 @@ pure module subroutine gauss_legendre_lobatto_fp64 (x, w, interval)
                 w(n+1) =  2._dp/(n*(n+1._dp))
 
                 do i = 1, (n+1)/2 - 1
+                    ! initial guess from an approximate form given by SV Parter (1999)
                     x(i+1) = -cos( (i+0.25_dp)*pi/n  - 3/(8*n*pi*(i+0.25_dp)))
                     do j = 0, newton_iters-1
                         leg  = legendre(n+1,x(i+1)) - legendre(n-1,x(i+1))

src/stdlib_functions.f90 renamed to src/stdlib_specialfunctions.f90

@@ -1,4 +1,4 @@
-module stdlib_functions
+module stdlib_specialfunctions
     use stdlib_kinds, only: sp, dp, qp
 
     implicit none
@@ -31,4 +31,4 @@ pure elemental module function dlegendre_fp64(n,x) result(dleg)
         end function
     end interface
 
-end module stdlib_functions
+end module stdlib_specialfunctions

src/stdlib_functions_legendre.f90 renamed to src/stdlib_specialfunctions_legendre.f90

@@ -1,4 +1,4 @@
-submodule (stdlib_functions) stdlib_functions_legendre
+submodule (stdlib_specialfunctions) stdlib_specialfunctions_legendre
     use stdlib_kinds, only: sp, dp, qp
     implicit none