* <code>pdsaupd</code> work and also how to set the parameters
* appropriately please take a look into the PARPACK manual.
*
- * @author Denis Davydov, 2014.
+ * @author Denis Davydov, 2015.
*/
template <typename VECTOR>
class PArpackSolver : public Subscriptor
#include <fstream>
#include <iostream>
+// test Parpack on Step-36 with PETSc algebra
+
const unsigned int dim = 2;//run in 2d to save time
const double eps = 1e-10;
void test ()
{
+#ifdef DEAL_II_ARPACK_WITH_PARPACK
const unsigned int global_mesh_refinement_steps = 5;
const unsigned int number_of_eigenvalues = 5;
dof_handler.clear ();
+#else
+ // just output expected results:
+ deallog <<"4.93877"<<std::endl;
+ deallog <<"12.3707"<<std::endl;
+ deallog <<"12.3707"<<std::endl;
+ deallog <<"19.8027"<<std::endl;
+ deallog <<"24.8370"<<std::endl;
+#endif
dealii::deallog << "Ok"<<std::endl;
}
#include <fstream>
#include <iostream>
+// test Parpack on Step-36 with Trilinos algebra
+
const unsigned int dim = 2;//run in 2d to save time
using namespace dealii;
void test ()
{
+#ifdef DEAL_II_ARPACK_WITH_PARPACK
const unsigned int global_mesh_refinement_steps = 5;
const unsigned int number_of_eigenvalues = 5;
dof_handler.clear ();
+#else
+ // just output expected results:
+ deallog <<"4.93877"<<std::endl;
+ deallog <<"12.3707"<<std::endl;
+ deallog <<"12.3707"<<std::endl;
+ deallog <<"19.8027"<<std::endl;
+ deallog <<"24.8370"<<std::endl;
+#endif
deallog << "Ok"<<std::endl;
}