[Getfem-commits] [getfem-commits] branch devel-tetsuo-fix-aclocal-python deleted (was a105a274)

2022-09-29 Thread Konstantinos Poulios via Getfem-commits
logari81 pushed a change to branch devel-tetsuo-fix-aclocal-python.

 was a105a274 Fix m4/ac_python_devel.m4 for Python3.10

The revisions that were on this branch are still contained in
other references; therefore, this change does not discard any commits
from the repository.



[Getfem-commits] [getfem-commits] branch devel-tetsuo-fix-error-in-ubuntu20.04 deleted (was bc309aa0)

2022-09-29 Thread Konstantinos Poulios via Getfem-commits
logari81 pushed a change to branch devel-tetsuo-fix-error-in-ubuntu20.04.

 was bc309aa0 Fix the compile error of nbp in ubuntu20.04

The revisions that were on this branch are still contained in
other references; therefore, this change does not discard any commits
from the repository.



[Getfem-commits] (no subject)

2022-09-29 Thread Konstantinos Poulios via Getfem-commits
branch: master
commit f50e3cbab06ca59e87decb0aea6f4fdfa1cf1b4c
Author: Konstantinos Poulios 
AuthorDate: Thu Sep 29 14:15:54 2022 +0200

Allow integration point variables for GWFL assembly in the scripting 
interface
---
 interface/src/gf_asm.cc | 6 --
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/interface/src/gf_asm.cc b/interface/src/gf_asm.cc
index a223b15e..1824fe01 100644
--- a/interface/src/gf_asm.cc
+++ b/interface/src/gf_asm.cc
@@ -488,7 +488,9 @@ static void do_high_level_generic_assembly(mexargs_in& in, 
mexargs_out& out) {
   nbdof += mf->nb_dof();
   workspace.add_fem_variable(varname, *mf, I, vectors[varname]);
 }  else if (mimd) {
-  THROW_BADARG("Data defined on integration points can not be a 
variable");
+  gmm::sub_interval I(nbdof, mimd->nb_filtered_index());
+  nbdof += mimd->nb_filtered_index();
+  workspace.add_im_variable(varname, *mimd, I, vectors[varname]);
 }  else {
   gmm::sub_interval I(nbdof, U.size());
   nbdof += U.size();
@@ -604,7 +606,7 @@ static void do_expression_analysis(mexargs_in& in, 
mexargs_out& out) {
   if (mf)
 workspace.add_fem_variable(varname, *mf, dummy_I, dummy_V);
   else if (mimd) {
-THROW_BADARG("Data defined on integration points can not be a 
variable");
+workspace.add_im_variable(varname, *mimd, dummy_I, dummy_V);
   } else
 workspace.add_fixed_size_variable(varname, dummy_I, dummy_V);
 }



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

2022-09-29 Thread Konstantinos Poulios via Getfem-commits
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 
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) {



[Getfem-commits] (no subject)

2022-09-29 Thread Konstantinos Poulios via Getfem-commits
branch: master
commit a499492ecdb32c9b73e83148541d79297b9cdbae
Author: Konstantinos Poulios 
AuthorDate: Thu Sep 29 13:16:47 2022 +0200

Add const attribute to class member function
---
 src/getfem/getfem_generic_assembly_tree.h | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/getfem/getfem_generic_assembly_tree.h 
b/src/getfem/getfem_generic_assembly_tree.h
index 9c9ba7ac..4384b9fd 100644
--- a/src/getfem/getfem_generic_assembly_tree.h
+++ b/src/getfem/getfem_generic_assembly_tree.h
@@ -235,7 +235,7 @@ namespace getfem {
 void set_sparsity(int sp, size_type q)
 { sparsity_ = sp; qdim_ = q; }
 
-size_type qdim() { return is_copied ? tensor_copied->qdim() : qdim_; }
+size_type qdim() const { return is_copied ? tensor_copied->qdim() : qdim_; 
}
 
 int sparsity() const
 { return is_copied ? tensor_copied->sparsity() : sparsity_; }



[Getfem-commits] [getfem-commits] master updated (49aedfae -> f50e3cba)

2022-09-29 Thread Konstantinos Poulios via Getfem-commits
logari81 pushed a change to branch master.

from 49aedfae fix a problem with a non adaptation of mim to mesh change
 new a499492e Add const attribute to class member function
 new 8cfa3b3b Typo fix and documentation improvements
 new f50e3cba Allow integration point variables for GWFL assembly in the 
scripting interface


Summary of changes:
 doc/sphinx/source/replaces.txt|  2 +-
 doc/sphinx/source/userdoc/bfem.rst| 50 +++
 doc/sphinx/source/userdoc/bmesh.rst   |  2 +-
 doc/web/index.php |  4 +--
 interface/src/gf_asm.cc   |  6 ++--
 src/getfem/getfem_generic_assembly_tree.h |  2 +-
 src/gmm/gmm_dense_matrix_functions.h  |  2 +-
 7 files changed, 35 insertions(+), 33 deletions(-)