Return to increment_and_fix_continuation.ads CVS log

Up to [local] / OpenXM_contrib / PHC / Ada / Continuation

Annotation of OpenXM_contrib/PHC/Ada/Continuation/increment_and_fix_continuation.ads, Revision 1.1.1.1

1.1       maekawa     1: with text_io;                            use text_io;
                      2: with Standard_Floating_Numbers;          use Standard_Floating_Numbers;
                      3: with Standard_Complex_Numbers;           use Standard_Complex_Numbers;
                      4: with Standard_Floating_Vectors;
                      5: with Standard_Floating_VecVecs;          use Standard_Floating_VecVecs;
                      6: with Standard_Complex_Vectors;           use Standard_Complex_Vectors;
                      7: with Standard_Complex_Matrices;          use Standard_Complex_Matrices;
                      8: with Standard_Complex_Solutions;         use Standard_Complex_Solutions;
                      9:
                     10: package Increment_and_Fix_Continuation is
                     11:
                     12: -- DESCRIPTION :
                     13: --   This package provides a general implementation of an increment-and-fix
                     14: --   continuation method.  The generic parameters are a norm, an evaluator and
                     15: --   a differentiator of the homotopy.
                     16: --   There are two basic versions: a silent and a reporting one.
                     17: --   The silent continuation simply performs its calculations without output
                     18: --   of intermediate results.  The reporting continuation routine allows to
                     19: --   put various kinds of intermediate results on a file.
                     20: --   It is assumed that the continuation parameters are already determined
                     21: --   before calling these routines (see Continuation_Parameters).
                     22: --   For both the silent and the reporting version, the facility is added
                     23: --   to estimate the directions of the solution paths, useful in a polyhedral
                     24: --   end game.
                     25:
                     26:   generic
                     27:
                     28:     with function Norm ( x : Vector ) return double_float;
                     29:     with function H  ( x : Vector; t : Complex_Number ) return Vector;
                     30:     with function dH ( x : Vector; t : Complex_Number ) return Vector;
                     31:     with function dH ( x : Vector; t : Complex_Number ) return Matrix;
                     32:
                     33:   procedure Silent_Continue
                     34:                ( sols : in out Solution_List; proj : in boolean;
                     35:                  target : in Complex_Number := Create(1.0) );
                     36:
                     37:   generic
                     38:
                     39:     with function Norm ( x : Vector ) return double_float;
                     40:     with function H  ( x : Vector; t : Complex_Number ) return Vector;
                     41:     with function dH ( x : Vector; t : Complex_Number ) return Vector;
                     42:     with function dH ( x : Vector; t : Complex_Number ) return Matrix;
                     43:
                     44:   procedure Reporting_Continue
                     45:                ( file : in file_type; sols : in out Solution_List;
                     46:                  proj : in boolean;
                     47:                  target : in Complex_Number := Create(1.0) );
                     48:
                     49:   -- DESCRIPTION :
                     50:   --   This routine implements the continuation strategy.
                     51:
                     52:   -- ON ENTRY :
                     53:   --   file      to write intermediate results on (if Reporting_);
                     54:   --   sols      the start solutions;
                     55:   --   proj      for projective-perpendicular path following;
                     56:   --   target    value for the continuation parameter at the end.
                     57:
                     58:   -- ON RETURN :
                     59:   --   sols      the computed solutions.
                     60:
                     61: -- WITH THE ESTIMATION OF THE PATH DIRECTIONS :
                     62:
                     63:   generic
                     64:
                     65:     with function Norm ( x : Vector ) return double_float;
                     66:     with function H  ( x : Vector; t : Complex_Number ) return Vector;
                     67:     with function dH ( x : Vector; t : Complex_Number ) return Vector;
                     68:     with function dH ( x : Vector; t : Complex_Number ) return Matrix;
                     69:
                     70:   procedure Silent_Toric_Continue
                     71:                ( sols : in out Solution_List; proj : in boolean;
                     72:                  v : in out VecVec;
                     73:                  errv : in out Standard_Floating_Vectors.Vector;
                     74:                  target : in Complex_Number := Create(1.0) );
                     75:
                     76:   generic
                     77:
                     78:     with function Norm ( x : Vector ) return double_float;
                     79:     with function H  ( x : Vector; t : Complex_Number ) return Vector;
                     80:     with function dH ( x : Vector; t : Complex_Number ) return Vector;
                     81:     with function dH ( x : Vector; t : Complex_Number ) return Matrix;
                     82:
                     83:   procedure Reporting_Toric_Continue
                     84:                ( file : in file_type; sols : in out Solution_List;
                     85:                  proj : in boolean; v : in out VecVec;
                     86:                  errv : in out Standard_Floating_Vectors.Vector;
                     87:                  target : in Complex_Number := Create(1.0) );
                     88:
                     89:   -- DESCRIPTION :
                     90:   --   This routine implements the continuation strategy with the estimation
                     91:   --   of the directions of the solution paths at the end.
                     92:
                     93:   -- ON ENTRY :
                     94:   --   file      to write intermediate results on (if Reporting_);
                     95:   --   sols      the start solutions;
                     96:   --   proj      for projective-perpendicular path following;
                     97:   --   v         v must be initialized with zero vectors
                     98:   --             and v'range is 1..Length_Of(sols);
                     99:   --   errv      errors on the computed directions;
                    100:   --   target    value for the continuation parameter at the end.
                    101:
                    102:   -- ON RETURN :
                    103:   --   sols      the computed solutions;
                    104:   --   v         directions of the solution paths;
                    105:   --   errv      errors on the computed directions.
                    106:
                    107: end Increment_and_Fix_Continuation;