Skip to content
Snippets Groups Projects
Commit e666fa32 authored by Elias Pipping's avatar Elias Pipping Committed by Elias Pipping
Browse files

Revert "Externalise some more"

This reverts commit 3c2fc4080e8af5ed4b17512914e6a7db866a11e6.

Conflicts:

	src/one-body-sample.org
parent 0e817230
No related branches found
No related tags found
No related merge requests found
...@@ -165,26 +165,6 @@ ...@@ -165,26 +165,6 @@
if (parset.get<bool>("printProgress")) if (parset.get<bool>("printProgress"))
std::cerr << std::endl; std::cerr << std::endl;
#+end_src #+end_src
#+name: setupAndSolveProblem
#+begin_src c++
{
auto myGlobalNonlinearity
= assemble_nonlinearity<MatrixType, VectorType>
(parset.sub("boundary.friction"), *nodalIntegrals, alpha,
surfaceNormalStress);
MyConvexProblemType const myConvexProblem
(problem_A, *myGlobalNonlinearity, problem_rhs);
MyBlockProblemType myBlockProblem(parset, myConvexProblem);
auto multigridStep = mySolver.getSolver();
multigridStep->setProblem(problem_iterate, myBlockProblem);
LoopSolver<VectorType> overallSolver
(multigridStep, parset.get<size_t>("solver.tnnmg.maxiterations"),
solver_tolerance, &energyNorm, verbosity, false); // absolute error
overallSolver.solve();
}
#+end_src
* Functions * Functions
#+name: defineInitTimeStepper #+name: defineInitTimeStepper
#+begin_src c++ #+begin_src c++
...@@ -414,7 +394,23 @@ ...@@ -414,7 +394,23 @@
// displacement as a convergence criterion // displacement as a convergence criterion
VectorType u_saved; VectorType u_saved;
for (size_t state_fpi = 1; state_fpi <= state_fpi_max; ++state_fpi) { for (size_t state_fpi = 1; state_fpi <= state_fpi_max; ++state_fpi) {
<<setupAndSolveProblem>>; {
auto myGlobalNonlinearity
= assemble_nonlinearity<MatrixType, VectorType>
(parset.sub("boundary.friction"), *nodalIntegrals, alpha,
surfaceNormalStress);
MyConvexProblemType const myConvexProblem
(problem_A, *myGlobalNonlinearity, problem_rhs);
MyBlockProblemType myBlockProblem(parset, myConvexProblem);
auto multigridStep = mySolver.getSolver();
multigridStep->setProblem(problem_iterate, myBlockProblem);
LoopSolver<VectorType> overallSolver
(multigridStep, parset.get<size_t>("solver.tnnmg.maxiterations"),
solver_tolerance, &energyNorm, verbosity, false); // absolute error
overallSolver.solve();
}
timeSteppingScheme->postProcess(problem_iterate); timeSteppingScheme->postProcess(problem_iterate);
timeSteppingScheme->extractDisplacement(u); timeSteppingScheme->extractDisplacement(u);
......
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