File: [local] / OpenXM_contrib / PHC / Ada / Continuation / vlprs_algorithm.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 vLpRs_Tables; use vLpRs_Tables;
with Standard_Floating_Numbers_io; use Standard_Floating_Numbers_io;
with Standard_Floating_Matrices_io; use Standard_Floating_Matrices_io;
package body vLpRs_Algorithm is
-- OUTPUT ROUTINES OF ERROR :
procedure Write_Init ( file : in file_type; s,l,v : in Vector ) is
-- DESCRIPTION :
-- Writes the beginning of the error table.
tmp,err : double_float;
w : integer;
begin
tmp := v(1)/l(1);
w := integer(tmp);
err := abs(tmp - double_float(w));
put(file,s(0),2,3,3); new_line(file);
put(file,s(1),2,3,3); put(file,err,2,3,3); new_line(file);
end Write_Init;
procedure Write ( file : in file_type; k : in natural;
s : in double_float; l,v : in Vector ) is
-- DESCRIPTION :
-- Writes an additional row of k columns of the error table.
-- Note that when m is wrong, the outcome is no longer integer.
tmp,err : double_float;
w : integer;
begin
put(file,s,2,3,3);
-- w := v(k)/l(k); --
w := integer(v(k)/l(k));
for i in 1..k loop
tmp := v(i)/l(i);
-- err := abs(tmp-w); --
err := abs(tmp - double_float(w));
put(file,err,2,3,3);
end loop;
new_line(file);
end Write;
-- TARGET ROUTINES :
procedure vLpRs_full
( r : in natural; s,logs,logx : in Vector;
srp,dsp,p,l,v : in out Vector; rt1,rt2 : in out Matrix ) is
begin
vL_full(s(0..r),logs(0..r),logx(0..r),srp,dsp,p,l,v,rt1,rt2);
rt1 := rt2;
for k in r+1..s'last loop
vlprs_pipe(s(k),logs(k),logx(k),srp,dsp,p,l,v,rt1,rt2);
end loop;
end vLpRs_full;
procedure vLpRs_pipe
( file : in file_type; r : in natural;
s,logs,logx : in Vector; srp,dsp,p,l,v : in out Vector;
rt1,rt2 : in out Matrix ) is
begin
p(0) := 1.0; -- initialization
v(0..1) := logx(0..1);
l(0..1) := logs(0..1);
L_pipe(l(0..1),p(0..0),logs(1));
v_pipe(v(0..1),p(0..0),logx(1));
-- Write_Init(file,s,l,v); -- write the error table
for k in 2..r loop
p_full(s(0..k),srp(1..k-1),dsp(1..k-1),p(0..k-1),rt1,rt2);
L_pipe(l(0..k),p(0..k-1),logs(k)); -- extrapolate
v_pipe(v(0..k),p(0..k-1),logx(k));
-- Write(file,k,s(k),l,v); -- write the error table
-- put_line(file,"rt2 :"); put(file,rt2,3,3,3);
end loop;
rt1 := rt2;
for k in r+1..s'last loop
vlprs_pipe(file,s(k),logs(k),logx(k),srp,dsp,p,l,v,rt1,rt2);
-- put_line(file,"rt2 : "); put(file,rt2,3,3,3);
end loop;
end vLpRs_pipe;
procedure vLpRs_pipe
( s,logs,logx : in double_float;
srp,dsp,p,l,v : in out Vector; rt1,rt2 : in out Matrix ) is
begin
s_pipe(srp,s,dsp);
RR_pipe(rt1,dsp,p,rt2);
p_pipe(rt1,rt2,p); rt1 := rt2;
L_pipe(l,p,logs);
v_pipe(v,p,logx);
end vLpRs_pipe;
procedure vLpRs_pipe
( file : in file_type; s,logs,logx : in double_float;
srp,dsp,p,l,v : in out Vector; rt1,rt2 : in out Matrix ) is
begin
vLpRs_pipe(s,logs,logx,srp,dsp,p,l,v,rt1,rt2);
-- Write(file,l'last,s,l,v);
end vLpRs_pipe;
end vLpRs_Algorithm;