Replaced deprecated function.

Author: albertjimenezramos

Solver/src/libs/jacobian/CSRMatrixClass.f90

@@ -69,8 +69,7 @@ subroutine mkl_dcsrgemv(transa, m, a, ia, ja, x, y)
          integer        :: ia(*), ja(*)
          real(kind=RP)  :: x(*), y(*)
       end subroutine mkl_dcsrgemv
-   end interface
-   
+   end interface   
 !
 !========
  CONTAINS
@@ -768,14 +767,17 @@ end subroutine CSR_MatMatMul
    !  -> v needs to be allocated beforehand
    !  ----------------------------------------------------
    function CSR_MatVecMul( A,u, trans) result(v)
+#ifdef HAS_MKL
+      use MKL_SPBLAS, only: SPARSE_OPERATION_NON_TRANSPOSE
+#endif
       implicit none
       !-arguments--------------------------------------------------------------------
       class(csrMat_t)  , intent(inout) :: A  !< Structure holding matrix
       real(kind=RP)    , intent(in)    :: u(A % num_of_Cols)  !< Vector to be multiplied
       logical, optional, intent(in)    :: trans   !< A matrix is transposed?
       real(kind=RP)                    :: v(A % num_of_Rows)  !> Result vector 
       !------------------------------------------------------------------------------
-      integer           :: i,j
+      integer           :: i,j,stat
       REAL(KIND=RP)     :: rsum
       character(len=1)  :: transInfo
       !------------------------------------------------------------------------------
@@ -796,6 +798,10 @@ function CSR_MatVecMul( A,u, trans) result(v)
 
 #ifdef HAS_MKL
 !      CALL mkl_dcsrgemv(transInfo, A % num_of_Rows, A % Values, A % Rows, A % Cols, u, v)
+      ! This call was changed because the function mkl_dcsrgemv was deprecated.
+      ! See: https://www.intel.com/content/www/us/en/docs/onemkl/developer-reference-fortran/2023-1/mkl-csrgemv.html
+      ! TODO: Check that the new function performs as expected.
+      stat = mkl_sparse_d_mv(SPARSE_OPERATION_NON_TRANSPOSE, 1.0_rp, A % mkl_options % csrA, A % mkl_options % descrA, u, 0.0_rp, v)
 #else
       if (transInfo == 't') error stop "CSR_MatVecMul with 't' only with MKL"
 !$omp parallel do private(j,rsum)