70 |
real_t wk0 = static_cast<real_t>(k) / static_cast<real_t>(nz_1); |
real_t wk0 = static_cast<real_t>(k) / static_cast<real_t>(nz_1); |
71 |
real_t wk1 = static_cast<real_t>(nz_1-k)/static_cast<real_t>(nz_1); |
real_t wk1 = static_cast<real_t>(nz_1-k)/static_cast<real_t>(nz_1); |
72 |
Vertex& v = vertex(i,j,k); |
Vertex& v = vertex(i,j,k); |
73 |
v.X() = wi0 * __verticesShape.b(0) + wi1 * __verticesShape.a(0); |
v.x() = wi0 * __verticesShape.b(0) + wi1 * __verticesShape.a(0); |
74 |
v.Y() = wj0 * __verticesShape.b(1) + wj1 * __verticesShape.a(1); |
v.y() = wj0 * __verticesShape.b(1) + wj1 * __verticesShape.a(1); |
75 |
v.Z() = wk0 * __verticesShape.b(2) + wk1 * __verticesShape.a(2); |
v.z() = wk0 * __verticesShape.b(2) + wk1 * __verticesShape.a(2); |
76 |
} |
} |
77 |
} |
} |
78 |
} |
} |
222 |
for (size_t i=0; i<this->numberOfVertices(); ++i) { |
for (size_t i=0; i<this->numberOfVertices(); ++i) { |
223 |
Vertex& V = vertex(i); |
Vertex& V = vertex(i); |
224 |
|
|
225 |
V.Ref() = S.inside(V); |
V.reference() = S.inside(V); |
226 |
} |
} |
227 |
} |
} |
228 |
|
|