Annotation of OpenXM_contrib/PHC/Ada/Main/tableau_formats.adb, Revision 1.1
1.1 ! maekawa 1: with integer_io; use integer_io;
! 2: with Communications_with_User; use Communications_with_User;
! 3: with Standard_Floating_Numbers; use Standard_Floating_Numbers;
! 4: with Standard_Floating_Numbers_io; use Standard_Floating_Numbers_io;
! 5: with Standard_Complex_Numbers; use Standard_Complex_Numbers;
! 6: with Standard_Complex_Numbers_io; use Standard_Complex_Numbers_io;
! 7: with Standard_Natural_Vectors;
! 8: with Standard_Integer_Vectors;
! 9: with Standard_Complex_Vectors;
! 10: with Standard_Complex_VecVecs; use Standard_Complex_VecVecs;
! 11: with Lists_of_Integer_Vectors; use Lists_of_Integer_Vectors;
! 12: with Lists_of_Integer_Vectors_io; use Lists_of_Integer_Vectors_io;
! 13: with Arrays_of_Integer_Vector_Lists; use Arrays_of_Integer_Vector_Lists;
! 14: with Arrays_of_Integer_Vector_Lists_io; use Arrays_of_Integer_Vector_Lists_io;
! 15: with Symbol_Table,Symbol_Table_io; use Symbol_Table;
! 16: with Standard_Complex_Polynomials; use Standard_Complex_Polynomials;
! 17: with Standard_Complex_Poly_Systems_io; use Standard_Complex_Poly_Systems_io;
! 18: with Power_Lists; use Power_Lists;
! 19:
! 20: package body Tableau_Formats is
! 21:
! 22: -- AUXILIARIES :
! 23:
! 24: procedure Read_Symbols ( file : in file_type; n : in natural ) is
! 25:
! 26: -- DESCRIPTION :
! 27: -- Reads n symbols from file.
! 28:
! 29: begin
! 30: Symbol_Table.Init(n);
! 31: for i in 1..n loop
! 32: declare
! 33: s : Symbol;
! 34: begin
! 35: Symbol_Table_io.get(file,s);
! 36: Symbol_Table.Add(s);
! 37: end;
! 38: end loop;
! 39: end Read_Symbols;
! 40:
! 41: procedure Write_Symbols ( file : in file_type ) is
! 42:
! 43: -- DESCRIPTION :
! 44: -- Writes the content of the symbol table on file.
! 45:
! 46: n : constant natural := Symbol_Table.Number;
! 47:
! 48: begin
! 49: put(file,n,1); new_line(file);
! 50: for i in 1..n loop
! 51: declare
! 52: s : Symbol := Symbol_Table.Get(i);
! 53: begin
! 54: put(file," ");
! 55: Symbol_Table_io.put(file,s);
! 56: end;
! 57: end loop;
! 58: new_line(file);
! 59: end Write_Symbols;
! 60:
! 61: procedure Read_Exponents ( file : in file_type;
! 62: exp : in out Array_of_Lists ) is
! 63:
! 64: -- DESCRIPTION :
! 65: -- Reads exponent vectors from file.
! 66: -- Each list of exponent vectors should be preceded by the number
! 67: -- of elements in the list.
! 68:
! 69: n : constant natural := exp'length;
! 70: m : natural;
! 71:
! 72: begin
! 73: for i in exp'range loop
! 74: get(file,m);
! 75: get(file,n,m,exp(i));
! 76: end loop;
! 77: end Read_Exponents;
! 78:
! 79: procedure Write_Exponents ( file : in file_type; p : in Poly_Sys ) is
! 80:
! 81: -- DESCRIPTION :
! 82: -- Writes the supports of p on file.
! 83:
! 84: exp : Array_of_Lists(p'range) := Create(p);
! 85:
! 86: begin
! 87: for i in exp'range loop
! 88: put(file,Length_Of(exp(i)),1); new_line(file);
! 89: put(file,exp(i));
! 90: end loop;
! 91: Deep_Clear(exp);
! 92: end Write_Exponents;
! 93:
! 94: procedure Read_Coefficients
! 95: ( file : in file_type;
! 96: exp : in Array_of_Lists; flt : in boolean;
! 97: cff : in out VecVec ) is
! 98:
! 99: -- DESCRIPTION :
! 100: -- Reads the coefficient from file, given the supports.
! 101:
! 102: f : double_float;
! 103:
! 104: begin
! 105: for i in cff'range loop
! 106: cff(i) := new Standard_Complex_Vectors.Vector(1..Length_Of(exp(i)));
! 107: for j in cff(i)'range loop
! 108: if flt
! 109: then get(file,f); cff(i)(j) := Create(f);
! 110: else get(file,cff(i)(j));
! 111: end if;
! 112: end loop;
! 113: end loop;
! 114: end Read_Coefficients;
! 115:
! 116: procedure Write_Coefficients ( file : in file_type; flt : in boolean;
! 117: p : in Poly ) is
! 118:
! 119: -- DESCRIPTION :
! 120: -- Writes the coefficients of the polynomial on file.
! 121: -- If flt, then the imaginary parts will be omitted.
! 122:
! 123: procedure Write_Coefficient_of_Term ( t : in Term; cont : out boolean ) is
! 124: begin
! 125: if flt
! 126: then put(file,REAL_PART(t.cf));
! 127: else put(file,t.cf);
! 128: end if;
! 129: new_line(file);
! 130: cont := true;
! 131: end Write_Coefficient_of_Term;
! 132:
! 133: procedure Write_Coefficients_of_Terms is
! 134: new Visiting_Iterator(Write_Coefficient_of_Term);
! 135:
! 136: begin
! 137: Write_Coefficients_of_Terms(p);
! 138: end Write_Coefficients;
! 139:
! 140: procedure Write_Coefficients ( file : in file_type; flt : in boolean;
! 141: p : in Poly_Sys ) is
! 142:
! 143: -- DESCRIPTION :
! 144: -- Writes the coefficients of the polynomial system on file.
! 145: -- If flt, then the imaginary parts will be omitted.
! 146:
! 147: begin
! 148: for i in p'range loop
! 149: Write_Coefficients(file,flt,p(i));
! 150: end loop;
! 151: end Write_Coefficients;
! 152:
! 153: function Create ( exp : Standard_Integer_Vectors.Vector;
! 154: cff : Complex_Number ) return Term is
! 155:
! 156: -- DESCRIPTION :
! 157: -- Creates a term from exponent vector and coefficient.
! 158:
! 159: res : Term;
! 160:
! 161: begin
! 162: res.cf := cff;
! 163: res.dg := new Standard_Natural_Vectors.Vector(exp'range);
! 164: for i in exp'range loop
! 165: res.dg(i) := exp(i);
! 166: end loop;
! 167: return res;
! 168: end Create;
! 169:
! 170: function Create ( exp : List; cff : Standard_Complex_Vectors.Vector )
! 171: return Poly is
! 172:
! 173: -- DESCRIPTION :
! 174: -- Creates a polynomial from the exponents and the coefficients.
! 175:
! 176: res : Poly;
! 177: tmp : List := exp;
! 178: lpt : Standard_Integer_Vectors.Link_to_Vector;
! 179:
! 180: begin
! 181: for i in cff'range loop
! 182: lpt := Head_Of(tmp);
! 183: declare
! 184: t : Term := Create(lpt.all,cff(i));
! 185: begin
! 186: Add(res,t);
! 187: end;
! 188: tmp := Tail_Of(tmp);
! 189: end loop;
! 190: return res;
! 191: end Create;
! 192:
! 193: function Create ( exp : Array_of_Lists; cff : VecVec ) return Poly_Sys is
! 194:
! 195: -- DESCRIPTION :
! 196: -- Creates the polynomial system from the supports and coefficients.
! 197:
! 198: res : Poly_Sys(exp'range);
! 199:
! 200: begin
! 201: for i in res'range loop
! 202: res(i) := Create(exp(i),cff(i).all);
! 203: end loop;
! 204: return res;
! 205: end Create;
! 206:
! 207: -- TARGET ROUTINES :
! 208:
! 209: procedure get ( file : in file_type; realcoeff : in boolean;
! 210: p : out Link_to_Poly_Sys ) is
! 211:
! 212: n : natural;
! 213:
! 214: begin
! 215: get(file,n);
! 216: new_line;
! 217: put("Dimension : "); put(n,1); new_line;
! 218: Read_Symbols(file,n);
! 219: put("Symbols :");
! 220: for i in 1..n loop
! 221: put(" "); put(Symbol_Table.get(i));
! 222: end loop;
! 223: new_line;
! 224: declare
! 225: exp : Array_of_Lists(1..n);
! 226: cff : VecVec(1..n);
! 227: sys : Poly_Sys(1..n);
! 228: begin
! 229: Read_Exponents(file,exp);
! 230: put_line("The exponents : ");
! 231: put(exp);
! 232: Read_Coefficients(file,exp,realcoeff,cff);
! 233: put_line("The coefficients : ");
! 234: for i in cff'range loop
! 235: for j in cff(i)'range loop
! 236: put(cff(i)(j)); new_line;
! 237: end loop;
! 238: new_line;
! 239: end loop;
! 240: sys := Create(exp,cff);
! 241: put_line("The polynomial system : "); put(sys);
! 242: p := new Poly_Sys'(sys);
! 243: end;
! 244: end get;
! 245:
! 246: procedure put ( file : in file_type; realcoeff : in boolean;
! 247: p : in Poly_Sys ) is
! 248:
! 249: begin
! 250: if not Symbol_Table.Empty
! 251: then Write_Symbols(file);
! 252: end if;
! 253: Write_Exponents(file,p);
! 254: Write_Coefficients(file,realcoeff,p);
! 255: end put;
! 256:
! 257: end Tableau_Formats;
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to tableau_formats.adb CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / PHC / Ada / Main