PrimaryCommunityPopulator.cpp

Go to the documentation of this file.

/*
 *  This is free software: you can redistribute it and/or modify it
 *  under the terms of the GNU General Public License as published by
 *  the Free Software Foundation, either version 3 of the License, or
 *  any later version.
 *  The software is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *  You should have received a copy of the GNU General Public License
 *  along with the software. If not, see <http://www.gnu.org/licenses/>.
 *
 *  Copyright 2018, 2019, Jan Broeckhove and Bistromatics group.
 */

#include "Populator.h"

#include "contact/ContactPool.h"
#include "geopop/GeoGrid.h"
#include "geopop/Location.h"
#include "pop/Person.h"

using namespace std;
using namespace stride::ContactType;

namespace geopop {

template<>
void Populator<stride::ContactType::Id::PrimaryCommunity>::Apply(GeoGrid& geoGrid, const GeoGridConfig&)
{
        m_logger->trace("Starting to populate Primary Communities");

        for (const shared_ptr<Location>& loc : geoGrid) {
                if (loc->GetPopCount() == 0) {
                        continue;
                }

                // 1. find all communities in an area of 10-k*10 km
                const auto nearbyPools = geoGrid.GetNearbyPools(Id::PrimaryCommunity, *loc);

                if (nearbyPools.empty()) { // apparently no commumity pools present, so nothing to populate
                        return;
                }

                // 2. for every household assign a community
                const auto dist = m_rn_man.GetUniformIntGenerator(0, static_cast<int>(nearbyPools.size()), 0U);
                for (auto& hhPool : loc->RefPools(Id::Household)) {
                        for (auto p : *hhPool) {
                                auto& pcPool = nearbyPools[dist()];
                                pcPool->AddMember(p);
                                p->SetPoolId(Id::PrimaryCommunity, pcPool->GetId());
                        }
                }
        }

        m_logger->trace("Done populating Primary Communities");
}

} // namespace geopop