fix wrong pardiso matrix type in complex case

Author: j-fu

Project.toml

@@ -1,7 +1,7 @@
 name = "ExtendableSparse"
 uuid = "95c220a8-a1cf-11e9-0c77-dbfce5f500b3"
 authors = ["Juergen Fuhrmann <juergen.fuhrmann@wias-berlin.de>"]
-version = "0.6.5"
+version = "0.6.6"
 
 [deps]
 DocStringExtensions = "ffbed154-4ef7-542d-bbb7-c09d3a79fcae"

src/pardiso_lu.jl

@@ -88,7 +88,7 @@ function default_initialize!(fact::AbstractPardisoLU{Tv,Ti}, iparm,dparm,mtype)
         my_mtype=Pardiso.REAL_NONSYM
     end
 
-    Pardiso.set_matrixtype!(fact.ps,Pardiso.REAL_NONSYM)
+    Pardiso.set_matrixtype!(fact.ps,my_mtype)
 
     if !isnothing(iparm)
         for i=1:min(length(iparm),length(fact.ps.iparm))