Molecule.h

#include <random>
#include <chrono>
#include <vector>
#include <tuple>
#include <algorithm>

#include "EspeceMoleculaire.h"


extern int diametre;
extern std::vector<EspeceMoleculaire*> especes;

int maxTaille();

class Environnement;

class Molecule
{
    private:
        double x, y, z; // Coordonnées actuelles
        double lastX, lastY, lastZ; // Coordonées avant le déplacement
        EspeceMoleculaire *espece; // Espèce de la molecule
        static std::minstd_rand generator; // Moteur de génération de nombres aléatoires

    public:
        bool traitee; // Si la molécule a été traitée
        Molecule(EspeceMoleculaire *e);
        Molecule(EspeceMoleculaire *e, double x, double y, double z);
        ~Molecule(){}
        double getX(){return x;}
        double getY(){return y;}
        double getZ(){return z;}
        double getR(){double sqR = x*x + y*y + z*z; return sqrt(sqR);} // Renvoie la distance par rapport au centre de l'environnement
        EspeceMoleculaire *getEspece(){return espece;}
        bool move(Environnement *env);
        void unmove(Environnement *env);

};

class Environnement
{
    private: 
        std::vector<std::vector<std::vector<std::vector<Molecule *>>>> env3D;
        int tailleCube; // Largeur du cube en nombre de petits cubes

    public:
        Environnement();
        ~Environnement(){}
        std::vector<Molecule *>& getListCoords(double x, double y, double z);
        std::vector<Molecule *>& getListIndices(int i, int j, int k);
        std::tuple<int, int, int> coords2Indices(double x, double y, double z);
        std::vector<std::vector<std::vector<Molecule *>>>& operator[](size_t i){ return env3D[i];}
        int cubeSize(){return env3D.size();} // Renvoie la taille d'un coté du cube
        void ajoutMolecule(Molecule *m);
        std::vector<Molecule *>& findMolecule(Molecule *m);
        void removeMolecule(Molecule *m);
};