/************************************************************/
/* Author : Gonzalo Tornaria <mailto:tornaria@math.utexas.edu>
 * Date   : December 2003
 *
 *
 * Verification of the formula for real twists
 * Set precision, bound, and epsilon, and run 
 *
 *    gp -q verify_real.gp                          
 *
 */

/************************************************************/

read("data_real.gp")
\p 28
bound=1000; epsilon=1e-20;
read("verify_default.gp");

{

for(i=1,length(tw_real),
  print1(tw_name(tw_real[i]));
  verify_twist_r(tw_real[i],bound,epsilon);
);

}

##
\q

/************************************************************/

 /* Example output (with default options):
  *
  * 11A     997     ok      (odd=140, N/A=0, good=162/162)
  * 19A     997     ok      (odd=145, N/A=0, good=158/158)
  * 43A     997     ok      (odd=144, N/A=0, good=159/159)
  * 67A     997     ok      (odd=152, N/A=0, good=151/151)
  * 79A     997     ok      (odd=143, N/A=0, good=158/158)
  * 83A     997     ok      (odd=141, N/A=0, good=161/161)
  * 131A    997     ok      (odd=136, N/A=0, ZERO=166/166)
  * 139A    997     ok      (odd=157, N/A=0, good=146/146)
  * 163A    997     ok      (odd=144, N/A=0, good=159/159)
  * 179A    997     ok      (odd=157, N/A=0, good=143/143)
  *
  * This took about 50 minutes in a Pentium 4 2.60GHz !!
  * Note that practically all the time is spent in elllseries
  * 
  */

/************************************************************/