Skip to content
Snippets Groups Projects
Commit 0e6db92c authored by Christian Engwer's avatar Christian Engwer
Browse files

pmgsolver had duplicate dump-names

[[Imported from SVN: r82]]
parent 1bcd40b3
No related branches found
No related tags found
No related merge requests found
......@@ -386,8 +386,6 @@ namespace Dune {
int cycle=0;
if (rank==0) std::cout << "MGC-Cycle " << cycle
<< " " << mydefect << " " << 0 << std::endl;
cycle ++;
double lastdefect=mydefect;
TIME_SMOOTHER = 0;
......@@ -396,10 +394,16 @@ namespace Dune {
TIME_DEFECT = 0;
#ifdef SOLVER_DUMPDX
dumpdx(g,g.smoothest(),d,"defect","pressure in the end");
dumpdx(g,g.smoothest(),x,"calpres","pressure in the end");
char dname[156];
char cname[256];
sprintf(dname,"defect-t%04i",cycle);
sprintf(cname,"calpres-t%04i",cycle);
dumpdx(g,g.smoothest(),d,dname,"pressure in the end");
dumpdx(g,g.smoothest(),x,cname,"pressure in the end");
#endif
cycle ++;
while (mydefect > maxdefect)
{
if (cycle > maxCycles) break;
......@@ -430,13 +434,18 @@ namespace Dune {
<< " " << mydefect/lastdefect
<< std::endl;
lastdefect = mydefect;
cycle ++;
#ifdef SOLVER_DUMPDX
dumpdx(g,g.smoothest(),d,"defect","pressure in the end");
dumpdx(g,g.smoothest(),x,"calpres","pressure in the end");
char dname[156];
char cname[256];
sprintf(dname,"defect-t%04i",cycle);
sprintf(cname,"calpres-t%04i",cycle);
dumpdx(g,g.smoothest(),d,dname,"pressure in the end");
dumpdx(g,g.smoothest(),x,cname,"pressure in the end");
#endif
cycle ++;
};
if (rank==0)
std::cout << "Time in smoother:" << TIME_SMOOTHER << std::endl
......
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