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 <logar...@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++; }