{simps,trapz}_x_*: error stop -> check; drop pure

Author: nshaffer

src/stdlib_experimental_quadrature.fypp

@@ -25,7 +25,7 @@ module stdlib_experimental_quadrature
         end function trapz_dx_${k1}$      
         #:endfor
         #:for k1, t1 in REAL_KINDS_TYPES
-        pure module function trapz_x_${k1}$(y, x) result(integral)
+        module function trapz_x_${k1}$(y, x) result(integral)
             ${t1}$, dimension(:), intent(in) :: y
             ${t1}$, dimension(:), intent(in) :: x
             ${t1}$ :: integral
@@ -59,7 +59,7 @@ module stdlib_experimental_quadrature
         end function simps_dx_${k1}$
         #:endfor
         #:for k1, t1 in REAL_KINDS_TYPES
-        pure recursive module function simps_x_${k1}$(y, x, even) result(integral)
+        recursive module function simps_x_${k1}$(y, x, even) result(integral)
             ${t1}$, dimension(:), intent(in) :: y
             ${t1}$, dimension(:), intent(in) :: x
             integer, intent(in), optional :: even

src/stdlib_experimental_quadrature_simps.fypp

@@ -1,7 +1,7 @@
 #:include "common.fypp"
 
 submodule (stdlib_experimental_quadrature) stdlib_experimental_quadrature_simps
-
+    use stdlib_experimental_error, only: check
     implicit none
 
     ! internal use only
@@ -88,7 +88,7 @@ contains
 #:endfor
 #:for k1, t1 in REAL_KINDS_TYPES
 
-    pure recursive module function simps_x_${k1}$(y, x, even) result(integral)
+    recursive module function simps_x_${k1}$(y, x, even) result(integral)
         ${t1}$, dimension(:), intent(in) :: y
         ${t1}$, dimension(:), intent(in) :: x
         integer, intent(in), optional :: even
@@ -101,7 +101,7 @@ contains
         ${t1}$ :: a, b, c
 
         n = size(y)
-        if (size(x) /= n) error stop "simps: Arguments `x` and `y` must be the same size."
+        call check(size(x) == n, "simps: Arguments `x` and `y` must be the same size.")
 
         select case (n)
         case (0:1)

src/stdlib_experimental_quadrature_trapz.fypp

@@ -1,7 +1,7 @@
 #:include "common.fypp"
 
 submodule (stdlib_experimental_quadrature) stdlib_experimental_quadrature_trapz
-
+    use stdlib_experimental_error, only: check
     implicit none
 
 contains
@@ -29,8 +29,8 @@ contains
 
     #:endfor
     #:for KIND in REAL_KINDS
-    
-    pure module function trapz_x_${KIND}$(y, x) result(integral)
+
+    module function trapz_x_${KIND}$(y, x) result(integral)
         real(${KIND}$), dimension(:), intent(in) :: y
         real(${KIND}$), dimension(:), intent(in) :: x
         real(${KIND}$) :: integral
@@ -39,7 +39,7 @@ contains
         integer :: n
 
         n = size(y)
-        if (size(x) /= n) error stop "trapz: Arguments `x` and `y` must be the same size."
+        call check(size(x) == n, "trapz: Arguments `x` and `y` must be the same size.")
 
         select case (n)
         case (0:1)