![]() ![]() | ![]() |
File: [local] / OpenXM_contrib / PHC / Ada / Continuation / increment_power_continuation.ads (download)
Revision 1.1.1.1 (vendor branch), Sun Oct 29 17:45:22 2000 UTC (23 years, 10 months ago) by maekawa
Import the second public release of PHCpack. OKed by Jan Verschelde. |
with text_io; use text_io; with Standard_Floating_Numbers; use Standard_Floating_Numbers; with Standard_Floating_Vectors; with Standard_Complex_Numbers; use Standard_Complex_Numbers; with Standard_Complex_Vectors; use Standard_Complex_Vectors; with Standard_Complex_Matrices; use Standard_Complex_Matrices; with Standard_Complex_Solutions; use Standard_Complex_Solutions; package Increment_Power_Continuation is -- DESCRIPTION : -- This package implements routines for tracing solution paths of a -- polynomial homotopy that is nonlinear in the continuation parameter. -- In particular, the following assumptions are made: -- 1) paths do not turn back, i.e.: increment-and-fix continuation; -- 2) the nonlinearity is due to high powers of the continuation parameter. -- The reporting version writes intermediate results on file. generic with function Norm ( x : Vector ) return double_float; with function H ( x : Vector; t : Complex_Number ) return Vector; with function dH ( x : Vector; t : Complex_Number ) return Matrix; procedure Silent_Power_Continue ( sols : in out Solution_List; powers : in Standard_Floating_Vectors.Vector ); generic with function Norm ( x : Vector ) return double_float; with function H ( x : Vector; t : Complex_Number ) return Vector; with function dH ( x : Vector; t : Complex_Number ) return Matrix; procedure Reporting_Power_Continue ( file : in file_type; sols : in out Solution_List; powers : in Standard_Floating_Vectors.Vector ); -- DESCRIPTION : -- Increment-and-fix continuation with power predictor. -- ON ENTRY : -- file to write intermediate results on (if Reporting_); -- sols the start solutions; -- powers vector of powers sorted in ascending order. -- ON RETURN : -- sols the computed solutions. end Increment_Power_Continuation;