Annotation of OpenXM_contrib/PHC/Ada/Continuation/ts_rootrefi.adb, Revision 1.1
1.1 ! maekawa 1: with text_io,integer_io; use text_io,integer_io;
! 2: with Communications_with_User; use Communications_with_User;
! 3: with Standard_Floating_Numbers; use Standard_Floating_Numbers;
! 4: with Standard_Floating_Numbers_io; use Standard_Floating_Numbers_io;
! 5: with Standard_Complex_Poly_Systems; use Standard_Complex_Poly_Systems;
! 6: with Standard_Complex_Poly_Systems_io; use Standard_Complex_Poly_Systems_io;
! 7: with Standard_Complex_Solutions; use Standard_Complex_Solutions;
! 8: with Standard_Complex_Solutions_io; use Standard_Complex_Solutions_io;
! 9: with Standard_Root_Refiners; use Standard_Root_Refiners;
! 10: with Standard_to_Multprec_Convertors; use Standard_to_Multprec_Convertors;
! 11: with Multprec_Floating_Numbers; use Multprec_Floating_Numbers;
! 12: with Multprec_Floating_Numbers_io; use Multprec_Floating_Numbers_io;
! 13: with Multprec_Complex_Poly_Systems; use Multprec_Complex_Poly_Systems;
! 14: with Multprec_Complex_Poly_SysFun; use Multprec_Complex_Poly_SysFun;
! 15: with Multprec_Complex_Solutions; use Multprec_Complex_Solutions;
! 16: with Multprec_Complex_Solutions_io; use Multprec_Complex_Solutions_io;
! 17: with Multprec_Residual_Evaluations; use Multprec_Residual_Evaluations;
! 18: with Multprec_Root_Refiners; use Multprec_Root_Refiners;
! 19:
! 20: procedure ts_rootrefi is
! 21:
! 22: -- DESCRIPTION :
! 23: -- This routine facilitates interactive testing of the root refiners.
! 24:
! 25: procedure Call_Standard_Root_Refiner
! 26: ( file : in file_type;
! 27: p : in Standard_Complex_Poly_Systems.Poly_Sys;
! 28: sols : in out Standard_Complex_Solutions.Solution_List ) is
! 29:
! 30: epsxa,epsfa,tolsing : double_float;
! 31: numit : natural;
! 32: max : constant natural := 5;
! 33:
! 34: begin
! 35: epsxa := 1.0E-14;
! 36: epsfa := 1.0E-14;
! 37: tolsing := 1.0E-08;
! 38: Reporting_Root_Refiner(file,p,sols,epsxa,epsfa,tolsing,numit,max,true);
! 39: end Call_Standard_Root_Refiner;
! 40:
! 41: procedure Call_Multprec_Root_Refiner
! 42: ( file : in file_type;
! 43: p : in Multprec_Complex_Poly_Systems.Poly_Sys;
! 44: sols : in out Multprec_Complex_Solutions.Solution_List ) is
! 45:
! 46: epsxa,epsfa,tolsing : Floating_Number;
! 47: numit,deci,size : natural;
! 48: max : constant natural := 5;
! 49:
! 50: begin
! 51: put("Give the number of decimal places : "); get(deci);
! 52: size := Decimal_to_Size(deci);
! 53: put("The size of the numbers : "); put(size,1); new_line;
! 54: Set_Size(sols,size);
! 55: put("Give tolerance for error : "); get(epsxa);
! 56: put("Give tolerance for residual : "); get(epsfa);
! 57: tolsing := Create(1.0E-08);
! 58: Reporting_Root_Refiner(file,p,sols,epsxa,epsfa,tolsing,numit,max,true);
! 59: end Call_Multprec_Root_Refiner;
! 60:
! 61: procedure Test_Standard_Root_Refiner is
! 62:
! 63: -- DESCRIPTION :
! 64: -- Test of root refining on list of solutions as standard vectors.
! 65:
! 66: file : file_type;
! 67: lp : Standard_Complex_Poly_Systems.Link_to_Poly_Sys;
! 68: sols : Standard_Complex_Solutions.Solution_List;
! 69:
! 70: begin
! 71: new_line;
! 72: put_line("Test on refining roots as standard complex vectors.");
! 73: new_line;
! 74: get(lp);
! 75: put_line("The system : "); put(lp.all);
! 76: new_line;
! 77: put_line("Reading the name of the output file.");
! 78: Read_Name_and_Create_File(file);
! 79: put(file,lp'last,lp.all);
! 80: new_line;
! 81: Read(sols);
! 82: new_line;
! 83: new_line(file);
! 84: put_line(file,"THE INITIAL SOLUTIONS : ");
! 85: put(file,sols);
! 86: Call_Standard_Root_Refiner(file,lp.all,sols);
! 87: end Test_Standard_Root_Refiner;
! 88:
! 89: procedure Test_Multprec_Root_Refiner is
! 90:
! 91: -- DESCRIPTION :
! 92: -- Test of root refining on list of solutions as standard vectors.
! 93:
! 94: file : file_type;
! 95: lp : Standard_Complex_Poly_Systems.Link_to_Poly_Sys;
! 96: stsols : Standard_Complex_Solutions.Solution_List;
! 97: mpsols : Multprec_Complex_Solutions.Solution_List;
! 98:
! 99: begin
! 100: new_line;
! 101: put_line("Test on refining roots as multi-precision complex vectors.");
! 102: new_line;
! 103: get(lp);
! 104: put_line("The system : "); put(lp.all);
! 105: new_line;
! 106: put_line("Reading the name of the output file.");
! 107: Read_Name_and_Create_File(file);
! 108: put(file,lp'last,lp.all);
! 109: new_line;
! 110: -- Read(stsols);
! 111: -- new_line(file);
! 112: -- put_line(file,"THE SOLUTION IN STANDARD PRECISION : ");
! 113: -- put(file,stsols);
! 114: -- mpsols := Create(stsols);
! 115: Read(mpsols);
! 116: new_line;
! 117: new_line(file);
! 118: put_line(file,"THE INITIAL SOLUTIONS : ");
! 119: put(file,Length_Of(mpsols),lp'last,mpsols);
! 120: declare
! 121: mp : Multprec_Complex_Poly_Systems.Poly_Sys(lp'range) := Convert(lp.all);
! 122: begin
! 123: Call_Multprec_Root_Refiner(file,mp,mpsols);
! 124: end;
! 125: end Test_Multprec_Root_Refiner;
! 126:
! 127: procedure Test_Multprec_Residual_Evaluator is
! 128:
! 129: file : file_type;
! 130: lp : Standard_Complex_Poly_Systems.Link_to_Poly_Sys;
! 131: mpsols : Multprec_Complex_Solutions.Solution_List;
! 132:
! 133: begin
! 134: new_line;
! 135: put_line("Test on evaluating residuals with multi-precision arithmetic.");
! 136: new_line;
! 137: get(lp);
! 138: put_line("The system : "); put(lp.all);
! 139: new_line;
! 140: put_line("Reading the name of the output file.");
! 141: Read_Name_and_Create_File(file);
! 142: put(file,lp'last,lp.all);
! 143: new_line;
! 144: Read(mpsols);
! 145: new_line(file);
! 146: put_line(file,"THE SOLUTIONS :");
! 147: put(file,Length_Of(mpsols),lp'last,mpsols);
! 148: declare
! 149: mp : Multprec_Complex_Poly_Systems.Poly_Sys(lp'range) := Convert(lp.all);
! 150: mp_eval : Multprec_Complex_Poly_SysFun.Eval_Poly_Sys(mp'range)
! 151: := Create(mp);
! 152: deci,size : natural;
! 153: begin
! 154: put("Give the number of decimal places : "); get(deci);
! 155: size := Decimal_to_Size(deci);
! 156: put("The size of the numbers : "); put(size,1); new_line;
! 157: Set_Size(mpsols,size);
! 158: put_line(file,"THE RESIDUALS :");
! 159: Residuals(file,mp_eval,mpsols);
! 160: end;
! 161: end Test_Multprec_Residual_Evaluator;
! 162:
! 163: procedure Main is
! 164:
! 165: ans : character;
! 166:
! 167: begin
! 168: new_line;
! 169: put_line("Interactive testing of root refiners.");
! 170: new_line;
! 171: put_line("Choose one of the following : ");
! 172: put_line(" 1. Test root refiner for standard complex numbers. ");
! 173: put_line(" 2. Test root refiner for multi-precision complex numbers.");
! 174: put_line(" 3. Evaluate residuals with multi-precision arithmetic.");
! 175: put("Type 0, 1, 2 or 3 to select : "); get(ans);
! 176: case ans is
! 177: when '1' => Test_Standard_Root_Refiner;
! 178: when '2' => Test_Multprec_Root_Refiner;
! 179: when '3' => Test_Multprec_Residual_Evaluator;
! 180: end case;
! 181: end Main;
! 182:
! 183: begin
! 184: Main;
! 185: end ts_rootrefi;
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to ts_rootrefi.adb CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / PHC / Ada / Continuation