CAutils.h

Go to the documentation of this file.

 
#pragma once
#include <utility> // pair
#include <fstream>
 
template <typename T>
void swap_states(T *vector, T *next_vector, int axis1_dim)
{
    for (int i = 0; i < axis1_dim; i++)
    {
        std::swap(vector[i], next_vector[i]);
        next_vector[i] = T();
    }
}
 
template <typename T>
void swap_states(T **matrix, T **next_matrix, int axis1_dim, int axis2_dim)
{
    for (int i = 0; i < axis1_dim; i++)
    {
        for (int j = 0; j < axis2_dim; j++)
        {
            std::swap(matrix[i][j], next_matrix[i][j]);
            next_matrix[i][j] = T();
        }
    }
}
 
template <typename T>
void swap_states(T ***tensor, T ***next_tensor, int axis1_dim, int axis2_dim, int axis3_dim)
{
    for (int i = 0; i < axis1_dim; i++)
    {
        for (int j = 0; j < axis2_dim; j++)
        {
            for (int k = 0; k < axis3_dim; k++)
            {
                std::swap(tensor[i][j][k], next_tensor[i][j][k]);
                next_tensor[i][j][k] = T();
            }
        }
    }
}
 
bool is_diagonal_neighboring_cell_2d(int i, int j);
 
bool is_diagonal_neighboring_cell_3d(int i, int j, int k);
 
bool less_than_votes(const std::pair<int, int> &a, const std::pair<int, int> &b);
 
int get_periodic_index(int i, int di, int axis_dim);
 
void get_periodic_moore_neighbor_index(int rank, int radius, int neighborhood_array_index, int *neighbor_index);
 
void get_periodic_von_neumann_neighbor_index(int rank, int radius, int neighborhood_array_index, int *neighbor_index);
 
void get_density(std::ifstream& data, std::ofstream& result);