code : Fortran

Bonjour tout le monde ;
Je me suis interesse par l'algorithme GCD : coordinate descent ,et j'ai un travail a faire ,j'ai un article :The cluster elastic net

dans l'algorithme ,il ya des iterations a faire jusqu'au convergence .Il ya une partie qui est faisable car , on a un code

source (package gcdnet ),mais il est ya un probleme aussi pour ajouter les correlations dans le code

et l'autre probleme est la dexieme partie qui est a chaque fois ,il faut choisir un groupe pour minimiser les betas dans ce dernier

et ainsi de suite ,until convergence .

pour l'algorithme ,il se trouve sur le net :pour la partie b de l'algorithme ,la solution est donnee par la formule 7 dans l'article ,l'initialisation peut

importe .Tout ce qui m'interesse est comment combiner les deux parties a et b .Merci.