/usr/include/liggghts/modified_andrew.h is in libliggghts-dev 2.3.8-1build1.
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 | /* ----------------------------------------------------------------------
LIGGGHTS - LAMMPS Improved for General Granular and Granular Heat
Transfer Simulations
LIGGGHTS is part of the CFDEMproject
www.liggghts.com | www.cfdem.com
Christoph Kloss, christoph.kloss@cfdem.com
Copyright 2009-2012 JKU Linz
Copyright 2012- DCS Computing GmbH, Linz
LIGGGHTS is based on LAMMPS
LAMMPS - Large-scale Atomic/Molecular Massively Parallel Simulator
http://lammps.sandia.gov, Sandia National Laboratories
Steve Plimpton, sjplimp@sandia.gov
This software is distributed under the GNU General Public License.
See the README file in the top-level directory.
------------------------------------------------------------------------- */
/* ----------------------------------------------------------------------
Contributing authors:
Stefan Amberger (JKU Linz)
------------------------------------------------------------------------- */
#ifndef LMP_MODIFIED_ANDREW_H
#define LMP_MODIFIED_ANDREW_H
#include <vector>
#include "pointers.h"
using namespace std;
namespace MODIFIED_ANDREW_AUX{
struct Point {
double x, y;
bool operator <(const Point &p) const {
return x < p.x || (x == p.x && y < p.y);
}
};
struct Circle {
double x, y, r;
};
}
using MODIFIED_ANDREW_AUX::Point;
using MODIFIED_ANDREW_AUX::Circle;
namespace LAMMPS_NS{
class ModifiedAndrew : protected Pointers {
public:
ModifiedAndrew(LAMMPS *lmp);
~ModifiedAndrew();
double area();
void add_contact(Circle c);
inline void clear_contacts()
{ contacts_.clear(); }
private:
double area(vector<Point> H);
double area(Point &p, Point &m, Point &q);
double cross(Point O, Point A, Point B);
Point mean_point(vector<Point> P);
vector<Point> convex_hull(vector<Point> P);
vector<Point> construct_hull_c_all(double *data0, int ndata0);
int construct_data(vector<Point> hull_c, double *&data);
// container for contacts: x, y
vector<Point> contacts_;
int npoints_per_circle_;
double **directions_;
};
}
#endif
|