/************************************************************/
/* Author : Gonzalo Tornaria <mailto:tornaria@math.utexas.edu>
 * Date   : December 2003
 */

/* Verification of the formula for imaginary twists
 * Set precision, bound, and epsilon, and run 
 *
 *    gp -q verify_imag.gp                          
 *
 */

 /* Example output (with default options):
  *
  * 11A     -996    ok (max diff = 1.837730365 E-26)
  * 17A     -996    ok (max diff = 9.84498410 E-27)
  * 19A     -996    ok (max diff = 3.251369107 E-26)
  * 37B     -996    ok (max diff = 4.503449136 E-26)
  * 67A     -996    ok (max diff = 3.441200596 E-24)
  * 73A     -996    ok (max diff = 9.39060022 E-27)
  * 89B     -996    ok (max diff = 1.837730365 E-26)
  * 109A    -996    ok (max diff = 1.841163488 E-24)
  * 113A    -996    ok (max diff = 2.837374904 E-26)
  * 139A    -996    ok (max diff = 3.796629765 E-26)
  * 179A    -996    ok (max diff = 1.288430740 E-25)
  * 233A    -996    ok (max diff = 1.187456544 E-25)
  *
  * This took about 1 hour 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 data1.gp
\p 28

{

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

}

\q

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