From e74b60314237e6bce790e6849dffe88d406a4658 Mon Sep 17 00:00:00 2001
From: Anthony Royer <anthony.royer@uliege.be>
Date: Mon, 16 May 2022 10:21:31 +0200
Subject: [PATCH] New post-pro
---
examples/helmholtz/crossPoints/ddm2D.cpp | 8 ++++++++
1 file changed, 8 insertions(+)
diff --git a/examples/helmholtz/crossPoints/ddm2D.cpp b/examples/helmholtz/crossPoints/ddm2D.cpp
index 1e984dea..ebaf8f8c 100755
--- a/examples/helmholtz/crossPoints/ddm2D.cpp
+++ b/examples/helmholtz/crossPoints/ddm2D.cpp
@@ -136,6 +136,8 @@ namespace D2 {
gmshDdm->userDefinedParameter(scalingPlot, "scaling");
bool saveMesh = false;
gmshDdm->userDefinedParameter(saveMesh, "saveMesh");
+ bool saveIteration = false;
+ gmshDdm->userDefinedParameter(saveIteration, "saveIteration");
if(benchmark == "scattering") {
checkerboard(nDomX, nDomY, sizeX, sizeY, R, lc, (boundary == "pml"), pmlSize, (boundaryExt == "pml"), pmlSizeExt, meshOrder, 0, 0);
@@ -788,6 +790,12 @@ namespace D2 {
gmshfem::common::Timer pre = formulation.pre();
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 j = 0; j < nDomY; ++j) {
unsigned int index = i * nDomY + j;
--
GitLab