Annotation of OpenXM_contrib/PHC/Ada/Root_Counts/Stalift/integer_lifting_functions.adb, Revision 1.1
1.1 ! maekawa 1: with Standard_Floating_Numbers; use Standard_Floating_Numbers;
! 2: with Standard_Complex_Numbers; use Standard_Complex_Numbers;
! 3: with Standard_Complex_Vectors;
! 4: with Standard_Random_Numbers;
! 5:
! 6: package body Integer_Lifting_Functions is
! 7:
! 8: -- AUXILIARIES :
! 9:
! 10: function Convert ( v : Standard_Integer_Vectors.Vector )
! 11: return Standard_Complex_Vectors.Vector is
! 12:
! 13: -- DESCRIPTION :
! 14: -- Converts the given vector into a vector with complex entries.
! 15:
! 16: res : Standard_Complex_Vectors.Vector(v'range);
! 17:
! 18: begin
! 19: for i in v'range loop
! 20: res(i) := Create(double_float(v(i)));
! 21: end loop;
! 22: return res;
! 23: end Convert;
! 24:
! 25: function Random_Vector ( m : natural; low,upp : integer ) return Vector is
! 26:
! 27: -- DESCRIPTION :
! 28: -- Returns a random vector of range 1..m of randomly generated integer
! 29: -- values between low and upp.
! 30:
! 31: res : Vector(1..m);
! 32:
! 33: begin
! 34: for k in res'range loop
! 35: res(k) := Standard_Random_Numbers.Random(low,upp);
! 36: end loop;
! 37: return res;
! 38: end Random_Vector;
! 39:
! 40: function Random ( vec : Standard_Integer_Vectors.Vector ) return integer is
! 41:
! 42: -- DESCRIPTION :
! 43: -- Returns a random number from the given vector.
! 44:
! 45: index : integer := Standard_Random_Numbers.Random(vec'first,vec'last);
! 46:
! 47: begin
! 48: return vec(index);
! 49: end Random;
! 50:
! 51: -- LINEAR LIFTING :
! 52:
! 53: function Linear_Lift ( lf,v : Vector ) return Vector is
! 54:
! 55: res : Vector(v'first..v'last+1);
! 56:
! 57: begin
! 58: res(v'range) := v;
! 59: res(res'last) := lf*v;
! 60: return res;
! 61: end Linear_Lift;
! 62:
! 63: function Linear_Lift ( lf : Vector; l : List ) return List is
! 64:
! 65: res,res_last : List;
! 66: tmp : List := l;
! 67:
! 68: begin
! 69: while not Is_Null(tmp) loop
! 70: Append(res,res_last,Linear_Lift(lf,Head_Of(tmp).all));
! 71: tmp := Tail_Of(tmp);
! 72: end loop;
! 73: return res;
! 74: end Linear_Lift;
! 75:
! 76: function Linear_Lift ( lf : VecVec;
! 77: l : Array_of_Lists ) return Array_of_Lists is
! 78:
! 79: res : Array_of_Lists(l'range);
! 80:
! 81: begin
! 82: for i in res'range loop
! 83: res(i) := Linear_Lift(lf(i).all,l(i));
! 84: end loop;
! 85: return res;
! 86: end Linear_Lift;
! 87:
! 88: -- POLYNOMIAL LIFTING :
! 89:
! 90: function Polynomial_Lift ( lf : Poly; v : Vector ) return Vector is
! 91:
! 92: res : Vector(v'first..v'last+1);
! 93:
! 94: begin
! 95: res(v'range) := v;
! 96: res(res'last) := integer(REAL_PART(Eval(lf,Convert(v))));
! 97: return res;
! 98: end Polynomial_Lift;
! 99:
! 100: function Polynomial_Lift ( lf : Eval_Poly; v : Vector ) return Vector is
! 101:
! 102: res : Vector(v'first..v'last+1);
! 103:
! 104: begin
! 105: res(v'range) := v;
! 106: res(res'last) := integer(REAL_PART(Eval(lf,Convert(v))));
! 107: return res;
! 108: end Polynomial_Lift;
! 109:
! 110: function Polynomial_Lift ( lf : Poly; l : List ) return List is
! 111:
! 112: res : List;
! 113: elf : Eval_Poly := Create(lf);
! 114:
! 115: begin
! 116: res := Polynomial_Lift(elf,l);
! 117: Clear(elf);
! 118: return res;
! 119: end Polynomial_Lift;
! 120:
! 121: function Polynomial_Lift ( lf : Eval_Poly; l : List ) return List is
! 122:
! 123: res,res_last : List;
! 124: tmp : List := l;
! 125:
! 126: begin
! 127: while not Is_Null(tmp) loop
! 128: Append(res,res_last,Polynomial_Lift(lf,Head_Of(tmp).all));
! 129: tmp := Tail_Of(tmp);
! 130: end loop;
! 131: return res;
! 132: end Polynomial_Lift;
! 133:
! 134: function Polynomial_Lift ( lf : Poly_Sys; l : Array_of_Lists )
! 135: return Array_of_Lists is
! 136:
! 137: res : Array_of_Lists(l'range);
! 138: elf : Eval_Poly_Sys(lf'range) := Create(lf);
! 139:
! 140: begin
! 141: res := Polynomial_Lift(elf,l);
! 142: Clear(elf);
! 143: return res;
! 144: end Polynomial_Lift;
! 145:
! 146: function Polynomial_Lift ( lf : Eval_Poly_Sys; l : Array_of_Lists )
! 147: return Array_of_Lists is
! 148:
! 149: res : Array_of_Lists(l'range);
! 150:
! 151: begin
! 152: for i in res'range loop
! 153: res(i) := Polynomial_Lift(lf(i),l(i));
! 154: end loop;
! 155: return res;
! 156: end Polynomial_Lift;
! 157:
! 158: -- RANDOM LIFTING :
! 159:
! 160: function Random_Lift ( lflow,lfupp : integer; v : Vector ) return Vector is
! 161:
! 162: res : Vector(v'first..v'last+1);
! 163:
! 164: begin
! 165: res(v'range) := v;
! 166: res(res'last) := Standard_Random_Numbers.Random(lflow,lfupp);
! 167: return res;
! 168: end Random_Lift;
! 169:
! 170: function Random_Lift ( lflow,lfupp : integer; l : List ) return List is
! 171:
! 172: res,res_last : List;
! 173: tmp : List := l;
! 174:
! 175: begin
! 176: while not Is_Null(tmp) loop
! 177: Append(res,res_last,Random_Lift(lflow,lfupp,Head_Of(tmp).all));
! 178: tmp := Tail_Of(tmp);
! 179: end loop;
! 180: return res;
! 181: end Random_Lift;
! 182:
! 183: function Random_Lift ( lflow,lfupp : Vector; l : Array_of_Lists )
! 184: return Array_of_Lists is
! 185:
! 186: res : Array_of_Lists(l'range);
! 187:
! 188: begin
! 189: for i in res'range loop
! 190: res(i) := Random_Lift(lflow(i),lfupp(i),l(i));
! 191: end loop;
! 192: return res;
! 193: end Random_Lift;
! 194:
! 195: -- RANDOM LINEAR LIFTING :
! 196:
! 197: function Random_Linear_Lift ( lflow,lfupp : integer; v : Vector )
! 198: return Vector is
! 199:
! 200: lf : Vector(v'range) := Random_Vector(v'last,lflow,lfupp);
! 201:
! 202: begin
! 203: return Linear_Lift(lf,v);
! 204: end Random_Linear_Lift;
! 205:
! 206: function Random_Linear_Lift ( lflow,lfupp : integer; l : List )
! 207: return List is
! 208: begin
! 209: if Is_Null(l)
! 210: then return l;
! 211: else declare
! 212: n : constant natural := Head_Of(l)'length;
! 213: lf : Vector(Head_Of(l)'range) := Random_Vector(n,lflow,lfupp);
! 214: begin
! 215: return Linear_Lift(lf,l);
! 216: end;
! 217: end if;
! 218: end Random_Linear_Lift;
! 219:
! 220: function Random_Linear_Lift ( lflow,lfupp : Vector; l : Array_of_Lists )
! 221: return Array_of_Lists is
! 222:
! 223: res : Array_of_Lists(l'range);
! 224:
! 225: begin
! 226: for i in res'range loop
! 227: res(i) := Random_Linear_Lift(lflow(i),lfupp(i),l(i));
! 228: end loop;
! 229: return res;
! 230: end Random_Linear_Lift;
! 231:
! 232: -- POINT-WISE LIFTING :
! 233:
! 234: function Point_Lift ( lf : integer; v : Vector ) return Vector is
! 235:
! 236: res : Vector(v'first..v'last+1);
! 237:
! 238: begin
! 239: res(v'range) := v;
! 240: res(res'last) := lf;
! 241: return res;
! 242: end Point_Lift;
! 243:
! 244: function Point_Lift ( lf : Vector; l : List ) return List is
! 245:
! 246: res,res_last : List;
! 247: tmp : List := l;
! 248: ind : integer := lf'first;
! 249:
! 250: begin
! 251: while not Is_Null(tmp) loop
! 252: Append(res,res_last,Point_Lift(lf(ind),Head_Of(tmp).all));
! 253: ind := ind + 1;
! 254: tmp := Tail_Of(tmp);
! 255: end loop;
! 256: return res;
! 257: end Point_Lift;
! 258:
! 259: function Point_Lift ( lf : VecVec;
! 260: l : Array_of_Lists ) return Array_of_Lists is
! 261:
! 262: res : Array_of_Lists(l'range);
! 263:
! 264: begin
! 265: for i in res'range loop
! 266: res(i) := Point_Lift(lf(i).all,l(i));
! 267: end loop;
! 268: return res;
! 269: end Point_Lift;
! 270:
! 271: end Integer_Lifting_Functions;
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to integer_lifting_functions.adb CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / PHC / Ada / Root_Counts / Stalift