File: [local] / OpenXM_contrib / PHC / Ada / Continuation / correctors.adb (download)
Revision 1.1.1.1 (vendor branch), Sun Oct 29 17:45:22 2000 UTC (23 years, 10 months ago) by maekawa
Branch: PHC, MAIN
CVS Tags: v2, maekawa-ipv6, RELEASE_1_2_3, RELEASE_1_2_2_KNOPPIX_b, RELEASE_1_2_2_KNOPPIX, RELEASE_1_2_2, RELEASE_1_2_1, HEAD Changes since 1.1: +0 -0
lines
Import the second public release of PHCpack.
OKed by Jan Verschelde.
|
with Standard_Natural_Vectors;
with Standard_Complex_Norms_Equals; use Standard_Complex_Norms_Equals;
with Standard_Complex_Linear_Solvers; use Standard_Complex_Linear_Solvers;
with Standard_Complex_QR_Decomposition; use Standard_Complex_QR_Decomposition;
with Standard_Complex_Least_Squares; use Standard_Complex_Least_Squares;
with Process_io; use Process_io;
package body Correctors is
-- AUXILIARIES FOR IMPLEMENTING MULTIPLE CORRECTORS :
procedure Equals ( s : in Solu_Info_Array; x : in Vector; i : in natural;
d : in double_float; j : in out natural ) is
eq : boolean;
begin
while j < i loop
eq := true;
for k in x'range loop
eq := Equal(s(j).sol.v(k),x(k),d);
exit when not eq;
end loop;
exit when eq;
j := j + 1;
end loop;
end Equals;
generic
with procedure Corrector ( s : in out Solu_Info; c : in Corr_Pars );
procedure Multiple_Silent_Corrector
( s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean );
-- DESCRIPTION :
-- Allows to build any multiple silent corrector,
-- depending on the corrector supplied as generic parameter.
generic
with procedure Corrector ( file : in file_type;
s : in out Solu_Info; c : in Corr_Pars );
procedure Multiple_Reporting_Corrector
( file : in file_type; s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean );
-- DESCRIPTION :
-- Allows to build any multiple reporting corrector,
-- depending on the corrector supplied as generic parameter.
procedure Multiple_Silent_Corrector
( s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
i,j : natural;
ffail : boolean := false;
begin
i := pivot;
loop
Corrector(s(i),c);
ffail := (((s(i).resa > c.epsaf) and then (s(i).cora > c.epsax))
or else ((s(i).resr > c.epsrf) and then (s(i).corr > c.epsrx)));
if not ffail
then j := 1;
Equals(s,s(i).sol.v,i,dist_sols,j);
if j /= i
then ffail := true;
end if;
end if;
exit when ffail;
i := i + 1;
if i > s'last
then i := s'first;
end if;
exit when (i = pivot);
end loop;
if ffail
then pivot := i;
end if;
fail := ffail;
end Multiple_Silent_Corrector;
procedure Multiple_Reporting_Corrector
( file : in file_type; s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
i,j : natural;
ffail : boolean := false;
begin
i := pivot;
loop
Write_path(file,i);
Corrector(file,s(i),c);
sWrite(file,s(i).sol.all);
ffail := (((s(i).resa > c.epsaf) and then (s(i).cora > c.epsax))
or else ((s(i).resr > c.epsrf) and then (s(i).corr > c.epsrx)));
if not ffail
then j := 1;
Equals(s,s(i).sol.v,i,dist_sols,j);
if j /= i
then ffail := true;
end if;
end if;
exit when ffail;
i := i + 1;
if i > s'last
then i := s'first;
end if;
exit when (i = pivot);
end loop;
if ffail
then pivot := i;
end if;
fail := ffail;
end Multiple_Reporting_Corrector;
-- TARGET PROCEDURES :
procedure Affine_Single_Loose_Normal_Silent_Corrector
( s : in out Solu_Info; c : in Corr_Pars ) is
info : integer;
y : Vector(s.sol.v'range) := H(s.sol.v,s.sol.t);
j : Matrix(y'range,y'range);
ipvt : Standard_Natural_Vectors.Vector(y'range);
nit : natural := 0;
normv : double_float;
begin
s.resa := Norm(y); s.cora := 1.0; -- INITIALIZATION
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
while nit < c.maxit loop -- STOP WHEN MAXIMUM #ITERATIONS REACHED
j := dH(s.sol.v,s.sol.t);
lufac(j,y'last,ipvt,info);
exit when info /= 0; -- STOP WHEN SINGULAR JACOBIAN
Min(y);
lusolve(j,y'last,ipvt,y); s.cora := Norm(y);
Add(s.sol.v,y); s.length_path := s.length_path + s.cora;
y := H(s.sol.v,s.sol.t); s.resa := Norm(y);
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
nit := nit + 1;
exit when (((s.cora <= c.epsax) or else (s.resa <= c.epsaf))
and then ((s.corr <= c.epsrx) or else (s.resr <= c.epsrf)));
end loop; -- STOP WHEN DESIRED PRECISION REACHED
s.niter := s.niter + nit; s.nsyst := s.nsyst + nit;
exception
when numeric_error => return;
end Affine_Single_Loose_Normal_Silent_Corrector;
procedure Affine_Single_Loose_Normal_Reporting_Corrector
( file : in file_type;
s : in out Solu_Info; c : in Corr_Pars ) is
info : integer;
y : Vector(s.sol.v'range) := H(s.sol.v,s.sol.t);
j : Matrix(y'range,y'range);
ipvt : Standard_Natural_Vectors.Vector(y'range);
nit : natural := 0;
normv : double_float;
begin
s.resa := Norm(y); s.cora := 1.0; -- INITIALIZATION
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv;
s.resr := s.resa / normv;
end if;
while nit < c.maxit loop -- STOP WHEN MAXIMUM #ITERATIONS REACHED
j := dH(s.sol.v,s.sol.t);
lufac(j,y'last,ipvt,info);
exit when info /= 0; -- STOP WHEN SINGULAR JACOBIAN
Min(y);
lusolve(j,y'last,ipvt,y); s.cora := Norm(y);
Add(s.sol.v,y); s.length_path := s.length_path + s.cora;
y := H(s.sol.v,s.sol.t); s.resa := Norm(y);
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
cWrite(file,s.cora,s.corr,s.resa,s.resr); -- WRITE STATISTICS
nit := nit + 1;
exit when (((s.cora <= c.epsax) or else (s.resa <= c.epsaf))
and then ((s.corr <= c.epsrx) or else (s.resr <= c.epsrf)));
end loop; -- STOP WHEN DESIRED PRECISION REACHED
s.niter := s.niter + nit; s.nsyst := s.nsyst + nit;
exception
when numeric_error => return;
end Affine_Single_Loose_Normal_Reporting_Corrector;
procedure Affine_Single_Loose_Conditioned_Silent_Corrector
( s : in out Solu_Info; c : in Corr_Pars ) is
y : Vector(s.sol.v'range) := H(s.sol.v,s.sol.t);
n : constant natural := y'last;
j : Matrix(y'range,y'range);
qraux : Standard_Complex_Vectors.Vector(y'range)
:= (y'range => Create(0.0));
dum : Standard_Complex_Vectors.Vector(y'range);
ipvt : Standard_Natural_Vectors.Vector(y'range);
info : integer;
nit : natural := 0;
normv : double_float;
begin
s.resa := Norm(y); s.cora := 1.0; -- INITIALIZATION
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
while nit < c.maxit loop -- STOP WHEN MAXIMUM #ITERATIONS REACHED
j := dH(s.sol.v,s.sol.t);
QRD(j,qraux,ipvt,false); -- QR WITHOUT PIVOTING
Min(y); -- FOLLOWED BY LEAST SQUARES
QRLS(j,n,n,n,qraux,y,dum,dum,y,dum,dum,100,info);
s.cora := Norm(y);
Add(s.sol.v,y); s.length_path := s.length_path + s.cora;
y := H(s.sol.v,s.sol.t); s.resa := Norm(y);
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
nit := nit + 1;
exit when (((s.cora <= c.epsax) or else (s.resa <= c.epsaf))
and then ((s.corr <= c.epsrx) or else (s.resr <= c.epsrf)));
end loop; -- STOP WHEN DESIRED PRECISION REACHED
s.niter := s.niter + nit; s.nsyst := s.nsyst + nit;
j := dH(s.sol.v,s.sol.t);
lufco(j,y'last,ipvt,s.rcond); -- ESTIMATE CONDITION NUMBER
exception
when numeric_error => return;
end Affine_Single_Loose_Conditioned_Silent_Corrector;
procedure Affine_Single_Loose_Conditioned_Reporting_Corrector
( file : in file_type;
s : in out Solu_Info; c : in Corr_Pars ) is
y : Vector(s.sol.v'range) := H(s.sol.v,s.sol.t);
n : constant natural := y'last;
j : Matrix(y'range,y'range);
qraux : Standard_Complex_Vectors.Vector(y'range)
:= (y'range => Create(0.0));
dum : Standard_Complex_Vectors.Vector(y'range);
ipvt : Standard_Natural_Vectors.Vector(y'range);
info : integer;
nit : natural := 0;
normv : double_float;
begin
s.resa := Norm(y); s.cora := 1.0; -- INITIALIZATION
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
while nit < c.maxit loop -- STOP WHEN MAXIMUM #ITERATIONS REACHED
j := dH(s.sol.v,s.sol.t);
QRD(j,qraux,ipvt,false); -- QR WITHOUT PIVOTING
Min(y); -- FOLLOWED BY LEAST SQUARES
QRLS(j,n,n,n,qraux,y,dum,dum,y,dum,dum,100,info);
s.cora := Norm(y);
Add(s.sol.v,y); s.length_path := s.length_path + s.cora;
y := H(s.sol.v,s.sol.t); s.resa := Norm(y);
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
cWrite(file,s.cora,s.corr,s.resa,s.resr); -- WRITE STATISTICS
cWrite(file,s.rcond,s.sol.m);
nit := nit + 1;
exit when (((s.cora <= c.epsax) or else (s.resa <= c.epsaf))
and then ((s.corr <= c.epsrx) or else (s.resr <= c.epsrf)));
end loop; -- STOP WHEN DESIRED PRECISION REACHED
s.niter := s.niter + nit; s.nsyst := s.nsyst + nit;
j := dH(s.sol.v,s.sol.t);
lufco(j,y'last,ipvt,s.rcond); -- ESTIMATE CONDITION NUMBER
exception
when numeric_error => return;
end Affine_Single_Loose_Conditioned_Reporting_Corrector;
procedure Affine_Single_Severe_Normal_Silent_Corrector
( s : in out Solu_Info; c : in Corr_Pars ) is
info : integer;
y : Vector(s.sol.v'range) := H(s.sol.v,s.sol.t);
j : Matrix(y'range,y'range);
ipvt : Standard_Natural_Vectors.Vector(y'range);
nit : natural := 0;
normv,ncora,nresa,ncorr,nresr : double_float;
stop : boolean;
begin
s.resa := Norm(y); s.cora := 1.0; -- INITIALIZATION
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
while nit < c.maxit loop -- STOP WHEN MAXIMUM #ITERATIONS REACHED
j := dH(s.sol.v,s.sol.t);
lufac(j,y'last,ipvt,info);
exit when info /= 0; -- STOP WHEN SINGULAR JACOBIAN
Min(y);
lusolve(j,y'last,ipvt,y); ncora := Norm(y);
Add(s.sol.v,y); s.length_path := s.length_path + s.cora;
y := H(s.sol.v,s.sol.t); nresa := Norm(y);
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then ncorr := ncora / normv; nresr := nresa / normv;
end if;
nit := nit + 1;
stop := (((ncora > s.cora) and then (nresa > s.resa))
or else ((ncorr > s.corr) and then (ncorr > s.corr)));
-- STOP WHEN DIVERGENCE
s.cora := ncora; s.resa := nresa; s.corr := ncorr; s.resr := nresr;
exit when stop;
exit when (((s.cora <= c.epsax) or else (s.resa <= c.epsaf))
and then ((s.corr <= c.epsrx) or else (s.resr <= c.epsrf)));
end loop; -- STOP WHEN DESIRED PRECISION REACHED
s.niter := s.niter + nit; s.nsyst := s.nsyst + nit;
exception
when numeric_error => return;
end Affine_Single_Severe_Normal_Silent_Corrector;
procedure Affine_Single_Severe_Normal_Reporting_Corrector
( file : in file_type;
s : in out Solu_Info; c : in Corr_Pars ) is
info : integer;
y : Vector(s.sol.v'range) := H(s.sol.v,s.sol.t);
j : Matrix(y'range,y'range);
ipvt : Standard_Natural_Vectors.Vector(y'range);
nit : natural := 0;
normv,ncora,nresa,ncorr,nresr : double_float;
stop : boolean;
begin
s.resa := Norm(y); s.cora := 1.0; -- INITIALIZATION
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
while nit < c.maxit loop -- STOP WHEN MAXIMUM #ITERATIONS REACHED
j := dH(s.sol.v,s.sol.t);
lufac(j,y'last,ipvt,info);
exit when info /= 0; -- STOP WHEN SINGULAR JACOBIAN
Min(y);
lusolve(j,y'last,ipvt,y); ncora := Norm(y);
Add(s.sol.v,y); s.length_path := s.length_path + s.cora;
y := H(s.sol.v,s.sol.t); nresa := Norm(y);
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then ncorr := ncora / normv; nresr := nresa / normv;
end if;
cWrite(file,ncora,ncorr,nresa,nresr); -- WRITE STATISTICS
nit := nit + 1;
stop := (((ncora > s.cora) and then (nresa > s.resa))
or else ((ncorr > s.corr) and then (ncorr > s.corr)));
-- STOP WHEN DIVERGENCE
s.cora := ncora; s.resa := nresa; s.corr := ncorr; s.resr := nresr;
exit when stop;
exit when (((s.cora <= c.epsax) or else (s.resa <= c.epsaf))
and then ((s.corr <= c.epsrx) or else (s.resr <= c.epsrf)));
end loop; -- STOP WHEN DESIRED PRECISION REACHED
s.niter := s.niter + nit; s.nsyst := s.nsyst + nit;
exception
when numeric_error => return;
end Affine_Single_Severe_Normal_Reporting_Corrector;
procedure Affine_Single_Severe_Conditioned_Silent_Corrector
( s : in out Solu_Info; c : in Corr_Pars ) is
y : Vector(s.sol.v'range) := H(s.sol.v,s.sol.t);
n : constant natural := y'last;
j : Matrix(y'range,y'range);
qraux : Standard_Complex_Vectors.Vector(y'range)
:= (y'range => Create(0.0));
dum : Standard_Complex_Vectors.Vector(y'range);
ipvt : Standard_Natural_Vectors.Vector(y'range);
info : integer;
nit : natural := 0;
normv,ncora,nresa,ncorr,nresr : double_float;
stop : boolean;
begin
s.resa := Norm(y); s.cora := 1.0; -- INITIALIZATION
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
while nit < c.maxit loop -- STOP WHEN MAXIMUM #ITERATIONS REACHED
j := dH(s.sol.v,s.sol.t);
QRD(j,qraux,ipvt,false); -- QR WITHOUT PIVOTING
Min(y); -- FOLLOWED BY LEAST SQUARES
QRLS(j,n,n,n,qraux,y,dum,dum,y,dum,dum,100,info);
ncora := Norm(y);
Add(s.sol.v,y); s.length_path := s.length_path + s.cora;
y := H(s.sol.v,s.sol.t); nresa := Norm(y);
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then ncorr := ncora / normv; nresr := nresa / normv;
end if;
nit := nit + 1;
stop := (((ncora > s.cora) and then (nresa > s.resa))
or else ((ncorr > s.corr) and then (ncorr > s.corr)));
-- STOP WHEN DIVERGENCE
s.cora := ncora; s.resa := nresa; s.corr := ncorr; s.resr := nresr;
exit when stop;
exit when (((s.cora <= c.epsax) or else (s.resa <= c.epsaf))
and then ((s.corr <= c.epsrx) or else (s.resr <= c.epsrf)));
end loop; -- STOP WHEN DESIRED PRECISION REACHED
s.niter := s.niter + nit; s.nsyst := s.nsyst + nit;
j := dH(s.sol.v,s.sol.t);
lufco(j,y'last,ipvt,s.rcond); -- ESTIMATE CONDITION NUMBER
exception
when numeric_error => return;
end Affine_Single_Severe_Conditioned_Silent_Corrector;
procedure Affine_Single_Severe_Conditioned_Reporting_Corrector
( file : in file_type;
s : in out Solu_Info; c : in Corr_Pars ) is
y : Vector(s.sol.v'range) := H(s.sol.v,s.sol.t);
n : constant natural := y'last;
j : Matrix(y'range,y'range);
qraux : Standard_Complex_Vectors.Vector(y'range)
:= (y'range => Create(0.0));
dum : Standard_Complex_Vectors.Vector(y'range);
ipvt : Standard_Natural_Vectors.Vector(y'range);
info : integer;
nit : natural := 0;
normv,ncora,nresa,ncorr,nresr : double_float;
stop : boolean;
begin
s.resa := Norm(y); s.cora := 1.0; -- INITIALIZATION
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
while nit < c.maxit loop -- STOP WHEN MAXIMUM #ITERATIONS REACHED
j := dH(s.sol.v,s.sol.t);
QRD(j,qraux,ipvt,false); -- QR WITHOUT PIVOTING
Min(y); -- FOLLOWED BY LEAST SQUARES
QRLS(j,n,n,n,qraux,y,dum,dum,y,dum,dum,100,info);
ncora := Norm(y);
Add(s.sol.v,y); s.length_path := s.length_path + s.cora;
y := H(s.sol.v,s.sol.t); nresa := Norm(y);
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then ncorr := ncora / normv; nresr := nresa / normv;
end if;
cWrite(file,ncora,ncorr,nresa,nresr); -- WRITE STATISTICS
cWrite(file,s.rcond,s.sol.m);
nit := nit + 1;
stop := (((ncora > s.cora) and then (nresa > s.resa))
or else ((ncorr > s.corr) and then (ncorr > s.corr)));
-- STOP WHEN DIVERGENCE
s.cora := ncora; s.resa := nresa; s.corr := ncorr; s.resr := nresr;
exit when stop;
exit when (((s.cora <= c.epsax) or else (s.resa <= c.epsaf))
and then ((s.corr <= c.epsrx) or else (s.resr <= c.epsrf)));
end loop; -- STOP WHEN DESIRED PRECISION REACHED
s.niter := s.niter + nit; s.nsyst := s.nsyst + nit;
j := dH(s.sol.v,s.sol.t);
lufco(j,y'last,ipvt,s.rcond); -- ESTIMATE CONDITION NUMBER
exception
when numeric_error => return;
end Affine_Single_Severe_Conditioned_Reporting_Corrector;
procedure Affine_Multiple_Loose_Normal_Silent_Corrector
( s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
procedure Single_Corrector is
new Affine_Single_Loose_Normal_Silent_Corrector(Norm,H,dH);
procedure Corrector is new Multiple_Silent_Corrector(Single_Corrector);
begin
Corrector(s,pivot,dist_sols,c,fail);
end Affine_Multiple_Loose_Normal_Silent_Corrector;
procedure Affine_Multiple_Loose_Normal_Reporting_Corrector
( file : in file_type; s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
procedure Single_Corrector is
new Affine_Single_Loose_Normal_Reporting_Corrector(Norm,H,dH);
procedure Corrector is new Multiple_Reporting_Corrector(Single_Corrector);
begin
Corrector(file,s,pivot,dist_sols,c,fail);
end Affine_Multiple_Loose_Normal_Reporting_Corrector;
procedure Affine_Multiple_Loose_Conditioned_Silent_Corrector
( s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
procedure Single_Corrector is
new Affine_Single_Loose_Conditioned_Silent_Corrector(Norm,H,dH);
procedure Corrector is new Multiple_Silent_Corrector(Single_Corrector);
begin
Corrector(s,pivot,dist_sols,c,fail);
end Affine_Multiple_Loose_Conditioned_Silent_Corrector;
procedure Affine_Multiple_Loose_Conditioned_Reporting_Corrector
( file : in file_type; s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
procedure Single_Corrector is
new Affine_Single_Loose_Conditioned_Reporting_Corrector(Norm,H,dH);
procedure Corrector is new Multiple_Reporting_Corrector(Single_Corrector);
begin
Corrector(file,s,pivot,dist_sols,c,fail);
end Affine_Multiple_Loose_Conditioned_Reporting_Corrector;
procedure Affine_Multiple_Severe_Normal_Silent_Corrector
( s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
procedure Single_Corrector is
new Affine_Single_Severe_Normal_Silent_Corrector(Norm,H,dH);
procedure Corrector is new Multiple_Silent_Corrector(Single_Corrector);
begin
Corrector(s,pivot,dist_sols,c,fail);
end Affine_Multiple_Severe_Normal_Silent_Corrector;
procedure Affine_Multiple_Severe_Normal_Reporting_Corrector
( file : in file_type; s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
procedure Single_Corrector is
new Affine_Single_Severe_Normal_Reporting_Corrector(Norm,H,dH);
procedure Corrector is new Multiple_Reporting_Corrector(Single_Corrector);
begin
Corrector(file,s,pivot,dist_sols,c,fail);
end Affine_Multiple_Severe_Normal_Reporting_Corrector;
procedure Affine_Multiple_Severe_Conditioned_Silent_Corrector
( s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
procedure Single_Corrector is
new Affine_Single_Severe_Conditioned_Silent_Corrector(Norm,H,dH);
procedure Corrector is new Multiple_Silent_Corrector(Single_Corrector);
begin
Corrector(s,pivot,dist_sols,c,fail);
end Affine_Multiple_Severe_Conditioned_Silent_Corrector;
procedure Affine_Multiple_Severe_Conditioned_Reporting_Corrector
( file : in file_type; s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
procedure Single_Corrector is
new Affine_Single_Severe_Conditioned_Reporting_Corrector(Norm,H,dH);
procedure Corrector is new Multiple_Reporting_Corrector(Single_Corrector);
begin
Corrector(file,s,pivot,dist_sols,c,fail);
end Affine_Multiple_Severe_Conditioned_Reporting_Corrector;
procedure Projective_Single_Loose_Normal_Silent_Corrector
( s : in out Solu_Info; c : in Corr_Pars ) is
info : integer;
y : Vector(s.sol.v'range) := H(s.sol.v,s.sol.t);
j : Matrix(y'range,y'range);
ipvt : Standard_Natural_Vectors.Vector(y'range);
nit : natural := 0;
normv : double_float;
begin
s.resa := Norm(y); s.cora := 1.0; -- INITIALIZATION
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
while nit < c.maxit loop -- STOP WHEN MAXIMUM #ITERATIONS REACHED
j := dH(s.sol.v,s.sol.t);
for jj in y'range loop
j(j'last(1),jj) := s.sol.v(jj); -- CORRECT PERPENDICULAR
end loop;
lufac(j,y'last,ipvt,info);
exit when info /= 0; -- STOP WHEN SINGULAR JACOBIAN
Min(y);
y(y'last) := Create(0.0); -- IGNORE SCALING EQUATION
lusolve(j,y'last,ipvt,y); s.cora := Norm(y);
Add(s.sol.v,y); s.length_path := s.length_path + s.cora;
y := H(s.sol.v,s.sol.t);
y(y'last) := Create(0.0); s.resa := Norm(y);
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
for jj in s.sol.v'range loop -- SCALE THE SOLUTION
s.sol.v(jj) := s.sol.v(jj)/Create(normv);
end loop;
end if;
nit := nit + 1;
exit when (((s.cora <= c.epsax) or else (s.resa <= c.epsaf))
and then ((s.corr <= c.epsrx) or else (s.resr <= c.epsrf)));
end loop; -- STOP WHEN DESIRED PRECISION REACHED
s.niter := s.niter + nit; s.nsyst := s.nsyst + nit;
exception
when numeric_error => return;
end Projective_Single_Loose_Normal_Silent_Corrector;
procedure Projective_Single_Loose_Normal_Reporting_Corrector
( file : in file_type;
s : in out Solu_Info; c : in Corr_Pars ) is
info : integer;
y : Vector(s.sol.v'range) := H(s.sol.v,s.sol.t);
j : Matrix(y'range,y'range);
ipvt : Standard_Natural_Vectors.Vector(y'range);
nit : natural := 0;
normv : double_float;
begin
s.resa := Norm(y); s.cora := 1.0; -- INITIALIZATION
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
while nit < c.maxit loop -- STOP WHEN MAXIMUM #ITERATIONS REACHED
j := dH(s.sol.v,s.sol.t);
for jj in y'range loop
j(j'last(1),jj) := s.sol.v(jj); -- CORRECT PERPENDICULAR
end loop;
lufac(j,y'last,ipvt,info);
exit when info /= 0; -- STOP WHEN SINGULAR JACOBIAN
Min(y);
y(y'last) := Create(0.0); -- IGNORE SCALING EQUATION
lusolve(j,y'last,ipvt,y); s.cora := Norm(y);
Add(s.sol.v,y); s.length_path := s.length_path + s.cora;
y := H(s.sol.v,s.sol.t);
y(y'last) := Create(0.0); s.resa := Norm(y);
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
for jj in s.sol.v'range loop -- SCALE THE SOLUTION
s.sol.v(jj) := s.sol.v(jj)/Create(normv);
end loop;
end if;
cWrite(file,s.cora,s.corr,s.resa,s.resr); -- WRITE STATISTICS
nit := nit + 1;
exit when (((s.cora <= c.epsax) or else (s.resa <= c.epsaf))
and then ((s.corr <= c.epsrx) or else (s.resr <= c.epsrf)));
end loop; -- STOP WHEN DESIRED PRECISION REACHED
s.niter := s.niter + nit; s.nsyst := s.nsyst + nit;
exception
when numeric_error => return;
end Projective_Single_Loose_Normal_Reporting_Corrector;
procedure Projective_Single_Loose_Conditioned_Silent_Corrector
( s : in out Solu_Info; c : in Corr_Pars ) is
y : Vector(s.sol.v'range) := H(s.sol.v,s.sol.t);
n : constant natural := y'last;
j : Matrix(y'range,y'range);
qraux : Standard_Complex_Vectors.Vector(y'range)
:= (y'range => Create(0.0));
dum : Standard_Complex_Vectors.Vector(y'range);
ipvt : Standard_Natural_Vectors.Vector(y'range);
info : integer;
nit : natural := 0;
normv : double_float;
begin
s.resa := Norm(y); s.cora := 1.0; -- INITIALIZATION
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
while nit < c.maxit loop -- STOP WHEN MAXIMUM #ITERATIONS REACHED
j := dH(s.sol.v,s.sol.t);
for jj in y'range loop
j(j'last(1),jj) := s.sol.v(jj); -- CORRECT PERPENDICULAR
end loop;
QRD(j,qraux,ipvt,false); -- QR WITHOUT PIVOTING
Min(y); -- FOLLOWED BY LEAST SQUARES
y(y'last) := Create(0.0); -- IGNORE SCALING EQUATION
QRLS(j,n,n,n,qraux,y,dum,dum,y,dum,dum,100,info);
s.cora := Norm(y);
Add(s.sol.v,y); s.length_path := s.length_path + s.cora;
y := H(s.sol.v,s.sol.t);
y(y'last) := Create(0.0); s.resa := Norm(y);
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
for jj in s.sol.v'range loop -- SCALE THE SOLUTION
s.sol.v(jj) := s.sol.v(jj)/Create(normv);
end loop;
end if;
nit := nit + 1;
exit when (((s.cora <= c.epsax) or else (s.resa <= c.epsaf))
and then ((s.corr <= c.epsrx) or else (s.resr <= c.epsrf)));
end loop; -- STOP WHEN DESIRED PRECISION REACHED
s.niter := s.niter + nit; s.nsyst := s.nsyst + nit;
j := dH(s.sol.v,s.sol.t);
lufco(j,y'last,ipvt,s.rcond); -- ESTIMATE CONDITION NUMBER
exception
when numeric_error => return;
end Projective_Single_Loose_Conditioned_Silent_Corrector;
procedure Projective_Single_Loose_Conditioned_Reporting_Corrector
( file : in file_type;
s : in out Solu_Info; c : in Corr_Pars ) is
y : Vector(s.sol.v'range) := H(s.sol.v,s.sol.t);
n : constant natural := y'last;
j : Matrix(y'range,y'range);
qraux : Standard_Complex_Vectors.Vector(y'range)
:= (y'range => Create(0.0));
dum : Standard_Complex_Vectors.Vector(y'range);
ipvt : Standard_Natural_Vectors.Vector(y'range);
info : integer;
nit : natural := 0;
normv : double_float;
begin
s.resa := Norm(y); s.cora := 1.0; -- INITIALIZATION
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
while nit < c.maxit loop -- STOP WHEN MAXIMUM #ITERATIONS REACHED
j := dH(s.sol.v,s.sol.t);
for jj in y'range loop
j(j'last(1),jj) := s.sol.v(jj); -- CORRECT PERPENDICULAR
end loop;
QRD(j,qraux,ipvt,false); -- QR WITHOUT PIVOTING
Min(y); -- FOLLOWED BY LEAST SQUARES
y(y'last) := Create(0.0); -- IGNORE SCALING EQUATION
QRLS(j,n,n,n,qraux,y,dum,dum,y,dum,dum,100,info);
s.cora := Norm(y);
Add(s.sol.v,y); s.length_path := s.length_path + s.cora;
y := H(s.sol.v,s.sol.t);
y(y'last) := Create(0.0); s.resa := Norm(y);
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
for jj in s.sol.v'range loop -- SCALE THE SOLUTION
s.sol.v(jj) := s.sol.v(jj)/Create(normv);
end loop;
end if;
cWrite(file,s.cora,s.corr,s.resa,s.resr); -- WRITE STATISTICS
cWrite(file,s.rcond,s.sol.m);
nit := nit + 1;
exit when (((s.cora <= c.epsax) or else (s.resa <= c.epsaf))
and then ((s.corr <= c.epsrx) or else (s.resr <= c.epsrf)));
end loop; -- STOP WHEN DESIRED PRECISION REACHED
s.niter := s.niter + nit; s.nsyst := s.nsyst + nit;
j := dH(s.sol.v,s.sol.t);
lufco(j,y'last,ipvt,s.rcond); -- ESTIMATE CONDITION NUMBER
exception
when numeric_error => return;
end Projective_Single_Loose_Conditioned_Reporting_Corrector;
procedure Projective_Single_Severe_Normal_Silent_Corrector
( s : in out Solu_Info; c : in Corr_Pars ) is
info : integer;
y : Vector(s.sol.v'range) := H(s.sol.v,s.sol.t);
j : Matrix(y'range,y'range);
ipvt : Standard_Natural_Vectors.Vector(y'range);
nit : natural := 0;
normv,ncora,nresa,ncorr,nresr : double_float;
stop : boolean;
begin
s.resa := Norm(y); s.cora := 1.0; -- INITIALIZATION
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
while nit < c.maxit loop -- STOP WHEN MAXIMUM #ITERATIONS REACHED
j := dH(s.sol.v,s.sol.t);
for jj in y'range loop
j(j'last(1),jj) := s.sol.v(jj); -- CORRECT PERPENDICULAR
end loop;
lufac(j,y'last,ipvt,info);
exit when info /= 0; -- STOP WHEN SINGULAR JACOBIAN
Min(y);
y(y'last) := Create(0.0); -- IGNORE SCALING EQUATION
lusolve(j,y'last,ipvt,y); ncora := Norm(y);
Add(s.sol.v,y); s.length_path := s.length_path + s.cora;
y := H(s.sol.v,s.sol.t);
y(y'last) := Create(0.0); nresa := Norm(y);
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then ncorr := ncora / normv; nresr := nresa / normv;
for jj in s.sol.v'range loop -- SCALE THE SOLUTION
s.sol.v(jj) := s.sol.v(jj)/Create(normv);
end loop;
end if;
nit := nit + 1;
stop := (((ncora > s.cora) and then (nresa > s.resa))
or else ((ncorr > s.corr) and then (ncorr > s.corr)));
-- STOP WHEN DIVERGENCE
s.cora := ncora; s.resa := nresa; s.corr := ncorr; s.resr := nresr;
exit when stop;
exit when (((s.cora <= c.epsax) or else (s.resa <= c.epsaf))
and then ((s.corr <= c.epsrx) or else (s.resr <= c.epsrf)));
end loop; -- STOP WHEN DESIRED PRECISION REACHED
s.niter := s.niter + nit; s.nsyst := s.nsyst + nit;
exception
when numeric_error => return;
end Projective_Single_Severe_Normal_Silent_Corrector;
procedure Projective_Single_Severe_Normal_Reporting_Corrector
( file : in file_type;
s : in out Solu_Info; c : in Corr_Pars ) is
info : integer;
y : Vector(s.sol.v'range) := H(s.sol.v,s.sol.t);
j : Matrix(y'range,y'range);
ipvt : Standard_Natural_Vectors.Vector(y'range);
nit : natural := 0;
normv,ncora,nresa,ncorr,nresr : double_float;
stop : boolean;
begin
s.resa := Norm(y); s.cora := 1.0; -- INITIALIZATION
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
while nit < c.maxit loop -- STOP WHEN MAXIMUM #ITERATIONS REACHED
j := dH(s.sol.v,s.sol.t);
for jj in y'range loop
j(j'last(1),jj) := s.sol.v(jj); -- CORRECT PERPENDICULAR
end loop;
lufac(j,y'last,ipvt,info);
exit when info /= 0; -- STOP WHEN SINGULAR JACOBIAN
Min(y);
y(y'last) := Create(0.0); -- IGNORE SCALING EQUATION
lusolve(j,y'last,ipvt,y); ncora := Norm(y);
Add(s.sol.v,y); s.length_path := s.length_path + s.cora;
y := H(s.sol.v,s.sol.t);
y(y'last) := Create(0.0); nresa := Norm(y);
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then ncorr := ncora / normv; nresr := nresa / normv;
for jj in s.sol.v'range loop -- SCALE THE SOLUTION
s.sol.v(jj) := s.sol.v(jj)/Create(normv);
end loop;
end if;
cWrite(file,ncora,ncorr,nresa,nresr); -- WRITE STATISTICS
nit := nit + 1;
stop := (((ncora > s.cora) and then (nresa > s.resa))
or else ((ncorr > s.corr) and then (ncorr > s.corr)));
-- STOP WHEN DIVERGENCE
s.cora := ncora; s.resa := nresa; s.corr := ncorr; s.resr := nresr;
exit when stop;
exit when (((s.cora <= c.epsax) or else (s.resa <= c.epsaf))
and then ((s.corr <= c.epsrx) or else (s.resr <= c.epsrf)));
end loop; -- STOP WHEN DESIRED PRECISION REACHED
s.niter := s.niter + nit; s.nsyst := s.nsyst + nit;
exception
when numeric_error => return;
end Projective_Single_Severe_Normal_Reporting_Corrector;
procedure Projective_Single_Severe_Conditioned_Silent_Corrector
( s : in out Solu_Info; c : in Corr_Pars ) is
y : Vector(s.sol.v'range) := H(s.sol.v,s.sol.t);
n : constant natural := y'last;
j : Matrix(y'range,y'range);
qraux : Standard_Complex_Vectors.Vector(y'range)
:= (y'range => Create(0.0));
dum : Standard_Complex_Vectors.Vector(y'range);
ipvt : Standard_Natural_Vectors.Vector(y'range);
info : integer;
nit : natural := 0;
normv,ncora,nresa,ncorr,nresr : double_float;
stop : boolean;
begin
s.resa := Norm(y); s.cora := 1.0; -- INITIALIZATION
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
while nit < c.maxit loop -- STOP WHEN MAXIMUM #ITERATIONS REACHED
j := dH(s.sol.v,s.sol.t);
for jj in y'range loop
j(j'last(1),jj) := s.sol.v(jj); -- CORRECT PERPENDICULAR
end loop;
QRD(j,qraux,ipvt,false); -- QR WITHOUT PIVOTING
Min(y); -- FOLLOWED BY LEAST SQUARES
y(y'last) := Create(0.0); -- IGNORE SCALING EQUATION
QRLS(j,n,n,n,qraux,y,dum,dum,y,dum,dum,100,info);
ncora := Norm(y);
Add(s.sol.v,y); s.length_path := s.length_path + s.cora;
y := H(s.sol.v,s.sol.t);
y(y'last) := Create(0.0); nresa := Norm(y);
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then ncorr := ncora / normv; nresr := nresa / normv;
for jj in s.sol.v'range loop -- SCALE THE SOLUTION
s.sol.v(jj) := s.sol.v(jj)/Create(normv);
end loop;
end if;
nit := nit + 1;
stop := (((ncora > s.cora) and then (nresa > s.resa))
or else ((ncorr > s.corr) and then (ncorr > s.corr)));
-- STOP WHEN DIVERGENCE
s.cora := ncora; s.resa := nresa; s.corr := ncorr; s.resr := nresr;
exit when stop;
exit when (((s.cora <= c.epsax) or else (s.resa <= c.epsaf))
and then ((s.corr <= c.epsrx) or else (s.resr <= c.epsrf)));
end loop; -- STOP WHEN DESIRED PRECISION REACHED
s.niter := s.niter + nit; s.nsyst := s.nsyst + nit;
j := dH(s.sol.v,s.sol.t);
lufco(j,y'last,ipvt,s.rcond); -- ESTIMATE CONDITION NUMBER
exception
when numeric_error => return;
end Projective_Single_Severe_Conditioned_Silent_Corrector;
procedure Projective_Single_Severe_Conditioned_Reporting_Corrector
( file : in file_type;
s : in out Solu_Info; c : in Corr_Pars ) is
y : Vector(s.sol.v'range) := H(s.sol.v,s.sol.t);
n : constant natural := y'last;
j : Matrix(y'range,y'range);
qraux : Standard_Complex_Vectors.Vector(y'range)
:= (y'range => Create(0.0));
dum : Standard_Complex_Vectors.Vector(y'range);
ipvt : Standard_Natural_Vectors.Vector(y'range);
info : integer;
nit : natural := 0;
normv,ncora,nresa,ncorr,nresr : double_float;
stop : boolean;
begin
s.resa := Norm(y); s.cora := 1.0; -- INITIALIZATION
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then s.corr := s.cora / normv; s.resr := s.resa / normv;
end if;
while nit < c.maxit loop -- STOP WHEN MAXIMUM #ITERATIONS REACHED
j := dH(s.sol.v,s.sol.t);
for jj in y'range loop
j(j'last(1),jj) := s.sol.v(jj); -- CORRECT PERPENDICULAR
end loop;
QRD(j,qraux,ipvt,false); -- QR WITHOUT PIVOTING
Min(y); -- FOLLOWED BY LEAST SQUARES
y(y'last) := Create(0.0); -- IGNORE SCALING EQUATION
QRLS(j,n,n,n,qraux,y,dum,dum,y,dum,dum,100,info);
ncora := Norm(y);
Add(s.sol.v,y); s.length_path := s.length_path + s.cora;
y := H(s.sol.v,s.sol.t);
y(y'last) := Create(0.0); nresa := Norm(y);
normv := Norm(s.sol.v);
if normv + 1.0 /= 1.0
then ncorr := ncora / normv; nresr := nresa / normv;
for jj in s.sol.v'range loop -- SCALE THE SOLUTION
s.sol.v(jj) := s.sol.v(jj)/Create(normv);
end loop;
end if;
cWrite(file,ncora,ncorr,nresa,nresr); -- WRITE STATISTICS
cWrite(file,s.rcond,s.sol.m);
nit := nit + 1;
stop := (((ncora > s.cora) and then (nresa > s.resa))
or else ((ncorr > s.corr) and then (ncorr > s.corr)));
-- STOP WHEN DIVERGENCE
s.cora := ncora; s.resa := nresa; s.corr := ncorr; s.resr := nresr;
exit when stop;
exit when (((s.cora <= c.epsax) or else (s.resa <= c.epsaf))
and then ((s.corr <= c.epsrx) or else (s.resr <= c.epsrf)));
end loop; -- STOP WHEN DESIRED PRECISION REACHED
s.niter := s.niter + nit; s.nsyst := s.nsyst + nit;
j := dH(s.sol.v,s.sol.t);
lufco(j,y'last,ipvt,s.rcond); -- ESTIMATE CONDITION NUMBER
exception
when numeric_error => return;
end Projective_Single_Severe_Conditioned_Reporting_Corrector;
procedure Projective_Multiple_Loose_Normal_Silent_Corrector
( s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
procedure Single_Corrector is
new Projective_Single_Loose_Normal_Silent_Corrector(Norm,H,dH);
procedure Corrector is new Multiple_Silent_Corrector(Single_Corrector);
begin
Corrector(s,pivot,dist_sols,c,fail);
end Projective_Multiple_Loose_Normal_Silent_Corrector;
procedure Projective_Multiple_Loose_Normal_Reporting_Corrector
( file : in file_type; s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
procedure Single_Corrector is
new Projective_Single_Loose_Normal_Reporting_Corrector(Norm,H,dH);
procedure Corrector is new Multiple_Reporting_Corrector(Single_Corrector);
begin
Corrector(file,s,pivot,dist_sols,c,fail);
end Projective_Multiple_Loose_Normal_Reporting_Corrector;
procedure Projective_Multiple_Loose_Conditioned_Silent_Corrector
( s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
procedure Single_Corrector is
new Projective_Single_Loose_Conditioned_Silent_Corrector(Norm,H,dH);
procedure Corrector is new Multiple_Silent_Corrector(Single_Corrector);
begin
Corrector(s,pivot,dist_sols,c,fail);
end Projective_Multiple_Loose_Conditioned_Silent_Corrector;
procedure Projective_Multiple_Loose_Conditioned_Reporting_Corrector
( file : in file_type; s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
procedure Single_Corrector is
new Projective_Single_Loose_Normal_Reporting_Corrector(Norm,H,dH);
procedure Corrector is new Multiple_Reporting_Corrector(Single_Corrector);
begin
Corrector(file,s,pivot,dist_sols,c,fail);
end Projective_Multiple_Loose_Conditioned_Reporting_Corrector;
procedure Projective_Multiple_Severe_Normal_Silent_Corrector
( s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
procedure Single_Corrector is
new Projective_Single_Severe_Normal_Silent_Corrector(Norm,H,dH);
procedure Corrector is new Multiple_Silent_Corrector(Single_Corrector);
begin
Corrector(s,pivot,dist_sols,c,fail);
end Projective_Multiple_Severe_Normal_Silent_Corrector;
procedure Projective_Multiple_Severe_Normal_Reporting_Corrector
( file : in file_type; s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
procedure Single_Corrector is
new Projective_Single_Loose_Normal_Reporting_Corrector(Norm,H,dH);
procedure Corrector is new Multiple_Reporting_Corrector(Single_Corrector);
begin
Corrector(file,s,pivot,dist_sols,c,fail);
end Projective_Multiple_Severe_Normal_Reporting_Corrector;
procedure Projective_Multiple_Severe_Conditioned_Silent_Corrector
( s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
procedure Single_Corrector is
new Projective_Single_Severe_Conditioned_Silent_Corrector(Norm,H,dH);
procedure Corrector is new Multiple_Silent_Corrector(Single_Corrector);
begin
Corrector(s,pivot,dist_sols,c,fail);
end Projective_Multiple_Severe_Conditioned_Silent_Corrector;
procedure Projective_Multiple_Severe_Conditioned_Reporting_Corrector
( file : in file_type; s : in out Solu_Info_Array;
pivot : in out natural; dist_sols : in double_float;
c : in Corr_Pars; fail : out boolean ) is
procedure Single_Corrector is
new Projective_Single_Severe_Conditioned_Reporting_Corrector(Norm,H,dH);
procedure Corrector is new Multiple_Reporting_Corrector(Single_Corrector);
begin
Corrector(file,s,pivot,dist_sols,c,fail);
end Projective_Multiple_Severe_Conditioned_Reporting_Corrector;
end Correctors;