#include <stddef.h>
#include <float.h>

#define SIMD_SIZE 8

inline double sqr(double x)
{
    return x*x;
}

size_t foo(double tourx[], double toury[], size_t ncities, double ThisX, double ThisY, size_t j)
{
  size_t ClosePt=j;
  double CloseDist = DBL_MAX;
  for (; j<ncities; j++) {
      double ThisDist = sqr(tourx[j]-ThisX)+sqr(toury[j]-ThisY);
      if (ThisDist < CloseDist) {
        CloseDist = ThisDist;
        ClosePt = j;
      }
  }
  return ClosePt;
}