Mathc complexes/h10f

Bibliothèque

Installer ce fichier dans votre répertoire de travail.

	wgpivot4.h

/* ------------------------------------ */
/*  Save as :   wgpivot4.h              */
/* ------------------------------------ */
/* ------------------------------------
   pivotbest = |a+bI|^2 =  = a^2+b^2
   ------------------------------------ */
double pivotbest4Ab_r0_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);
                        swapcoef_R0_mZ(Ab,(((pivot_r-C1)*C2)+C1),
                                             (((best_r -C1)*C2)+C1) );
                     }

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

Déclaration des fichiers h.