Mathemagix

[Joris van der Hoeven <address@concealed>]

Hello Johann,

Please notice that the official language of for our mailing is English.

On Mon, Aug 29, 2005 at 04:41:01PM +0200, address@concealed wrote:
> je suis étudiant a l'ensimag et je suis en stage a l'université de limoges.
> j'ai travaillé avec olivier ruatta sur l'algorithme de strassen pour le 
> produit
> matriciel (olivier devrait bientot mettre a jour les sources). Je suis assez
> interressé par ce que vous faites et j'aimerais bien apporter d'autres 
> petites
> contributions. Pour finir l'algorithme de Strassen, il reste à mettre en 
> place
> un algorithme de decision quand les tailles ne sont pas des puissances de 2
> (pour l'instant, j'agrandis la matrice a la puissance de 2 supérieure avec
> l'identité puis j'extrait la sous matrice significative). Si quelqu'un 
> connait
> l'algorithme de choix, je serais heureux de poursuivre ou de lui expliquer 
> mon
> code si besoin.

For thresholds between different algorithms, I currently encode them
in the type for the "variant of the implementation"; see for instance
the implementation of polynomials.

In the future, we will implement a special variant which heuristically
chooses a good threshold as a function of the actual entries (or type
of the coefficients). This threshold may depend on the architecture.

Other note: it would be good to use a pair/odd recurrence instead of
the standard block decomposition: this has a better numerical stability.

Best wishes, Joris