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
110
111
112
|
diff --git a/V_HexMetric.cpp b/V_HexMetric.cpp
index fda4771..74eecb2 100644
--- a/V_HexMetric.cpp
+++ b/V_HexMetric.cpp
@@ -2974,10 +2974,10 @@ double hex_distortion(int num_nodes, const double coordinates[][3])
double weight[maxTotalNumberGaussPoints];
// create an object of GaussIntegration
- GaussIntegration gint{};
- gint.initialize(number_of_gauss_points, num_nodes, number_dimension);
- gint.calculate_shape_function_3d_hex();
- gint.get_shape_func(shape_function[0], dndy1[0], dndy2[0], dndy3[0], weight);
+ GaussIntegration gaussint;
+ gaussint.initialize(number_of_gauss_points, num_nodes, number_dimension);
+ gaussint.calculate_shape_function_3d_hex();
+ gaussint.get_shape_func(shape_function[0], dndy1[0], dndy2[0], dndy3[0], weight);
VerdictVector xxi, xet, xze, xin;
@@ -3014,7 +3014,7 @@ double hex_distortion(int num_nodes, const double coordinates[][3])
double dndy2_at_node[maxNumberNodes][maxNumberNodes];
double dndy3_at_node[maxNumberNodes][maxNumberNodes];
- gint.calculate_derivative_at_nodes_3d(dndy1_at_node, dndy2_at_node, dndy3_at_node);
+ gaussint.calculate_derivative_at_nodes_3d(dndy1_at_node, dndy2_at_node, dndy3_at_node);
int node_id;
for (node_id = 0; node_id < num_nodes; node_id++)
{
diff --git a/V_QuadMetric.cpp b/V_QuadMetric.cpp
index 2486146..68af002 100644
--- a/V_QuadMetric.cpp
+++ b/V_QuadMetric.cpp
@@ -1409,10 +1409,10 @@ double quad_distortion(int num_nodes, const double coordinates[][3])
double weight[maxTotalNumberGaussPoints];
// create an object of GaussIntegration
- GaussIntegration gint{};
- gint.initialize(number_of_gauss_points, num_nodes);
- gint.calculate_shape_function_2d_quad();
- gint.get_shape_func(shape_function[0], dndy1[0], dndy2[0], weight);
+ GaussIntegration gaussint;
+ gaussint.initialize(number_of_gauss_points, num_nodes);
+ gaussint.calculate_shape_function_2d_quad();
+ gaussint.get_shape_func(shape_function[0], dndy1[0], dndy2[0], weight);
// calculate element area
int ife, ja;
@@ -1435,7 +1435,7 @@ double quad_distortion(int num_nodes, const double coordinates[][3])
double dndy1_at_node[maxNumberNodes][maxNumberNodes];
double dndy2_at_node[maxNumberNodes][maxNumberNodes];
- gint.calculate_derivative_at_nodes(dndy1_at_node, dndy2_at_node);
+ gaussint.calculate_derivative_at_nodes(dndy1_at_node, dndy2_at_node);
VerdictVector normal_at_nodes[9];
diff --git a/V_TetMetric.cpp b/V_TetMetric.cpp
index db3b7c3..b227006 100644
--- a/V_TetMetric.cpp
+++ b/V_TetMetric.cpp
@@ -1360,10 +1360,10 @@ double tet_distortion(int num_nodes, const double coordinates[][3])
double weight[maxTotalNumberGaussPoints];
// create an object of GaussIntegration for tet
- GaussIntegration gint{};
- gint.initialize(number_of_gauss_points, num_nodes, number_dims, is_tri);
- gint.calculate_shape_function_3d_tet();
- gint.get_shape_func(shape_function[0], dndy1[0], dndy2[0], dndy3[0], weight);
+ GaussIntegration gaussint;
+ gaussint.initialize(number_of_gauss_points, num_nodes, number_dims, is_tri);
+ gaussint.calculate_shape_function_3d_tet();
+ gaussint.get_shape_func(shape_function[0], dndy1[0], dndy2[0], dndy3[0], weight);
// vector xxi is the derivative vector of coordinates w.r.t local xi coordinate in the
// computation space
@@ -1408,7 +1408,7 @@ double tet_distortion(int num_nodes, const double coordinates[][3])
double dndy2_at_node[maxNumberNodes][maxNumberNodes];
double dndy3_at_node[maxNumberNodes][maxNumberNodes];
- gint.calculate_derivative_at_nodes_3d_tet(dndy1_at_node, dndy2_at_node, dndy3_at_node);
+ gaussint.calculate_derivative_at_nodes_3d_tet(dndy1_at_node, dndy2_at_node, dndy3_at_node);
int node_id;
for (node_id = 0; node_id < num_nodes; node_id++)
{
diff --git a/V_TriMetric.cpp b/V_TriMetric.cpp
index 3a2a16a..a2caba5 100644
--- a/V_TriMetric.cpp
+++ b/V_TriMetric.cpp
@@ -778,10 +778,10 @@ double tri_distortion(int num_nodes, const double coordinates[][3])
// create an object of GaussIntegration
int number_dims = 2;
int is_tri = 1;
- GaussIntegration gint{};
- gint.initialize(number_of_gauss_points, num_nodes, number_dims, is_tri);
- gint.calculate_shape_function_2d_tri();
- gint.get_shape_func(shape_function[0], dndy1[0], dndy2[0], weight);
+ GaussIntegration gaussint;
+ gaussint.initialize(number_of_gauss_points, num_nodes, number_dims, is_tri);
+ gaussint.calculate_shape_function_2d_tri();
+ gaussint.get_shape_func(shape_function[0], dndy1[0], dndy2[0], weight);
// calculate element area
int ife, ja;
@@ -805,7 +805,7 @@ double tri_distortion(int num_nodes, const double coordinates[][3])
double dndy1_at_node[maxNumberNodes][maxNumberNodes];
double dndy2_at_node[maxNumberNodes][maxNumberNodes];
- gint.calculate_derivative_at_nodes_2d_tri(dndy1_at_node, dndy2_at_node);
+ gaussint.calculate_derivative_at_nodes_2d_tri(dndy1_at_node, dndy2_at_node);
VerdictVector normal_at_nodes[7];
|
