rad_option.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_RAD_OPTION_HH
#define VOROPP_RAD_OPTION_HH
 
#include <cmath>
 
namespace voro {
 
class radius_mono {
    public:
        struct radius_info // intentionally empty; no special radius info needs to be passed in for radius_mono case
        {
        };
    protected:
        inline void r_init(int ijk,int s,radius_info& radius_out) const {}
        inline void r_prime(double rv,radius_info& radius_out) const {}
        inline bool r_ctest(double crs,double mrs,const radius_info& rad_info) const {return crs>mrs;}
        inline double r_cutoff(double lrs, const radius_info& rad_info) const {return lrs;}
        inline double r_max_add(double rs) const {return rs;}
        inline double r_current_sub(double rs,int ijk,int q) const {return rs;}
        inline double r_scale(double rs,int ijk,int q, const radius_info& rad_info) const {return rs;}
        inline bool r_scale_check(double &rs,double mrs,int ijk,int q, const radius_info& rad_info) const {return rs<mrs;}
};
 
class radius_poly {
    public:
        struct radius_info
        {
            double r_rad, r_mul, r_val;
        };
        double **ppr;
        double max_radius;
        radius_poly() : max_radius(0) {}
    protected:
        inline void r_init(int ijk, int s, radius_info& radius_out) const {
            radius_out.r_rad=ppr[ijk][4*s+3]*ppr[ijk][4*s+3];
            radius_out.r_mul= radius_out.r_rad-max_radius*max_radius;
        }
        inline void r_prime(double rv, radius_info& radius_out) const {
            radius_out.r_val = 1 + radius_out.r_mul / rv;
        }
        inline bool r_ctest(double crs, double mrs, const radius_info& rad_info) const {return crs+ rad_info.r_mul>sqrt(mrs*crs);}
        inline double r_cutoff(double lrs, const radius_info& rad_info) const {return lrs* rad_info.r_val;}
        inline double r_max_add(double rs) const {return rs+max_radius*max_radius;}
        inline double r_current_sub(double rs,int ijk,int q) const {
            return rs-ppr[ijk][4*q+3]*ppr[ijk][4*q+3];
        }
        inline double r_scale(double rs,int ijk,int q, const radius_info& rad_info) const {
            return rs+ rad_info.r_rad-ppr[ijk][4*q+3]*ppr[ijk][4*q+3];
        }
        inline bool r_scale_check(double &rs,double mrs,int ijk,int q, const radius_info& rad_info) const {
            double trs=rs;
            rs+= rad_info.r_rad-ppr[ijk][4*q+3]*ppr[ijk][4*q+3];
            return rs<sqrt(mrs*trs);
        }
        
};
 
}
#endif