Return to face_cardinalities.adb CVS log

Up to [local] / OpenXM_contrib / PHC / Ada / Math_Lib / Supports

Annotation of OpenXM_contrib/PHC/Ada/Math_Lib/Supports/face_cardinalities.adb, Revision 1.1.1.1

1.1       maekawa     1: with Standard_Floating_Numbers;          use Standard_Floating_Numbers;
                      2: with Integer_Face_Enumerators;
                      3: with Floating_Face_Enumerators;
                      4:
                      5: package body Face_Cardinalities is
                      6:
                      7:   function fvector ( pts : in Standard_Integer_VecVecs.VecVec )
                      8:                    return Vector is
                      9:
                     10:   -- ALGORITHM : plain enumeration of vertices, edges, k-faces...
                     11:
                     12:     use Integer_Face_Enumerators;
                     13:
                     14:     n : constant natural := pts(pts'first).all'length;
                     15:     f : Vector(-1..n);
                     16:
                     17:     procedure Count_Vertex ( i : in integer; cont : out boolean ) is
                     18:     begin
                     19:       f(0) := f(0) + 1;
                     20:       cont := true;
                     21:     end Count_Vertex;
                     22:     procedure Count_Vertices is new Enumerate_Vertices(Count_Vertex);
                     23:
                     24:     procedure Count_Edge ( i,j : in integer; cont : out boolean ) is
                     25:     begin
                     26:       f(1) := f(1) + 1;
                     27:       cont := true;
                     28:     end Count_Edge;
                     29:     procedure Count_Edges is new Enumerate_Edges(Count_Edge);
                     30:
                     31:     procedure Count_Face ( face : in Vector; cont : out boolean ) is
                     32:     begin
                     33:       f(face'length-1) := f(face'length-1) + 1;
                     34:       cont := true;
                     35:     end Count_Face;
                     36:     procedure Count_Faces is new Enumerate_Faces(Count_Face);
                     37:
                     38:   begin
                     39:     f(-1) := 1;
                     40:     f(0..n) := (0..n => 0);
                     41:     Count_Vertices(pts);
                     42:     Count_Edges(pts);
                     43:     for i in 2..(n-1) loop
                     44:       Count_Faces(i,pts);
                     45:       exit when (f(i) = 0);
                     46:     end loop;
                     47:     if f(n-1) > 1
                     48:      then f(n) := 1;
                     49:     end if;
                     50:     return f;
                     51:   end fvector;
                     52:
                     53:   function fvector ( pts : in Standard_Floating_VecVecs.VecVec )
                     54:                    return Vector is
                     55:
                     56:   -- ALGORITHM : plain enumeration of vertices, edges, k-faces...
                     57:
                     58:     use Floating_Face_Enumerators;
                     59:
                     60:     n : constant natural := pts(pts'first).all'length;
                     61:     f : Vector(-1..n);
                     62:     tol : constant double_float := 10.0**(-8); --10.0**(-12);
                     63:
                     64:     procedure Count_Vertex ( i : in integer; cont : out boolean ) is
                     65:     begin
                     66:       f(0) := f(0) + 1;
                     67:       cont := true;
                     68:     end Count_Vertex;
                     69:     procedure Count_Vertices is new Enumerate_Vertices(Count_Vertex);
                     70:
                     71:     procedure Count_Edge ( i,j : in integer; cont : out boolean ) is
                     72:     begin
                     73:       f(1) := f(1) + 1;
                     74:       cont := true;
                     75:     end Count_Edge;
                     76:     procedure Count_Edges is new Enumerate_Edges(Count_Edge);
                     77:
                     78:     procedure Count_Face ( face : in Vector; cont : out boolean ) is
                     79:     begin
                     80:       f(face'length-1) := f(face'length-1) + 1;
                     81:       cont := true;
                     82:     end Count_Face;
                     83:     procedure Count_Faces is new Enumerate_Faces(Count_Face);
                     84:
                     85:   begin
                     86:     f(-1) := 1;
                     87:     f(0..n) := (0..n => 0);
                     88:     Count_Vertices(pts,tol);
                     89:     Count_Edges(pts,tol);
                     90:     for i in 2..(n-1) loop
                     91:       Count_Faces(i,pts,tol);
                     92:       exit when (f(i) = 0);
                     93:     end loop;
                     94:     if f(n-1) > 1
                     95:      then f(n) := 1;
                     96:     end if;
                     97:     return f;
                     98:   end fvector;
                     99:
                    100: end Face_Cardinalities;