Return to ts_matinv.adb CVS log

Up to [local] / OpenXM_contrib / PHC / Ada / Schubert

Annotation of OpenXM_contrib/PHC/Ada/Schubert/ts_matinv.adb, Revision 1.1

1.1     ! maekawa     1: with text_io,integer_io;                 use text_io,integer_io;
        !             2: with Standard_Floating_Matrices;         use Standard_Floating_Matrices;
        !             3: with Standard_Floating_Matrices_io;      use Standard_Floating_Matrices_io;
        !             4: with Standard_Complex_Matrices;          use Standard_Complex_Matrices;
        !             5: with Standard_Complex_Matrices_io;       use Standard_Complex_Matrices_io;
        !             6: with Standard_Random_Matrices;           use Standard_Random_Matrices;
        !             7: with Standard_Matrix_Inversion;          use Standard_Matrix_Inversion;
        !             8:
        !             9: procedure ts_matinv is
        !            10:
        !            11: -- DESCRIPTION :
        !            12: --   Test on matrix inversion.
        !            13:
        !            14:   procedure Main is
        !            15:
        !            16:     n : natural;
        !            17:
        !            18:   begin
        !            19:     put("Give the dimension : "); get(n);
        !            20:     declare
        !            21:       cmpmat : Standard_Complex_Matrices.Matrix(1..n,1..n)
        !            22:              := Random_Matrix(n,n);
        !            23:       invcmp : Standard_Complex_Matrices.Matrix(1..n,1..n) := Inverse(cmpmat);
        !            24:       fltmat : Standard_Floating_Matrices.Matrix(1..n,1..n)
        !            25:              := Random_Matrix(n,n);
        !            26:       invflt : Standard_Floating_Matrices.Matrix(1..n,1..n) := Inverse(fltmat);
        !            27:     begin
        !            28:       put_line("Test on complex matrix inversion : ");
        !            29:       put(cmpmat*invcmp,3);
        !            30:       put_line("Test on commutativity with complex matrix inversion : ");
        !            31:       put(invcmp*cmpmat,3);
        !            32:       put_line("Test on floating matrix inversion : ");
        !            33:       put(fltmat*invflt,3);
        !            34:       put_line("Test on commutativity with floating matrix inversion : ");
        !            35:       put(invflt*fltmat,3);
        !            36:     end;
        !            37:   end Main;
        !            38:
        !            39: begin
        !            40:   Main;
        !            41: end ts_matinv;