Merge pull request #116 from goulart-paul/panua_license_check

add panua license check

Author: j-fu

src/Pardiso.jl

@@ -51,7 +51,9 @@ export solve, solve!
 export get_matrix
 export schur_complement, pardisogetschur
 export fix_iparm!
-export  mkl_is_available, panua_is_available
+export mkl_is_available, panua_is_available
+
+VERSION >= v"1.11.0-DEV.469" && eval(Meta.parse("public panua_is_loaded, panua_is_licensed"))
 
 struct PardisoException <: Exception
     info::String
@@ -112,9 +114,36 @@ const pardiso_chkvec = Ref{Ptr}()
 const pardiso_chkvec_z = Ref{Ptr}()
 const pardiso_get_schur_f = Ref{Ptr}()
 const PARDISO_LOADED = Ref(false)
+const PARDISO_LICENSED = Ref(false)
+
+function panua_is_licensed()
+
+    if !PARDISO_LOADED[]
+        return false
+    elseif PARDISO_LICENSED[]
+        return true
+    end
+    # Suppress unwanted output from pardisoinit, which prints license info to stdout
+    redirect_stdout(devnull) do
+        try
+            ps = PardisoSolver(;loadchecks = false)
+            pardisoinit(ps)   # errors if unlicensed
+            PARDISO_LICENSED[] = true
+            return true
+        catch e
+            if isa(e, PardisoException)
+                return false
+            else
+                rethrow(e)
+            end
+        end
+    end
+end
 
-panua_is_available() = PARDISO_LOADED[]
+panua_is_loaded() = PARDISO_LOADED[]
+panua_is_available() = panua_is_loaded() && panua_is_licensed()
 
+    
 function __init__()
     global MKL_LOAD_FAILED
     if LOCAL_MKL_FOUND

src/panua_pardiso.jl

@@ -13,9 +13,13 @@ mutable struct PardisoSolver <: AbstractPardisoSolver
     rowval::Vector{Int32}
 end
 
-function PardisoSolver()
-    if !panua_is_available()
-      error("Panua pardiso library was not loaded")
+function PardisoSolver(; loadchecks::Bool = true)
+    if loadchecks
+        if !panua_is_loaded()
+        error("Panua pardiso library was not loaded")
+        elseif !panua_is_licensed()
+        error("Panua pardiso library requires license")
+        end
     end
 
     pt = zeros(Int, 64)