physics.h

#ifndef PHYSICS_H
#define PHYSICS_H
#include <math.h>
#include "particle.h"
#include "constants.h"
#include <boost/random.hpp>

#define _USE_MATH_DEFINES

class Physics
{
public:
    Physics();
    ~Physics();

    Particle* getCenterOfGravity();

    bool isOrbital();

    virtual void updatePositionAndVelocity(Particle* p, int timestep) = 0;

    bool particlesOverlap(Particle* a, Particle* b);

    void resolveCollision(Particle* lP, Particle* rP, int timestep);

    std::valarray<double> getRelativeVelocity(Particle* a, Particle* b);

    double getMagnitude(std::valarray<double> vel);

    double calculateKE(Particle* p);

    virtual double calculatePE(Particle* p);

    virtual void initializeAcceleration(Particle* p){};

    double randDouble(double min, double max);

    double randGamma(double mean, double var);

    int randBit();

protected:
    bool orbital;
    Particle* centerOfGravity;

    boost::mt19937 rng;

};

#endif // PHYSICS_H