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

procedure OUT ( var sa : real ; lda, n , m : integer ) ;
     {                                                     }
     { General purpose output routine for the matrix A     }
     { which has leading dimension lda and is n by m.      }
     {                                                     }
     { J.J. Dongarra, SICE, ... (?)                        }
     { Adam Fritz, TURBO Pascal, 2/22/85                   }
     {                                                     }

var
     i, j, k        : integer ;
     ic, icb, ice   : integer ;
     a              : RowPointer ;

begin
   if n > 0 then begin
      a := Ptr(Addr(sa)) ;
      if m > 1 then begin
         ic :=  (m + 4) div 5 ;
         icb := 1 ;
         ice := 5 ;
         for k := 1 to ic do begin
            if ice > m then
               ice := m ;
            for i := 1 to n do begin
               for j := icb to ice do
                  write (a^.s[(i-1)*lda+j]:12, '  ') ;
               writeln
            end ;
            icb := icb + 5 ;
            ice := ice + 5 ;
            writeln
         end
      end
      else begin
         ic := (n + 4) div 5 ;
         icb := 1 ;
         ice := 5 ;
         for k := 1 to ic do begin
            if ice > n then
               ice := n ;
            for i := icb to ice do
               write (a^.s[i]:12, '  ') ;
            writeln ;
            icb := icb + 5 ;
            ice := ice + 5
         end ;
         writeln
      end
   end
end ;

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