/* Copyright (C) 1985,86 Adam Fritz, 133 Main St., Afton, NY 13730 */

OUT ( A, lda, n , m )
   float A[] ;
   int lda, n, m ;

     /* General purpose output routine for the matrix A     */
     /* which has leading dimension lda and is n by m.      */
     /*                                                     */
     /* J.J. Dongarra, SICE, ... (?)                        */
     /* Adam Fritz, C, 2/22/85                              */

{
   int          i, j, k ;
   int          ic, icb, ice ;

   if ( n > 0 )
   {
      if ( m > 1 )
      {
         ic =  (m+4)/5 ;
         icb = 0 ;
         ice = 5 ;
         for ( k = 0; k < ic; k++ )
         {
            if (ice > m )
               ice = m ;
            for ( i = 0; i < n; i++ )
            {
               for ( j = icb; j < ice; j++ )
                  printf("%5e  ", A[i*lda+j]) ;
               printf("\n") ;
            } ;
            icb = icb + 5 ;
            ice = ice + 5 ;
            printf("\n") ;
         }
      }
      else {
         ic = (n+4)/5 ;
         icb = 0 ;
         ice = 5 ;
         for ( k = 0; k < ic; k++ )
         {
            if (ice > n )
               ice = n ;
            for ( i = icb; i < ice; i++ )
               printf("%5e  ", A[i]) ;
            printf("\n") ;
            icb = icb + 5 ;
            ice = ice + 5 ;
         } ;
         printf("\n") ;
      }
   }
}

/* Copyright (C) 1985,86 Adam Fritz, 133 Main St., Afton, NY 13730 */