34 void Foam::block::blockPoints()
37 const label ni = blockDef_.
n().
x();
38 const label nj = blockDef_.
n().
y();
39 const label nk = blockDef_.
n().
z();
57 for (label
k = 0;
k <= nk;
k++)
59 for (label j = 0; j <= nj; j++)
61 for (label i = 0; i <= ni; i++)
63 label vertexNo = vtxLabel(i, j,
k);
66 vector edgex1 = p000 + (p100 - p000)*w[0][i];
67 vector edgex2 = p010 + (p110 - p010)*w[1][i];
68 vector edgex3 = p011 + (p111 - p011)*w[2][i];
69 vector edgex4 = p001 + (p101 - p001)*w[3][i];
71 vector edgey1 = p000 + (p010 - p000)*w[4][j];
72 vector edgey2 = p100 + (p110 - p100)*w[5][j];
73 vector edgey3 = p101 + (p111 - p101)*w[6][j];
74 vector edgey4 = p001 + (p011 - p001)*w[7][j];
76 vector edgez1 = p000 + (p001 - p000)*w[8][
k];
77 vector edgez2 = p100 + (p101 - p100)*w[9][
k];
78 vector edgez3 = p110 + (p111 - p110)*w[10][
k];
79 vector edgez4 = p010 + (p011 - p010)*w[11][
k];
85 (1.0 - w[0][i])*(1.0 - w[4][j])*(1.0 - w[8][
k])
86 + w[0][i]*(1.0 - w[5][j])*(1.0 - w[9][
k])
91 (1.0 - w[1][i])*w[4][j]*(1.0 - w[11][
k])
92 + w[1][i]*w[5][j]*(1.0 - w[10][
k])
97 (1.0 - w[2][i])*w[7][j]*w[11][k]
98 + w[2][i]*w[6][j]*w[10][k]
104 (1.0 - w[3][i])*(1.0 - w[7][j])*w[8][
k]
105 + w[3][i]*(1.0 - w[6][j])*w[9][k]
108 scalar magImpx = impx1 + impx2 + impx3 + impx4;
118 (1.0 - w[4][j])*(1.0 - w[0][i])*(1.0 - w[8][
k])
119 + w[4][j]*(1.0 - w[1][i])*(1.0 - w[11][
k])
124 (1.0 - w[5][j])*w[0][i]*(1.0 - w[9][k])
125 + w[5][j]*w[1][i]*(1.0 - w[10][
k])
130 (1.0 - w[6][j])*w[3][i]*w[9][k]
131 + w[6][j]*w[2][i]*w[10][k]
136 (1.0 - w[7][j])*(1.0 - w[3][i])*w[8][
k]
137 + w[7][j]*(1.0 - w[2][i])*w[11][k]
140 scalar magImpy = impy1 + impy2 + impy3 + impy4;
150 (1.0 - w[8][
k])*(1.0 - w[0][i])*(1.0 - w[4][j])
151 + w[8][k]*(1.0 - w[3][i])*(1.0 - w[7][j])
156 (1.0 - w[9][
k])*w[0][i]*(1.0 - w[5][j])
157 + w[9][
k]*w[3][i]*(1.0 - w[6][j])
162 (1.0 - w[10][
k])*w[1][i]*w[5][j]
163 + w[10][k]*w[2][i]*w[6][j]
168 (1.0 - w[11][
k])*(1.0 - w[1][i])*w[4][j]
169 + w[11][
k]*(1.0 - w[2][i])*w[7][j]
172 scalar magImpz = impz1 + impz2 + impz3 + impz4;
179 vector corx1 = impx1*(p[0][i] - edgex1);
180 vector corx2 = impx2*(p[1][i] - edgex2);
181 vector corx3 = impx3*(p[2][i] - edgex3);
182 vector corx4 = impx4*(p[3][i] - edgex4);
184 vector cory1 = impy1*(p[4][j] - edgey1);
185 vector cory2 = impy2*(p[5][j] - edgey2);
186 vector cory3 = impy3*(p[6][j] - edgey3);
187 vector cory4 = impy4*(p[7][j] - edgey4);
189 vector corz1 = impz1*(p[8][
k] - edgez1);
190 vector corz2 = impz2*(p[9][
k] - edgez2);
191 vector corz3 = impz3*(p[10][
k] - edgez3);
192 vector corz4 = impz4*(p[11][
k] - edgez4);
217 vertices_[vertexNo] = edgex1 + edgex2 + edgex3 + edgex4;
218 vertices_[vertexNo] += edgey1 + edgey2 + edgey3 + edgey4;
219 vertices_[vertexNo] += edgez1 + edgez2 + edgez3 + edgez4;
221 vertices_[vertexNo] /= 3.0;
223 vertices_[vertexNo] += corx1 + corx2 + corx3 + corx4;
224 vertices_[vertexNo] += cory1 + cory2 + cory3 + cory4;
225 vertices_[vertexNo] += corz1 + corz2 + corz3 + corz4;