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     ! 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;