/************************************************************/
/* 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                          
 *
 */

 /* Example output (with default options):
  *
  * 11A     997     ok (max diff = 2.403185862 E-26)
  * 19A     997     ok (max diff = 2.342601344 E-26)
  * 43A     997     ok (max diff = 1.615587134 E-25)
  * 67A     997     ok (max diff = 1.431410201 E-23)
  * 79A     997     ok (max diff = 2.342601344 E-26)
  * 83A     997     ok (max diff = 1.171300672 E-26)
  * 139A    997     ok (max diff = 3.453317499 E-26)
  * 163A    997     ok (max diff = 3.190784590 E-26)
  *
  * This took about 40 minutes in a Pentium 4 2.60GHz !! (5 min each)
  * Note that practically all the time is spent in elllseries
  * 
  */

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

\r twist.gp
\r data_ell.gp
\r data2.gp
\p 28

{

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

}

\q

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