Add a Tpetra-based version of Trilinos SparseMatrix and SparsityPattern.

dealii · Jan 9, 2024 · f0eaee6 · f0eaee6
1 parent 239f11e
commit f0eaee6
Show file tree

Hide file tree

Showing 12 changed files with 3,955 additions and 1 deletion.
diff --git a/cmake/config/template-arguments.in b/cmake/config/template-arguments.in
@@ -242,6 +242,8 @@ TRIANGULATIONS := { Triangulation<deal_II_dimension, deal_II_space_dimension>;
                     parallel::distributed::Triangulation<deal_II_dimension, deal_II_space_dimension>;
                     parallel::fullydistributed::Triangulation<deal_II_dimension, deal_II_space_dimension>; }
 
+TRILINOS_SCALARS := { @DEAL_II_EXPAND_TPETRA_TYPES@; }
+
 // all supported logical dimensions
 DIMENSIONS := { 1; 2; 3 }
 

diff --git a/cmake/configure/configure_20_trilinos.cmake b/cmake/configure/configure_20_trilinos.cmake
@@ -481,6 +481,7 @@ macro(feature_trilinos_configure_external)
 
   if(${DEAL_II_TRILINOS_WITH_TPETRA})
     if(DEAL_II_HAVE_TPETRA_INST_DOUBLE)
+      set(DEAL_II_EXPAND_TPETRA_TYPES "double")
       set(DEAL_II_EXPAND_TPETRA_VECTOR_DOUBLE "LinearAlgebra::TpetraWrappers::Vector<double>")
     endif()
     if(DEAL_II_HAVE_TPETRA_INST_FLOAT)

diff --git a/doc/news/changes/major/20231122SebastianKinnewig b/doc/news/changes/major/20231122SebastianKinnewig
@@ -0,0 +1,5 @@
+New: LinearAlgebra::TpetraWrappers::SparseMatrix class
+that implements a wrapper for Tpetra::CrsMatrix.
+<br>
+(Sebastian Kinnewig, 2023/11/22)
+
diff --git a/include/deal.II/base/template_constraints.h b/include/deal.II/base/template_constraints.h
@@ -670,7 +670,10 @@ namespace LinearAlgebra
   {
     template <typename Number>
     class Vector;
-  }
+
+    template <typename Number, typename NodeType>
+    class SparseMatrix;
+  } // namespace TpetraWrappers
 #  endif
 } // namespace LinearAlgebra
 #endif

diff --git a/include/deal.II/lac/affine_constraints.templates.h b/include/deal.II/lac/affine_constraints.templates.h
@@ -47,6 +47,8 @@
 #include <deal.II/lac/trilinos_block_sparse_matrix.h>
 #include <deal.II/lac/trilinos_parallel_block_vector.h>
 #include <deal.II/lac/trilinos_sparse_matrix.h>
+#include <deal.II/lac/trilinos_tpetra_sparse_matrix.h>
+#include <deal.II/lac/trilinos_tpetra_vector.h>
 #include <deal.II/lac/trilinos_vector.h>
 
 #include <boost/serialization/complex.hpp>