//#define XERR
#include "tumorinfo.ih"

    // see void calcTotalRisk(int carrierIndex)

void TumorInfo::cumTotalRisk(DoubleVect const &radiationRisk)
{
    double sum = 0;
    d_cumRisk.clear();

xerr("0");
    uint16_t idx = d_incidence.index();

    auto iter = radiationRisk.begin();

    for (double risk: d_incidence.tumorRisk(idx))
    {
        //g_log << risk << ' ' << *iter << '\n';
        d_cumRisk.push_back(sum += risk * *iter++);
    }

xerr("0");
}


//void calcTotalRisk(int carrierIndex) {
//    float totRisk;
//    for (int T = 0; T < 101; T++) {
//        totRisk = radRisk[T] * tumorRisk[carrierIndex][T];  <<------ !!!
//        if (T == 0) {
//            cumTotRisk[T] = totRisk;
//        }
//        else {
//            cumTotRisk[T] = cumTotRisk[T-1] + totRisk;
//        }
//    }
//
