Context Navigation

tesselationhelpers.hpp@ 60cbe5

Visit:

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

Various fixes and attempt to get convex hull working.

Tesselation::FlipBaseline: Now returns pointer to new baseline or NULL if not flipped
Tesselation::PickFarthestofTwoBaselines() returns volume change due to flipping
WriteTecplotFile(): ZONE name now contains nodes of last triangle, incorporates upper two changes
ConvexizeNonconvexEnvelope(): additional dummy string to add single step output of intermediate envelopes
TesselStruct::RemovePointFromTesselatedSurface(): Now returns volume lost due to point removal
ConvexizeNonconvexEnvelope() evaluates change in volume due to added/removed points
FindNonConvexBorder(): last two parameters have changed position, we check for NULL filename before writing, we add straddling points by InsertStraddlingPoints() (is better than AddBoundaryPointByDegeneratedTriangle()), IsInnerAtom() test removed
new function AddBoundaryPointByDegeneratedTriangle() to add outside points
StoreTrianglesinFile() contains the triangle storing
boundary.cpp: typos, Verbose() changes and Verbosity reduced
new function BoundaryTriangleSet::IsPresentTupel(), checks whether two triangles are the same, used in Tesselation::FindTriangles()
new attribute Tesselation::LastTriangle contains pointer to last added triangle
FIX: Tesselation::FindNextSuitableTriangle() checks for existentTriangleCount between 1 and 3 (as is the number of already present lines) and checks for zero NormalVector of created triangle
Tesselation::FindSecondPointForTesselation(): removed parameter Candidate
Tesselation::FindClosestTrianglesToPoint() now has a list of triangles from the nearest point which we look through for nearest triangle
BUGFIX: Tesselation::IsInnerPoint() now takes center of triangle with respect to point to decide whether in- or outside
FIX: Tesselation::GetPathsOfConnectedPoints() notes down touched triangles as well in order move along the surrounding triangles sequentially
FIX: Tesselation::RemovePointFromTesselatedSurface(): Huge changes, from the closed path of connected points we remove one by one by adding the respective triangle and try to do it in such an order, that the volume is maximized

new function Tesselation::FindAllDegeneratedTriangles(), pendant to Tesselation::FindAllDegeneratedLines()
FIX: new function Tesselation::RemoveDegeneratedTriangles()
new function Tesselation::AddBoundaryPointByDegeneratedTriangle()

Moved tesselation writing to tesselation.cpp

WriteVrmlFile(), WriteRaster3dFile(), WriteTecplotFile() moved to tesselation.cpp
these also don't use molecule anymore, but the PointCloud structure (hence, bonds are no more visible)
new function TesselStruct::Output() which performs the writing
new function InsertSphereInRaster3D(), which places a sphere at the position of the last triangle

Property mode set to 100644

File size: 2.4 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
13	using namespace std;
14
15	// include config.h
16	#ifdef HAVE_CONFIG_H
17	#include <config.h>
18	#endif
19
20	#include <gsl/gsl_linalg.h>
21	#include <gsl/gsl_matrix.h>
22	#include <gsl/gsl_multimin.h>
23	#include <gsl/gsl_permutation.h>
24	#include <gsl/gsl_vector.h>
25
26	#include "defs.hpp"
27	#include "tesselation.hpp"
28	#include "vector.hpp"
29
30	#define HULLEPSILON 1e-10
31
32	double DetGet(gsl_matrix *A, int inPlace);
33	void GetSphere(Vector *center, Vector &a, Vector &b, Vector &c, double RADIUS);
34	void 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);
35	void GetCenterofCircumcircle(Vector Center, Vector a, Vector b, Vector c);
36	double GetPathLengthonCircumCircle(Vector &CircleCenter, Vector &CirclePlaneNormal, double CircleRadius, Vector &NewSphereCenter, Vector &OldSphereCenter, Vector &NormalVector, Vector &SearchDirection);
37	double MinIntersectDistance(const gsl_vector * x, void *params);
38	bool existsIntersection(Vector point1, Vector point2, Vector point3, Vector point4);
39	double CalculateVolumeofGeneralTetraeder(Vector a, Vector b, Vector c, Vector d);
40	double GetAngle(const Vector &point, const Vector &reference, const Vector OrthogonalVector);
41
42	bool CheckLineCriteriaForDegeneratedTriangle(class BoundaryPointSet *nodes[3]);
43	bool SortCandidates(class CandidateForTesselation* candidate1, class CandidateForTesselation* candidate2);
44	TesselPoint* FindClosestPoint(const Vector* Point, TesselPoint &SecondPoint, LinkedCell LC);
45	TesselPoint* FindSecondClosestPoint(const Vector, LinkedCell);
46	Vector * GetClosestPointBetweenLine(ofstream out, class BoundaryLineSet Base, class BoundaryLineSet *OtherBase);
47
48	void WriteTecplotFile(ofstream out, ofstream tecplot, class Tesselation TesselStruct, PointCloud cloud, int N);
49	void WriteRaster3dFile(ofstream out, ofstream rasterfile, class Tesselation Tess, PointCloud cloud);
50	void IncludeSphereinRaster3D(ofstream out, ofstream rasterfile, class Tesselation Tess, PointCloud cloud);
51	void WriteVrmlFile(ofstream out, ofstream vrmlfile, class Tesselation Tess, PointCloud cloud);
52
53
54	#endif /* TESSELATIONHELPERS_HPP_ */

Note: See TracBrowser for help on using the repository browser.

Context Navigation

source: molecuilder/src/tesselationhelpers.hpp@ 60cbe5

Download in other formats: