GridFire 0.0.1a
General Purpose Nuclear Network
Loading...
Searching...
No Matches
screening_weak.cpp
Go to the documentation of this file.
2
3#include "fourdst/composition/atomicSpecies.h"
4
5#include "cppad/cppad.hpp"
6
7#include <vector>
8
9
10namespace gridfire::screening {
11 using ADDouble = CppAD::AD<double>;
13 const reaction::LogicalReactionSet &reactions,
14 const std::vector<fourdst::atomic::Species>& species,
15 const std::vector<ADDouble> &Y,
16 const ADDouble T9,
17 const ADDouble rho
18 ) const {
19 return calculateFactors_impl<ADDouble>(reactions, species, Y, T9, rho);
20 }
21
23 const reaction::LogicalReactionSet &reactions,
24 const std::vector<fourdst::atomic::Species>& species,
25 const std::vector<double> &Y,
26 const double T9,
27 const double rho
28 ) const {
29 return calculateFactors_impl<double>(reactions, species, Y, T9, rho);
30 }
31}
CppAD::AD< double > ADDouble
Alias for CppAD Automatic Differentiation type for double precision.
std::vector< T > calculateFactors_impl(const reaction::LogicalReactionSet &reactions, const std::vector< fourdst::atomic::Species > &species, const std::vector< T > &Y, const T T9, const T rho) const
Template implementation for calculating weak screening factors.
std::vector< double > calculateScreeningFactors(const reaction::LogicalReactionSet &reactions, const std::vector< fourdst::atomic::Species > &species, const std::vector< double > &Y, const double T9, const double rho) const override
Calculates weak screening factors for a set of reactions.
TemplatedReactionSet< LogicalReaction > LogicalReactionSet
A set of logical reactions.
Definition reaction.h:557
CppAD::AD< double > ADDouble