pre_container.hh

Go to the documentation of this file.

// Copyright Epic Games, Inc. All Rights Reserved.
// Modified version of Voro++'s source file
 
// Voro++, a 3D cell-based Voronoi library
//
// Author   : Chris H. Rycroft (LBL / UC Berkeley)
// Email    : chr@alum.mit.edu
// Date     : August 30th 2011
 
#ifndef VOROPP_PRE_CONTAINER_HH
#define VOROPP_PRE_CONTAINER_HH
 
#include <cstdio>
 
#include "c_loops.hh"
#include "container.hh"
 
#include "HAL/PlatformMisc.h"
 
namespace voro {
 
void guess_optimal(int siteCount, double sizeX, double sizeY, double sizeZ, int &nx, int &ny, int &nz);
 
class pre_container_base {
    public:
        const double ax;
        const double bx;
        const double ay;
        const double by;
        const double az;
        const double bz;
        const bool xperiodic;
        const bool yperiodic;
        const bool zperiodic;
        void guess_optimal(int &nx,int &ny,int &nz);
        pre_container_base(double ax_,double bx_,double ay_,double by_,double az_,double bz_,bool xperiodic_,bool yperiodic_,bool zperiodic_,int ps_);
        ~pre_container_base();
PRAGMA_DISABLE_UNSAFE_TYPECAST_WARNINGS
        inline int total_particles() {
            return (end_id-pre_id)*pre_container_chunk_size+(ch_id-*end_id);
        }
PRAGMA_RESTORE_UNSAFE_TYPECAST_WARNINGS
    protected:
        const int ps;
        void new_chunk();
        void extend_chunk_index();
        int index_sz;
        int **pre_id;
        int **end_id;
        int **l_id;
        int *ch_id;
        int *e_id;
        double **pre_p;
        double **end_p;
        double *ch_p;
};
 
class pre_container : public pre_container_base {
    public:
        pre_container(double ax_,double bx_,double ay_,double by_,double az_,double bz_,
                bool xperiodic_,bool yperiodic_,bool zperiodic_)
            : pre_container_base(ax_,bx_,ay_,by_,az_,bz_,xperiodic_,yperiodic_,zperiodic_,3) {};
        void put(int n,double x,double y,double z);
        void setup(container &con);
        void setup(particle_order &vo,container &con);
};
 
class pre_container_poly : public pre_container_base {
    public:
        pre_container_poly(double ax_,double bx_,double ay_,double by_,double az_,double bz_,
                bool xperiodic_,bool yperiodic_,bool zperiodic_)
            : pre_container_base(ax_,bx_,ay_,by_,az_,bz_,xperiodic_,yperiodic_,zperiodic_,4) {};
        void put(int n,double x,double y,double z,double r);
        void setup(container_poly &con);
        void setup(particle_order &vo,container_poly &con);
};
 
}
 
#endif