https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/How-to-reorganize-sparse-matrix-columns-dss-reorder/m-p/819487#M4603 Hi,<DIV></DIV><DIV>I use MKL to solve FEM linear systems using the DSS library of commands. My question is related to ordering columns.</DIV><DIV></DIV><DIV>Currently, when I am currently setting up a matrix, my column numbers (due to connectivity) will not be in ascending order. Therefore I created a simple subroutine to reorganize them, however I feel that there is gotta be a better way. What happens if I send to the solver a unorganized matrix? And, is there a more efficient way to reorganize these matrices?</DIV><DIV></DIV><DIV>I know DSS_Reorder takes care of permutation issues. Does it has anything to do with my issue?</DIV><DIV></DIV><DIV>I appreciate any suggestion,</DIV><DIV></DIV><DIV>Roberto</DIV><DIV></DIV><DIV>I use Windows XP, 32 bits, Intel Fortran 11, VS2008.</DIV> Thu, 20 Jan 2011 21:55:37 GMT Roberto_Soares 2011-01-20T21:55:37Z https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/How-to-reorganize-sparse-matrix-columns-dss-reorder/m-p/819487#M4603 Hi,<DIV></DIV><DIV>I use MKL to solve FEM linear systems using the DSS library of commands. My question is related to ordering columns.</DIV><DIV></DIV><DIV>Currently, when I am currently setting up a matrix, my column numbers (due to connectivity) will not be in ascending order. Therefore I created a simple subroutine to reorganize them, however I feel that there is gotta be a better way. What happens if I send to the solver a unorganized matrix? And, is there a more efficient way to reorganize these matrices?</DIV><DIV></DIV><DIV>I know DSS_Reorder takes care of permutation issues. Does it has anything to do with my issue?</DIV><DIV></DIV><DIV>I appreciate any suggestion,</DIV><DIV></DIV><DIV>Roberto</DIV><DIV></DIV><DIV>I use Windows XP, 32 bits, Intel Fortran 11, VS2008.</DIV> Thu, 20 Jan 2011 21:55:37 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/How-to-reorganize-sparse-matrix-columns-dss-reorder/m-p/819487#M4603 Roberto_Soares 2011-01-20T21:55:37Z https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/How-to-reorganize-sparse-matrix-columns-dss-reorder/m-p/819488#M4604 <P>Hi, Roberto,</P><P>The non-zero elements of each row of the matrix <VAR>A</VAR> must be stored in increasing order (its note from MKL PARDISO manual). Reordering step permutes matrix in order to decrease the time of computations on multi-core computers. This permutation stored in array perm(). Let me note, that input matrix (arrays rowIndex() and columns() ) arent changed during all phases (reordering, factorization and solve).<BR />Sergey</P> Mon, 24 Jan 2011 05:26:15 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/How-to-reorganize-sparse-matrix-columns-dss-reorder/m-p/819488#M4604 Sergey_Solovev__Inte 2011-01-24T05:26:15Z