UnstructuredMeshAnalyzer/GMGW_geom_utils.hxx at develop · cfog/UnstructuredMeshAnalyzer · GitHub

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
/*
 * geom_utils.hxx
 *
 *  Created on: Oct 13, 2017
 *      Author: cfog
 */

/*
 Copyright (C) The University of British Columbia, 2018.

 This file is part of UnstructuredMeshAnalyzer.

 UnstructuredMeshAnalyzer is free software: you can redistribute it
 and/or modify it under the terms of the GNU General Public License
 as published by the Free Software Foundation, either version 3 of
 the License, or (at your option) any later version.

 UnstructuredMeshAnalyzer is distributed in the hope that it will be
 useful, but WITHOUT ANY WARRANTY; without even the implied warranty
 of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 General Public License for more details.

 You should have received a copy of the GNU General Public License
 along with UnstructuredMeshAnalyzer.  If not, see
 <https://www.gnu.org/licenses/>.
 */

#ifndef GMGW_GEOM_UTILS_HXX_
#define GMGW_GEOM_UTILS_HXX_

#include "GMGW_FileWrapper.hxx"

#define DOT(a,b) (a[0]*b[0] + a[1]*b[1] + a[2]*b[2])
#define MAG(a) sqrt(DOT(a,a))
#define NORMALIZE(a) do {double dInvMag = 1./MAG(a); \
a[0]*=dInvMag; a[1]*=dInvMag; a[2]*=dInvMag;} while(0)

#define SUB(a,b) {a[0]-b[0], a[1]-b[1], a[2]-b[2]}
#define CROSS(a,b) {a[1]*b[2]-a[2]*b[1],\
                    a[2]*b[0]-a[0]*b[2],\
                    a[0]*b[1]-a[1]*b[0]}

void
analyzeBdryFace(const double coords[][3], const GMGW_int nConn,
		const GMGW_int connect[], GMGW_int triFaceAngles[],
		GMGW_int quadFaceAngles[], GMGW_int quadDistortion[]);

void
analyzeCellQuality(const double coords[][3], const GMGW_int nConn,
		   const GMGW_int connect[], GMGW_int triFaceAngles[],
		   GMGW_int quadFaceAngles[], GMGW_int quadDistortion[],
		   GMGW_int dihedralsQuadQuad[], GMGW_int dihedralsQuadTri[],
		   GMGW_int dihedralsTriTri[]);

double
cellVolume(const double coords[][3], const GMGW_int nConn,
	   const GMGW_int connect[]);

double
distance(const double a[3], const double b[3]);

void
findClosestPoint(const GMGW_int nConn, const GMGW_int connect[], bool onBdry[],
		 const double coords[][3], double minDist[]);

void
findOffWallSpacing(const FileWrapper* wrapper, const double coords[][3],
		   const GMGW_int nConn, const GMGW_int connect[],
		   const GMGW_int newIndex[], double volSpacing[]);

void
findOnWallSpacing(const double coords[][3], const GMGW_int nConn,
		  const GMGW_int connect[], const GMGW_int newIndex[],
		  double skinSpacing[]);

void
outputAngleHistograms(FILE* angleFile, const GMGW_int quadFaceAngles[],
		      const GMGW_int triFaceAngles[],
		      const GMGW_int dihedralsQuadQuad[],
		      const GMGW_int dihedralsQuadTri[],
		      const GMGW_int dihedralsTriTri[]);

void
outputDistortionHistogram(const char* distortFileName,
			  const GMGW_int quadDistortion[]);

#endif /* GMGW_GEOM_UTILS_HXX_ */