// primesg2.c  generate a file of primes

#include <stdio.h>

static long int primes[1000001];
static int np = 0;    // initially, first prime is 2 at index zero
static int numb = 1000000;

int main(void)
{
  FILE * outp;
  int i, j, k, kk, ix, icount, sum;
  long fct;
  static int kp = 10;

  printf("primesg2.c running \n");
  primes[0] = 2;
  printf("%d th prime is %ld \n", 1, primes[0]);
  for(i=3; np<numb; i+=2) // only odd after 2
  {
    for(j=0; j<=np; j++)
      if(i%primes[j] == 0) goto non;

    np++;
    primes[np] = i; // build list of primes
    if(np+1==kp)
    {
       printf("the %dth  prime is %d \n", kp, i);
       if(kp<100000) kp = kp*10;
       else          kp = kp+100000;
    }
    non:;
  }
  printf("%d th prime is %ld \n", np+1, primes[np]);

  outp = fopen("primef2.h","w");
  fprintf(outp, "// primef.h\n");
  fprintf(outp, "static int np = %d; \n", np);
  fprintf(outp, "static long int primes[]={ \n");
  for(k=0; k<numb-1; k++)
  {
    fprintf(outp,"%ldl, ", primes[k]);
    if(k%7==0) fprintf(outp,"\n");
  }
  fprintf(outp,"%ldl};\n", primes[k]);

  printf("primesg2.c finished, primef2.h written \n");
  return 0;
} // end psimeq2