allow only two matrix types in Pardiso: real nonsym and complex nonsym until

symmetric case is fixed

Author: j-fu

src/pardiso_lu.jl

@@ -80,9 +80,11 @@ MKLPardisoLU(;valuetype::Type=Float64, indextype::Type=Int64,kwargs...)=MKLPardi
 
 ##########################################################################################
 function default_initialize!(fact::AbstractPardisoLU{Tv,Ti}, iparm,dparm,mtype) where {Tv, Ti}
-    if !isnothing(mtype)
-        my_mtype=mtype
-    elseif Tv<:Complex
+    # if !isnothing(mtype)
+    #     my_mtype=mtype fix this!
+    # else
+
+    if Tv<:Complex
         my_mtype=Pardiso.COMPLEX_NONSYM
     else
         my_mtype=Pardiso.REAL_NONSYM

test/runtests.jl

@@ -372,9 +372,9 @@ end
     @test test_lu2(25,40,1,lufac=MKLPardisoLU())
     @test test_lu2(1000,1,1,lufac=MKLPardisoLU())
 
-    @test test_lu2(10,10,10,lufac=MKLPardisoLU(mtype=2))
-    @test test_lu2(25,40,1,lufac=MKLPardisoLU(mtype=2))
-    @test test_lu2(1000,1,1,lufac=MKLPardisoLU(mtype=2))
+    # @test test_lu2(10,10,10,lufac=MKLPardisoLU(mtype=2))
+    # @test test_lu2(25,40,1,lufac=MKLPardisoLU(mtype=2))
+    # @test test_lu2(1000,1,1,lufac=MKLPardisoLU(mtype=2))
 end
 
 
@@ -388,9 +388,9 @@ if Pardiso.PARDISO_LOADED[]
         @test test_lu2(25,40,1,lufac=PardisoLU())
         @test test_lu2(1000,1,1,lufac=PardisoLU())
 
-        @test test_lu2(10,10,10,lufac=PardisoLU(mtype=2))
-        @test test_lu2(25,40,1,lufac=PardisoLU(mtype=2))
-        @test test_lu2(1000,1,1,lufac=PardisoLU(mtype=2))
+        # @test test_lu2(10,10,10,lufac=PardisoLU(mtype=2))
+        # @test test_lu2(25,40,1,lufac=PardisoLU(mtype=2))
+        # @test test_lu2(1000,1,1,lufac=PardisoLU(mtype=2))
 
     end
 end