/* Output from p2c, the Pascal-to-C translator */
/* From input file "ilink.p" */
/* This file contains some of the old nuclear family update routines */
/* used in a modified version of the ILINK program*/

#include "commondefs.h"
#include "ildefs.h"

Local Void getapprox(LINK)
struct LOC_seg *LINK;
{
  int first;
  double maxval;
  thisarray *WITH;
  int FORLIM;
  approxrec *WITH1;

  maxval = (*LINK->p)->gen->genarray[0];
  WITH = (*LINK->p)->gen;
  FORLIM = fgeno;
  for (first = 0; first < FORLIM; first++) {
    if (WITH->genarray[first] > maxval)
      maxval = WITH->genarray[first];
  }
  WITH1 = approxstruct;
  WITH = (*LINK->p)->gen;
  FORLIM = fgeno;
  for (first = 0; first < FORLIM; first++) {
    WITH1->approxarray[LINK->LINK->thisped - 1][first] =
      (WITH->genarray[first] > maxval * epsilon);
  }
  if (lasttime)
    return;
  WITH1 = approxstruct;
  WITH = (*LINK->p)->gen;
  FORLIM = fgeno;
  for (first = 0; first < FORLIM; first++) {
    if (!(WITH1->approxarray[LINK->LINK->thisped - 1][first]))
      WITH->genarray[first] = 0.0;
  }
}  /*getapprox*/


#include "comnuclear.c"