netgen/libsrc/meshing/improve3.hpp

147 lines
5.0 KiB
C++
Raw Normal View History

2009-01-13 04:40:13 +05:00
#ifndef FILE_IMPROVE3
#define FILE_IMPROVE3
2011-07-25 17:33:19 +06:00
extern double CalcTotalBad (const Mesh::T_POINTS & points,
2019-08-09 03:23:12 +05:00
const Array<Element> & elements,
2011-07-25 17:33:19 +06:00
const MeshingParameters & mp);
2009-01-13 04:40:13 +05:00
///
class MeshOptimize3d
{
2011-07-25 17:33:19 +06:00
const MeshingParameters & mp;
2009-01-13 04:40:13 +05:00
public:
2011-07-25 17:33:19 +06:00
MeshOptimize3d (const MeshingParameters & amp) : mp(amp) { ; }
2019-08-13 21:47:35 +05:00
double CombineImproveEdge (Mesh & mesh, const MeshingParameters & mp,
Table<ElementIndex, PointIndex> & elements_of_point,
2019-08-13 21:47:35 +05:00
Array<double> & elerrs, PointIndex pi0, PointIndex pi1,
FlatArray<bool, PointIndex> is_point_removed, bool check_only=false);
2009-01-13 04:40:13 +05:00
void CombineImprove (Mesh & mesh, OPTIMIZEGOAL goal = OPT_QUALITY);
2019-08-13 21:47:35 +05:00
void CombineImproveSequential (Mesh & mesh, OPTIMIZEGOAL goal = OPT_QUALITY);
2009-01-13 04:40:13 +05:00
void SplitImprove (Mesh & mesh, OPTIMIZEGOAL goal = OPT_QUALITY);
2019-10-25 17:57:00 +05:00
void SplitImproveSequential (Mesh & mesh, OPTIMIZEGOAL goal = OPT_QUALITY);
double SplitImproveEdge (Mesh & mesh, OPTIMIZEGOAL goal, Table<ElementIndex,PointIndex> & elementsonnode, Array<double> &elerrs, NgArray<INDEX_3> &locfaces, double badmax, PointIndex pi1, PointIndex pi2, PointIndex ptmp, bool check_only=false);
void SplitImprove2 (Mesh & mesh);
double SplitImprove2Element (Mesh & mesh, ElementIndex ei, const Table<ElementIndex, PointIndex> & elements_of_point, const Array<double> & elerrs, bool check_only);
2019-10-25 17:57:00 +05:00
2019-09-24 18:12:39 +05:00
double SwapImproveEdge (Mesh & mesh, OPTIMIZEGOAL goal, const NgBitArray * working_elements, Table<ElementIndex,PointIndex> & elementsonnode, INDEX_3_HASHTABLE<int> & faces, PointIndex pi1, PointIndex pi2, bool check_only=false);
2009-01-13 04:40:13 +05:00
void SwapImprove (Mesh & mesh, OPTIMIZEGOAL goal = OPT_QUALITY,
2019-08-28 17:00:49 +05:00
const NgBitArray * working_elements = NULL);
2019-09-24 18:12:39 +05:00
void SwapImproveSequential (Mesh & mesh, OPTIMIZEGOAL goal = OPT_QUALITY,
const NgBitArray * working_elements = NULL);
2009-01-13 04:40:13 +05:00
void SwapImproveSurface (Mesh & mesh, OPTIMIZEGOAL goal = OPT_QUALITY,
2019-08-28 17:00:49 +05:00
const NgBitArray * working_elements = NULL,
2019-07-09 13:39:16 +05:00
const NgArray< NgArray<int,PointIndex::BASE>* > * idmaps = NULL);
2019-09-06 19:51:34 +05:00
void SwapImprove2Sequential (Mesh & mesh, OPTIMIZEGOAL goal = OPT_QUALITY);
2009-01-13 04:40:13 +05:00
void SwapImprove2 (Mesh & mesh, OPTIMIZEGOAL goal = OPT_QUALITY);
double SwapImprove2 ( Mesh & mesh, OPTIMIZEGOAL goal, ElementIndex eli1, int face, Table<ElementIndex, PointIndex> & elementsonnode, TABLE<SurfaceElementIndex, PointIndex::BASE> & belementsonnode, bool check_only=false );
2009-01-13 04:40:13 +05:00
2011-07-25 17:33:19 +06:00
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, mp);
return 0;
}
2011-07-25 17:33:19 +06:00
2011-07-25 17:33:19 +06:00
double CalcTotalBad (const Mesh::T_POINTS & points,
2019-08-09 03:23:12 +05:00
const Array<Element> & elements)
2011-07-25 17:33:19 +06:00
{
return netgen::CalcTotalBad (points, elements, mp);
}
};
2009-01-13 04:40:13 +05:00
inline double
CalcBad (const Mesh::T_POINTS & points, const Element & elem, double h, const MeshingParameters & mp)
{
if (elem.GetType() == TET)
return CalcTetBadness (points[elem[0]], points[elem[1]],
points[elem[2]], points[elem[3]], h, mp);
return 0;
}
2009-01-13 04:40:13 +05:00
extern int WrongOrientation (const Mesh::T_POINTS & points, const Element & el);
/* Functional depending of inner point inside triangular surface */
class MinFunctionSum : public MinFunction
{
protected:
2019-07-09 13:39:16 +05:00
NgArray<MinFunction*> functions;
2009-01-13 04:40:13 +05:00
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;
2019-07-09 13:39:16 +05:00
const NgArray<INDEX_3> & faces;
2011-07-25 17:33:19 +06:00
const MeshingParameters & mp;
2009-01-13 04:40:13 +05:00
double h;
public:
PointFunction1 (Mesh::T_POINTS & apoints,
2019-07-09 13:39:16 +05:00
const NgArray<INDEX_3> & afaces,
2011-07-25 17:33:19 +06:00
const MeshingParameters & amp,
2009-01-13 04:40:13 +05:00
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;
2019-08-09 03:23:12 +05:00
const Array<Element> & elements;
2009-01-13 04:40:13 +05:00
TABLE<INDEX> elementsonpoint;
PointIndex actpind;
bool onplane;
Vec<3> nv;
public:
JacobianPointFunction (Mesh::T_POINTS & apoints,
2019-08-09 03:23:12 +05:00
const Array<Element> & aelements);
2014-08-19 17:01:15 +06:00
virtual ~JacobianPointFunction () { ; }
2009-01-13 04:40:13 +05:00
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