From 3fe0ad6eec242e98d3365e44522190495543b76e Mon Sep 17 00:00:00 2001
From: Simon Praetorius <simon.praetorius@tu-dresden.de>
Date: Tue, 4 Sep 2018 17:09:07 +0200
Subject: [PATCH] polygongrid writer with legacy vtk function
---
src/polygongrid.cc | 15 ++++++++++-----
1 file changed, 10 insertions(+), 5 deletions(-)
diff --git a/src/polygongrid.cc b/src/polygongrid.cc
index d3ad10d..d497706 100644
--- a/src/polygongrid.cc
+++ b/src/polygongrid.cc
@@ -58,12 +58,17 @@ int main(int argc, char** argv)
using Writer = VtkUnstructuredGridWriter<GridView>;
Writer vtkWriter(gridView, Vtk::ASCII);
- auto p1Analytic = makeAnalyticGridViewFunction([](auto const& x) {
- return std::sin(10*x[0])*std::cos(10*x[1]);
- }, gridView);
- vtkWriter.addPointData(p1Analytic, "p1");
- vtkWriter.addCellData(p1Analytic, "p0");
+ std::vector<double> p1function(gridView.size(dim), 1.0);
+ using P1Function = P1VTKFunction<GridView,std::vector<double>>;
+ std::shared_ptr<VTKFunction<GridView> const> p1FctWrapped(new P1Function(gridView, p1function, "p1"));
+
+ std::vector<double> p0function(gridView.size(0), 1.0);
+ using P0Function = P0VTKFunction<GridView,std::vector<double>>;
+ std::shared_ptr<VTKFunction<GridView> const> p0FctWrapped(new P0Function(gridView, p0function, "p0"));
+
+ vtkWriter.addPointData(p1FctWrapped, "p1");
+ vtkWriter.addCellData(p0FctWrapped, "p0");
vtkWriter.write("poly_ascii_float32.vtu");
}
--
GitLab