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


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

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

{

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

}

##
\q

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

 /* Example output (with default options, for conductor less than 100):
  *
  * 14A     -996    ok      (odd=92, N/A=100, good=113/113)
  * 15A     -996    ok      (odd=95, N/A=100, good=109/109)
  * 21A     -996    ok      (odd=101, N/A=96, good=105/105)
  * 26A     -996    ok      (odd=96, N/A=102, good=107/107)
  * 26B     -996    ok      (odd=94, N/A=99, good=109/109)
  * 30A     -996    ok      (odd=63, N/A=167, good=75/75)
  * 33A     -996    ok      (odd=105, N/A=97, good=102/102)
  * 34A     -996    ok      (odd=96, N/A=101, good=107/107)
  * 35A     -996    ok      (odd=116, N/A=90, good=99/99)
  * 38A     -996    ok      (odd=97, N/A=96, good=104/104)
  * 38B     -996    ok      (odd=100, N/A=98, good=107/107)
  * 39A     -996    ok      (odd=108, N/A=96, good=101/101)
  * 42A     -996    ok      (odd=67, N/A=166, good=69/69)
  * 46A     -996    ok      (odd=96, N/A=104, good=103/103)
  * 51A     -996    ok      (odd=105, N/A=94, good=103/103)
  * 55A     -996    ok      (odd=116, N/A=90, good=96/96)
  * 57B     -996    ok      (odd=106, N/A=98, good=101/101)
  * 57C     -996    ok      (odd=106, N/A=98, good=101/101)
  * 58B     -996    ok      (odd=95, N/A=102, good=107/107)
  * 62A     -996    ok      (odd=104, N/A=93, good=107/107)
  * 66A     -996    ok      (odd=69, N/A=166, good=70/70)
  * 66B     -996    ok      (odd=69, N/A=167, good=69/69)
  * 66C     -996    ok      (odd=67, N/A=167, good=67/67)
  * 69A     -996    ok      (odd=110, N/A=94, good=101/101)
  * 70A     -996    ok      (odd=71, N/A=160, good=68/68)
  * 77B     -996    ok      (odd=128, N/A=87, good=90/90)
  * 77C     -996    ok      (odd=125, N/A=85, good=90/90)
  * 78A     -996    ok      (odd=71, N/A=167, good=67/67)
  * 85A     -996    ok      (odd=120, N/A=90, good=95/95)
  * 94A     -996    ok      (odd=95, N/A=98, good=109/109)
  *
  * This took about 100 minutes on a Pentium 4 2.60GHz
  * Note that practically all the time is spent in elllseries
  *
  */

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