Return to symmetric_randomize.adb CVS log

Up to [local] / OpenXM_contrib / PHC / Ada / Root_Counts / Symmetry

Annotation of OpenXM_contrib/PHC/Ada/Root_Counts/Symmetry/symmetric_randomize.adb, Revision 1.1

1.1     ! maekawa     1: with Standard_Complex_Numbers;           use Standard_Complex_Numbers;
        !             2: with Standard_Random_Numbers;            use Standard_Random_Numbers;
        !             3: with Standard_Integer_Vectors;           use Standard_Integer_Vectors;
        !             4: with Standard_Complex_Laur_Polys;        use Standard_Complex_Laur_Polys;
        !             5: with Permutations,Permute_Operations;    use Permutations,Permute_Operations;
        !             6:
        !             7: function Symmetric_Randomize ( p : Laur_Sys; v,w : List_of_Permutations )
        !             8:                              return Laur_Sys is
        !             9:
        !            10:   res : Laur_Sys(p'range);
        !            11:
        !            12:   procedure Symmetric_Randomize_Terms ( index : in natural; py : in Poly ) is
        !            13:
        !            14:     tpy : Term;
        !            15:
        !            16:     procedure Permute_and_Randomize ( t : in Term ) is
        !            17:
        !            18:       tmpv,tmpw : List_of_Permutations;
        !            19:
        !            20:     begin
        !            21:       tmpv := v;  tmpw := w;
        !            22:       while not Is_Null(tmpv) loop
        !            23:         declare
        !            24:           permt : Term := Permutation(Head_Of(tmpv).all)*t;
        !            25:           indw : natural := Head_Of(tmpw)(index);
        !            26:         begin
        !            27:           if Coeff(res(indw),permt.dg) = Create(0.0)
        !            28:            then Add(res(indw),permt);
        !            29:           end if;
        !            30:           Clear(permt);
        !            31:         end;
        !            32:         tmpv := Tail_Of(tmpv);
        !            33:         tmpw := Tail_Of(tmpw);
        !            34:       end loop;
        !            35:     end Permute_and_Randomize;
        !            36:
        !            37:     procedure Pick_Term ( t : in Term; cont : out boolean ) is
        !            38:     begin
        !            39:       if Coeff(res(index),t.dg) = Create(0.0)
        !            40:        then Copy(t,tpy);
        !            41:             tpy.cf := Random1;
        !            42:             cont := false;
        !            43:        else cont := true;
        !            44:       end if;
        !            45:     end Pick_Term;
        !            46:     procedure Pick_A_Term is new Visiting_Iterator(Pick_Term);
        !            47:
        !            48:   begin
        !            49:     tpy.cf := Create(0.0);
        !            50:     Pick_A_Term(py);
        !            51:     if tpy.cf /= Create(0.0)
        !            52:      then Permute_and_Randomize(tpy);
        !            53:           Clear(tpy);
        !            54:     end if;
        !            55:   end Symmetric_Randomize_Terms;
        !            56:
        !            57: begin
        !            58:   res := (res'range => Null_Poly);
        !            59:   for k in res'range loop
        !            60:     while Number_of_Terms(res(k)) < Number_of_Terms(p(k)) loop
        !            61:       Symmetric_Randomize_Terms(k,p(k));
        !            62:     end loop;
        !            63:   end loop;
        !            64:   return res;
        !            65: end Symmetric_Randomize;