/usr/include/liggghts/region_wedge.h is in libliggghts-dev 3.3.1+repack1-1ubuntu3.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 | /* ----------------------------------------------------------------------
This is the
██╗ ██╗ ██████╗ ██████╗ ██████╗ ██╗ ██╗████████╗███████╗
██║ ██║██╔════╝ ██╔════╝ ██╔════╝ ██║ ██║╚══██╔══╝██╔════╝
██║ ██║██║ ███╗██║ ███╗██║ ███╗███████║ ██║ ███████╗
██║ ██║██║ ██║██║ ██║██║ ██║██╔══██║ ██║ ╚════██║
███████╗██║╚██████╔╝╚██████╔╝╚██████╔╝██║ ██║ ██║ ███████║
╚══════╝╚═╝ ╚═════╝ ╚═════╝ ╚═════╝ ╚═╝ ╚═╝ ╚═╝ ╚══════╝®
DEM simulation engine, released by
DCS Computing Gmbh, Linz, Austria
http://www.dcs-computing.com, office@dcs-computing.com
LIGGGHTS® is part of CFDEM®project:
http://www.liggghts.com | http://www.cfdem.com
Core developer and main author:
Christoph Kloss, christoph.kloss@dcs-computing.com
LIGGGHTS® is open-source, distributed under the terms of the GNU Public
License, version 2 or later. It is distributed in the hope that it will
be useful, but WITHOUT ANY WARRANTY; without even the implied warranty
of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. You should have
received a copy of the GNU General Public License along with LIGGGHTS®.
If not, see http://www.gnu.org/licenses . See also top-level README
and LICENSE files.
LIGGGHTS® and CFDEM® are registered trade marks of DCS Computing GmbH,
the producer of the LIGGGHTS® software and the CFDEM®coupling software
See http://www.cfdem.com/terms-trademark-policy for details.
-------------------------------------------------------------------------
Contributing author and copyright for this file:
Stefan Amberger (JKU Linz)
Christoph Kloss (DCS Computing GmbH, Linz)
Christoph Kloss (JKU Linz)
Copyright 2012- DCS Computing GmbH, Linz
Copyright 2009-2012 JKU Linz
------------------------------------------------------------------------- */
#ifdef REGION_CLASS
RegionStyle(wedge,RegWedge)
#else
#ifndef domain_region_wedge_h
#define domain_region_wedge_h
#include "region.h"
// command: region id wedge dim c1 c2 radius lo hi angle1 angle2
// dim ... x or y or z
// c1 ... if dim == x then y-coord
// if dim == y then z-coord
// if dim == z then x-coord
// c2 ... if dim == x then z-coord
// if dim == y then x-coord
// if dim == z then y-coord
// radius ... radius of cylinder
// lo ... dim-coord of lower flat face of cylinder of wedge
// hi ... dim-coord of higher flat face of cylinder of wedge
// angle1 ... mathematically positive angle of the wedge's starting face
// if axis == x then starting from y-axis
// if axis == y then starting from z-axis
// if axis == z then starting from x-axis
// dang ... angle between the wedge's ending face
namespace LAMMPS_NS
{
class RegWedge : public Region
{
friend class DomainWedge;
public:
RegWedge(class LAMMPS *, int, char **);
~RegWedge();
int inside(double,double,double);
int surface_exterior(double *, double);
int surface_interior(double *, double);
private:
// internal data
char axis;
double c1,c2;
double radius;
double lo,hi;
double angle1, cosang1, sinang1;
double dang, cosdang, sindang, cosmdang, sinmdang;
double angle2, cosang2, sinang2;
// helper data
double pi_half;
double onedivr;
// normal vectors
// the normal vectors point outside the wedge and are normalized
double normal1[2]; // normal vector to 1st plane of section (angle1)
double normal2[2]; // normal vector to 2nd plane of section (angle2)
// helper functions
inline void snormalize2(const double length, const double *v, double *ans);
void printRegion();
void printProperty(const char *name, double val);
void printContacts(double *, int);
};
}
#endif
#endif
|