From: David Wells <drwells@email.unc.edu>
Date: Tue, 21 Mar 2023 21:17:48 +0000 (-0400)
Subject: Unconditionally instantiate complex FullMatrix.
X-Git-Tag: v9.5.0-rc1~436^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F14939%2Fhead;p=dealii.git

Unconditionally instantiate complex FullMatrix.

Like the previous commit - we need it for complex eigenvalues.
---

diff --git a/source/lac/full_matrix.cc b/source/lac/full_matrix.cc
index b0c19ccb38..63f50b2631 100644
--- a/source/lac/full_matrix.cc
+++ b/source/lac/full_matrix.cc
@@ -23,11 +23,14 @@ DEAL_II_NAMESPACE_OPEN
 
 #ifndef DOXYGEN
 
-#  ifndef DEAL_II_WITH_COMPLEX_VALUES
-// instantiate for std::complex<double> because we use it internally in
-// FESeries.
+// We use FullMatrix<std::complex<T>> for complex eigenvalues so ignore the
+// value of DEAL_II_WITH_COMPLEX_VALUES and always instantiate. As a consequence
+// we cannot use REAL_AND_COMPLEX_SCALARS without getting duplicate
+// instantiations so do float and double here too.
+template class FullMatrix<float>;
+template class FullMatrix<double>;
+template class FullMatrix<std::complex<float>>;
 template class FullMatrix<std::complex<double>>;
-#  endif
 
 // instantiate for long double manually because we use it in a few places
 // inside the library
@@ -74,12 +77,10 @@ FullMatrix<long double>::add<long double>(const long double,
 TEMPL_OP_EQ(double, float);
 TEMPL_OP_EQ(float, double);
 
-#  ifdef DEAL_II_WITH_COMPLEX_VALUES
 TEMPL_OP_EQ(std::complex<double>, std::complex<float>);
 TEMPL_OP_EQ(std::complex<float>, std::complex<double>);
 TEMPL_OP_EQ(std::complex<double>, double);
 TEMPL_OP_EQ(std::complex<float>, float);
-#  endif
 
 #  undef TEMPL_OP_EQ
 
diff --git a/source/lac/full_matrix.inst.in b/source/lac/full_matrix.inst.in
index 90b48df135..23e334aeb8 100644
--- a/source/lac/full_matrix.inst.in
+++ b/source/lac/full_matrix.inst.in
@@ -17,8 +17,6 @@
 
 for (S : REAL_AND_COMPLEX_SCALARS)
   {
-    template class FullMatrix<S>;
-
     template void FullMatrix<S>::print(LogStream &,
                                        const unsigned int,
                                        const unsigned int) const;