* dans les ai-reml-logit-* : on peut optimiser un peu en évitant de stocker W(n,n)
  -> utiliser un vecteur W(n) avec W.asDiagonal() ...

* regarder dans score_lmm.cpp si il y a des optimisations (diagonalisation...)

* dans m4_stats.cpp on passe trois vecteurs de Logical pour les chr... -> passer directement le numero de chr...?

** si on utilise RcppArmadillo
   il faut installer sous linux le package libomp-dev pour que ça compile
   
   
