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

Parset-control printing of frictional boundary

parent f5eb05e4
Branches
No related tags found
No related merge requests found
......@@ -435,15 +435,17 @@ int main(int argc, char *argv[]) {
std::cout << "sup |u1 - u4| = " << diff4.infinity_norm() << ", "
<< "|u1 - u4| = " << diff4.two_norm() << std::endl;
// Print displacement on frictional boundary
boost::format const formatter("u1[%02d] = %+3e, "
"%|40t|u2[%02d] = %+3e, "
"%|80t|u3[%02d] = %+3e, "
"%|120t|u4[%02d] = %+3e");
for (size_t i = 0; i < frictionalNodes.size(); ++i)
if (frictionalNodes[i][0])
std::cout << boost::format(formatter) % i % u1[i] % i % u2[i] % i %
u3[i] % i % u4[i] << std::endl;
if (parset.get<bool>("printFrictionalBoundary")) {
// Print displacement on frictional boundary
boost::format const formatter("u1[%02d] = %+3e, "
"%|40t|u2[%02d] = %+3e, "
"%|80t|u3[%02d] = %+3e, "
"%|120t|u4[%02d] = %+3e");
for (size_t i = 0; i < frictionalNodes.size(); ++i)
if (frictionalNodes[i][0])
std::cout << boost::format(formatter) % i % u1[i] % i % u2[i] % i %
u3[i] % i % u4[i] << std::endl;
}
}
catch (Dune::Exception &e) {
Dune::derr << "Dune reported error: " << e << std::endl;
......
......@@ -2,6 +2,7 @@
timesteps = 100
verbose = false
printFrictionalBoundary = false
useNonlinearGS = false # true
useGS = false
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment