> > so 'a' is of type Sacado::Fad::DFad<double>. It then needs to call > > numbers::is_finite (Sacado::Fad::DFad<double>), which doesn't exist. > It probably just tries to go through all of the overloads of > numbers::is_finite() and wants to see whether it can convert the > argument Sacado::Fad::DFad<double> to the argument type of these > overloads. The error message you show then would just explain why this > one possibility (namely, numbers::is_finite(std::complex<long double>&)) > does not work. But that doesn't mean that this is the right overload > anyway -- I suspect that your compiler produces similar error messages > above or below the one you show for all of the other overloads, right? > > True, the error message gets long pretty quickly, but the undefined is_finite was another issue. Even if is_finite exists, the complex constructor is still an issue.
> I *think* that the solution is to simply provide an overload for > numbers::is_finite (const Sacado::Fad::DFad<double> &x) > Can you try this? You could declare it in your own .cc file before you > #include <deal.II/lac/la_parallel_vector.templates.h> > Attached is a tiny .cc file that simply instantiates the vector. Copy pasted here here convenience. #include <Sacado.hpp> namespace dealii{ namespace numbers{ bool is_finite(const Sacado::Fad::DFad<double> &x) { ¦ (void) x; ¦ return true; } }} #include <deal.II/lac/vector.h> #include <deal.II/lac/la_parallel_vector.templates.h> int main (int /*argc*/, char * /*argv*/[]){ using namespace dealii; dealii::LinearAlgebra::distributed::Vector<double> vector_double; using ADtype = Sacado::Fad::DFad<double>; dealii::LinearAlgebra::distributed::Vector<ADtype> vector_ad; return 0; } I also provide an function for numbers::is_finite (const Sacado::Fad::DFad<double> &x) to avoid the first set of errors. However, I still get the error below In file included from /home/ddong/Libraries/dealii/install/include/deal.II/base/aligned_vector.h:22:0, from /home/ddong/Libraries/dealii/install/include/deal.II/lac/vector.h:22, from /home/ddong/Codes/PHiLiP/src/instantiate_vector_ad.cpp:13: /home/ddong/Libraries/dealii/install/include/deal.II/lac/la_parallel_vector.templates.h: In instantiation of ‘dealii::LinearAlgebra::distributed::Vector<Number, MemorySpace>& dealii::LinearAlgebra::distributed::Vector<Number, MemorySpace>::operator*=(Number) [with Number = Sacado::Fad::DFad<double>; MemorySpace = dealii::MemorySpace::Host]’: /home/ddong/Codes/PHiLiP/src/instantiate_vector_ad.cpp:26:1: required from here /home/ddong/Libraries/dealii/install/include/deal.II/lac/la_parallel_vector.templates.h:1651:7: error: no matching function for call to ‘std::complex<double>::complex(const Sacado::Fad::DFad<double>&)’ AssertIsFinite(factor); In file included from /usr/include/trilinos/Teuchos_ConfigDefs.hpp:94:0, from /usr/include/trilinos/Teuchos_PromotionTraits.hpp:45, from /usr/include/trilinos/Sacado_Fad_Exp_GeneralFadTraits.hpp:139, from /usr/include/trilinos/Sacado.hpp:52, from /home/ddong/Codes/PHiLiP/src/instantiate_vector_ad.cpp:1: /usr/include/c++/7/complex:1512:3: note: candidate: constexpr std::complex<double>::complex(const std::complex<long double>&) complex<double>::complex(const complex<long double>& __z) ^~~~~~~~~~~~~~~ /usr/include/c++/7/complex:1512:3: note: no known conversion for argument 1 from ‘const Sacado::Fad::DFad<double>’ to ‘const std::complex<long double>&’ /usr/include/c++/7/complex:1219:26: note: candidate: constexpr std::complex<double>::complex(const std::complex<float>&) _GLIBCXX_CONSTEXPR complex(const complex<float>& __z) ^~~~~~~ /usr/include/c++/7/complex:1219:26: note: no known conversion for argument 1 from ‘const Sacado::Fad::DFad<double>’ to ‘const std::complex<float>&’ /usr/include/c++/7/complex:1209:26: note: candidate: constexpr std::complex<double>::complex(double, double) _GLIBCXX_CONSTEXPR complex(double __r = 0.0, double __i = 0.0) ^~~~~~~ /usr/include/c++/7/complex:1209:26: note: no known conversion for argument 1 from ‘const Sacado::Fad::DFad<double>’ to ‘double’ /usr/include/c++/7/complex:1207:26: note: candidate: constexpr std::complex<double>::complex(std::complex<double>::_ComplexT) _GLIBCXX_CONSTEXPR complex(_ComplexT __z) : _M_value(__z) { } ^~~~~~~ /usr/include/c++/7/complex:1207:26: note: no known conversion for argument 1 from ‘const Sacado::Fad::DFad<double>’ to ‘std::complex<double>::_ComplexT {aka __complex__ double}’ /usr/include/c++/7/complex:1202:12: note: candidate: constexpr std::complex<double>::complex(const std::complex<double>&) struct complex<double> ^~~~~~~~~~~~~~~ /usr/include/c++/7/complex:1202:12: note: no known conversion for argument 1 from ‘const Sacado::Fad::DFad<double>’ to ‘const std::complex<double>&’ /usr/include/c++/7/complex:1202:12: note: candidate: constexpr std::complex<double>::complex(std::complex<double>&&) /usr/include/c++/7/complex:1202:12: note: no known conversion for argument 1 from ‘const Sacado::Fad::DFad<double>’ to ‘std::complex<double>&&’ In file included from /home/ddong/Libraries/dealii/install/include/deal.II/base/aligned_vector.h:22:0, from /home/ddong/Libraries/dealii/install/include/deal.II/lac/vector.h:22, from /home/ddong/Codes/PHiLiP/src/instantiate_vector_ad.cpp:13: So, even if the is_finite function exists, the exception to be thrown still requires the complex constructor std::complex<double>::complex(const Sacado::Fad::DFad<double>&) to exist. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/7c3d83a7-154c-4bd2-a736-6c5c83f6f829%40googlegroups.com.
#include <Sacado.hpp> namespace dealii{ namespace numbers{ bool is_finite(const Sacado::Fad::DFad<double> &x) { (void) x; return true; } } } #include <deal.II/lac/vector.h> #include <deal.II/lac/la_parallel_vector.templates.h> int main (int /*argc*/, char * /*argv*/[]) { using namespace dealii; dealii::LinearAlgebra::distributed::Vector<double> vector_double; using ADtype = Sacado::Fad::DFad<double>; dealii::LinearAlgebra::distributed::Vector<ADtype> vector_ad; return 0; }