[Getfem-commits] [getfem-commits] branch master updated: Fix non-UTF8 file encodings

[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 e41de9db Fix non-UTF8 file encodings
e41de9db is described below

commit e41de9dbc9a1172ad6b2025ecfee860c0e90e327
Author: Konstantinos Poulios <logari81@gmail.com>
AuthorDate: Wed Jun 14 20:55:13 2023 +0200

    Fix non-UTF8 file encodings
---
 src/bgeot_geotrans_inv.cc                |  4 ++--
 src/getfem/bgeot_sparse_tensors.h        | 38 ++++++++++++++++----------------
 src/getfem/getfem_assembling_tensors.h   |  2 +-
 src/getfem/getfem_mesher.h               |  6 ++---
 src/getfem_contact_and_friction_nodal.cc |  2 +-
 src/getfem_error_estimate.cc             |  4 ++--
 src/getfem_mat_elem.cc                   |  4 ++--
 src/getfem_mesh_slicers.cc               |  2 +-
 src/gmm/gmm_solver_idgmres.h             | 32 +++++++++++++--------------
 9 files changed, 47 insertions(+), 47 deletions(-)

diff --git a/src/bgeot_geotrans_inv.cc b/src/bgeot_geotrans_inv.cc
index f6301fdb..175fbfb0 100644
--- a/src/bgeot_geotrans_inv.cc
+++ b/src/bgeot_geotrans_inv.cc
@@ -99,8 +99,8 @@ namespace bgeot
       gmm::mult(K, CS, B);
     }
     else {
-      // L'inversion peut �tre optimis�e par le non calcul global de B
-      // et la resolution d'un syst�me lin�aire.
+      // L'inversion peut être optimisée par le non calcul global de B
+      // et la resolution d'un système linéaire.
       base_matrix KT(K.nrows(), K.ncols());
       pgt->compute_K_matrix(G, pc, KT);
       gmm::copy(gmm::transposed(KT), K);
diff --git a/src/getfem/bgeot_sparse_tensors.h 
b/src/getfem/bgeot_sparse_tensors.h
index d2283811..768f9b92 100644
--- a/src/getfem/bgeot_sparse_tensors.h
+++ b/src/getfem/bgeot_sparse_tensors.h
@@ -76,7 +76,7 @@
 
 namespace bgeot {
   typedef gmm::uint32_type index_type;
-  typedef gmm::int32_type stride_type; /* sign�! */
+  typedef gmm::int32_type stride_type; /* signé! */
 
   //  typedef std::vector<index_type> tensor_ranges;
   class tensor_ranges : public std::vector<index_type> {
@@ -313,7 +313,7 @@ namespace bgeot {
     tensor_mask_container masks_;
 
     /* verifie si un masque est completement vide,
-       si c'est le cas alors tous les autres masques sont vid�s
+       si c'est le cas alors tous les autres masques sont vidés
        (le tenseur est identiquement nul) */
     void check_empty_mask() {
       if (card() == 0) {
@@ -506,13 +506,13 @@ namespace bgeot {
       update_idx2mask();
     }
 
-    /* forme d'une tranche (c'est la forme qu'on applique � un tenseur pour
+    /* forme d'une tranche (c'est la forme qu'on applique à un tenseur pour
        en extraire la tranche) */
     tensor_shape slice_shape(tensor_mask::Slice slice) const {
       assert(slice.dim < ndim() && slice.i0 < dim(slice.dim));
       tensor_shape ts(ndim());
       ts.push_mask(tensor_mask(dim(slice.dim), slice));
-      ts.merge(*this); /* le masque peut se retrouver brutalement vid� si on a 
tranch� au mauvais endroit! */
+      ts.merge(*this); /* le masque peut se retrouver brutalement vidé si on a 
tranché au mauvais endroit! */
       return ts;
     }
 
@@ -548,8 +548,8 @@ namespace bgeot {
   */
   class tensor_ref : public tensor_shape {
     std::vector< tensor_strides > strides_;
-    TDIter *pbase_; /* pointeur sur un pointeur qui designe les donn�es
-                      �a permet de changer la base pour toute une serie
+    TDIter *pbase_; /* pointeur sur un pointeur qui designe les données
+                      ça permet de changer la base pour toute une serie
                       de tensor_ref en un coup */
 
     stride_type base_shift_;
@@ -578,7 +578,7 @@ namespace bgeot {
 
     
 
-    /* s'assure que le stride du premier indice est toujours bien �gal � z�ro 
*/
+    /* s'assure que le stride du premier indice est toujours bien égal à zéro 
*/
     void  ensure_0_stride() {
       for (index_type i=0; i < strides_.size(); ++i) {
        if (strides_[i].size() >= 1 && strides_[i][0] != 0) {
@@ -589,8 +589,8 @@ namespace bgeot {
       }
     }
 
-    /* constructeur � partir d'une forme : ATTENTION ce constructeur n'alloue 
pas la
-       m�moire n�cessaire pour les donn�es !! */
+    /* constructeur à partir d'une forme : ATTENTION ce constructeur n'alloue 
pas la
+       mémoire nécessaire pour les données !! */
     explicit tensor_ref(const tensor_shape& ts) : tensor_shape(ts), pbase_(0), 
base_shift_(0) {
       strides_.reserve(16);
       init_strides();
@@ -614,8 +614,8 @@ namespace bgeot {
 
     void set_sub_tensor(const tensor_ref& tr, const tensor_shape& sub);
 
-    /* constructeur � partir d'un sous-tenseur � partir d'un tenseur et d'une 
forme 
-       hypothese: la forme 'sub' doit �tre un sous-ensemble de la forme du 
tenseur
+    /* constructeur à partir d'un sous-tenseur à partir d'un tenseur et d'une 
forme
+       hypothese: la forme 'sub' doit être un sous-ensemble de la forme du 
tenseur
     */
     explicit tensor_ref(const tensor_ref& tr, const tensor_shape& sub) {
       set_sub_tensor(tr,sub);
@@ -658,7 +658,7 @@ namespace bgeot {
       if (n < pi.n) return true;
       else return false;
     }
-    stride_type mean_increm; /* valeur moyenne de l'increment (utilis� pour le 
tri) */
+    stride_type mean_increm; /* valeur moyenne de l'increment (utilisé pour le 
tri) */
     tensor_strides inc; /* not strides but increments to the next index value,
                                     with inc[range-1] == -sum(inc[0..range-2]) 
(automatic rewinding!) 
                                     of course, stride_type MUST be signed
@@ -680,13 +680,13 @@ namespace bgeot {
     tensor_strides itbase;
     struct  index_value_data {
       dim_type cnt_num;
-      const stride_type **ppinc; /* pointe vers pr[cnt_num].pinc, initialis� 
par rewind()
-                                 et pas avant (� cause de pbs lors de la copie 
de multi_tensor_iterator sinon) 
-                                 permet de d�duire la valeur du compteur: 
(*ppinc - pincbase) (� diviser par nn=(pri[cnt_num].n-N))
+      const stride_type **ppinc; /* pointe vers pr[cnt_num].pinc, initialisé 
par rewind()
+                                 et pas avant (à cause de pbs lors de la copie 
de multi_tensor_iterator sinon)
+                                 permet de déduire la valeur du compteur: 
(*ppinc - pincbase) (à diviser par nn=(pri[cnt_num].n-N))
                               */
       const stride_type *pincbase;
       const stride_type *pposbase; /* pointe dans pri[cnt_num].mask_pos, 
retrouve la position dans le masque en fonction
-                                 du compteur d�duit ci-dessus et des champs 
div et mod ci-dessous */
+                                 du compteur déduit ci-dessus et des champs 
div et mod ci-dessous */
       index_type div, mod, nn;
       stride_type pos_; /* stores the position when the indexe is not part of 
the pri array
                          (hence the index only has 1 value, and ppos == &pos_, 
and pcnt = &zero */
@@ -753,8 +753,8 @@ namespace bgeot {
     bool vnext() { return next(unsigned(-1), vectorized_pr_dim); }
     bool bnext(dim_type b) { return next(bloc_rank[b]-1, bloc_rank[b+1]-1); }
     bool bnext_useful(dim_type b) { return bloc_rank[b] != bloc_rank[b+1]; }
-    /* version speciale pour it�rer sur des tenseurs de m�me dimensions
-       (doit �tre un poil plus rapide) */    
+    /* version speciale pour itérer sur des tenseurs de même dimensions
+       (doit être un poil plus rapide) */
     bool qnext1() {
       if (pr.size() == 0) return false;
       std::vector<packed_range>::reverse_iterator p_ = pr.rbegin();
@@ -870,7 +870,7 @@ namespace bgeot {
     virtual ~tensor_reduction() { clear(); }
     void clear();
 
-    /* renvoie les formes diagonalis�es 
+    /* renvoie les formes diagonalisées
        pour bien faire, il faudrait que cette fonction prenne en argument
        le required_shape de l'objet ATN_reducted_tensor, et fasse le merge
        avec ce qu'elle renvoie... non trivial
diff --git a/src/getfem/getfem_assembling_tensors.h 
b/src/getfem/getfem_assembling_tensors.h
index 89559fec..a2490a48 100644
--- a/src/getfem/getfem_assembling_tensors.h
+++ b/src/getfem/getfem_assembling_tensors.h
@@ -589,7 +589,7 @@ namespace getfem {
       int n=int(tok_number_dval());
       if (n != curr_tok_dval) ASM_THROW_PARSE_ERROR("not an integer"); 
       if (n > maxval) ASM_THROW_PARSE_ERROR("out of bound integer");
-      return n-1; /* -1 pour un indicage qui commence � 1! */ 
+      return n-1; /* -1 pour un indicage qui commence à 1! */
     }
     size_type tok_mfref_num()
     { assert(tok_type()==MFREF); return curr_tok_ival; }
diff --git a/src/getfem/getfem_mesher.h b/src/getfem/getfem_mesher.h
index 7e41bf22..4a072298 100644
--- a/src/getfem/getfem_mesher.h
+++ b/src/getfem/getfem_mesher.h
@@ -480,7 +480,7 @@ namespace getfem {
        for (size_type k = 1; k < dists.size(); ++k)
          d = std::min(d, (*(dists[k]))(P));
       }
-      else { // essai rat� ...
+      else { // essai raté ...
        isin = false;
        for (size_type k = 0; k < dists.size(); ++k) {
          vd[k] = (*(dists[k]))(P);
@@ -506,7 +506,7 @@ namespace getfem {
        }
        return d;
       }
-      else { // essai rat� ...
+      else { // essai raté ...
        vd[0] = (*(dists[0]))(P);
        bool ok = (vd[0] > -SEPS);
        for (size_type k = 1; k < dists.size(); ++k) {
@@ -534,7 +534,7 @@ namespace getfem {
        }
        return dists[i]->grad(P, G);
       }
-      else { // essai rat� ...
+      else { // essai raté ...
        d = operator()(P);
        base_small_vector Gloc;
        for (size_type k = 0; k < dists.size(); ++k) {
diff --git a/src/getfem_contact_and_friction_nodal.cc 
b/src/getfem_contact_and_friction_nodal.cc
index 7776ec2f..44efc994 100644
--- a/src/getfem_contact_and_friction_nodal.cc
+++ b/src/getfem_contact_and_friction_nodal.cc
@@ -528,7 +528,7 @@ namespace getfem {
     int augmentation_version; // 0 for non-symmetric Alart-Curnier version
                               // 1 for symmetric Alart-Curnier version
                               // 2 for new version (augmented multipliers)
-                              // 3 for new version with De Saxc� projection
+                              // 3 for new version with De Saxcé projection
 
     void init_BBN_BBT(void) const {
       gmm::resize(BBN1, gmm::mat_nrows(BN1), gmm::mat_ncols(BN1));
diff --git a/src/getfem_error_estimate.cc b/src/getfem_error_estimate.cc
index 21673bb4..ec759f44 100644
--- a/src/getfem_error_estimate.cc
+++ b/src/getfem_error_estimate.cc
@@ -93,11 +93,11 @@ namespace getfem {
           for (size_type j = 0; j < N; ++j)
             res[i] += (lambda + mu) * hess1(j, i*N+j) + mu * hess1(i, 
j*N+j)+F[i];
        
-        ERR[cv] += 
radius*radius*ctx1.J()*pai1->coeff(ii)*gmm::vect_norm2_sqr(res); //norme carr�e
+        ERR[cv] += 
radius*radius*ctx1.J()*pai1->coeff(ii)*gmm::vect_norm2_sqr(res); //norme carrée
         eta1 += 
(radius*radius*ctx1.J()*pai1->coeff(ii)*gmm::vect_norm2_sqr(res));
       }    
         //if (ERR[cv] > 100)
-        //cout << "Erreur en r�sidu sur element " << cv << " : " << ERR[cv] << 
endl;    
+        //cout << "Erreur en résidu sur element " << cv << " : " << ERR[cv] << 
endl;
         
         
        // jump of the stress between the element ant its neighbors.
diff --git a/src/getfem_mat_elem.cc b/src/getfem_mat_elem.cc
index 68236326..f99b628b 100644
--- a/src/getfem_mat_elem.cc
+++ b/src/getfem_mat_elem.cc
@@ -138,7 +138,7 @@ namespace getfem {
                       "Exact integration not allowed in this context");
 
           if ((*it).t != GETFEM_NONLINEAR_ && !((*it).pfi->is_equivalent())) {
-            // TODO : le numero d'indice � reduire peut changer ...
+            // TODO : le numero d'indice à reduire peut changer ...
             trans_reduction.push_back(k);
             trans_reduction_pfi.push_back((*it).pfi);
           }
@@ -445,7 +445,7 @@ namespace getfem {
         std::fill(mref.begin()+1, mref.end(), aux);
       }
 
-      if (is_ppi) // pour accelerer, il faudrait pr�calculer les d�riv�es
+      if (is_ppi) // pour accelerer, il faudrait précalculer les dérivées
       {
         base_poly P(dim, 0), Q(dim, 0), R(dim, 0);
         size_type old_ind = size_type(-1), ind;
diff --git a/src/getfem_mesh_slicers.cc b/src/getfem_mesh_slicers.cc
index f2f12ea7..cdc6a409 100644
--- a/src/getfem_mesh_slicers.cc
+++ b/src/getfem_mesh_slicers.cc
@@ -672,7 +672,7 @@ namespace getfem {
       gmm::mult(G,pc,K);
       scalar_type J = bgeot::lu_det(&(*(K.begin())), pgt->dim());
       // bgeot::geotrans_interpolation_context ctx(pgp,0,G);
-      // scalar_type J = gmm::lu_det(ctx.B()); // pb car inverse K m�me
+      // scalar_type J = gmm::lu_det(ctx.B()); // pb car inverse K même
       if (J < 0) return true;
       //cout << "cv = " << cv << ", J = " << J << "\n";
     }
diff --git a/src/gmm/gmm_solver_idgmres.h b/src/gmm/gmm_solver_idgmres.h
index d6cc0460..d675dc81 100644
--- a/src/gmm/gmm_solver_idgmres.h
+++ b/src/gmm/gmm_solver_idgmres.h
@@ -179,7 +179,7 @@ namespace gmm {
        combine(KS, y, x, i);
        mult(A, gmm::scaled(x, T(-1)), b, w);
        mult(M, w, r);
-       beta = gmm::vect_norm2(r); // + verif sur beta ... � faire
+       beta = gmm::vect_norm2(r); // + verif sur beta ... à faire
        break;
       }
 
@@ -192,13 +192,13 @@ namespace gmm {
       beta = gmm::vect_norm2(r);
       
       mult(Hess, scaled(y, T(-1)), gamma, ztest);
-      // En fait, d'apr�s Caroline qui s'y connait ztest et gam devrait
-      // �tre confondus
-      // Quand on aura v�rifi� que �a marche, il faudra utiliser gam � la 
+      // En fait, d'après Caroline qui s'y connait ztest et gam devrait
+      // être confondus
+      // Quand on aura vérifié que ça marche, il faudra utiliser gam à la
       // place de ztest.
       if (st.tb_def < p) {
         T nss = H(m,m-1) / ztest[m];
-       nss /= gmm::abs(nss); // ns � calculer plus tard aussi
+       nss /= gmm::abs(nss); // ns à calculer plus tard aussi
        gmm::copy(KS.mat(), W); gmm::copy(scaled(r, nss /beta), mat_col(W, m));
        
        // Computation of the oblique matrix
@@ -226,8 +226,8 @@ namespace gmm {
                         sub_matrix(YB,  sub_interval(0, m-st.tb_def)),
                         sub_interval(st.tb_def, m-st.tb_def), 
                         (st.tb_defwant < p)); 
-         // ns *= alpha; // � calculer plus tard ??
-         //  V(:,m+1) = alpha*V(:, m+1); �a devait servir � qlq chose ...
+         // ns *= alpha; // à calculer plus tard ??
+         //  V(:,m+1) = alpha*V(:, m+1); ça devait servir à qlq chose ...
 
 
          //       Clean the portions below the diagonal corresponding
@@ -451,8 +451,8 @@ namespace gmm {
          for (i = k1; i <= k2; ++i) {
             if (i > k1) Givens_rotation(H(i, i-1), H(i+1, i-1), c, s);
             
-           // Ne pas oublier de nettoyer H(i+1,i-1) (le mettre � z�ro).
-           // V�rifier qu'au final H(i+1,i) est bien un r�el positif.
+           // Ne pas oublier de nettoyer H(i+1,i-1) (le mettre à zéro).
+           // Vérifier qu'au final H(i+1,i) est bien un réel positif.
 
             // apply rotation from left to rows of H
            row_rot(sub_matrix(H, sub_interval(i,2), sub_interval(i, kend-i)),
@@ -468,7 +468,7 @@ namespace gmm {
             
             // accumulate e'  Q so residual can be updated k+p
            Apply_Givens_rotation_left(q(0,i), q(0,i+1), c, s);
-           // peut �tre que nous utilisons G au lieu de G* et que
+           // peut être que nous utilisons G au lieu de G* et que
            // nous allons trop loin en k2.
          }
        }
@@ -486,7 +486,7 @@ namespace gmm {
          num = num + 2;    // mark that a complex conjugate
          mark = true;      // pair has been applied
 
-         // Indices de fin de boucle � surveiller... de pr�s !
+         // Indices de fin de boucle à surveiller... de près !
          for (size_type k1 = 0, k3 = 0; k3 != kend-2; k1 = k3+1) {
            k3 = k1;
            while (h(k3+1, k3) != T(0) && k3 < kend-2) ++k3;
@@ -504,13 +504,13 @@ namespace gmm {
                y = H(i+1, i-1);
                z = H(i+2, i-1);
                // Ne pas oublier de nettoyer H(i+1,i-1) et H(i+2,i-1) 
-               // (les mettre � z�ro).
+               // (les mettre à zéro).
              }
 
              hv[0] = x; hv[1] = y; hv[2] = z;
              house_vector(v);
 
-             // V�rifier qu'au final H(i+1,i) est bien un r�el positif
+             // Vérifier qu'au final H(i+1,i) est bien un réel positif
 
              // apply transformation from left to rows of H
              w.resize(kend-i);
@@ -545,8 +545,8 @@ namespace gmm {
          c = x; s = y;
          if (i > k1) Givens_rotation(H(i, i-1), H(i+1, i-1), c, s);
             
-         // Ne pas oublier de nettoyer H(i+1,i-1) (le mettre � z�ro).
-         // V�rifier qu'au final H(i+1,i) est bien un r�el positif.
+         // Ne pas oublier de nettoyer H(i+1,i-1) (le mettre à zéro).
+         // Vérifier qu'au final H(i+1,i) est bien un réel positif.
 
          // apply rotation from left to rows of H
          row_rot(sub_matrix(H, sub_interval(i,2), sub_interval(i, kend-i)),
@@ -653,7 +653,7 @@ namespace gmm {
     for (j = 0, ind = 0; j < m-p; ++j) {
       if (ritznew[j] == R(-1)) {
        if (std::imag(eval[j]) != R(0)) {
-         st.nb_nolong += 2; ++j; //  � adapter dans le cas complexe ...
+         st.nb_nolong += 2; ++j; //  à adapter dans le cas complexe ...
        } 
        else st.nb_nolong++;
       }