Pardiso option

CMakeLists.txt

@@ -48,11 +48,15 @@ if ( LEMMA_VTK6_SUPPORT OR LEMMA_VTK7_SUPPORT OR LEMMA_VTK8_SUPPORT AND LEMMA_MO
     	CXX_STANDARD_REQUIRED_ON
 	)
 
-        # No OpenMP
+	option ( EMSCHUR3D_ON_KINGSPEAK FALSE )
+	if ( EMSCHUR3D_ON_KINGSPEAK )
+    # No OpenMP
+		add_compile_options(-DHAVE_PARDISO)
 	#set (MKLLINK  "-L$MKLROOT/lib/intel64 -lmkl_intel_lp64 -lmkl_sequential -lmkl_core -lpthread -Wl,-rpath=$MKLROOT/lib/intel64")
-        # OpenMP 
+        # OpenMP R3
         set (MKLLINK  "-L$MKLROOT/lib/intel64 -Wl,-rpath=$MKLROOT/lib/intel64 -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_core -liomp5 -lpthread")
-	target_link_libraries(emschur3d "lemmacore" ${MKLLINK})
+		target_link_libraries(emschur3d "lemmacore" ${MKLLINK})
+	endif()
 
 	# Linking
 	target_link_libraries(emschur3d "lemmacore" "fdem1d" )

examples/EMSchur3D-vtk.cpp

@@ -114,11 +114,11 @@ int main( int argc, char** argv ) {
     // And solve
 
     // Use BiCGSTAB Diagonal preconditioner
-    //auto EM3D = EMSchur3D< Eigen::BiCGSTAB<Eigen::SparseMatrix<Complex, Eigen::ColMajor>, Eigen::IncompleteLUT<Complex> > >::NewSP();
+    auto EM3D = EMSchur3D< Eigen::BiCGSTAB<Eigen::SparseMatrix<Complex, Eigen::ColMajor>, Eigen::IncompleteLUT<Complex> > >::NewSP();
     //auto EM3D = EMSchur3D< Eigen::BiCGSTAB<Eigen::SparseMatrix<Complex, Eigen::ColMajor> > >::NewSP();
-    
+
     //auto EM3D = EMSchur3D< Eigen::PardisoLU<Eigen::SparseMatrix<Complex, Eigen::ColMajor> > >::NewSP();
-    auto EM3D = EMSchur3D< Eigen::PardisoLDLT<Eigen::SparseMatrix<Complex, Eigen::ColMajor>, Eigen::Symmetric > >::NewSP();
+    //auto EM3D = EMSchur3D< Eigen::PardisoLDLT<Eigen::SparseMatrix<Complex, Eigen::ColMajor>, Eigen::Symmetric > >::NewSP();
 
     // LS CG
     //auto EM3D = EMSchur3D< Eigen::LeastSquaresConjugateGradient<Eigen::SparseMatrix<Complex, Eigen::ColMajor> > >::NewSP();

examples/EMSchur3D.cpp

@@ -34,12 +34,12 @@ int main( int argc, char** argv ) {
 
     // SUPERLU
     //auto EM3D = EMSchur3D< Eigen::SuperLU<Eigen::SparseMatrix<Complex, Eigen::ColMajor> > >::NewSP();
-    
+
     //auto EM3D = EMSchur3D< Eigen::PardisoLU<Eigen::SparseMatrix<Complex, Eigen::ColMajor> > >::NewSP();
-    auto EM3D = EMSchur3D< Eigen::PardisoLDLT<Eigen::SparseMatrix<Complex, Eigen::ColMajor>, Eigen::Symmetric&Eigen::Lower > >::NewSP();
+    //auto EM3D = EMSchur3D< Eigen::PardisoLDLT<Eigen::SparseMatrix<Complex, Eigen::ColMajor>, Eigen::Symmetric&Eigen::Lower > >::NewSP();
 
     // Eigen built-in
-    //auto EM3D = EMSchur3D< Eigen::SparseLU<Eigen::SparseMatrix<Complex, Eigen::ColMajor> > >::NewSP();
+    auto EM3D = EMSchur3D< Eigen::SparseLU<Eigen::SparseMatrix<Complex, Eigen::ColMajor> > >::NewSP();
 
     // UmfPack
     //auto EM3D = EMSchur3D< Eigen::UmfPackLU<Eigen::SparseMatrix<Complex, Eigen::ColMajor> > >::NewSP();

include/EMSchur3D.h

@@ -36,7 +36,9 @@
 #include <Eigen/PaStiXSupport>
 #endif
 
+#ifdef HAVE_PARDISO
 #include <Eigen/PardisoSupport>
+#endif
 
 //#include "CSymSimplicialCholesky.h"
 
@@ -449,17 +451,17 @@ namespace Lemma {
             std::cout << " done in " << timer.end() / 60. << " [m]" << std::endl;
         }
     }
-    
-// Pardiso 
+
+#ifdef HAVE_PARDISO
     template<>
     void EMSchur3D< Eigen::PardisoLU<Eigen::SparseMatrix<Complex, Eigen::ColMajor> > >::BuildCDirectSolver() {
         CSolver = new Eigen::PardisoLU<Eigen::SparseMatrix<Complex, Eigen::ColMajor> > [Omegas.size()];
         for (int iw=0; iw<Omegas.size(); ++iw) {
             jsw_timer timer;
             timer.begin();
-        
+
             //CSolver[iw].pardisoParameterArray()[2] = 28; // OMP_NUM_THREADS?
-            
+
             /*  Complex system */
             std::cout << "PardisoLU pattern analyzing C_" << iw << ",";
             std::cout.flush();
@@ -481,10 +483,10 @@ namespace Lemma {
         for (int iw=0; iw<Omegas.size(); ++iw) {
             jsw_timer timer;
             timer.begin();
-        
+
             //CSolver[iw].pardisoParameterArray()[2] = 28; // OMP_NUM_THREADS?
             //Eigen::SparseMatrix<Complex> Csym = Cvec[iw].selfadjointView< Eigen::Lower >();
-            
+
             /*  Complex system */
             std::cout << "PardisoLDLT pattern analyzing C_" << iw << ",";
             std::cout.flush();
@@ -501,6 +503,7 @@ namespace Lemma {
             std::cout << " done in " << timer.end() / 60. << " [m]" << std::endl;
         }
     }
+#endif
 
 #ifdef HAVE_UMFPACK
     // Umfpack only seems to work when LOWER and UPPER are set to 1. Workarounds this have not been found.