From 9510df634ba8b569ddac8fa6d3229e9ad757e26d Mon Sep 17 00:00:00 2001 From: Rainer Backofen <rainer.backofen@tu-dresden.de> Date: Mon, 5 Jul 2010 07:46:00 +0000 Subject: [PATCH] changed heat demo to have some working example --- demo/init/heat.dat.2d | 18 ++++++++++-------- demo/src/heat.cc | 2 +- 2 files changed, 11 insertions(+), 9 deletions(-) diff --git a/demo/init/heat.dat.2d b/demo/init/heat.dat.2d index b0affbd9..6cd362ed 100644 --- a/demo/init/heat.dat.2d +++ b/demo/init/heat.dat.2d @@ -15,16 +15,17 @@ heat->space->solver->left precon: no heat->space->solver->right precon: no heat->space->estimator: residual -heat->space->estimator->C0: 1.0 -heat->space->estimator->C1: 1.0 +heat->space->estimator->C0: 0.1 +heat->space->estimator->C1: 0.1 heat->theta: 1.0 -heat->adapt->timestep: 0.1 +heat->adapt->timestep: 0.01 +heat->adapt->min timestep: 0.01 heat->adapt->start time: 0.0 -heat->adapt->end time: 1.0 +heat->adapt->end time: 10.0 -heat->adapt->tolerance: 0.1 +heat->adapt->tolerance: 0.0001 heat->adapt->rel space error: 0.5 heat->adapt->rel time error: 0.5 heat->adapt->time theta 1: 1.0 @@ -36,13 +37,14 @@ heat->adapt->time delta 2: 1.4142 heat->adapt->max iteration: 100 heat->adapt->coarsen allowed: 1 % 0|1 -heat->adapt->refine bisections: 2 -heat->adapt->coarsen bisections: 2 +heat->adapt->refine bisections: 1 +heat->adapt->coarsen bisections: 1 -heat->space->marker->strategy: 2 % 0=none, 1=GR, 2=MS, 3=ES, 4=GERS +heat->space->marker->strategy: 4 % 0=none, 1=GR, 2=MS, 3=ES, 4=GERS heat->space->marker->MSGamma: 0.5 heat->space->marker->MSGammaC: 0.1 +heat->space->output->ParaView animation: 1 heat->space->output->filename: output/heat heat->space->output->ParaView format: 1 heat->space->output->write every i-th timestep: 1 diff --git a/demo/src/heat.cc b/demo/src/heat.cc index a3a42d05..80a60ecb 100644 --- a/demo/src/heat.cc +++ b/demo/src/heat.cc @@ -230,7 +230,7 @@ int main(int argc, char** argv) boundaryFct->setTimePtr(heat.getBoundaryTimePtr()); heat.setExactSolution(boundaryFct); - heatSpace.addDirichletBC(DIRICHLET, boundaryFct); + // heatSpace.addDirichletBC(DIRICHLET, boundaryFct); // ===== start adaption loop ===== -- GitLab