make x arrays assumed-shape w/ runtime size check

Author: nshaffer

src/stdlib_experimental_quadrature.fypp

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

src/stdlib_experimental_quadrature_simps.fypp

@@ -90,7 +90,7 @@ contains
 
     pure recursive module function simps_x_${k1}$(y, x, even) result(integral)
         ${t1}$, dimension(:), intent(in) :: y
-        ${t1}$, dimension(size(y)), intent(in) :: x
+        ${t1}$, dimension(:), intent(in) :: x
         integer, intent(in), optional :: even
         ${t1}$ :: integral
 
@@ -101,6 +101,7 @@ contains
         ${t1}$ :: a, b, c, d
 
         n = size(y)
+        if (size(x) /= n) error stop "simps: Arguments `x` and `y` must be the same size."
 
         select case (n)
         case (0:1)

src/stdlib_experimental_quadrature_trapz.fypp

@@ -32,13 +32,14 @@ contains
 
     pure module function trapz_x_${KIND}$(y, x) result(integral)
         real(${KIND}$), dimension(:), intent(in) :: y
-        real(${KIND}$), dimension(size(y)), intent(in) :: x
+        real(${KIND}$), dimension(:), intent(in) :: x
         real(${KIND}$) :: integral
 
         integer :: i
         integer :: n
 
         n = size(y)
+        if (size(x) /= n) error stop "trapz: Arguments `x` and `y` must be the same size."
 
         select case (n)
         case (0:1)