54 #ifndef AMESOS2_SHYLUBASKER_DECL_HPP
55 #define AMESOS2_SHYLUBASKER_DECL_HPP
58 #include "Amesos2_SolverCore.hpp"
59 #include "Amesos2_ShyLUBasker_FunctionMap.hpp"
61 #include "shylubasker_decl.hpp"
62 #include "shylubasker_def.hpp"
63 #include "shylubasker_trilinos_decl.hpp"
75 template <
class Matrix,
class Vector>
92 typedef typename super_type::scalar_type scalar_type;
93 typedef typename super_type::local_ordinal_type local_ordinal_type;
94 typedef typename super_type::global_ordinal_type global_ordinal_type;
95 typedef typename super_type::global_size_type global_size_type;
96 typedef typename super_type::node_type node_type;
99 typedef typename type_map::type shylubasker_type;
100 typedef typename type_map::dtype shylubasker_dtype;
102 typedef typename type_map::type slu_type;
106 typedef Matrix matrix_type;
109 typedef Kokkos::DefaultHostExecutionSpace HostExecSpaceType;
111 typedef Kokkos::View<local_ordinal_type*, HostExecSpaceType> host_ordinal_type_array;
112 typedef Kokkos::View<shylubasker_type*, HostExecSpaceType> host_value_type_array;
116 Teuchos::RCP<Vector> X,
117 Teuchos::RCP<const Vector> B);
138 int symbolicFactorization_impl();
160 int
solve_impl(const Teuchos::Ptr<MultiVecAdapter<Vector> > X,
161 const Teuchos::Ptr<const MultiVecAdapter<Vector> > B) const;
170 void setParameters_impl(
171 const Teuchos::RCP<Teuchos::ParameterList> & parameterList );
208 typedef typename Kokkos::View<shylubasker_type**, Kokkos::LayoutLeft,
209 typename HostExecSpaceType::memory_space> host_solve_array_t;
221 #if defined( HAVE_AMESOS2_KOKKOS ) && defined( KOKKOS_ENABLE_OPENMP )
228 typedef Kokkos::OpenMP Exe_Space;
229 ::BaskerNS::BaskerTrilinosInterface<local_ordinal_type, shylubasker_dtype, Exe_Space> *ShyLUbasker;
231 #pragma message("Amesos_ShyLUBasker_decl Error: ENABLED SHYLU_NODEBASKER BUT NOT KOKKOS or NOT OPENMP!")
242 #ifdef HAVE_TEUCHOS_COMPLEX
249 typedef Meta::make_list6<float,
251 Kokkos::complex<float>,
252 Kokkos::complex<double>,
254 std::complex<double> > supported_scalars;
256 typedef Meta::make_list2<float, double> supported_scalars;
260 template <
typename Scalar,
typename LocalOrdinal,
typename ExecutionSpace>
262 KokkosSparse::CrsMatrix<Scalar, LocalOrdinal, ExecutionSpace>> {
263 static const bool value =
true;
268 #endif // AMESOS2_SHYLUBASKER_DECL_HPP
Amesos2::SolverCore: A templated interface for interaction with third-party direct sparse solvers...
Definition: Amesos2_SolverCore_decl.hpp:105
int solve_impl(const Teuchos::Ptr< MultiVecAdapter< Vector > > X, const Teuchos::Ptr< const MultiVecAdapter< Vector > > B) const
ShyLUBasker specific solve.
Definition: Amesos2_ShyLUBasker_def.hpp:325
host_ordinal_type_array rowind_view_
Stores the location in Ai_ and Aval_ that starts row j.
Definition: Amesos2_ShyLUBasker_decl.hpp:201
host_value_type_array nzvals_view_
Stores the values of the nonzero entries for Umfpack.
Definition: Amesos2_ShyLUBasker_decl.hpp:199
Map types to solver-specific data-types and enums.
Definition: Amesos2_TypeMap.hpp:82
host_solve_array_t bValues_
Persisting 1D store for B.
Definition: Amesos2_ShyLUBasker_decl.hpp:216
Amesos2 interface to the Baker package.
Definition: Amesos2_ShyLUBasker_decl.hpp:76
Teuchos::RCP< const Teuchos::ParameterList > getValidParameters_impl() const
Definition: Amesos2_ShyLUBasker_def.hpp:526
int preOrdering_impl()
Performs pre-ordering on the matrix to increase efficiency.
Definition: Amesos2_ShyLUBasker_def.hpp:141
Provides traits about solvers.
Definition: Amesos2_SolverTraits.hpp:71
host_ordinal_type_array colptr_view_
Stores the row indices of the nonzero entries.
Definition: Amesos2_ShyLUBasker_decl.hpp:203
bool matrixShapeOK_impl() const
Determines whether the shape of the matrix is OK for this solver.
Definition: Amesos2_ShyLUBasker_def.hpp:400
bool loadA_impl(EPhase current_phase)
Reads matrix data into internal structures.
Definition: Amesos2_ShyLUBasker_def.hpp:591
A Matrix adapter interface for Amesos2.
Definition: Amesos2_MatrixAdapter_decl.hpp:76
std::string name() const override
Return the name of this solver.
Definition: Amesos2_SolverCore_def.hpp:759
int numericFactorization_impl()
ShyLUBasker specific numeric factorization.
Definition: Amesos2_ShyLUBasker_def.hpp:225
Passes functions to TPL functions based on type.
Definition: Amesos2_FunctionMap.hpp:76
host_solve_array_t xValues_
Persisting 1D store for X.
Definition: Amesos2_ShyLUBasker_decl.hpp:212
Provides access to interesting solver traits.
bool single_proc_optimization() const
can we optimize size_type and ordinal_type for straight pass through, also check that is_contiguous_ ...
Definition: Amesos2_ShyLUBasker_def.hpp:135