SubConfiguration.cpp

Go to the documentation of this file.

/*
 * SubConfiguration.cpp
 *
 *  Created on: Dec 16, 2019
 *      Author: Nikhil
 */

#include "SubConfiguration.h"

SubConfiguration::SubConfiguration(const Stencil& stencil) :
    Configuration(stencil.particleContributingMap.size(), stencil.parent.coordinates.at(Reference).rows()>0),
    parent(stencil.parent)
{
    int numberOfParticlesInParent= parent.numberOfParticles;

    particleContributing.resize(numberOfParticles);
    particleContributing.setZero();


    int referenceAndFinal= (coordinates.at(Reference).rows()>0);

    int i_localParticle=0;
    for(const auto& particleContributingPair : stencil.particleContributingMap)
    {
        int particle= particleContributingPair.first;
        int contributing= particleContributingPair.second;

        auto result= globalLocalMap.insert({particle,i_localParticle});
        assert(result.second && "Insertion failed in the creation of a glocalLocalMap");
        particleContributing[i_localParticle]= contributing;

        this->coordinates.at(Current).row(i_localParticle)= parent.coordinates.at(Current).row(particle);
        if (referenceAndFinal) this->coordinates.at(Reference).row(i_localParticle)= parent.coordinates.at(Reference).row(particle);
        this->velocities.row(i_localParticle)= parent.velocities.row(particle);
        this->species.push_back(parent.species[particle]);
        i_localParticle++;
    }

}

SubConfiguration::~SubConfiguration() {
    // TODO Auto-generated destructor stub
}