Return to increment_power_continuation.ads CVS log

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

Annotation of OpenXM_contrib/PHC/Ada/Continuation/increment_power_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_Floating_Vectors;
                      4: with Standard_Complex_Numbers;           use Standard_Complex_Numbers;
                      5: with Standard_Complex_Vectors;           use Standard_Complex_Vectors;
                      6: with Standard_Complex_Matrices;          use Standard_Complex_Matrices;
                      7: with Standard_Complex_Solutions;         use Standard_Complex_Solutions;
                      8:
                      9: package Increment_Power_Continuation is
                     10:
                     11: -- DESCRIPTION :
                     12: --   This package implements routines for tracing solution paths of a
                     13: --   polynomial homotopy that is nonlinear in the continuation parameter.
                     14: --   In particular, the following assumptions are made:
                     15: --     1) paths do not turn back, i.e.: increment-and-fix continuation;
                     16: --     2) the nonlinearity is due to high powers of the continuation parameter.
                     17: --   The reporting version writes intermediate results on file.
                     18:
                     19:   generic
                     20:
                     21:     with function Norm ( x : Vector ) return double_float;
                     22:     with function H  ( x : Vector; t : Complex_Number ) return Vector;
                     23:     with function dH ( x : Vector; t : Complex_Number ) return Matrix;
                     24:
                     25:   procedure Silent_Power_Continue
                     26:                ( sols : in out Solution_List;
                     27:                  powers : in Standard_Floating_Vectors.Vector );
                     28:
                     29:   generic
                     30:
                     31:     with function Norm ( x : Vector ) return double_float;
                     32:     with function H  ( x : Vector; t : Complex_Number ) return Vector;
                     33:     with function dH ( x : Vector; t : Complex_Number ) return Matrix;
                     34:
                     35:   procedure Reporting_Power_Continue
                     36:                ( file : in file_type; sols : in out Solution_List;
                     37:                  powers : in Standard_Floating_Vectors.Vector );
                     38:
                     39:   -- DESCRIPTION :
                     40:   --   Increment-and-fix continuation with power predictor.
                     41:
                     42:   -- ON ENTRY :
                     43:   --   file      to write intermediate results on (if Reporting_);
                     44:   --   sols      the start solutions;
                     45:   --   powers    vector of powers sorted in ascending order.
                     46:
                     47:   -- ON RETURN :
                     48:   --   sols      the computed solutions.
                     49:
                     50: end Increment_Power_Continuation;