This is an automated email from the git hooks/post-receive script. logari81 pushed a commit to branch master in repository getfem.
The following commit(s) were added to refs/heads/master by this push: new 2a1dd105 Relax convergence threshold for inverse geometric transformation 2a1dd105 is described below commit 2a1dd1059ed54a3dfceb5662f7e542132a774c99 Author: Konstantinos Poulios <logar...@gmail.com> AuthorDate: Thu Sep 29 14:50:01 2022 +0200 Relax convergence threshold for inverse geometric transformation - convergence criterion depending on size and position of the real element --- src/bgeot_geotrans_inv.cc | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/src/bgeot_geotrans_inv.cc b/src/bgeot_geotrans_inv.cc index 04b0e007..f6301fdb 100644 --- a/src/bgeot_geotrans_inv.cc +++ b/src/bgeot_geotrans_inv.cc @@ -220,7 +220,9 @@ namespace bgeot base_node x0_real(N); while (res > IN_EPS/100.) { if ((gmm::abs(res - res0) < IN_EPS/100.) || (factor < IN_EPS)) { - converged = false; + // relaxed convergence criterion depending on the size and position + // of the real element + converged = (res < gmm::mat_maxnorm(G) * IN_EPS/100.); return (pgt->convex_ref()->is_in(x) < IN_EPS) && (res < IN_EPS); } if (res > res0) {