Anomalie #288
Messages loi de comportement sur V6.983
Description
Bonjour Gérard,
En relançant un calcul complexe avec V6.983 le fichier de sortie standard a rempli mon disque alors que le calcul passe bien avec l’ancienne version, même si'l peut apparaitre parfois certaines WARNING dans le loi CP.
J'ai essayé de mettre des permet_affichage_ 0 partout où je pouvais sans changement...
Voici les messages lors de l'activation de la loi CP :
- pb sur le calcul de la contrainte: le resultat est soit infini soit un nan mail: 1 ele: 8844 npti: 3
-- deformation: ---
nb noeud concerne = 4 numInteg= 3
-- information concernant les noeuds:
--noeud --
Numero : 286, du maillage: 1
...................Puis 4 millions de lignes :
--- LoiContraintesPlanes::Calcul_SigmaHH --- mail: 1 ele: 3553 npti: 2
..............Après 200 itérations et avant convergence le calcul s'arrête avec le message :
cas de pilotage de fin d'iteration inconnues ? ou alors qui suit une erreur signalee precedemment
Fichiers
Mis à jour par Gérard Rio il y a plus de 3 ans
Bonjour Frank,
je vais regarder...
pourrais-tu joindre le fichier qui manque semble-t-il ??
cas d'un fichier decrit par la commande <
erreur en ouverture du fichier = fonctions_1d_R28.template
Mis à jour par Frank Petitjean il y a plus de 3 ans
- Fichier fonctions_1d_R28.template fonctions_1d_R28.template ajouté
Pardon, je l'avais oublié celui-là puisqu'il est effacé automatiquement à la fin des scripts Python. Le voici
Mis à jour par Gérard Rio il y a plus de 3 ans
- Statut changé de Nouveau à En cours
- % réalisé changé de 0 à 10
Au niveau des commentaires, cela provient du fait que j'ai abaissé le seuil de sortie de message ...
Je vais relever le seuil.
En fait, si on met un niveau local de commentaire à 0, c'est équivalent à ne pas indiquer de niveau. C'est alors le niveau global qui est utilisé: dans ton cas 3, et il se trouve que j'avais abaissé le seuil à 2 !!, je vais remettre à > 3.
Si on veut vraiment, aucun commentaire localement à une loi, il faut indiquer par exemple un nombre négatif (ex: -1). Dans ce cas, il n'y a plus de commentaire du tout (sauf les erreurs détectés).
Pour l'autre pb, au niveau de la fonction nD, je continue mes investigations !!
Mis à jour par Gérard Rio il y a plus de 3 ans
- % réalisé changé de 10 à 90
bon... je crois que j'ai trouvé quelque chose:
en fait dans le message d'erreur il est indiqué que le premier argument vaut une valeur légèrement négative : -9.567e-16
or dans la fonction f_R_bulle qui pose pb, tu utilises une élévation à une puissance non entière, qui donc doit avoir un argument positif (because log...)
J'ai regardé l'intégrale de la fonction F_V qui est calculée et effectivement, au début elle passe par cette valeur légèrement négative ... due peut-être à des erreurs d'arrondi (?)
Donc en résumé, je crois (mais à confirmer) que le fonctionnement est normal...
Maintenant, pourquoi dans une version précédente il n'y avait pas de valeur négative ????
là je ne sais pas ...
Mis à jour par Gérard Rio il y a plus de 3 ans
Je viens de mettre à jour : version V 6.984, avec un seuil de déclenchement de commentaire plus haut.
Peux-tu me dire si cela te convient pour que je ferme éventuellement le ticket ?
Mis à jour par Frank Petitjean il y a plus de 3 ans
Merci Gérard pour ton intervention express.
Je suis passé avec un niveau_commentaire 2 et j'ai plus les premiers messages mais j'ai toujours le message pour tous les pti :
--- LoiContraintesPlanes::Calcul_SigmaHH ---
Avec permet_affichage_ -1 c'est messages n'apparaissent plus. Je n'ai pas vu dans la doc l'option -1 ?
Concernant le calcul du rayon de la bulle qui utilise le volume du ballon j'ai éviter le problème du volume nulle en n'activant l'appel de la fonction f_R_bulle qu'après l'itération 500 (au premier incrément), voir la fonction f_ps qui l'utilise. Pour une meilleur garanti j'ai mis min(1, volume).
Tu peux clore le ticket. Merci.
Frank
Mis à jour par Gérard Rio il y a plus de 3 ans
- Statut changé de En cours à Résolu
- % réalisé changé de 90 à 100
Bonjour Frank,
effectivement le niveau -1 n'est pas indiqué dans la doc mais je pensais que c'était sous-jacent compte tenu du fait que seuls les niveaux supérieur à 0 sont effectif.
bon... je vais l'indiquer pour que cela soit clair (et c'est vrai pour tous les niveaux)
Pour le volume, effectivement il suffit de s'assurer que les valeurs utilisées soient toujours > 0
La fonction analytique appelle la fonction en dur d'élévation à la puissance qui ramène nan dès que c'est un chouilla négatif ...
Je crois que j'ai mis récemment des tests supplémentaires au cas où on avait un nan, et c'est sans doute cela qui a changé (mais je ne suis pas sûr) ...
Je ferme le ticket,
à +
Mis à jour par Frank Petitjean il y a plus de 3 ans
- Fichier Extrait fichier log.txt Extrait fichier log.txt ajouté
- Fichier modele.info modele.info ajouté
Entre la 982 et la 984 le format numérique de la sortie standard a changé : il y a 17 décimales au lieu de 6.
Le niveau de commentaire général à 2 me va bien il n'y a pas le n° des itérations, c'est dommage...
Comment supprimer les infos concernant les forces de stabilisation de membranes ?
limitation stabilisation due aux forces externes, F_stab gene=
Avec le niveau de commentaire global à 3 j'ai mis des permet_affichage_ -1 partout dans la loi de comportement et j'ai encore des messages trop détaillés pour mon usage courant. Est-il possible de revenir comme avec la V6.982 ? Ci-joint un extrait de sortie.
Frank
Mis à jour par Gérard Rio il y a plus de 3 ans
- Statut changé de Résolu à En cours
- % réalisé changé de 100 à 90
bonjour Frank,
il me manque le fichier : plan_sym.her, pour faire le test.
Concernant les forces de stabilisation de membranes,
- je ne vois pas de sortie dans ton fichier exemple
- si le niveau global est > à 0, il y a systématiquement affichage de la force... effectivement c'est un peu beaucoup!
- pour l'affichage dans la loi: normalement c'est uniquement les messages d'erreurs qu'il devrait avoir : c'est d'ailleurs le cas dans ton exemple.
Ceci étant, effectivement avec -1, on ne devrait plus rien avoir... il manque sur ton cas, le niveau -1 pour les lois élastique, cependant j'attends ton fichier pour faire le test.
- pour les 17 décimales... c'est rigolo, car je n'ai rien changé à ce niveau... je pense que cela vient des sorties précédentes ??? bon, je vais spécifier le nb de déci (qui n'est pas fait à ce niveau) ...
à suivre !
Mis à jour par Frank Petitjean il y a plus de 3 ans
- Fichier plan_sym.her plan_sym.her ajouté
Décidément j'oublie toujours un fichier, désolé.
L'affichage des forces de stabilisation n'est pas trop gênante mais il pourrait se produire avec un niveau global supérieur à 3 ou 4 par exemple et non simplement >0.
Je vais ajouter le niveau -1 pour les lois élastique, sait-on jamais !
Mis à jour par Gérard Rio il y a plus de 3 ans
en fait il manque également les ref:
E_biel_bord1 E_biel_bord2 E_biel_centre ...
peut-être que ce n'est pas le bon fichier de maillage ??
Mis à jour par Frank Petitjean il y a plus de 3 ans
- Fichier fonctions_1d_R28.template fonctions_1d_R28.template ajouté
- Fichier fonctionsAtmo.dat fonctionsAtmo.dat ajouté
- Fichier maillage.her maillage.her ajouté
- Fichier modele.CVisu modele.CVisu ajouté
- Fichier modele.info modele.info ajouté
Je n'avais pas identifié clairement le dossier correspondant aux fichiers envoyés, j'ai essayé de le retrouver avec le .info mais sans succès...
Je relancé un cas en ajoutant un permet_affichage_ -1 pour la loi élastique et je n'ai plus les mêmes messages d'erreur. Il reste celui-ci
l'application du critere conduit a une contrainte nulle...
Manque t'il un permet_affichage_ -1 quelque part ?
Cette fois j'ai mis tous les fichiers !
Frank
Mis à jour par Frank Petitjean il y a plus de 3 ans
Bonjour Gérard,
Je relance la demande : est-il possible de supprimer de la sortie standard ce massage :
**>> l'application du critere conduit a une contrainte nulle....
C'est très pénalisant en terme d'affichage.
Merci
Frank
Mis à jour par Frank Petitjean il y a plus de 3 ans
Bonjour Gérard,
J'ai résolu mon problème d'affichage non désiré dans les lois de comportement. Il restait un dernier endroit où je pouvais placer l'instruction permet_affichage_ -1 :
parametres_controle_pour_contraintes_planes_double_
NEWTON_LOCAL avec_parametres_de_reglage_
permet_affichage_ -1
fin_parametres_reglage_Algo_Newton_
fin_parametres_controle_pour_contraintes_planes_double_
Tu peux fermer le ticket, merci
Frank
Mis à jour par Gérard Rio il y a plus de 3 ans
merci Frank pour ton retour.
De mon coté voici ce que j'ai fait:
- pour la loi critère:
. par défaut, l'affichage des pb, s'effectue si le niveau de commentaire est supérieur à 4
. si l'affichage est supérieur à 5 toutes les informations intermédiaires sont affichées sauf les composantes des opérateurs tangent (qui sont nombreuses ...)
. si l'affichage est supérieur à 6 : affichage des composantes des opérateurs tangents
Un niveau élevé de commentaire est seulement intéressant en ciblant un pti particulier ou un élément, via un critère particulier suite un pb que l'on recherche, ceci à l'aide d'une fonction nD pour piloter le niveau de commentaire (cf. doc)
Je vais l'indiquer dans la documentation,
- pour l'affichage du nombre de décimal, en fait il y a un paramètre de contrôle (cf. paramètres de contrôle de l'affichage) qui par défaut conduit à 17 décimales ceci pour l'écran et les fichiers de sauvegarde. On peut changer l'affichage, mais cela change aussi la sauvegarde sur le .BI ce qui peut influencer les restart.
à tester... avec la version 6.985
Mis à jour par Frank Petitjean il y a plus de 3 ans
Bonjour Gérard,
Merci pour ces évolutions, je vais pouvoir supprimer la dizaines de permet_affichage_ -1 !
Je n'ai pas trouvé dans la doc (paramètres de contrôle de l'affichage) le moyen de contrôler le nombre de décimale mais maintenant il est revenu à la normale.
Frank
Mis à jour par Gérard Rio il y a plus de 3 ans
Bonjour Frank,
est-ce que c'est ok pour toi maintenant ?
Mis à jour par Frank Petitjean il y a plus de 3 ans
Oui c'est ok. Je ne suis plus embêté par ces messages avec la commande permet_affichage_ 0 aux bons endroits. Merci
Frank
Mis à jour par Gérard Rio il y a plus de 3 ans
- Statut changé de En cours à Résolu
- % réalisé changé de 90 à 100
Parfait !