LinkedCells.h

Go to the documentation of this file.

//
// Created by daniel on 22.05.24.
//
#pragma once
 
#include "../Model.h"
#include "../../particleRepresentation/container/linkedCellsContainer/LinkedCellsContainer.h"
#include "moleculeSimulator/forceCalculation/harmonic/HarmonicForce.h"
 
class LinkedCells final : public Model {
private:
    bool membraneSetting;
    bool pull;
    int pullingActiveUntil;
    std::array<double, 3> pullingForce{};
    std::unique_ptr<Force> forceBetweenDiagonalNeighborsInMembrane;
    std::unique_ptr<Force> forceBetweenDirectNeighborsInMembrane;
    enumsStructs::ParallelizationStrategy parallelizationStrategy;
 
    LinkedCellsContainer particles;
    std::vector<std::pair<Side, enumsStructs::BoundaryCondition> > boundarySettings;
 
    void processBoundaryForces();
 
    void processHaloCells();
 
    void pullMarkedParticles();
 
    void updateForcesMembrane();
 
public:
    LinkedCells(Force &force, double deltaT, std::array<double, 3> domainSize, double rCutOff,
                FileHandler::outputFormat outputFormat, BoundarySet boundaryConditions, bool gravityOn,
                std::array<double, 3> g = {}, MembraneParameters membraneParameters = MembraneParameters{},
                enumsStructs::ParallelizationStrategy parallelizationStrategy =
                        enumsStructs::ParallelizationStrategy::none);
 
    void step(int iteration) override;
 
    void updateForcesOptimized();
 
    void initializeForces() override;
 
 //The following methods are only needed for parallelization.
 
#ifdef _OPENMP
    void updateForcesParallelReduction();
 
    void updateForcesParallelSkipping();
 
    void updateForcesParallelLinear();
 
    void updateVelocitiesParallel();
 
    void updatePositionsParallel();
 
    void processBoundaryForcesParallel();
 
    void initializeReductionVectors(int maxNumThreads);
#endif
};