52#ifndef AMESOS2_MUMPS_DECL_HPP
53#define AMESOS2_MUMPS_DECL_HPP
56#include "Amesos2_SolverCore.hpp"
83template <
class Matrix,
class Vector>
94 typedef
MUMPS<Matrix,Vector> type;
98 typedef typename super_type::scalar_type scalar_type;
99 typedef typename super_type::local_ordinal_type local_ordinal_type;
100 typedef typename super_type::global_ordinal_type global_ordinal_type;
101 typedef typename super_type::global_size_type global_size_type;
103 typedef
TypeMap<Amesos2::MUMPS,scalar_type> type_map;
105 typedef typename type_map::type slu_type;
106 typedef typename type_map::magnitude_type magnitude_type;
107 typedef typename type_map::MUMPS_STRUC_C MUMPS_STRUC_C;
109 typedef
FunctionMap<Amesos2::MUMPS,slu_type> function_map;
111 MUMPS(Teuchos::RCP<const Matrix> A,
112 Teuchos::RCP<Vector> X,
113 Teuchos::RCP<const Vector> B);
127 int symbolicFactorization_impl();
149 int
solve_impl(const Teuchos::Ptr<MultiVecAdapter<Vector> > X,
150 const Teuchos::Ptr<const MultiVecAdapter<Vector> > B) const;
174 const Teuchos::RCP<Teuchos::ParameterList> & parameterList );
197 int ConvertToTriplet();
199 void MUMPS_ERROR() const;
206 bool MUMPS_MATRIX_LOAD;
208 mutable bool MUMPS_MATRIX_LOAD_PREORDERING;
219 mutable Teuchos::Array<slu_type>
xvals_; local_ordinal_type ldx_;
221 mutable Teuchos::Array<slu_type>
bvals_; local_ordinal_type ldb_;
223 mutable MUMPS_STRUC_C mumps_par;
230template <class Matrix,class Vector>
231class MUMPSNS :
public SolverCore<Amesos2::MUMPS, Matrix, Vector>
239 MUMPSNS(Teuchos::RCP<const Matrix> A,
240 Teuchos::RCP<Vector> X,
241 Teuchos::RCP<const Vector> B);
256 int preOrdering_impl()
259 TEUCHOS_TEST_FOR_EXCEPTION(
true, std::logic_error,
260 "This solver is not support for these types.");
264 int symbolicFactorization_impl()
267 TEUCHOS_TEST_FOR_EXCEPTION(
true, std::logic_error,
268 "This solver is not support for these types.");
272 int numericFactorization_impl()
275 TEUCHOS_TEST_FOR_EXCEPTION(
true, std::logic_error,
276 "This solver is not support for these types.");
280 int solve_impl(
const Teuchos::Ptr<MultiVecAdapter<Vector> > X,
281 const Teuchos::Ptr<
const MultiVecAdapter<Vector> > B)
const
284 TEUCHOS_TEST_FOR_EXCEPTION(
true, std::logic_error,
285 "This solver is not support for these types.");
289 bool matrixShapeOK_impl()
const
292 TEUCHOS_TEST_FOR_EXCEPTION(
true, std::logic_error,
293 "This solver is not support for these types.");
297 void setParameters_impl(
298 const Teuchos::RCP<Teuchos::ParameterList> & parameterList )
301 TEUCHOS_TEST_FOR_EXCEPTION(
true, std::logic_error,
302 "This solver is not support for these types.");
306 Teuchos::RCP<const Teuchos::ParameterList> getValidParameters_impl()
const
309 TEUCHOS_TEST_FOR_EXCEPTION(
true, std::logic_error,
310 "This solver is not support for these types.");
315 bool loadA_impl(
EPhase current_phase)
318 TEUCHOS_TEST_FOR_EXCEPTION(
true, std::logic_error,
319 "This solver is not support for these types.");
324 int ConvertToTriplet()
327 TEUCHOS_TEST_FOR_EXCEPTION(
true, std::logic_error,
328 "This solver is not support for these types.");
333 void MUMPS_ERROR()
const
336 TEUCHOS_TEST_FOR_EXCEPTION(
true, std::logic_error,
337 "This solver is not support for these types.");
356#ifdef HAVE_TEUCHOS_COMPLEX
357 typedef Meta::make_list4<float,
360 std::complex<double> > supported_scalars;
362 typedef Meta::make_list2<float, double> supported_scalars;
Provides a mechanism to map function calls to the correct Solver function based on the scalar type of...
Provides access to interesting solver traits.
Amesos2 interface to the MUMPS package.
Definition Amesos2_MUMPS_decl.hpp:85
bool loadA_impl(EPhase current_phase)
Reads matrix data into internal structures.
Definition Amesos2_MUMPS_def.hpp:387
bool matrixShapeOK_impl() const
Determines whether the shape of the matrix is OK for this solver.
Definition Amesos2_MUMPS_def.hpp:307
int solve_impl(const Teuchos::Ptr< MultiVecAdapter< Vector > > X, const Teuchos::Ptr< const MultiVecAdapter< Vector > > B) const
MUMPS specific solve.
Definition Amesos2_MUMPS_def.hpp:233
Teuchos::Array< slu_type > xvals_
Persisting 1D store for X.
Definition Amesos2_MUMPS_decl.hpp:219
static const char * name
Name of this solver interface.
Definition Amesos2_MUMPS_decl.hpp:92
Teuchos::Array< local_ordinal_type > colptr_
Stores the row indices of the nonzero entries.
Definition Amesos2_MUMPS_decl.hpp:216
int numericFactorization_impl()
MUMPS specific numeric factorization.
Definition Amesos2_MUMPS_def.hpp:203
Teuchos::Array< slu_type > nzvals_
Stores the values of the nonzero entries for MUMPS.
Definition Amesos2_MUMPS_decl.hpp:212
int preOrdering_impl()
Performs pre-ordering on the matrix to increase efficiency.
Definition Amesos2_MUMPS_def.hpp:174
Teuchos::Array< local_ordinal_type > rowind_
Stores the location in Ai_ and Aval_ that starts row j.
Definition Amesos2_MUMPS_decl.hpp:214
Teuchos::RCP< const Teuchos::ParameterList > getValidParameters_impl() const
Definition Amesos2_MUMPS_def.hpp:359
void setParameters_impl(const Teuchos::RCP< Teuchos::ParameterList > ¶meterList)
Definition Amesos2_MUMPS_def.hpp:316
Teuchos::Array< slu_type > bvals_
Persisting 1D store for B.
Definition Amesos2_MUMPS_decl.hpp:221
Amesos2::SolverCore: A templated interface for interaction with third-party direct sparse solvers.
Definition Amesos2_SolverCore_decl.hpp:106
EPhase
Used to indicate a phase in the direct solution.
Definition Amesos2_TypeDecl.hpp:65
Passes functions to TPL functions based on type.
Definition Amesos2_FunctionMap.hpp:77
Map types to solver-specific data-types and enums.
Definition Amesos2_TypeMap.hpp:82
Provides traits about solvers.
Definition Amesos2_SolverTraits.hpp:71