Exemple de fichier d'entrée :
32 Nb points entree (TRAN) 8 U 0. mu (sans le facteur U/2) 1000 nbloc(important si lecture diagonale) 20 C dimmension 1 B dimmension 1 KX Vecteur d'onde 1 KY Vecteur d'onde 301 Nb points sortie (nw) -15 wmin (0. pour les bosons) 15 wmax 30 Nb de valeurs de alpha (nalpha) 0.1 Alpha1 1.4 Alpha2 8 Nb d'iterations (niter) 0.25 Ski 4 Nb options (nopt) y Matrice de covariance n Densite d'etat n Lecture diagonale n Imposer la symetrie de G |
Notes :
Nous effectuons ensuite un développement au deuxième ordre
de la fonction Q, où nous vérifions facilement que
le terme quadratique est :
(1)
que nous pouvons diagonaliser sous la forme Q=Uf
UT. Ce changement de base est appliqué au poids spectrale,
de sorte que
(2)
Ce changement de base nous permet d'écrire . Nous pouvons donc minimiser la norme du vecteur y avec certaines contraintes de positivité.