CellProcessor.h

 
#ifndef CELLPROCESSOR_H_
#define CELLPROCESSOR_H_
 
#include <cstddef>
#include <cmath>
 
#include "molecules/MoleculeForwardDeclaration.h"
#include "particleContainer/ParticleCellForwardDeclaration.h"
 
class CellProcessor {
protected:
    double _cutoffRadiusSquare;
    double _LJCutoffRadiusSquare;
 
public:
    CellProcessor(const double cutoffRadius, const double LJCutoffRadius) : 
        _cutoffRadiusSquare(cutoffRadius * cutoffRadius), 
        _LJCutoffRadiusSquare(LJCutoffRadius * LJCutoffRadius) {}
    virtual ~CellProcessor() {}
 
    double getCutoffRadius() const {return sqrt(_cutoffRadiusSquare);}
    double getLJCutoffRadius() const {return sqrt(_LJCutoffRadiusSquare);}
    void setCutoffRadius(const double c) {_cutoffRadiusSquare = c * c;}
    void setLJCutoffRadius(const double ljc) {_LJCutoffRadiusSquare = ljc * ljc;}
 
 
    double getCutoffRadiusSquare() const {return _cutoffRadiusSquare;}
    double getLJCutoffRadiusSquare() const {return _LJCutoffRadiusSquare;}
    void setCutoffRadiusSquare(const double c) {_cutoffRadiusSquare = c;}
    void setLJCutoffRadiusSquare(const double ljc) {_LJCutoffRadiusSquare = ljc;}
 
 
    virtual void initTraversal() = 0;
 
    virtual void preprocessCell(ParticleCell& cell) = 0;
 
    virtual void processCellPair(ParticleCell& cell1, ParticleCell& cell2, bool sumAll = false) = 0;
 
    virtual void processCell(ParticleCell& cell) = 0;
 
    virtual double processSingleMolecule(Molecule* m1, ParticleCell& cell2) = 0;
 
    virtual void postprocessCell(ParticleCell& cell) = 0;
 
    virtual void endTraversal() = 0;
};
 
 
#endif