galaxydatatypes.h

Go to the documentation of this file.

 
#pragma once
#include "CAdatatypes.h"
#include <algorithm> // copy
#include <vector>
 
class GalaxyCell
{
public:
    int state; 
 
    // additional application specific member variables and static variables
    double velocity[3]; 
    double mass;        
 
    bool operator!=(const GalaxyCell &other);
 
    GalaxyCell &operator=(const GalaxyCell &other);
 
    GalaxyCell(const GalaxyCell &other);
 
    GalaxyCell();
};
 
class Galaxy
{
public:
    int min_mass;                           
    int max_mass;                           
    double density;                         
    int boundary_radius;                    
    int axis1_dim;                          
    int axis2_dim;                          
    int axis3_dim;                          
    static double time_step;                
    static CellularAutomata<GalaxyCell> CA; 
 
    Galaxy();
 
    Galaxy(double time_step, int min_mass, int max_mass, double density, int boundary_radius, int axis1_dim, int axis2_dim, int axis3_dim);
 
    ~Galaxy() = default;
 
    int init_galaxy();
 
    int simulation(int steps);
 
    static void galaxy_formation_rule(int *cell_index, const int index_size,
                                      GalaxyCell *neighborhood_cells, const int neighborhood_size,
                                      GalaxyCell &new_cell_state);
 
    static void set_new_position(int *cell_index, GalaxyCell &new_cell_state,
                                 const std::vector<double> &displacement_vector);
 
    static bool did_galaxies_collide(int *cell_index, const std::vector<int> &offset_index, GalaxyCell &new_cell_state);
 
    static void update_velocity_after_collision(GalaxyCell &new_cell, GalaxyCell collied_cell);
 
    static GalaxyCell &get_cell_state(int *cell_index, const std::vector<int> &offset_index);
 
    static std::vector<int> get_periodic_vector(int *cell_index, const std::vector<int> &offset_index);
 
    static int round_int(double dbl);
 
    static std::vector<double> compute_gravitational_force(const GalaxyCell &cell_of_interest, const GalaxyCell &neighbor_cell, const std::vector<int> &neighbor_index);
 
    static std::vector<double> compute_accel(const std::vector<double> &total_force, double mass);
 
    static std::vector<double> compute_velocity(const std::vector<double> &accel, const GalaxyCell &cell_of_interest, double time_step);
 
    static std::vector<double> compute_displacement(const std::vector<double> &velocity, const GalaxyCell &cell_of_interest, double time_step);
 
    static double compute_vector_norm(const std::vector<int> &vector);
 
    static std::vector<double> compute_vector_difference(const std::vector<double> &vec1, const std::vector<double> &vec2);
};