Return to multprec_residual_evaluations.adb CVS log

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

Annotation of OpenXM_contrib/PHC/Ada/Continuation/multprec_residual_evaluations.adb, Revision 1.1.1.1

1.1       maekawa     1: with integer_io;                         use integer_io;
                      2: with Multprec_Floating_Numbers_io;       use Multprec_Floating_Numbers_io;
                      3: with Multprec_Complex_Norms_Equals;      use Multprec_Complex_Norms_Equals;
                      4:
                      5: package body Multprec_Residual_Evaluations is
                      6:
                      7:   function Residual ( p_eval : Eval_Poly_Sys; zero : Vector )
                      8:                     return Floating_Number is
                      9:
                     10:     res : Floating_Number;
                     11:     eva : Vector(p_eval'range) := Eval(p_eval,zero);
                     12:
                     13:   begin
                     14:     res := Max_Norm(eva);
                     15:     Clear(eva);
                     16:     return res;
                     17:   end Residual;
                     18:
                     19:   procedure Residual ( file : in file_type;
                     20:                        p_eval : in Eval_Poly_Sys; sol : in Solution ) is
                     21:
                     22:     res : Floating_Number := Residual(p_eval,sol.v);
                     23:
                     24:   begin
                     25:     put(file,res);
                     26:     Clear(res);
                     27:   end Residual;
                     28:
                     29:   procedure Residuals ( file : in file_type;
                     30:                         p_eval : in Eval_Poly_Sys; sols : in Solution_List ) is
                     31:
                     32:     tmp : Solution_List := sols;
                     33:     ls : Link_to_Solution;
                     34:
                     35:   begin
                     36:     for i in 1..Length_Of(sols) loop
                     37:       ls := Head_Of(tmp);
                     38:       put(file,"residual "); put(file,i,1); put(file," : ");
                     39:       Residual(file,p_eval,ls.all);
                     40:       new_line(file);
                     41:       tmp := Tail_Of(tmp);
                     42:     end loop;
                     43:   end Residuals;
                     44:
                     45: end Multprec_Residual_Evaluations;