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

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

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

{

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

}

##
\q

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

 /* Example output (with default options, for conductor less than 500):
  *
  * 11A     -996    ok      (odd=141, N/A=0, good=163/163)
  * 17A     -996    ok      (odd=145, N/A=0, good=160/160)
  * 19A     -996    ok      (odd=146, N/A=0, good=158/158)
  * 37B     -996    ok      (odd=148, N/A=0, good=155/155)
  * 67A     -996    ok      (odd=147, N/A=0, good=155/155)
  * 73A     -996    ok      (odd=150, N/A=0, good=155/155)
  * 89B     -996    ok      (odd=152, N/A=0, good=151/151)
  * 109A    -996    ok      (odd=149, N/A=0, good=155/155)
  * 113A    -996    ok      (odd=152, N/A=0, good=153/153)
  * 139A    -996    ok      (odd=147, N/A=0, good=154/154)
  * 179A    -996    ok      (odd=156, N/A=0, good=149/149)
  * 233A    -996    ok      (odd=151, N/A=0, good=152/152)
  * 307A    -996    ok      (odd=146, N/A=0, good=157/157)
  * 307B    -996    ok      (odd=146, N/A=0, good=157/157)
  * 307C    -996    ok      (odd=146, N/A=0, good=157/157)
  * 307D    -996    ok      (odd=146, N/A=0, good=157/157)
  * 353A    -996    ok      (odd=148, N/A=0, good=157/157)
  * 389A    -996    ok      (odd=149, N/A=0, ZERO=155/155)
  * 431B    -996    ok      (odd=154, N/A=0, good=148/148)
  * 443C    -996    ok      (odd=148, N/A=0, good=156/156)
  *
  * This took about 165 minutes in a Pentium 4 2.60GHz !!
  * Note that practically all the time is spent in elllseries
  * 
  */

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