Annotation of OpenXM_contrib/PHC/Ada/Math_Lib/Matrices/ts_intmat.adb, Revision 1.1
1.1 ! maekawa 1: with text_io,integer_io; use text_io,integer_io;
! 2: with Standard_Integer_Vectors;
! 3: with Standard_Integer_Vectors_io;
! 4: with Standard_Integer_Matrices;
! 5: with Standard_Integer_Matrices_io;
! 6: with Standard_Integer_VecMats;
! 7: with Standard_Integer_VecMats_io;
! 8: with Standard_Integer_Linear_Solvers;
! 9: with Standard_Random_Matrices; use Standard_Random_Matrices;
! 10: with Multprec_Random_Matrices; use Multprec_Random_Matrices;
! 11: with Multprec_Integer_Vectors;
! 12: with Multprec_Integer_Vectors_io;
! 13: with Multprec_Integer_Matrices;
! 14: with Multprec_Integer_Matrices_io;
! 15: with Multprec_Integer_Linear_Solvers;
! 16:
! 17: procedure ts_intmat is
! 18:
! 19: -- DESCRIPTION :
! 20: -- Tests the matrix packages of standard and multi-precision integers.
! 21:
! 22: procedure Test_Standard_io is
! 23:
! 24: use Standard_Integer_Matrices,Standard_Integer_Matrices_io;
! 25:
! 26: n,m : natural;
! 27:
! 28: begin
! 29: put("Give the number of rows : "); get(n);
! 30: put("Give the number of columns : "); get(m);
! 31: declare
! 32: mat : Matrix(1..n,1..m);
! 33: begin
! 34: put("Give "); put(n,1); put("x"); put(m,1);
! 35: put_line(" integer matrix : "); get(mat);
! 36: put_line("Your matrix : "); put(mat);
! 37: end;
! 38: end Test_Standard_io;
! 39:
! 40: procedure Test_Standard_VecMat_io is
! 41:
! 42: use Standard_Integer_Matrices,Standard_Integer_Matrices_io;
! 43: use Standard_Integer_VecMats,Standard_Integer_VecMats_io;
! 44:
! 45: n,n1,n2 : natural;
! 46: lv : Link_to_VecMat;
! 47:
! 48: begin
! 49: put("Give the number of matrices : "); get(n);
! 50: put("Give #rows : "); get(n1);
! 51: put("Give #columns : "); get(n2);
! 52: put("Give "); put(n,1); put(" "); put(n1,1); put("-by-"); put(n2,1);
! 53: put_line(" integer matrices : ");
! 54: get(n,n1,n2,lv);
! 55: put_line("The vector of matrices :"); put(lv);
! 56: end Test_Standard_VecMat_io;
! 57:
! 58: procedure Random_Test_Standard_Solver
! 59: ( n,m : in natural; low,upp : in integer ) is
! 60:
! 61: use Standard_Integer_Vectors,Standard_Integer_Vectors_io;
! 62: use Standard_Integer_Matrices;
! 63: use Standard_Integer_Linear_Solvers;
! 64:
! 65: mat : Matrix(1..n,1..m) := Random_Matrix(n,m,low,upp);
! 66: wrk : Matrix(1..n,1..m) := mat;
! 67: sol : Vector(1..m);
! 68: res : Vector(1..n);
! 69:
! 70: begin
! 71: Upper_Triangulate(wrk);
! 72: Solve0(wrk,sol);
! 73: res := mat*sol;
! 74: put("The residual : "); put(res);
! 75: put(" of solution : "); put(sol); new_line;
! 76: end Random_Test_Standard_Solver;
! 77:
! 78: procedure Random_Test_Multprec_Solver ( n,m,sz : in natural ) is
! 79:
! 80: use Multprec_Integer_Vectors,Multprec_Integer_Vectors_io;
! 81: use Multprec_Integer_Matrices,Multprec_Integer_Matrices_io;
! 82: use Multprec_Integer_Linear_Solvers;
! 83:
! 84: mat : Matrix(1..n,1..m) := Random_Matrix(n,m,sz);
! 85: wrk : Matrix(1..n,1..m);
! 86: sol : Vector(1..m);
! 87: res : Vector(1..n);
! 88:
! 89: begin
! 90: Copy(mat,wrk);
! 91: Upper_Triangulate(wrk);
! 92: Solve0(wrk,sol);
! 93: res := mat*sol;
! 94: put("The residual : "); put(res);
! 95: put(" of solution : "); put(sol); new_line;
! 96: end Random_Test_Multprec_Solver;
! 97:
! 98: procedure Random_Test_Standard_Solvers is
! 99:
! 100: n,m,nb : natural;
! 101: low,upp : integer;
! 102:
! 103: begin
! 104: put("Give number of tests : "); get(nb);
! 105: put("Give number of rows : "); get(n);
! 106: put("Give number of columns : "); get(m);
! 107: put("Give lower bound for numbers : "); get(low);
! 108: put("Give upper bound for numbers : "); get(upp);
! 109: for i in 1..nb loop
! 110: Random_Test_Standard_Solver(n,m,low,upp);
! 111: end loop;
! 112: end Random_Test_Standard_Solvers;
! 113:
! 114: procedure Random_Test_Multprec_Solvers is
! 115:
! 116: n,m,nb,sz : natural;
! 117:
! 118: begin
! 119: put("Give number of tests : "); get(nb);
! 120: put("Give number of rows : "); get(n);
! 121: put("Give number of columns : "); get(m);
! 122: put("Give the size of the numbers : "); get(sz);
! 123: for i in 1..nb loop
! 124: Random_Test_Multprec_Solver(n,m,sz);
! 125: end loop;
! 126: end Random_Test_Multprec_Solvers;
! 127:
! 128: procedure Interactive_Test_Standard_Solvers is
! 129:
! 130: use Standard_Integer_Vectors,Standard_Integer_Vectors_io;
! 131: use Standard_Integer_Matrices,Standard_Integer_Matrices_io;
! 132: use Standard_Integer_Linear_Solvers;
! 133:
! 134: n,m : natural;
! 135:
! 136: begin
! 137: put("Give the number of rows : "); get(n);
! 138: put("Give the number of columns : "); get(m);
! 139: declare
! 140: sol : Vector(1..m) := (1..m => 0);
! 141: res : Vector(1..n);
! 142: mat,wrk : Matrix(1..n,1..m);
! 143: l : Matrix(1..n,1..n);
! 144: begin
! 145: put("Give "); put(n,1); put("x"); put(m,1);
! 146: put_line(" integer matrix : "); get(mat);
! 147: put_line("Your matrix : "); put(mat);
! 148: wrk := mat;
! 149: Upper_Triangulate(l,wrk);
! 150: put_line("The matrix in upper triangular form : "); put(wrk);
! 151: put_line("The transformation matrix T : "); put(l);
! 152: put_line("product of T and original matrix : "); put(l*mat);
! 153: Solve0(wrk,sol);
! 154: put_line("The solution of the homogeneous system : ");
! 155: put(sol); new_line;
! 156: res := mat*sol;
! 157: put("The residual : "); put(res); new_line;
! 158: end;
! 159: end Interactive_Test_Standard_Solvers;
! 160:
! 161: procedure Test_Multprec_io is
! 162:
! 163: use Multprec_Integer_Matrices,Multprec_Integer_Matrices_io;
! 164:
! 165: n,m : natural;
! 166:
! 167: begin
! 168: put("Give the number of rows : "); get(n);
! 169: put("Give the number of columns : "); get(m);
! 170: declare
! 171: mat : Matrix(1..n,1..m);
! 172: begin
! 173: put("Give "); put(n,1); put("x"); put(m,1);
! 174: put_line(" integer matrix : "); get(mat);
! 175: put_line("Your matrix : "); put(mat); new_line;
! 176: end;
! 177: end Test_Multprec_io;
! 178:
! 179: procedure Interactive_Test_Multprec_Solvers is
! 180:
! 181: use Multprec_Integer_Vectors,Multprec_Integer_Vectors_io;
! 182: use Multprec_Integer_Matrices,Multprec_Integer_Matrices_io;
! 183: use Multprec_Integer_Linear_Solvers;
! 184:
! 185: n,m : natural;
! 186:
! 187: begin
! 188: put("Give the number of rows : "); get(n);
! 189: put("Give the number of columns : "); get(m);
! 190: declare
! 191: sol : Vector(1..m);
! 192: res : Vector(1..n);
! 193: mat,wrk : Matrix(1..n,1..m);
! 194: l : Matrix(1..n,1..n);
! 195: begin
! 196: put("Give "); put(n,1); put("x"); put(m,1);
! 197: put_line(" integer matrix : "); get(mat);
! 198: put_line("Your matrix : "); put(mat);
! 199: Copy(mat,wrk);
! 200: Upper_Triangulate(l,wrk);
! 201: put_line("The matrix in upper triangular form : "); put(wrk);
! 202: put_line("The transformation matrix T : "); put(l);
! 203: put_line("product of T and original matrix : "); put(l*mat);
! 204: Solve0(wrk,sol);
! 205: put_line("The solution of the homogeneous system : ");
! 206: put(sol); new_line;
! 207: res := mat*sol;
! 208: put("The residual : "); put(res); new_line;
! 209: end;
! 210: end Interactive_Test_Multprec_Solvers;
! 211:
! 212: procedure Test_Multprec_Matrix_Vector_Product is
! 213:
! 214: use Multprec_Integer_Vectors,Multprec_Integer_Vectors_io;
! 215: use Multprec_Integer_Matrices,Multprec_Integer_Matrices_io;
! 216:
! 217: n,m : natural;
! 218:
! 219: begin
! 220: put("Give the number of rows : "); get(n);
! 221: put("Give the number of columns : "); get(m);
! 222: declare
! 223: mat : Matrix(1..n,1..m);
! 224: vec : Vector(1..m);
! 225: prod : Vector(1..n);
! 226: begin
! 227: put("Give "); put(n,1); put("x"); put(m,1);
! 228: put_line(" integer matrix : "); get(mat);
! 229: put_line("Your matrix : "); put(mat);
! 230: put("Give an " ); put(m,1); put("-vector : "); get(vec);
! 231: put("Your vector : "); put(vec); new_line;
! 232: prod := mat*vec;
! 233: put_line("The matrix-vector product : "); put(prod); new_line;
! 234: end;
! 235: end Test_Multprec_Matrix_Vector_Product;
! 236:
! 237: procedure Main is
! 238:
! 239: ans : character;
! 240:
! 241: begin
! 242: new_line;
! 243: put_line("Interactive testing of matrices of integer numbers");
! 244: new_line;
! 245: loop
! 246: put_line("Choose one of the following : ");
! 247: put_line(" 0. exit this program.");
! 248: put_line(" 1. io of matrices of standard numbers.");
! 249: put_line(" 2. io of vectors of matrices of standard numbers.");
! 250: put_line(" 3. solve given linear systems of standard numbers.");
! 251: put_line(" 4. solve random linear systems of standard numbers.");
! 252: put_line(" 5. io of matrices of multi-precision numbers.");
! 253: put_line(" 6. solve given linear systems of multi-precision numbers.");
! 254: put_line(" 7. solve random linear systems of multi-precision numbers.");
! 255: put_line(" 8. test multi-precision matrix-vector product.");
! 256: put("Make your choice (0,1,2,3,4,5,6,7 or 8) : "); get(ans);
! 257: exit when (ans = '0');
! 258: case ans is
! 259: when '1' => Test_Standard_io;
! 260: when '2' => Test_Standard_VecMat_io;
! 261: when '3' => Interactive_Test_Standard_Solvers;
! 262: when '4' => Random_Test_Standard_Solvers;
! 263: when '5' => Test_Multprec_io;
! 264: when '6' => Interactive_Test_Multprec_Solvers;
! 265: when '7' => Random_Test_Multprec_Solvers;
! 266: when '8' => Test_Multprec_Matrix_Vector_Product;
! 267: when others => null;
! 268: end case;
! 269: end loop;
! 270: end Main;
! 271:
! 272: begin
! 273: Main;
! 274: end ts_intmat;
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to ts_intmat.adb CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / PHC / Ada / Math_Lib / Matrices