Skip to content
Snippets Groups Projects
Commit c1e8f31d authored by oliver.sander_at_tu-dresden.de's avatar oliver.sander_at_tu-dresden.de
Browse files

Move linelast.cc into the 'src' directory

This is where such files should reside nowadays.
parent 30dd39a4
No related branches found
No related tags found
No related merge requests found
...@@ -16,7 +16,7 @@ add_subdirectory("dune") ...@@ -16,7 +16,7 @@ add_subdirectory("dune")
add_subdirectory("src") add_subdirectory("src")
if(HAVE_AMIRAMESH AND HAVE_IPOPT AND HAVE_PSURFACE AND HAVE_UG) if(HAVE_AMIRAMESH AND HAVE_IPOPT AND HAVE_PSURFACE AND HAVE_UG)
set(programs linelast viscoelast nonlinelast) set(programs viscoelast nonlinelast)
foreach(_program ${programs}) foreach(_program ${programs})
add_executable(${_program} ${_program}.cc) add_executable(${_program} ${_program}.cc)
add_dune_amiramesh_flags(${_program}) add_dune_amiramesh_flags(${_program})
......
if(ADOLC_FOUND AND IPOPT_FOUND AND PYTHONLIBS_FOUND AND dune-uggrid_FOUND) if(ADOLC_FOUND AND IPOPT_FOUND AND PYTHONLIBS_FOUND AND dune-uggrid_FOUND)
set(programs finite-strain-elasticity) set(programs finite-strain-elasticity
linear-elasticity)
foreach(_program ${programs}) foreach(_program ${programs})
add_executable(${_program} ${_program}.cc) add_executable(${_program} ${_program}.cc)
add_dune_adolc_flags(${_program}) add_dune_adolc_flags(${_program})
......
...@@ -89,24 +89,36 @@ int main (int argc, char *argv[]) try ...@@ -89,24 +89,36 @@ int main (int argc, char *argv[]) try
GridType* grid= new GridType; GridType* grid= new GridType;
grid->setRefinementType(GridType::COPY); grid->setRefinementType(GridType::COPY);
#if HAVE_AMIRAMESH
if (paramBoundaries) if (paramBoundaries)
grid = AmiraMeshReader<GridType>::read(path + gridFile, PSurfaceBoundary<dim-1>::read(path + parFile)); grid = AmiraMeshReader<GridType>::read(path + gridFile, PSurfaceBoundary<dim-1>::read(path + parFile));
else else
grid = AmiraMeshReader<GridType>::read(path + gridFile); grid = AmiraMeshReader<GridType>::read(path + gridFile);
#else
#warning You need libamiramesh for this code!
#endif
LevelBoundaryPatch coarseDirichletBoundary; LevelBoundaryPatch coarseDirichletBoundary;
coarseDirichletBoundary.setup(grid->levelGridView(0)); coarseDirichletBoundary.setup(grid->levelGridView(0));
readBoundaryPatch<GridType>(coarseDirichletBoundary, path + dirichletFile); readBoundaryPatch<GridType>(coarseDirichletBoundary, path + dirichletFile);
VectorType coarseDirichletValues(grid->size(0, dim)); VectorType coarseDirichletValues(grid->size(0, dim));
#if HAVE_AMIRAMESH
AmiraMeshReader<GridType>::readFunction(coarseDirichletValues, path + dirichletValuesFile); AmiraMeshReader<GridType>::readFunction(coarseDirichletValues, path + dirichletValuesFile);
#else
#warning You need libamiramesh for this code!
#endif
LevelBoundaryPatch coarseNeumannBoundary; LevelBoundaryPatch coarseNeumannBoundary;
coarseNeumannBoundary.setup(grid->levelGridView(0)); coarseNeumannBoundary.setup(grid->levelGridView(0));
readBoundaryPatch<GridType>(coarseNeumannBoundary, path + neumannFile); readBoundaryPatch<GridType>(coarseNeumannBoundary, path + neumannFile);
VectorType coarseNeumannValues(grid->size(0, dim)); VectorType coarseNeumannValues(grid->size(0, dim));
#if HAVE_AMIRAMESH
AmiraMeshReader<GridType>::readFunction(coarseNeumannValues, path + neumannValuesFile); AmiraMeshReader<GridType>::readFunction(coarseNeumannValues, path + neumannValuesFile);
#else
#warning You need libamiramesh for this code!
#endif
P1NodalBasis<GridType::LevelGridView> coarseBasis(grid->levelGridView(0)); P1NodalBasis<GridType::LevelGridView> coarseBasis(grid->levelGridView(0));
auto coarseDir = ::Functions::makeFunction(coarseBasis,coarseDirichletValues); auto coarseDir = ::Functions::makeFunction(coarseBasis,coarseDirichletValues);
...@@ -284,10 +296,14 @@ int main (int argc, char *argv[]) try ...@@ -284,10 +296,14 @@ int main (int argc, char *argv[]) try
} }
// Output result // Output result
#if HAVE_AMIRAMESH
LeafAmiraMeshWriter<GridType> amiramesh; LeafAmiraMeshWriter<GridType> amiramesh;
amiramesh.addLeafGrid(*grid,true); amiramesh.addLeafGrid(*grid,true);
amiramesh.addVertexData(x, grid->leafGridView(),true); amiramesh.addVertexData(x, grid->leafGridView(),true);
amiramesh.write("resultGrid",1); amiramesh.write("resultGrid",1);
#else
#warning You need libamiramesh for this code!
#endif
} catch (Exception e) { } catch (Exception e) {
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment