calculate_score.hpp

Go to the documentation of this file.

// SPDX-License-Identifier: BSD-2-Clause
 
#ifndef cMHN_COMMON_CALCULATE_SCORE_H
#define cMHN_COMMON_CALCULATE_SCORE_H
 
#include <map>
 
#include <cmhn/common/calculate_pTheta.hpp>
#include <cmhn/common/regulator.hpp>
#include <cmhn/common/score.hpp>
#include <cmhn/nontt/mhn_operator.hpp>
#include <cmhn/tt/mhn_operator.hpp>
 
#include <prc.hpp>
 
namespace cMHN
{
    template<class T, pRC::Size D, class S>
    T calculateScore(nonTT::MHNOperator<T, D> const &op,
        std::map<S, T> const &pD, cMHN::Score<T> const &Score,
        cMHN::Regulator<T, D> const &Regulator, T const &toleranceSolverP = 1e-4)
    {
        T score = pRC::zero();
 
        auto const pTheta = calculatePTheta(op, toleranceSolverP);
 
        for(auto const &[k, v] : pD)
        {
            score += Score.pointwiseScore(v, pTheta(k));
        }
 
        score -= Regulator.score(log(op.theta()));
 
        return score;
    }
 
    template<pRC::Size RP, class T, pRC::Size D, class S>
    T calculateScore(TT::MHNOperator<T, D> const &op,
        std::map<S, T> const &pD, cMHN::Score<T> const &Score,
        cMHN::Regulator<T, D> const &Regulator, T const &toleranceSolverP = 1e-4)
    {
        T score = pRC::zero();
 
        auto const pTheta = calculatePTheta<RP>(op, toleranceSolverP);
 
        for(auto const &[k, v] : pD)
        {
            score += Score.pointwiseScore(v, pTheta(k));
        }
 
        score -= Regulator.score(log(op.theta()));
 
        return score;
    }
}
 
#endif // cMHN_COMMON_CALCULATE_SCORE_H