56 #ifndef AMESOS2_PARDISOMKL_DECL_HPP
57 #define AMESOS2_PARDISOMKL_DECL_HPP
61 #include <Teuchos_StandardParameterEntryValidators.hpp>
64 #include "Amesos2_SolverCore.hpp"
65 #include "Amesos2_PardisoMKL_FunctionMap.hpp"
81 template <
class Matrix,
97 typedef typename super_type::scalar_type scalar_type;
98 typedef typename super_type::local_ordinal_type local_ordinal_type;
99 typedef typename super_type::global_ordinal_type global_ordinal_type;
100 typedef typename super_type::global_size_type global_size_type;
104 typedef typename type_map::type solver_scalar_type;
105 typedef typename type_map::magnitude_type solver_magnitude_type;
109 typedef typename
TypeMap<Amesos2::PardisoMKL,local_ordinal_type>
::type int_t;
118 typedef Kokkos::DefaultHostExecutionSpace HostExecSpaceType;
119 typedef Kokkos::View<int_t*, HostExecSpaceType> host_size_type_array;
120 typedef Kokkos::View<int_t*, HostExecSpaceType> host_ordinal_type_array;
121 typedef Kokkos::View<solver_scalar_type*, HostExecSpaceType> host_value_type_array;
133 Teuchos::RCP<Vector> X,
134 Teuchos::RCP<const Vector> B);
182 int
solve_impl(const Teuchos::Ptr<MultiVecAdapter<Vector> > X,
183 const Teuchos::Ptr<const MultiVecAdapter<Vector> > B) const;
284 mutable Teuchos::Array<solver_scalar_type>
xvals_;
286 mutable Teuchos::Array<solver_scalar_type>
bvals_;
307 static const int_t maxfct_;
308 static const int_t mnum_;
311 static const bool complex_
312 = Meta::or_<std::is_same_v<solver_scalar_type, PMKL::_MKL_Complex8>,
313 std::is_same_v<solver_scalar_type, PMKL::_DOUBLE_COMPLEX_t>>::value;
323 #ifdef HAVE_TEUCHOS_COMPLEX
324 typedef Meta::make_list6<float,
327 std::complex<double>,
329 PMKL::_DOUBLE_COMPLEX_t> supported_scalars;
331 typedef Meta::make_list2<float,
332 double> supported_scalars;
338 #endif // AMESOS2_PARDISOMKL_DECL_HPP
Amesos2::SolverCore: A templated interface for interaction with third-party direct sparse solvers...
Definition: Amesos2_SolverCore_decl.hpp:105
Map types to solver-specific data-types and enums.
Definition: Amesos2_TypeMap.hpp:82
int_t n_
Number of equations in the sparse linear system.
Definition: Amesos2_PardisoMKL_decl.hpp:293
int_t nrhs_
number of righthand-side vectors
Definition: Amesos2_PardisoMKL_decl.hpp:297
void set_pardiso_mkl_matrix_type(int_t mtype=0)
Definition: Amesos2_PardisoMKL_def.hpp:581
static const int_t msglvl_
The messaging level. Set to 1 if you wish for Pardiso MKL to print statistical info.
Definition: Amesos2_PardisoMKL_decl.hpp:304
Teuchos::Array< solver_scalar_type > bvals_
Persisting, contiguous, 1D store for B.
Definition: Amesos2_PardisoMKL_decl.hpp:286
void * pt_[64]
PardisoMKL internal data address pointer.
Definition: Amesos2_PardisoMKL_decl.hpp:289
Provides traits about solvers.
Definition: Amesos2_SolverTraits.hpp:71
Teuchos::Array< int_t > perm_
Permutation vector.
Definition: Amesos2_PardisoMKL_decl.hpp:295
Amesos2 interface to the PardisoMKL package.
Definition: Amesos2_PardisoMKL_decl.hpp:83
int numericFactorization_impl()
PardisoMKL specific numeric factorization.
Definition: Amesos2_PardisoMKL_def.hpp:176
int symbolicFactorization_impl()
Perform symbolic factorization of the matrix using PardisoMKL.
Definition: Amesos2_PardisoMKL_def.hpp:144
host_value_type_array nzvals_view_
Stores the values of the nonzero entries for PardisoMKL.
Definition: Amesos2_PardisoMKL_decl.hpp:278
bool loadA_impl(EPhase current_phase)
Reads matrix data into internal structures.
Definition: Amesos2_PardisoMKL_def.hpp:491
std::string name() const override
Return the name of this solver.
Definition: Amesos2_SolverCore_def.hpp:759
Teuchos::Array< solver_scalar_type > xvals_
Persisting, contiguous, 1D store for X.
Definition: Amesos2_PardisoMKL_decl.hpp:284
host_size_type_array rowptr_view_
Stores the row indices of the nonzero entries.
Definition: Amesos2_PardisoMKL_decl.hpp:282
void setParameters_impl(const Teuchos::RCP< Teuchos::ParameterList > ¶meterList)
Definition: Amesos2_PardisoMKL_def.hpp:286
int solve_impl(const Teuchos::Ptr< MultiVecAdapter< Vector > > X, const Teuchos::Ptr< const MultiVecAdapter< Vector > > B) const
PardisoMKL specific solve.
Definition: Amesos2_PardisoMKL_def.hpp:203
Passes functions to TPL functions based on type.
Definition: Amesos2_FunctionMap.hpp:76
void check_pardiso_mkl_error(EPhase phase, int_t error) const
Throws an appropriate runtime error in the event that error < 0 .
Definition: Amesos2_PardisoMKL_def.hpp:528
Provides access to interesting solver traits.
host_ordinal_type_array colind_view_
Stores the location in Ai_ and Aval_ that starts row j.
Definition: Amesos2_PardisoMKL_decl.hpp:280
Teuchos::RCP< const Teuchos::ParameterList > getValidParameters_impl() const
Definition: Amesos2_PardisoMKL_def.hpp:389
bool matrixShapeOK_impl() const
Determines whether the shape of the matrix is OK for this solver.
Definition: Amesos2_PardisoMKL_def.hpp:277
int_t mtype_
The matrix type. We deal only with unsymmetrix matrices.
Definition: Amesos2_PardisoMKL_decl.hpp:291
int_t iparm_[64]
Definition: Amesos2_PardisoMKL_decl.hpp:301
int preOrdering_impl()
Performs pre-ordering on the matrix to increase efficiency.
Definition: Amesos2_PardisoMKL_def.hpp:133