Python-Sage
Bonjour,
Est il possible d'importer Sage (ou l'équivalent, Sagemath ...) dans un EDI comme Pyzo par exemple( et pas le contraire) ?
Je voudrais pouvoir disposer dans Pyzo de fonction comme prime, nextprime pour commencer, puis pouvoir faire un peu de théorie des nombres, puis un peu d'algèbre théorique, déclarer des groupes, des anneaux, faire ce que font CQ et d'autres.
Bien sûr, si c'est pour me dire de me mettre à Sage lui-même, ce n'est pas la peine.
Cordialement,
Rescassol
Est il possible d'importer Sage (ou l'équivalent, Sagemath ...) dans un EDI comme Pyzo par exemple( et pas le contraire) ?
Je voudrais pouvoir disposer dans Pyzo de fonction comme prime, nextprime pour commencer, puis pouvoir faire un peu de théorie des nombres, puis un peu d'algèbre théorique, déclarer des groupes, des anneaux, faire ce que font CQ et d'autres.
Bien sûr, si c'est pour me dire de me mettre à Sage lui-même, ce n'est pas la peine.
Cordialement,
Rescassol
Connectez-vous ou Inscrivez-vous pour répondre.
Réponses
-- Schnoebelen, Philippe
Oui, j'ai essayé ça, mais rien.
Cordialement,
Rescassol
Ben, yaka lire ça pour voir qu'il faut être un geek, ce que je ne suis pas.
Cordialemennt,
Rescassol
Ce que j'ai mentionné dans l'autre fil est un peu compliqué car il s'agit de compiler SageMath, ceci car Dr. Geodingus a déjà essayé sans succès l'installation « normale ». Mais si cette installation « normale » marche bien dans ton cas, cela ne doit pas être très compliqué.
Si tu es sous windows, je ne pense pas que tu arrives à faire ce que tu souhaites car ton installation SageMath passe par Cygwin et n'est pas accessible en dehors de cette interface.
Cordialement,
Geodingus
-- Schnoebelen, Philippe
En utilisant l'interpréteur Python de SageMath, il est possible d'appeler des fonctions SageMath au sein d'un programme Python « normal ». La preuve avec cette adaptation directe du code de cette question à Python 3 (oui, SageMath est passé à Python 3 depuis janvier dernier avec la version 9.0) :
#! /chemin/vers/mon/SAGE_ROOT/sage -python import os import sys from sage.all import * if len(sys.argv) != 2: print("Usage: {} N".format(os.path.basename(sys.argv[0]))) print("Output the prime factorization of N.") sys.exit(1) print(factor(sage_eval(sys.argv[1])))Test : La première exécution est très lente ; après, ça va beaucoup plus vite chez moi (effet des divers niveaux de cache, sans doute).Tu as un installeur Windows pour Sage
https://wiki.sagemath.org/SageWindows
J'ai vu tout un tas de gens installer Sage sur leur système Windows.
Je ne pense pas que tu sois plus maladroit que le plupart de ces gens !
giacpy
J'envisage Linux dans quelque temps, mais pour le moment, je suis sous Windows 10.
Ce que je veux est très simple.
Dans Pyzo, et pas autre chose, faire tourner: ou:
Cordialement,
Rescassol
from giacpy import * def primes(n): l=[] k=1 for j in range(n): k=nextprime(k+1) l.append(k) return lensuite primes(20) affiche la liste des 20 premiers premiers (les experts Python sauront mieux que moi creer un generateur avec nextprime). La commande isprime fonctionne.Pour travailler avec des entiers de Gauss il suffit d'ecrire i=giac('i'), ensuite (2+3*i)*(1+i) fait ce qui est attendu. Pour calculer dans Z[sqrt(5)*i], il faut faire des commande normal(), par exemple normal((1+sqrt(5)*i)/(2-sqrt(5)*i))
Merci pour la référence à giacpy et cet exemple. Voici une variante de ce code avec primes() comme fonction génératrice (je traduis bêtement generator function, aucune idée s'il y a une traduction plus canonique) :
#! /usr/bin/env python3 import giacpy as giac import sys def primes(n): l = [] k = 1 for j in range(n): k = giac.nextprime(k+1) l.append(k) yield k for p in primes(int(sys.argv[1])): print(p)Pour ceux qui ne connaissent pas, le 'yield' suspend l'exécution de la fonction jusqu'à la prochaine fois où l'itérateur obtenu par évaluation de primes(int(sys.argv[1])) est sommé de fournir une valeur. C'est un moyen très économique d'itérer sur une liste d'objets Python qui « coûtent » un peu cher (calculs, entrées/sorties...) sans avoir besoin de tous les calculer et stocker à l'avance. C'est particulièrement intéressant si l'on décide de terminer l'itération de manière prématurée, en n'ayant utilisé que les premiers objets retournés par l'itérateur.Je conseille de ne pas utiliser 'from giacpy import *' car cela écrase au minimum la fonction int() de Python... ce qui n'est pas très agréable.
Exemple d'appel du programme ci-dessus : Ce serait bien si l'on pouvait éviter l'affichage des lignes suivantes :
Dans ton script, brian, j’enlèverais toute référence à l, on peut importer un seul objet d’un module et on peut envoyer j à la benne :
#! /usr/bin/env python3 from giacpy import nextprime import sys def primes(n): k = 1 for _ in range(n): k = nextprime(k+1) yield k for p in primes(int(sys.argv[1])): print(p)-- Schnoebelen, Philippe
-- Schnoebelen, Philippe
J'ai pu installer giacpy dans Pyzo facilement.
Pour les nombres premiers, ok, mais où peut on trouver une doc sur giacpy, de préférence en français ?
Je ne l'ai pas encore exploré, mais fait il des anneaux, des idéaux, des groupes etc ... ?
Cordialement,
Rescassol
#! /usr/bin/env python3 import giacpy as giac import sys def primes(n): k = 1 for _ in range(n): k = giac.nextprime(k) yield k for p in primes(int(sys.argv[1])): print(p)Test :-- Schnoebelen, Philippe
Tu peux aussi regarder ici, notamment à About Xcas... (Incomplete) documentation. parisse saura sans doute te recommander le ou les documents les plus appropriés.
En ce qui concerne l'algebre, Xcas est un logiciel de calcul symbolique qui a une philosophie proche de maple, i.e. il s'adresse a un public scientifique pas forcement tres avance en maths pures : contrairement a sage, on ne doit pas declarer au prealable qu'on va travailler dans l'anneau des polynomes a coefficients dans Q ou dans Z/5Z, ce sont les commandes qui s'adaptent a leurs arguments. Par exemple
a=giac('1 mod 7') quorem(a*(x**3+x+1),5*x**2+3*x)va calculer le quotient et le reste de la division euclidienne dans Z/7Z[X], alors que quorem(x**3+x+1,5*x**2+3*x) fait le calcul dans Q[X].Il n'y a pas actuellement d'interface giac/gap donc tres peu de commandes de groupes, cela se resume a quelques commandes sur les permutations (p1op2
Merci à tous.
Cordialement,
Rescassol