subroutine reac_diff !INCLUDE 'param_com.com' !INCLUDE 'param_out.com' !INCLUDE 'param_lg_out.com' !INCLUDE 'param_agr.com' !INCLUDE 'timers.com' INTEGER, PARAMETER :: nNonZeros=9 INTEGER, PARAMETER :: nRows=5,nCols=5 DOUBLE PRECISION, DIMENSION(nNonZeros) :: rValues DOUBLE PRECISION, DIMENSION(nRows) :: rRhsValues DOUBLE PRECISION, DIMENSION(nRows) :: rSolValues INTEGER, DIMENSION(nRows+1) :: rowIndex INTEGER, DIMENSION(nNonZeros) :: columns INTEGER :: error INTEGER*8 :: fhandle rValues=(/9.,1.5,6.,0.75,3.,0.5,-12.,0.625,16./) rRhsValues=(/1.,2.,3.,4.,5./) rowIndex=(/1,6,7,8,9,10/) columns=(/1,2,3,4,5,2,3,4,5/) error=DSS_CREATE(fhandle,MKL_DSS_DEFAULTS) error=DSS_DEFINE_STRUCTURE( fhandle, MKL_DSS_SYMMETRIC, & rowIndex, nRows, nCols, columns, nNonZeros ) error=DSS_REORDER(fhandle, MKL_DSS_DEFAULTS,perm) error=DSS_FACTOR_REAL(fhandle,MKL_DSS_INDEFINITE,rValues) error=DSS_SOLVE_REAL(fhandle, MKL_DSS_DEFAULTS, & rRhsValues, 1, rSolValues) WRITE(*,*) rSolValues end