/usr/share/netgen/libsrc/meshing/improve3.hpp is in netgen-headers 4.9.13.dfsg-8build2.
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 | #ifndef FILE_IMPROVE3
#define FILE_IMPROVE3
///
class MeshOptimize3d
{
public:
void CombineImprove (Mesh & mesh, OPTIMIZEGOAL goal = OPT_QUALITY);
void SplitImprove (Mesh & mesh, OPTIMIZEGOAL goal = OPT_QUALITY);
void SwapImprove (Mesh & mesh, OPTIMIZEGOAL goal = OPT_QUALITY,
const BitArray * working_elements = NULL);
void SwapImproveSurface (Mesh & mesh, OPTIMIZEGOAL goal = OPT_QUALITY,
const BitArray * working_elements = NULL,
const Array< Array<int,PointIndex::BASE>* > * idmaps = NULL);
void SwapImprove2 (Mesh & mesh, OPTIMIZEGOAL goal = OPT_QUALITY);
};
inline double
CalcBad (const Mesh::T_POINTS & points, const Element & elem,
double h)
{
if (elem.GetType() == TET)
return CalcTetBadness (points[elem[0]], points[elem[1]],
points[elem[2]], points[elem[3]], h);
return 0;
}
extern double CalcTotalBad (const Mesh::T_POINTS & points,
const Mesh::T_VOLELEMENTS & elements);
extern int WrongOrientation (const Mesh::T_POINTS & points, const Element & el);
/* Functional depending of inner point inside triangular surface */
class MinFunctionSum : public MinFunction
{
protected:
Array<MinFunction*> functions;
public:
virtual double Func (const Vector & x) const;
virtual void Grad (const Vector & x, Vector & g) const;
virtual double FuncGrad (const Vector & x, Vector & g) const;
virtual double FuncDeriv (const Vector & x, const Vector & dir, double & deriv) const;
virtual double GradStopping (const Vector & x) const;
void AddFunction(MinFunction & fun);
const MinFunction & Function(int i) const;
MinFunction & Function(int i);
};
class PointFunction1 : public MinFunction
{
Mesh::T_POINTS & points;
const Array<INDEX_3> & faces;
double h;
public:
PointFunction1 (Mesh::T_POINTS & apoints,
const Array<INDEX_3> & afaces,
double ah);
virtual double Func (const Vector & x) const;
virtual double FuncDeriv (const Vector & x, const Vector & dir, double & deriv) const;
virtual double FuncGrad (const Vector & x, Vector & g) const;
virtual double GradStopping (const Vector & x) const;
};
class JacobianPointFunction : public MinFunction
{
public:
Mesh::T_POINTS & points;
const Mesh::T_VOLELEMENTS & elements;
TABLE<INDEX> elementsonpoint;
PointIndex actpind;
bool onplane;
Vec<3> nv;
public:
JacobianPointFunction (Mesh::T_POINTS & apoints,
const Mesh::T_VOLELEMENTS & aelements);
virtual void SetPointIndex (PointIndex aactpind);
virtual double Func (const Vector & x) const;
virtual double FuncGrad (const Vector & x, Vector & g) const;
virtual double FuncDeriv (const Vector & x, const Vector & dir, double & deriv) const;
inline void SetNV(const Vec<3> & anv) {nv = anv; onplane = true;}
inline void UnSetNV(void) {onplane = false;}
};
#endif
|