Skip to content

[ODK] Meeting 2019 03 06

A. Breust edited this page Mar 7, 2019 · 1 revision

ODK Linbox meeting 2019-03-06 CR

Alexis:

refacto solve:

  • solveInPlace. Ok pour sparseElim. Les autres se rabattent sur solve classique
  • remplacement des MethodsWIP par des Methods en cours:
    • fini pour solve
    • ca casse sur minpoly, det, etc -> fix en cours
    • test-suite à finaliser Pb: segfault sur des solve
  • NumSymOverlap
  • BlockLanczos
  • Wiedemann avec BlasMatrix -> ajouter un commentaire devant la définition de la méthode pour dire que ça ne marche pas -> ajouter une issue pour chaque pb

Dixon refacto:

  • commentaires explicatifs du monolithic solve (rational-solver.inl) -> nettoyage à faire (découper en morceaux, factoriser le codes, tester
  • lifting container est bien factorisée
  • mais le rational solver fait des duplications inutiles pour chaque cas (sparse, BB, dense)
  • sortir le certif d'inconsistance
  • faire une passe de nettoyage/ mise à jour par rapport aux fonctionnalités dans ffpack. A garder en tete pour le design:
  • possibilité d'avoir une famille de vecteurs inconnus (matrix)
  • hybrid CRT-Dixon -> matrix modulo différents premiers pour chaque colonne.

Planning: merger la PR solveCRAMPI >> fixer solve refacto pour ouvrir une PR >> merger cette PR

MDRvector: supprimé AB a trouvé un autre type (VectorFraction) de donné équivalent en interne -> on les laisse pour l'instant. Ajout d'une surcharge de solve avec le type VectorFraction (dans algorithms/vector-fraction.h)

Zhu:

restructurer le benchmark-fgemv pour y introduire le parallélisme paladin -> done -> data de bench dans wiki fflas-ffpack -> refaire avec q=0 -> changer les field par défaut, enlever static, etc

tutoriel pour solveCRAMPI

  • ajouter comment lancer l'exemple benchmark-solveCRAMPI
  • ajouter l'option --enable-mpi

Bug ModularBalanced incomplet: n'a jamais existé: en attendant qu'il soit fait (refactor ModularBalanced) utiliser Modular

Notes diverses

Abandon:

  • pb du générateur heuristique: ne détecte pas qu'un nb tiré a déjà été pris car pas de progress. -> implementer un order_set pour vérifier si on l'a déjà tiré.
  • OMP: completement passé en taches paladin (branche solveCRA-paladin)

finir SV1-SV2 et SV3 minimaliste (JGD)

Comportement en cas de sol impossible -> exception

Voir avec Pascal sur le statut apply.h

Coding sprints les vendredi matins

TODOS: SG1: HZ + CP fgemm, PLUQ, etc avec paladin exposés dans sage: on commence vendredi 8/03 L1: AB + JGD D1: JGD+CP+AB+HZ plus tard

Clone this wiki locally