with text_io,integer_io;                 use text_io,integer_io;
with Standard_Floating_Matrices;         use Standard_Floating_Matrices;
with Standard_Floating_Matrices_io;      use Standard_Floating_Matrices_io;
with Standard_Complex_Matrices;          use Standard_Complex_Matrices;
with Standard_Complex_Matrices_io;       use Standard_Complex_Matrices_io;
with Standard_Random_Matrices;           use Standard_Random_Matrices;
with Standard_Matrix_Inversion;          use Standard_Matrix_Inversion;

procedure ts_matinv is

-- DESCRIPTION :
--   Test on matrix inversion.

  procedure Main is
    
    n : natural;

  begin
    put("Give the dimension : "); get(n);
    declare
      cmpmat : Standard_Complex_Matrices.Matrix(1..n,1..n)
             := Random_Matrix(n,n);
      invcmp : Standard_Complex_Matrices.Matrix(1..n,1..n) := Inverse(cmpmat);
      fltmat : Standard_Floating_Matrices.Matrix(1..n,1..n)
             := Random_Matrix(n,n);
      invflt : Standard_Floating_Matrices.Matrix(1..n,1..n) := Inverse(fltmat);
    begin
      put_line("Test on complex matrix inversion : ");
      put(cmpmat*invcmp,3);
      put_line("Test on commutativity with complex matrix inversion : ");
      put(invcmp*cmpmat,3);
      put_line("Test on floating matrix inversion : ");
      put(fltmat*invflt,3);
      put_line("Test on commutativity with floating matrix inversion : ");
      put(invflt*fltmat,3);
    end;
  end Main;

begin 
  Main;
end ts_matinv;