173 |
} |
} |
174 |
|
|
175 |
for (size_t i=0; i<mesh.numberOfVertices(); ++i) { |
for (size_t i=0; i<mesh.numberOfVertices(); ++i) { |
176 |
const TinyVector<3>& X = mesh.vertex(i).Vector3(); |
const TinyVector<3>& X = mesh.vertex(i); |
177 |
fout << (*__function).value(X[0], X[1], X[2]) << CR; |
fout << (*__function).value(X[0], X[1], X[2]) << CR; |
178 |
} |
} |
179 |
break; |
break; |
200 |
default: { |
default: { |
201 |
FunctionExpression& f = *(*__function).value(); |
FunctionExpression& f = *(*__function).value(); |
202 |
for (size_t i=0; i<mesh.numberOfVertices(); ++i) { |
for (size_t i=0; i<mesh.numberOfVertices(); ++i) { |
203 |
const TinyVector<3>& X = mesh.vertex(i).Vector3(); |
const TinyVector<3>& X = mesh.vertex(i); |
204 |
fsol << f.value(X[0], X[1], X[2]) << CR; |
fsol << f.value(X[0], X[1], X[2]) << CR; |
205 |
} |
} |
206 |
} |
} |
264 |
// } |
// } |
265 |
|
|
266 |
// for (size_t i=0; i<mesh.numberOfVertices(); ++i) { |
// for (size_t i=0; i<mesh.numberOfVertices(); ++i) { |
267 |
// const TinyVector<3>& X = mesh.vertex(i).Vector3(); |
// const TinyVector<3>& X = mesh.vertex(i); |
268 |
// fout << (*__function).value(X[0], X[1], X[2]) << CR; |
// fout << (*__function).value(X[0], X[1], X[2]) << CR; |
269 |
// } |
// } |
270 |
// break; |
// break; |
277 |
fsol << "3 3 " << mesh.numberOfVertices() << " 2" << CR; |
fsol << "3 3 " << mesh.numberOfVertices() << " 2" << CR; |
278 |
|
|
279 |
for (size_t i=0; i<mesh.numberOfVertices(); ++i) { |
for (size_t i=0; i<mesh.numberOfVertices(); ++i) { |
280 |
const TinyVector<3>& X = mesh.vertex(i).Vector3(); |
const TinyVector<3>& X = mesh.vertex(i); |
281 |
fsol << (*__field).f1().value(X[0], X[1], X[2]) << ' ' |
fsol << (*__field).f1().value(X[0], X[1], X[2]) << ' ' |
282 |
<< (*__field).f2().value(X[0], X[1], X[2]) << ' ' |
<< (*__field).f2().value(X[0], X[1], X[2]) << ' ' |
283 |
<< (*__field).f3().value(X[0], X[1], X[2]) << CR; |
<< (*__field).f3().value(X[0], X[1], X[2]) << CR; |