[Getfem-commits] [getfem-commits] branch master updated: Relax convergence threshold for inverse geometric transformation

[Konstantinos Poulios]

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 <logari81@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) {