Return to vlprs_algorithm.adb CVS log

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

Annotation of OpenXM_contrib/PHC/Ada/Continuation/vlprs_algorithm.adb, Revision 1.1

1.1     ! maekawa     1: with vLpRs_Tables;                       use vLpRs_Tables;
        !             2: with Standard_Floating_Numbers_io;       use Standard_Floating_Numbers_io;
        !             3: with Standard_Floating_Matrices_io;      use Standard_Floating_Matrices_io;
        !             4:
        !             5: package body vLpRs_Algorithm is
        !             6:
        !             7: -- OUTPUT ROUTINES OF ERROR :
        !             8:
        !             9:   procedure Write_Init ( file : in file_type; s,l,v : in Vector ) is
        !            10:
        !            11:   -- DESCRIPTION :
        !            12:   --   Writes the beginning of the error table.
        !            13:
        !            14:     tmp,err : double_float;
        !            15:     w : integer;
        !            16:
        !            17:   begin
        !            18:     tmp := v(1)/l(1);
        !            19:     w := integer(tmp);
        !            20:     err := abs(tmp - double_float(w));
        !            21:     put(file,s(0),2,3,3); new_line(file);
        !            22:     put(file,s(1),2,3,3); put(file,err,2,3,3); new_line(file);
        !            23:   end Write_Init;
        !            24:
        !            25:   procedure Write ( file : in file_type; k : in natural;
        !            26:                     s : in double_float; l,v : in Vector ) is
        !            27:
        !            28:   -- DESCRIPTION :
        !            29:   --   Writes an additional row of k columns of the error table.
        !            30:   --   Note that when m is wrong, the outcome is no longer integer.
        !            31:
        !            32:     tmp,err : double_float;
        !            33:     w : integer;
        !            34:
        !            35:   begin
        !            36:     put(file,s,2,3,3);
        !            37:    -- w := v(k)/l(k); --
        !            38:     w := integer(v(k)/l(k));
        !            39:     for i in 1..k loop
        !            40:       tmp := v(i)/l(i);
        !            41:      -- err := abs(tmp-w); --
        !            42:       err := abs(tmp - double_float(w));
        !            43:       put(file,err,2,3,3);
        !            44:     end loop;
        !            45:     new_line(file);
        !            46:   end Write;
        !            47:
        !            48: -- TARGET ROUTINES :
        !            49:
        !            50:   procedure vLpRs_full
        !            51:                 ( r : in natural; s,logs,logx : in Vector;
        !            52:                   srp,dsp,p,l,v : in out Vector; rt1,rt2 : in out Matrix ) is
        !            53:   begin
        !            54:     vL_full(s(0..r),logs(0..r),logx(0..r),srp,dsp,p,l,v,rt1,rt2);
        !            55:     rt1 := rt2;
        !            56:     for k in r+1..s'last loop
        !            57:       vlprs_pipe(s(k),logs(k),logx(k),srp,dsp,p,l,v,rt1,rt2);
        !            58:     end loop;
        !            59:   end vLpRs_full;
        !            60:
        !            61:   procedure vLpRs_pipe
        !            62:                 ( file : in file_type; r : in natural;
        !            63:                   s,logs,logx : in Vector; srp,dsp,p,l,v : in out Vector;
        !            64:                   rt1,rt2 : in out Matrix ) is
        !            65:   begin
        !            66:     p(0) := 1.0;                                             -- initialization
        !            67:     v(0..1) := logx(0..1);
        !            68:     l(0..1) := logs(0..1);
        !            69:     L_pipe(l(0..1),p(0..0),logs(1));
        !            70:     v_pipe(v(0..1),p(0..0),logx(1));
        !            71:    -- Write_Init(file,s,l,v);                           -- write the error table
        !            72:     for k in 2..r loop
        !            73:       p_full(s(0..k),srp(1..k-1),dsp(1..k-1),p(0..k-1),rt1,rt2);
        !            74:       L_pipe(l(0..k),p(0..k-1),logs(k));                        -- extrapolate
        !            75:       v_pipe(v(0..k),p(0..k-1),logx(k));
        !            76:      -- Write(file,k,s(k),l,v);                         -- write the error table
        !            77:      -- put_line(file,"rt2 :"); put(file,rt2,3,3,3);
        !            78:     end loop;
        !            79:     rt1 := rt2;
        !            80:     for k in r+1..s'last loop
        !            81:       vlprs_pipe(file,s(k),logs(k),logx(k),srp,dsp,p,l,v,rt1,rt2);
        !            82:      -- put_line(file,"rt2 : "); put(file,rt2,3,3,3);
        !            83:     end loop;
        !            84:   end vLpRs_pipe;
        !            85:
        !            86:   procedure vLpRs_pipe
        !            87:                  ( s,logs,logx : in double_float;
        !            88:                    srp,dsp,p,l,v : in out Vector; rt1,rt2 : in out Matrix ) is
        !            89:   begin
        !            90:     s_pipe(srp,s,dsp);
        !            91:     RR_pipe(rt1,dsp,p,rt2);
        !            92:     p_pipe(rt1,rt2,p);  rt1 := rt2;
        !            93:     L_pipe(l,p,logs);
        !            94:     v_pipe(v,p,logx);
        !            95:   end vLpRs_pipe;
        !            96:
        !            97:   procedure vLpRs_pipe
        !            98:                  ( file : in file_type; s,logs,logx : in double_float;
        !            99:                    srp,dsp,p,l,v : in out Vector; rt1,rt2 : in out Matrix ) is
        !           100:   begin
        !           101:     vLpRs_pipe(s,logs,logx,srp,dsp,p,l,v,rt1,rt2);
        !           102:    -- Write(file,l'last,s,l,v);
        !           103:   end vLpRs_pipe;
        !           104:
        !           105: end vLpRs_Algorithm;