1.5.0 built without glpk

Librairie C++ de calcul formel/ C++ symbolic computation library

Modérateur : xcasadmin

frederic han
Messages : 1113
Inscription : dim. mai 20, 2007 7:09 am
Localisation : Paris
Contact :

1.5.0 built without glpk

Message par frederic han » mer. sept. 19, 2018 5:33 am

Salut,
la compilation de giac sans glpk ne passe plus avec la 1.5.0 à cause d'un painter dans graphe.cc ligne 4860
qui n'est pas dans un

#ifdef HAVE_LIBGLPK

Aussi j'ai encore un faux doctest failure dans TP21 :

Code : Tout sélectionner

fred@localhost check]$ diff TP21-sol.cas.out1 TP21.tst
77c77
< [1,sqrt(3)*(x+3)/3,sqrt(3)*(x^2+3*x-3)/3,(-I)*sqrt(273)*(3*x^3+20*x^2+21*x-42)/273,273*(-I)*sqrt(32214)*(x^4+398/91*x^3+166/91*x^2+11/39*x+23/13)/32214],
---
> [1,sqrt(3)*(x+3)/3,sqrt(3)*(x^2+3*x-3)/3,(-I)*sqrt(273)*(3*x^3+20*x^2+21*x-42)/273,(-I)*sqrt(32214)*(273*x^4+1194*x^3+498*x^2+77*x+483)/32214],
Fred

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

Re: 1.5.0 built without glpk

Message par parisse » mer. sept. 19, 2018 8:53 am

Chez moi le test 21 passe. Tu as une idee pour qu'il soit moins sensible?

lukamar
Messages : 297
Inscription : ven. juin 30, 2017 9:55 am
Localisation : Zagreb, Croatia

Re: 1.5.0 built without glpk

Message par lukamar » mer. sept. 19, 2018 9:01 am

There was indeed a problem, one line more in graphe.cc needed to be excluded from the compilation without GLPK. I think it should work now.

Luka

frederic han
Messages : 1113
Inscription : dim. mai 20, 2007 7:09 am
Localisation : Paris
Contact :

Re: 1.5.0 built without glpk

Message par frederic han » mer. sept. 19, 2018 10:00 am

parisse a écrit :
mer. sept. 19, 2018 8:53 am
Chez moi le test 21 passe. Tu as une idee pour qu'il soit moins sensible?
Je propose de ne pas afficher le resultat mais de le ranger dans une variable: GS
j'ajoute la ligne matrix pour verifier que l'on a bien une base orthonormee
et je change la derniere ligne pour montrer que les polynomes de GS sont proportionnels à ceux de S. La liste des contantes devrait etre moins sensible aux versions des librairies et systèmes. Ca teste à peu pres la meme chose.

Code : Tout sélectionner

S:=factor(gramschmidt([1,x,x^2,x^3,x^4],(p,q)->scalU(p,q,4))):
matrix(5,5,(u,v)->simplify(scalU(GS[u],GS[v],4)));
S:=[seq(pade(U,x,2*i-1,i),i=1..4)];
/* On constate bien que les denominateurs des approximants de pade coincident a facteur�pres avec les polyn\^omes reciproques des orthonormali$
recip:=proc(P)
normal(x^degree(P)*subs(x=1/x,P))
end;
seq(simplify(recip(denom(S[i]))/GS[i+1]),i=1..4); #donne bien les polyn\^omes obtenus par gramschmidt
du coup chez moi cela donne:

Code : Tout sélectionner

"Done",
matrix[[1,0,0,0,0],[0,1,0,0,0],[0,0,1,0,0],[0,0,0,1,0],[0,0,0,0,1]],
[1/(3*x+1),-1/(3*x^2-3*x-1),(3*x^2-11*x-3)/(42*x^3-21*x^2-20*x-3),(626*x^3+192*x^2+375*x+273)/(483*x^4+77*x^3+498*x^2+1194*x+273)],
proc(P) 
  normal(x^(degree(P))*subs(x=(1/x),P));  
 
end;,
sqrt(3),-sqrt(3),(-I)*sqrt(273),I*sqrt(32214)

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

Re: 1.5.0 built without glpk

Message par parisse » jeu. sept. 20, 2018 9:19 am

C'est fait!

Répondre