promote use of assumed shape

Author: jalvesz

src/stdlib_sparse_spmv.fypp

@@ -513,23 +513,23 @@ contains
         #:for chunk in CHUNKS
         pure subroutine chunk_kernel_${chunk}$(n,a,col,x,y)
             integer, value      :: n
-            ${t1}$, intent(in)  :: a(${chunk}$,n), x(*)
+            ${t1}$, intent(in)  :: a(${chunk}$,n), x(:)
             integer(ilp), intent(in) :: col(${chunk}$,n)
             ${t1}$, intent(inout) :: y(${chunk}$)
             integer :: j
             do j = 1, n
-                y = y + alpha_ * a(:,j) * x(col(:,j))
+                y(:) = y(:) + alpha_ * a(:,j) * x(col(:,j))
             end do
         end subroutine
         pure subroutine chunk_kernel_trans_${chunk}$(n,a,col,x,y)
             integer, value      :: n
             ${t1}$, intent(in)  :: a(${chunk}$,n), x(${chunk}$)
             integer(ilp), intent(in) :: col(${chunk}$,n)
-            ${t1}$, intent(inout) :: y(*)
+            ${t1}$, intent(inout) :: y(:)
             integer :: j, k
             do j = 1, n
                 do k = 1, ${chunk}$
-                    y(col(k,j)) = y(col(k,j)) + alpha_ * a(k,j) * x(k)
+                    y(col(:,j)) = y(col(:,j)) + alpha_ * a(:,j) * x(:)
                 end do
             end do
         end subroutine
@@ -538,7 +538,7 @@ contains
             integer, value      :: n
             ${t1}$, intent(in)  :: a(${chunk}$,n), x(${chunk}$)
             integer(ilp), intent(in) :: col(${chunk}$,n)
-            ${t1}$, intent(inout) :: y(*)
+            ${t1}$, intent(inout) :: y(:)
             integer :: j, k
             do j = 1, n
                 do k = 1, ${chunk}$
@@ -551,19 +551,19 @@ contains
 
         pure subroutine chunk_kernel_rm(n,cs,r,a,col,x,y)
             integer, value      :: n, cs, r
-            ${t1}$, intent(in)  :: a(cs,n), x(*)
+            ${t1}$, intent(in)  :: a(cs,n), x(:)
             integer(ilp), intent(in) :: col(cs,n)
             ${t1}$, intent(inout) :: y(r)
             integer :: j
             do j = 1, n
-                y = y + alpha_ * a(1:r,j) * x(col(1:r,j))
+                y(1:r) = y(1:r) + alpha_ * a(1:r,j) * x(col(1:r,j))
             end do
         end subroutine
         pure subroutine chunk_kernel_rm_trans(n,cs,r,a,col,x,y)
             integer, value      :: n, cs, r
             ${t1}$, intent(in)  :: a(cs,n), x(r)
             integer(ilp), intent(in) :: col(cs,n)
-            ${t1}$, intent(inout) :: y(*)
+            ${t1}$, intent(inout) :: y(:)
             integer :: j, k
             do j = 1, n
                 do k = 1, r
@@ -576,7 +576,7 @@ contains
             integer, value      :: n, cs, r
             ${t1}$, intent(in)  :: a(cs,n), x(r)
             integer(ilp), intent(in) :: col(cs,n)
-            ${t1}$, intent(inout) :: y(*)
+            ${t1}$, intent(inout) :: y(:)
             integer :: j, k
             do j = 1, n
                 do k = 1, r