background.cpp

////////////////////////////////////////////////////////////////////////////////
// File: background.cpp                                                       //
// Project: respond                                                           //
// Created Date: 2026-02-05                                                   //
// Author: Matthew Carroll                                                    //
// -----                                                                      //
// Last Modified: 2026-02-12                                                  //
// Modified By: Matthew Carroll                                               //
// -----                                                                      //
// Copyright (c) 2026 Syndemics Lab at Boston Medical Center                  //
////////////////////////////////////////////////////////////////////////////////

#include "internals/background.hpp"

#include <memory>
#include <string>

namespace respond {
Eigen::VectorXd
BackgroundDeath::Execute(const Eigen::VectorXd &state,
                         std::map<std::string, History> &h) const {
    if (GetTransitionMatrices().size() != 1) {
        throw std::runtime_error(
            "Mortality Transitions must have 1 Transition Matrix.");
    }
    auto deaths =
        state.cwiseProduct(GetTransitionMatrices()[0]); // calculate the deaths
    if (h.find("background_death") != h.end()) {
        h["background_death"].AddState(deaths);
    }
    if (!(state.array() >= deaths.array()).all()) {
        std::runtime_error(
            "The state is not larger than the estimated background deaths!");
    }
    auto new_state = state - deaths; // remove deaths from state

    if (h.find("background_death") != h.end()) {
        h["state"].AddState(new_state);
    }
    return new_state;
}

std::unique_ptr<Transition>
BackgroundDeath::Create(const std::string &name, const std::string &log_name) {
    return std::make_unique<BackgroundDeath>(name, log_name);
}
} // namespace respond