Skip to content
Snippets Groups Projects
Commit e74b6031 authored by Anthony Royer's avatar Anthony Royer
Browse files

New post-pro

parent 6f659a83
Branches
Tags
No related merge requests found
Pipeline #9643 passed
...@@ -136,6 +136,8 @@ namespace D2 { ...@@ -136,6 +136,8 @@ namespace D2 {
gmshDdm->userDefinedParameter(scalingPlot, "scaling"); gmshDdm->userDefinedParameter(scalingPlot, "scaling");
bool saveMesh = false; bool saveMesh = false;
gmshDdm->userDefinedParameter(saveMesh, "saveMesh"); gmshDdm->userDefinedParameter(saveMesh, "saveMesh");
bool saveIteration = false;
gmshDdm->userDefinedParameter(saveIteration, "saveIteration");
if(benchmark == "scattering") { if(benchmark == "scattering") {
checkerboard(nDomX, nDomY, sizeX, sizeY, R, lc, (boundary == "pml"), pmlSize, (boundaryExt == "pml"), pmlSizeExt, meshOrder, 0, 0); checkerboard(nDomX, nDomY, sizeX, sizeY, R, lc, (boundary == "pml"), pmlSize, (boundaryExt == "pml"), pmlSizeExt, meshOrder, 0, 0);
...@@ -788,6 +790,12 @@ namespace D2 { ...@@ -788,6 +790,12 @@ namespace D2 {
gmshfem::common::Timer pre = formulation.pre(); gmshfem::common::Timer pre = formulation.pre();
gmshfem::common::Timer solve = formulation.solve("gmres", res, iterMax); gmshfem::common::Timer solve = formulation.solve("gmres", res, iterMax);
if(saveIteration) {
gmshfem::common::CSVio file(fileName, ';', gmshfem::common::OpeningMode::Append);
file << N << thetaPade << formulation.numberOfIterations() << gmshfem::csv::endl;
file.close();
}
for(unsigned int i = 0; i < nDomX; ++i) { for(unsigned int i = 0; i < nDomX; ++i) {
for(unsigned int j = 0; j < nDomY; ++j) { for(unsigned int j = 0; j < nDomY; ++j) {
unsigned int index = i * nDomY + j; unsigned int index = i * nDomY + j;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment