Skip to content
Snippets Groups Projects
Commit 1c213147 authored by Oliver Sander's avatar Oliver Sander Committed by sander
Browse files

Minor cleanup

[[Imported from SVN: r10108]]
parent 0856b5ab
No related branches found
No related tags found
No related merge requests found
...@@ -83,8 +83,6 @@ int main (int argc, char *argv[]) try ...@@ -83,8 +83,6 @@ int main (int argc, char *argv[]) try
// parse data file // parse data file
ParameterTree parameterSet; ParameterTree parameterSet;
// if (argc != 2)
// DUNE_THROW(Exception, "Usage: ./hencky-material <parameter file>");
ParameterTreeParser::readINITree(argv[1], parameterSet); ParameterTreeParser::readINITree(argv[1], parameterSet);
...@@ -138,7 +136,6 @@ int main (int argc, char *argv[]) try ...@@ -138,7 +136,6 @@ int main (int argc, char *argv[]) try
typedef GridType::LeafGridView GridView; typedef GridType::LeafGridView GridView;
GridView gridView = grid->leafGridView(); GridView gridView = grid->leafGridView();
// typedef P1NodalBasis<GridView,double> FEBasis;
typedef P1NodalBasis<GridView,double> FEBasis; typedef P1NodalBasis<GridView,double> FEBasis;
FEBasis feBasis(gridView); FEBasis feBasis(gridView);
...@@ -219,13 +216,10 @@ int main (int argc, char *argv[]) try ...@@ -219,13 +216,10 @@ int main (int argc, char *argv[]) try
{ {
size_t idx = grid->leafGridView().indexSet().index(*it); size_t idx = grid->leafGridView().indexSet().index(*it);
displacement[idx] = x[idx] - it->geometry().corner(0); displacement[idx] = x[idx] - it->geometry().corner(0);
//std::cout << "idx: " << idx << " coordinate: " << it->geometry().corner(0) << std::endl;
} }
Dune::shared_ptr<VTKBasisGridFunction<FEBasis,BlockVector<FieldVector<double,3> > > > vtkDisplacement auto vtkDisplacement = Dune::make_shared<VTKBasisGridFunction<FEBasis,BlockVector<FieldVector<double,3> > > >
= Dune::make_shared<VTKBasisGridFunction<FEBasis,BlockVector<FieldVector<double,3> > > > (feBasis, displacement, "Displacement");
(feBasis, displacement, "Displacement");
vtkWriter.addVertexData(vtkDisplacement); vtkWriter.addVertexData(vtkDisplacement);
vtkWriter.write(resultPath + "hencky_homotopy_0"); vtkWriter.write(resultPath + "hencky_homotopy_0");
...@@ -330,9 +324,8 @@ int main (int argc, char *argv[]) try ...@@ -330,9 +324,8 @@ int main (int argc, char *argv[]) try
// Output result // Output result
///////////////////////////////// /////////////////////////////////
Dune::shared_ptr<VTKBasisGridFunction<FEBasis,BlockVector<FieldVector<double,3> > > > vtkDisplacement auto vtkDisplacement = Dune::make_shared<VTKBasisGridFunction<FEBasis,BlockVector<FieldVector<double,3> > > >
= Dune::make_shared<VTKBasisGridFunction<FEBasis,BlockVector<FieldVector<double,3> > > > (feBasis, displacement, "Displacement");
(feBasis, displacement, "Displacement");
vtkWriter.addVertexData(vtkDisplacement); vtkWriter.addVertexData(vtkDisplacement);
vtkWriter.write(resultPath + "hencky_homotopy_" + std::to_string(i+1)); vtkWriter.write(resultPath + "hencky_homotopy_" + std::to_string(i+1));
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment