Return to driver_for_minkowski_polynomials.adb CVS log

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

Annotation of OpenXM_contrib/PHC/Ada/Root_Counts/Dynlift/driver_for_minkowski_polynomials.adb, Revision 1.1

1.1     ! maekawa     1: with integer_io;                         use integer_io;
        !             2: with Timing_Package;                     use Timing_Package;
        !             3: with Standard_Complex_Numbers;           use Standard_Complex_Numbers;
        !             4: with Standard_Integer_Vectors_io;        use Standard_Integer_Vectors_io;
        !             5: with Standard_Complex_Polynomials;       use Standard_Complex_Polynomials;
        !             6: with Integer_Mixed_Subdivisions_io;      use Integer_Mixed_Subdivisions_io;
        !             7: with Minkowski_Polynomials;              use Minkowski_Polynomials;
        !             8:
        !             9: procedure Driver_for_Minkowski_Polynomials
        !            10:                 ( file : in file_type;
        !            11:                   n : in natural; mix : in Vector; t : in Triangulation;
        !            12:                   alltri : in boolean; mixsub : out Mixed_Subdivision ) is
        !            13:
        !            14:   mp : Poly;
        !            15:   timer : timing_widget;
        !            16:
        !            17:   procedure Coefficient_Subdivision
        !            18:                  ( submix : in Vector; sub : in Mixed_Subdivision;
        !            19:                    vol : out natural ) is
        !            20:
        !            21:     len : constant natural := Length_Of(sub);
        !            22:     wrksub : Mixed_Subdivision := sub;
        !            23:     volu : natural;
        !            24:
        !            25:   begin
        !            26:     put(file,"#Cells of type "); put(file,submix); put(file," is " );
        !            27:     put(file,len,1); new_line(file);
        !            28:     if len = 0
        !            29:      then volu := 0;
        !            30:      else put(file,n,submix,wrksub,volu);
        !            31:     end if;
        !            32:     put(file,"The volume equals : "); put(file,volu,1); new_line(file);
        !            33:     new_line(file);
        !            34:     vol := volu;
        !            35:   end Coefficient_Subdivision;
        !            36:   procedure Coefficient_Subdivisions is
        !            37:     new Minkowski_Polynomial_Subdivisions(Coefficient_Subdivision);
        !            38:
        !            39:   procedure Write_Minkowski_Polynomial ( file : in file_type; p : in Poly ) is
        !            40:
        !            41:     first : boolean;
        !            42:     cnt : natural := 0;
        !            43:
        !            44:     procedure Minkowski_Term ( t : in Term; cont : out boolean ) is
        !            45:     begin
        !            46:       if first
        !            47:        then first := false;
        !            48:        else put(file," + "); cnt := cnt + 3;
        !            49:       end if;
        !            50:       put(file,integer(REAL_PART(t.cf)),1); cnt := cnt + 2;
        !            51:       for i in t.dg'range loop
        !            52:         if t.dg(i) /= 0
        !            53:          then put(file,"*l"); put(file,i,1); cnt := cnt + 3;
        !            54:               if t.dg(i) /= 1
        !            55:                then put(file,"^"); put(file,t.dg(i),1);
        !            56:                     cnt := cnt + 2;
        !            57:               end if;
        !            58:         end if;
        !            59:       end loop;
        !            60:       if cnt > 60
        !            61:        then new_line(file); cnt := 0;
        !            62:       end if;
        !            63:       cont := true;
        !            64:     end Minkowski_Term;
        !            65:     procedure Minkowski_Terms is new Visiting_Iterator(Minkowski_Term);
        !            66:
        !            67:   begin
        !            68:     Minkowski_Terms(p);
        !            69:   end Write_Minkowski_Polynomial;
        !            70:
        !            71: begin
        !            72:   tstart(timer);
        !            73:   mp := Minkowski_Polynomial(n,mix'last);
        !            74:   if not alltri
        !            75:    then Minkowski_Polynomial(mp,t,n,mix,mixsub);
        !            76:    else Coefficient_Subdivisions(mp,t,n,mix,mixsub);
        !            77:   end if;
        !            78:   tstop(timer);
        !            79:   put_line(file,"the Minkowski-polynomial : ");
        !            80:   Write_Minkowski_Polynomial(file,mp);
        !            81:   new_line(file);
        !            82:   new_line(file);
        !            83:   print_times(file,timer,"computing the Minkowski-polynomial");
        !            84:   Clear(mp);
        !            85: end Driver_for_Minkowski_Polynomials;