source: molecuilder/src/tesselationhelpers.hpp@ a9b2a0a

Last change on this file since a9b2a0a was a9b2a0a, checked in by Frederik Heber <heber@…>, 16 years ago

Huge refactoring to make const what is const (ticket #38), continued.

  • too many changes because of too many cross-references to be able to list them up here.
  • NOTE that "make check" runs fine and did catch several error.
  • note that we had to use const_iterator several times when the map, ... was declared const.
  • at times we changed an allocated LinkedCell LCList(...) into

const LinkedCell *LCList;
LCList = new LinkedCell(...);

  • also mutable (see ticket #5) was used, e.g. for molecule::InternalPointer (PointCloud changes are allowed, because they are just accounting).

Signed-off-by: Frederik Heber <heber@…>
  • Property mode set to 100644
File size: 3.2 KB
Line 
1/*
2 * TesselationHelpers.hpp
3 *
4 * Auxiliary functions for the tesselation.
5 *
6 * Created on: Aug 3, 2009
7 * Author: heber
8 */
9
10#ifndef TESSELATIONHELPERS_HPP_
11#define TESSELATIONHELPERS_HPP_
12
13using namespace std;
14
15/*********************************************** includes ***********************************/
16
17// include config.h
18#ifdef HAVE_CONFIG_H
19#include <config.h>
20#endif
21
22#include <gsl/gsl_linalg.h>
23#include <gsl/gsl_matrix.h>
24#include <gsl/gsl_multimin.h>
25#include <gsl/gsl_permutation.h>
26#include <gsl/gsl_vector.h>
27
28#include <iostream>
29
30#include "defs.hpp"
31
32/****************************************** forward declarations *****************************/
33
34class BoundaryPointSet;
35class BoundaryLineSet;
36class BoundaryTriangleSet;
37class LinkedCell;
38class TesselPoint;
39class PointCloud;
40class Tesselation;
41class Vector;
42
43/********************************************** definitions *********************************/
44
45#define HULLEPSILON 1e-10
46
47/********************************************** declarations *******************************/
48
49double DetGet(gsl_matrix *A, int inPlace);
50void GetSphere(Vector *center, Vector &a, Vector &b, Vector &c, double RADIUS);
51void GetCenterOfSphere(Vector* Center, Vector a, Vector b, Vector c, Vector *NewUmkreismittelpunkt, Vector* Direction, Vector* AlternativeDirection, double HalfplaneIndicator, double AlternativeIndicator, double alpha, double beta, double gamma, double RADIUS, double Umkreisradius);
52void GetCenterofCircumcircle(Vector *Center, Vector *a, Vector *b, Vector *c);
53double GetPathLengthonCircumCircle(Vector &CircleCenter, Vector &CirclePlaneNormal, double CircleRadius, Vector &NewSphereCenter, Vector &OldSphereCenter, Vector &NormalVector, Vector &SearchDirection);
54double MinIntersectDistance(const gsl_vector * x, void *params);
55bool existsIntersection(Vector point1, Vector point2, Vector point3, Vector point4);
56double CalculateVolumeofGeneralTetraeder(Vector *a, Vector *b, Vector *c, Vector *d);
57double GetAngle(const Vector &point, const Vector &reference, const Vector OrthogonalVector);
58
59bool CheckLineCriteriaForDegeneratedTriangle(const BoundaryPointSet *nodes[3]);
60bool SortCandidates(const CandidateForTesselation* candidate1, const CandidateForTesselation* candidate2);
61TesselPoint* FindClosestPoint(const Vector* Point, TesselPoint *&SecondPoint, const LinkedCell* const LC);
62TesselPoint* FindSecondClosestPoint(const Vector*, const LinkedCell*);
63Vector * GetClosestPointBetweenLine(ofstream *out, const BoundaryLineSet *Base, const BoundaryLineSet *OtherBase);
64
65void WriteTecplotFile(ofstream *out, ofstream *tecplot, const Tesselation *TesselStruct, const PointCloud *cloud, const int N);
66void WriteRaster3dFile(ofstream *out, ofstream *rasterfile, const Tesselation *Tess, const PointCloud *cloud);
67void IncludeSphereinRaster3D(ofstream *out, ofstream *rasterfile, const Tesselation *Tess, const PointCloud *cloud);
68void WriteVrmlFile(ofstream *out, ofstream *vrmlfile, const Tesselation *Tess, const PointCloud *cloud);
69void CalculateConcavityPerBoundaryPoint(ofstream *out, const Tesselation *TesselStruct);
70double DistanceToTrianglePlane(ofstream *out, const Vector *x, const BoundaryTriangleSet *triangle);
71
72bool CheckListOfBaselines(ofstream *out, const Tesselation *TesselStruct);
73
74
75#endif /* TESSELATIONHELPERS_HPP_ */
Note: See TracBrowser for help on using the repository browser.