#include <stdio.h>
#include <stdlib.h>
#include <assert.h>
#include <limits.h>
#include <math.h>
#include <complex.h>
#include <string.h>

#include "a7p2-0-base.c"

#define MaxNumElts 30000
Elt_t Elts[MaxNumElts];

int main(int argc,char *(argv[])) {
  int Ix,NumElts,NumPts;
  Pt_t ZeroPt={0,0},MPt;

  EltsFile=fopen(argv[1],"r");
  assert(EltsFile != NULL);

  NumElts=ReadElts(Elts,MaxNumElts);

  EltsFile=fopen("tmpPts","w");
  assert(EltsFile != NULL);

  NumPts=MIN(MaxNumPts,NumElts);
  for(Ix=0;Ix<NumPts;Ix++)
    memcpy(&Pts[Ix],&Elts[Ix].p,sizeof(Pt_t));

  for(Ix=0;Ix<NumElts;Ix++) {
    Midpoint(&MPt,ZeroPt,Elts[Ix].p);
    if(InDirichlet(MPt),NumPts)OutputPt(Elts[Ix].p);
  }
}