SHLocalParticle.h

/*
 * SHLocalParticle.h
 *
 *  Created on: Nov 27, 2014
 *      Author: tchipevn
 */
 
#ifndef SHLOCALPARTICLE_H_
#define SHLOCALPARTICLE_H_
 
#include "bhfmm/pseudoParticles/LocalParticle.h"
#include "bhfmm/expansions/SolidHarmonicsExpansion.h"
#include "bhfmm/utils/RotationParameter.h"
#include <map>
 
namespace bhfmm {
 
class SHLocalParticle: public LocalParticle {
public:
    SHLocalParticle(int order, bool initializeExpansionToZero = true);
    ~~SHLocalParticle();
 
    void addSource(const Vector3<double>& position, double charge);
 
    void addMultipoleParticle(const MultipoleParticle& multipole, Vector3<double> periodicShift);
 
    void addMultipoleParticle_Wigner(const MultipoleParticle& multipole, Vector3<double> periodicShift,
        double* cellWid,
        std::map<Vector3<int>, RotationParams, Vector3<int>::compare>& M2L_Wigner);
 
    void actOnLocalParticle(LocalParticle& small) const;
 
    void actOnLocalParticle_Wigner(LocalParticle& small,
        const WignerMatrix& W_pos, const WignerMatrix& W_neg, const double* CosSinPhi, const int negate, const double& magnitude) const;
 
    void actOnTarget(const Vector3<double>& position, double charge, double& potential, Vector3<double>& force) const;
 
    void clear();
    int getNumEntries() const;
 
 
    const SolidHarmonicsExpansion& getConstExpansion() const {
        return _expansionM;
    }
 
    SolidHarmonicsExpansion& getExpansion() {
        return _expansionM;
    }
 
    void writeValuesToMPIBuffer(double * buf, int& position) const {
        _expansionM.writeValuesToMPIBuffer(buf, position);
    }
 
    void readValuesFromMPIBuffer(double * buf, int& position) {
        _expansionM.readValuesFromMPIBuffer(buf,position);
    }
    void addValuesFromMPIBuffer(double * buf, int& position) {
        _expansionM.addValuesFromMPIBuffer(buf,position);
    }
 
private:
    SolidHarmonicsExpansion _expansionM;
};
 
} /* namespace bhfmm */
 
#endif /* SHLOCALPARTICLE_H_ */