bridge/class-master/depend/pari/src/basemath/alglin1.c
/*******************************************************************/
/*                                                                 */
/*                        EIGENVECTORS                             */
/*   (independent eigenvectors, sorted by increasing eigenvalue)   */
/*                                                                 */
/*******************************************************************/
/* assume x is square of dimension > 0 */
EIGENVECTORS
(independent eigenvectors, sorted by increasing eigenvalue)

assume x is square of dimension > 0
EIGENVECTORS ( independent eigenvectors sorted by increasing eigenvalue ) assume x is square of dimension > 0
static int
RgM_is_symmetric_cx(GEN x, long bit)
{
  pari_sp av = avma;
  long i, j, l = lg(x);
  for (i = 1; i < l; i++)
    for (j = 1; j < i; j++)
    {
      GEN a = gcoeff(x,i,j), b = gcoeff(x,j,i), c = gsub(a,b);
      if (!gequal0(c) && gexpo(c) - gexpo(a) > -bit) { avma = av; return 0; }
    }
  avma = av; return 1;
}
static int RgM_is_symmetric_cx ( GEN x , long bit ) { pari_sp av = avma ; long i , j , l = lg ( x ) ; for ( i = 1 ; i < l ; i ++ ) for ( j = 1 ; j < i ; j ++ ) { GEN a = gcoeff ( x , i , j ) , b = gcoeff ( x , j , i ) , c = gsub ( a , b ) ; if ( ! gequal0 ( c ) && gexpo ( c ) - gexpo ( a ) > - bit ) { avma = av ; return 0 ; } } avma = av ; return 1 ; }
EIGENVECTORS
(independent eigenvectors, sorted by increasing eigenvalue)
EIGENVECTORS ( independent eigenvectors sorted by increasing eigenvalue )