From 687f534bc1599900d6a5e6731ad415da9a407c83 Mon Sep 17 00:00:00 2001
From: Christophe Geuzaine <cgeuzaine@uliege.be>
Date: Sat, 18 Mar 2023 14:57:49 +0100
Subject: [PATCH] validation on marge_request would require keeping the gmshfem
dependency - removing until this is handled correctly
---
.gitlab-ci.yml | 1 -
examples/helmholtzflow/nacelle3D/nacelleFunctions.cpp | 2 +-
tutorials/helmholtz/waveguide/main.cpp | 10 ++++++++++
3 files changed, 11 insertions(+), 2 deletions(-)
diff --git a/.gitlab-ci.yml b/.gitlab-ci.yml
index 801e3440..695b12ac 100644
--- a/.gitlab-ci.yml
+++ b/.gitlab-ci.yml
@@ -81,7 +81,6 @@ linux64_validations:
- docker
only:
- master
- - merge_requests
#linux64_coverage:
# stage: test
diff --git a/examples/helmholtzflow/nacelle3D/nacelleFunctions.cpp b/examples/helmholtzflow/nacelle3D/nacelleFunctions.cpp
index 24799f05..3e196b96 100644
--- a/examples/helmholtzflow/nacelle3D/nacelleFunctions.cpp
+++ b/examples/helmholtzflow/nacelle3D/nacelleFunctions.cpp
@@ -25,7 +25,7 @@ using namespace gmshfem::equation;
void Nacelle(const double h_near, const double h_far, const double D, const double R, const double L_duct, const double L_pml, const double L_pml_a, const int nDom)
{
gmsh::model::add("Nacelle");
-
+ gmsh::option::setNumber("Geometry.OCCUnionUnify", 0);
double RefCoef = 2.; // refinement coefficient
// Nacelle (points)
gmsh::model::occ::addPoint(-0.3000000000000000, 1.5000000000000000, 0, h_near, 1);
diff --git a/tutorials/helmholtz/waveguide/main.cpp b/tutorials/helmholtz/waveguide/main.cpp
index b96b50ce..837e29d5 100644
--- a/tutorials/helmholtz/waveguide/main.cpp
+++ b/tutorials/helmholtz/waveguide/main.cpp
@@ -81,6 +81,7 @@ void waveguide(const unsigned int nDom, const double l, const double L, const do
gmsh::model::geo::synchronize();
gmsh::model::mesh::generate();
+ //gmsh::write("test.msh");
}
int main(int argc, char **argv)
@@ -207,5 +208,14 @@ int main(int argc, char **argv)
}
}
+ bool save_matrix = false;
+ gmshDdm.userDefinedParameter(save_matrix, "save_iteration_matrix");
+ if (save_matrix) {
+ auto mat = formulation.computeMatrix();
+ msg::info << "nnz = " << mat.numberOfNonZero() << ", symmetric = " << mat.isSymmetric() << ", hermitian = " << mat.isHermitian() << msg::endl;
+ mat.save("Mymatrix");
+ mat.saveSpyPlot("Mymatrix");
+ }
+
return 0;
}
--
GitLab