A non-linear system of 100 unknowns of type xi and 10 of
type wi is quite large, symbolically.
If I understand your French, you have solved for the xi in terms
of the wi, and gotten through 82 equations.
Maybe, instead of writing a Gauss-Jordan program, you should try
the built-in linsolve or solve programs, starting with much smaller
examples.
There may be bugs in the Maxima system, but it is also possible that,
since you are running at the limits of your memory, you have found
bugs in the underlying Lisp system, and that just running on a larger
machine, or using a different implementation could allow this to work.
It is entirely possible that the work needed to go from the 83eme to the
final equation is many many times larger than the work up to that time,
and that the answer, if you obtain it, will be so large as to be
uninteresting.
I think that if you posted all the equations, someone else might
try to solve it. Even better would be to post a subset (say 15)
that might be interesting, and see if there is a better way of
solving that system.
Rewriting Maxima data structures would not be high on the list
of things to try.
RJF
laurent couraud wrote:
> Bonjour a tous,
>
> Le système d'équation que j'essai de résoudre est de la forme suivante:
>
> Toute les équations sont de la forme :
> w1*x00+w2*x01+w1*x02+w2*x10+x11+w2*x12+w1*x20+w2*x21+w1*x22-p1
> ou les wi et les xi sont les inconnues et les pi des paramètres.
>
> Et le plus petit sytème que j'ai put trouver jusque là comporte 100
> inconnues de type xi et 10 inconnues de type wi. Ou plus présisement
> c'est le plus petit systeme me permetant de resoudre potentielement
> les xi en fonction des wi et des pi. J'ai essaye de resoudre par la
> methode de GAUS-JORDAN les xi en fonction des wi et des pi et meme
> ca je n'y arrive pas. Le mieux que j'ai pu faire jusqu'à present
> c'est d'éliminé les 82 premier xi mais arrivé au 83eme la capacité
> memoire utilisé par Maxima passe de 70 ou 80 Mo à plus de 450Mo.
> Comme je n'ai que 512Mo de RAM le systeme swap et donc je suis obligé
> d'arreter le calcul. De toute facon il me semble tres currieux que
> le volume de la matrice augmente autant en un seul passage.
> Je soupconne la presence d'un BUG dans Maxima. Si tel est le cas
> je soupconne qu'il ce trouve dans la fonction "ratsimp" etant donne
> que lorsque je ne l'utilise pas la memoire grossi moins vite. La
> fonction "gcd" aussi à l'air d'etre gourmande en memoire, et ne
> semble pas liberer la memoire temporaire quelle utilise. Ceci
> etant dit j'ai peut etre trouver un moyen de reduire un peut la
> taille du systeme. Mais je ne descendrais guere au dessous de
> 90 inconnus de type xi.
> Actuelement je suis en train de reecrire mon algo de resolution
> en remplacant les liste par des tableaux lorsque c'est possible.
> Mais j'ai quelque soucis.
>
> Merci à tous pour votre aide.
>
> _______________________________________________
> Maxima mailing list
> Maxima@www.math.utexas.edu
> http://www.math.utexas.edu/mailman/listinfo/maxima