I want to use the direct sparse solver with the intel fortran interface from the mkl library. The matrix I want to invert is symmetric complex non-hermitian. Can someone give me a hint what the correct option in dss_define_structure, dss_factor_complex would be. I use a 64bit linux cluster version of mkl10.
I used the dss_sym_f90 just as a template of how to use the solver. My problem is different in some characteristics. The matrix I want to solve is 210000 by 210000 with ca 23 * 10**6 non zero entries, symmetric with complex values on the diagonal(therefore non-hermitian). In the array "values" I did put all non zero values(not just the upper triangle). When I put all non zero elements in the "values" arrays do I have to choose MKL_DSS_NON_SYMMETRIC in the line with:
the solution if you have a symmetric non hermitian matrix is:
1) put all non zero elements in the input array (not only the upper triangle) 2) declare the structure as non symmetric (although it's symmetric) -> DSS_DEFINE_STRUCTURE( handle, MKL_DSS_NON_SYMMETRIC, rowindex, nrows, nrows, columns, totalnonzero) 3) choose the correct factorzation - in this case LDL^T -> DSS_FACTOR_COMPLEX( handle, MKL_DSS_POSITIVE_DEFINITE, values)