Aller au contenu

Mathc complexes/h10aa

Un livre de Wikilivres.


Bibliothèque


Installer ce fichier dans votre répertoire de travail.

wgpivot3.h
/* ------------------------------------ */
/*  Save as :   wgpivot3.h              */
/* ------------------------------------ */
/* ------------------------------------
   pivotbest = |a+bI|^2 =  = a^2+b^2
   ------------------------------------ */
double pivotbest_gj3Ab_mZ(
double **Ab,
int pivot_r,
int pivot_c
)
{
double pivotbest=Ab[pivot_r][pivot_c]   *Ab[pivot_r][pivot_c]
                                     +
                 Ab[pivot_r][pivot_c+C1]*Ab[pivot_r][pivot_c+C1];

double pivot;

int best_r = pivot_r;
int r;

  for( r=pivot_r; r<Ab[R_SIZE][C0]; r++)
         {
          pivot = Ab[r][pivot_c]*Ab[r][pivot_c]    + 
                  Ab[r][pivot_c+C1]*Ab[r][pivot_c+C1];
          
           if(pivot>pivotbest)
             {
               pivotbest=pivot;
               best_r = r;
             }
         }

 if(best_r!=pivot_r) swapR_mZ(Ab,pivot_r,best_r);

 return(pivotbest);
}
/* ------------------------------------ */
/* ------------------------------------ */


Déclaration des fichiers h.