Branch data Line data Source code
1 : : // *****************************************************************************
2 : : /*!
3 : : \file src/DiffEq/Beta/ConfigureMassFractionBeta.cpp
4 : : \copyright 2012-2015 J. Bakosi,
5 : : 2016-2018 Los Alamos National Security, LLC.,
6 : : 2019-2021 Triad National Security, LLC.
7 : : All rights reserved. See the LICENSE file for details.
8 : : \brief Register and compile configuration on the mass fraction beta SDE
9 : : \details Register and compile configuration on the mass fraction beta SDE.
10 : : */
11 : : // *****************************************************************************
12 : :
13 : : #include <set>
14 : : #include <map>
15 : : #include <vector>
16 : : #include <string>
17 : : #include <utility>
18 : :
19 : : #include <brigand/algorithms/for_each.hpp>
20 : :
21 : : #include "Tags.hpp"
22 : : #include "CartesianProduct.hpp"
23 : : #include "DiffEqFactory.hpp"
24 : : #include "Walker/Options/DiffEq.hpp"
25 : : #include "Walker/Options/InitPolicy.hpp"
26 : :
27 : : #include "ConfigureMassFractionBeta.hpp"
28 : : #include "MassFractionBeta.hpp"
29 : : #include "MassFractionBetaCoeffPolicy.hpp"
30 : :
31 : : namespace walker {
32 : :
33 : : void
34 : 374 : registerMassFractionBeta( DiffEqFactory& f, std::set< ctr::DiffEqType >& t )
35 : : // *****************************************************************************
36 : : // Register mass fraction beta SDE into DiffEq factory
37 : : //! \param[in,out] f Differential equation factory to register to
38 : : //! \param[in,out] t Counters for equation types registered
39 : : // *****************************************************************************
40 : : {
41 : : // Construct vector of vectors for all possible policies for SDE
42 : : using MassFractionBetaPolicies =
43 : : tk::cartesian_product< InitPolicies, MassFractionBetaCoeffPolicies >;
44 : : // Register SDE for all combinations of policies
45 : : brigand::for_each< MassFractionBetaPolicies >(
46 : 374 : registerDiffEq< MassFractionBeta >( f, ctr::DiffEqType::MASSFRACBETA, t ) );
47 : 374 : }
48 : :
49 : : std::vector< std::pair< std::string, std::string > >
50 : 4 : infoMassFractionBeta( std::map< ctr::DiffEqType, tk::ctr::ncomp_t >& cnt )
51 : : // *****************************************************************************
52 : : // Return information on the mass fraction beta SDE
53 : : //! \param[inout] cnt std::map of counters for all differential equation types
54 : : //! \return vector of string pairs describing the SDE configuration
55 : : // *****************************************************************************
56 : : {
57 [ + - ]: 4 : auto c = ++cnt[ ctr::DiffEqType::MASSFRACBETA ]; // count eqs
58 : : --c; // used to index vectors starting with 0
59 : :
60 : : std::vector< std::pair< std::string, std::string > > nfo;
61 : :
62 [ + - ][ + - ]: 8 : nfo.emplace_back( ctr::DiffEq().name( ctr::DiffEqType::MASSFRACBETA ), "" );
63 : :
64 [ + - ]: 4 : nfo.emplace_back( "start offset in particle array", std::to_string(
65 [ + - ][ + - ]: 4 : g_inputdeck.get< tag::component >().offset< tag::massfracbeta >(c) ) );
66 : : auto ncomp =
67 : 4 : g_inputdeck.get< tag::component >().get< tag::massfracbeta >()[c] / 3;
68 [ + - ][ + - ]: 4 : nfo.emplace_back( "number of components", std::to_string( ncomp ) );
69 : :
70 [ + - ]: 4 : nfo.emplace_back( "kind", "stochastic" );
71 : 4 : nfo.emplace_back( "dependent variable", std::string( 1,
72 [ + - ][ + - ]: 4 : g_inputdeck.get< tag::param, tag::massfracbeta, tag::depvar >()[c] ) );
73 [ + - ]: 4 : nfo.emplace_back( "initialization policy", ctr::InitPolicy().name(
74 [ + - ]: 8 : g_inputdeck.get< tag::param, tag::massfracbeta, tag::initpolicy >()[c] ) );
75 [ + - ]: 4 : nfo.emplace_back( "coefficients policy", ctr::CoeffPolicy().name(
76 [ + - ]: 8 : g_inputdeck.get< tag::param, tag::massfracbeta, tag::coeffpolicy >()[c] ) );
77 : :
78 [ + - ][ - + ]: 8 : nfo.emplace_back( "random number generator", tk::ctr::RNG().name(
79 [ + - ]: 8 : g_inputdeck.get< tag::param, tag::massfracbeta, tag::rng >()[c] ) );
80 : : nfo.emplace_back(
81 [ + - ][ + - ]: 8 : "coeff b [" + std::to_string( ncomp ) + "]",
[ + - ][ - + ]
[ - + ][ - + ]
[ - - ][ - - ]
[ - - ]
82 [ + - ][ - + ]: 8 : parameters(g_inputdeck.get< tag::param, tag::massfracbeta, tag::b >().at(c))
83 [ + - ]: 4 : );
84 : : nfo.emplace_back(
85 [ + - ][ + - ]: 8 : "coeff S [" + std::to_string( ncomp ) + "]",
[ + - ][ - + ]
[ - + ][ - + ]
[ - - ][ - - ]
[ - - ]
86 [ - + ]: 4 : parameters(
87 [ + - ][ + - ]: 8 : g_inputdeck.get< tag::param, tag::massfracbeta, tag::S >().at(c) ) );
88 : : nfo.emplace_back(
89 [ + - ][ + - ]: 8 : "coeff kappa [" + std::to_string( ncomp ) + "]",
[ + - ][ - + ]
[ - + ][ - + ]
[ - - ][ - - ]
[ - - ]
90 [ - + ]: 4 : parameters(
91 [ + - ][ + - ]: 8 : g_inputdeck.get< tag::param, tag::massfracbeta, tag::kappa >().at(c) ) );
92 : : nfo.emplace_back(
93 [ + - ][ + - ]: 8 : "coeff rho2 [" + std::to_string( ncomp ) + "]",
[ + - ][ - + ]
[ - + ][ - + ]
[ - - ][ - - ]
[ - - ]
94 [ - + ]: 4 : parameters(
95 [ + - ][ + - ]: 8 : g_inputdeck.get< tag::param, tag::massfracbeta, tag::rho2 >().at(c) ) );
96 : : nfo.emplace_back(
97 [ + - ][ + - ]: 8 : "coeff r [" + std::to_string( ncomp ) + "]",
[ + - ][ - + ]
[ - + ][ - + ]
[ - - ][ - - ]
[ - - ]
98 [ - + ]: 4 : parameters(
99 [ + - ][ + - ]: 8 : g_inputdeck.get< tag::param, tag::massfracbeta, tag::r >().at(c) ) );
100 [ + - ]: 4 : spikes( nfo, g_inputdeck.get< tag::param, tag::massfracbeta, tag::init,
101 : : tag::spike >().at(c) );
102 [ + - ]: 4 : betapdfs( nfo, g_inputdeck.get< tag::param, tag::massfracbeta, tag::init,
103 : : tag::betapdf >().at(c) );
104 : :
105 : 4 : return nfo;
106 : : }
107 : :
108 : : } // walker::
|
