ShieldedCoulombForce.h

/**
* @file  ShieldedCoulombForce.h
* @brief Defines the data and methods of the ShieldedCoulombForce class
*
* @license This file is distributed under the BSD Open Source License. 
*          See LICENSE.TXT for details. 
**/

#ifndef SHIELDEDCOULOMBFORCE_H
#define SHIELDEDCOULOMBFORCE_H

#include "Force.h"
#include "VectorCompatibility.h"

class ShieldedCoulombForce : public Force {
public:
	ShieldedCoulombForce(Cloud * const C, const double shieldingConstant);
	~ShieldedCoulombForce();

	void force1(const double currentTime); //rk substep 1
	void force2(const double currentTime); //rk substep 2
	void force3(const double currentTime); //rk substep 3
	void force4(const double currentTime); //rk substep 4

	void writeForce(fitsfile * const file, int * const error) const;
	void readForce(fitsfile * const file, int * const error);

private:
	double shielding; //<! Inverse of shielding distance [m^-1]
    SEMAPHORES
	
	static const double coulomb; //<! Coulomb constant: 8.987551787 [m/F]

	void force(const cloud_index currentParticle,
	           const double charges, const doubleV displacementV);
	void force(const cloud_index currentParticle, const cloud_index iParticle, 
	           const doubleV charges, const doubleV displacementX, const doubleV displacementY);
	void forcer(const cloud_index currentParticle, const cloud_index iParticle,
	            const doubleV charges, const doubleV displacementX, const doubleV displacementY);
    
	static doubleV exp_pd(const int mask, const doubleV a);
    static void plusEqualr_pd(double * const a, const doubleV b);
    static void minusEqualr_pd(double * const a, const doubleV b);
};

#endif // SHIELDEDCOULOMBFORCE_H