gbasis : texte "parasite" et réponse manquante

Utilisation de Xcas

Modérateur : xcasadmin

Dominique MICOLLET
Messages : 52
Inscription : jeu. janv. 11, 2007 10:39 am
Localisation : Université de Bourgogne - France

gbasis : texte "parasite" et réponse manquante

Message par Dominique MICOLLET » lun. mars 11, 2024 4:04 pm

Bonjour,

Sous icas (et j'imagine xcas aussi) la commande :

Code : Tout sélectionner

gbasis([X9^2*X3^2-Y10,X9^2*X2^2-Y11,X9^2*X1^2-Y12,X8*X9*X2*X3-Y16,X7*X9*X3^2-Y22,X7*X9*X2^2-Y23,X7*X9*X1^2-Y24,X6*X9*X1*X3-Y26,X5*X9*X1*X2-Y28,X4*X9*X3^2-Y30,X4*X9*X2^2-Y31,X4*X9*X1^2-Y32,X8^2*X3^2-Y39,X8^2*X2^2-Y40,X8^2*X1^2-Y41,X7*X8*X2*X3-Y43,X6*X8*X1*X2-Y45,X5*X8*X1*X3-Y46,X4*X8*X2*X3-Y47,X7^2*X3^2-Y53,X7^2*X2^2-Y54,X7^2*X1^2-Y55,X6*X7*X1*X3-Y57,X5*X7*X1*X2-Y59,X4*X7*X3^2-Y61,X4*X7*X2^2-Y62,X4*X7*X1^2-Y63,X6^4-Y64,X6^2*X3^2-Y67,X6^2*X2^2-Y68,X6^2*X1^2-Y69,X5*X6*X2*X3-Y70,X4*X6*X1*X3-Y71,X5^2*X3^2-Y74,X5^2*X2^2-Y75,X5^2*X1^2-Y76,X4*X5*X1*X2-Y77,X4^2*X3^2-Y79,X4^2*X2^2-Y80,X4^2*X1^2-Y81,X3^4-Y82,X2^2*X3^2-Y83,X1^2*X3^2-Y84,X2^4-Y85,X1^2*X2^2-Y86,X1^4-Y87,X9*X3^2-Y97,X9*X2^2-Y98,X9*X1^2-Y99,X5*X6*X8-Y102,X8*X2*X3-Y103,X7*X3^2-Y109,X7*X2^2-Y110,X7*X1^2-Y111,X6*X1*X3-Y113,X5*X1*X2-Y115,X4*X3^2-Y117,X4*X2^2-Y118,X4*X1^2-Y119,X8^2-Y123,X6^2-Y126,X5^2-Y127,X3^2-Y129,X2^2-Y130,X1^2-Y131,X9-Y132,X7-Y133,X4-Y134],[X1,X2,X3,X4,X5,X6,X7,X8,X9,Y10,Y11,Y12,Y16,Y22,Y23,Y24,Y26,Y28,Y30,Y31,Y32,Y39,Y40,Y41,Y43,Y45,Y46,Y47,Y53,Y54,Y55,Y57,Y59,Y61,Y62,Y63,Y64,Y67,Y68,Y69,Y70,Y71,Y74,Y75,Y76,Y77,Y79,Y80,Y81,Y82,Y83,Y84,Y85,Y86,Y87,Y97,Y98,Y99,Y102,Y103,Y109,Y110,Y111,Y113,Y115,Y117,Y118,Y119,Y123,Y126,Y127,Y129,Y130,Y131,Y132,Y133,Y134],plex);
me renvoie (au bout d'une petite dizaine de minute) une information intermédiaire que je souhaiterais effacer :

Code : Tout sélectionner

// Timing for 2nd run 39.7354 1st run 250.939 speed ratio 0.158347 [current reconstructed ratio for reinjection=0.2 speed_ratio for reinjection=0.166667 modifiable by gbasis_reinject(reconstructed_ratio,speed_ratio) command]
Running a probabilistic check for the reconstructed Groebner basis. If successful, error probability is less than 1e-15 and is estimated to be less than 10^-18. Use proba_epsilon:=0 to certify (this takes more time).
Par ailleurs, la réponse à la commande n'apparaît pas ; seul le temps de calcul est affiché

Code : Tout sélectionner

Evaluation time: 566.16
Done
// Time 566.16
J'ai un souci similaire lorsque j'emploie giac sous forme de bibliothèque dans un programme compilé : l'information intermédiare se mêle aux autres résultats. La réponse de l'appel à gbasis est correct.

J'utilise la version 1.9.9-93 compilée sans option particulière (1).

Auriez-vous une suggestion à me faire ?

(1) : j'ai tenté d'attacher le fichier config.log, sans succès : l'extension n'est pas acceptée.

Cordialement.
Cordialement

Dominique MICOLLET

parisse
Messages : 5741
Inscription : mar. déc. 20, 2005 4:02 pm
Contact :

Re: gbasis : texte "parasite" et réponse manquante

Message par parisse » jeu. mars 14, 2024 2:15 pm

Lorsque le résultat d'un calcul est trop long, icas ne l'affiche pas, et renvoie à la place "Done". La bonne méthode, c'est de stocker le résultat dans une variable puis de sauvegarder la variable dans un fichier, avec

Code : Tout sélectionner

G:=gbasis(...):; write("base",G);
et du coup on n'est pas géné par l'information intermédiaire qui indique le temps nécessaire pour faire le 2ème run avec 4 nombres premiers par rapport au 1er run. Ce temps est utilisé pour savoir s'il faut à un moment donné réinjecter des éléments déjà reconstruits dans la base ou pas. L'affichage est fait dans cocoa.cc ligne 16663 (dans la version actuelle, ça a sans doute un peu bougé), il suffit de mettre un 0 && au début du test if (...) pour le supprimer.
La sauvegarde peut se faire de manière plus efficace avec archive qu'avec write si les données sont de grande taille, mais c'est moins lisible.

Dominique MICOLLET
Messages : 52
Inscription : jeu. janv. 11, 2007 10:39 am
Localisation : Université de Bourgogne - France

Re: gbasis : texte "parasite" et réponse manquante

Message par Dominique MICOLLET » jeu. mars 21, 2024 9:16 am

Merci.
Cordialement

Dominique MICOLLET

Répondre