``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64``` ```// ***************************************************************************** /*! \file src/PDE/MultiMat/Problem/UserDefined.hpp \copyright 2012-2015 J. Bakosi, 2016-2018 Los Alamos National Security, LLC., 2019-2021 Triad National Security, LLC. All rights reserved. See the LICENSE file for details. \brief Problem configuration for the multi-material compressible flow equations \details This file defines a Problem policy class for the multi-material compressible flow equations, defined under PDE/MultiMat. See PDE/MultiMat/Problem.h for general requirements on Problem policy classes for MultiMat. */ // ***************************************************************************** #ifndef MultiMatProblemUserDefined_h #define MultiMatProblemUserDefined_h #include #include "Types.hpp" #include "Inciter/InputDeck/InputDeck.hpp" #include "FunctionPrototypes.hpp" #include "Inciter/Options/Problem.hpp" #include "MultiMat/MultiMatIndexing.hpp" namespace inciter { extern ctr::InputDeck g_inputdeck; //! MultiMat system of PDEs problem: user defined class MultiMatProblemUserDefined { private: using ncomp_t = tk::ctr::ncomp_t; using eq = tag::multimat; public: //! Initialize numerical solution static tk::InitializeFn::result_type initialize( ncomp_t system, ncomp_t ncomp, tk::real, tk::real, tk::real, tk::real ); //! Evaluate analytical solution at (x,y,z,t) for all components static std::vector< tk::real > analyticSolution( ncomp_t system, ncomp_t ncomp, tk::real x, tk::real y, tk::real z, tk::real t ) { return initialize( system, ncomp, x, y, z, t ); } //! Compute and return source term for Rayleigh-Taylor manufactured solution //! \details No-op for user-deefined problems. static tk::MultiMatSrcFn::result_type src( ncomp_t, ncomp_t ncomp, tk::real, tk::real, tk::real, tk::real ) { std::vector< tk::real > s( ncomp, 0.0 ); }<--- Variable 's' is assigned a value that is never used. //! Return names of integral variables to be output to diagnostics file static std::vector< std::string > names( ncomp_t ); static ctr::ProblemType type() noexcept { return ctr::ProblemType::USER_DEFINED; } }; } // inciter:: #endif // MultiMatProblemUserDefined_h ```