/*-
 * Written by H. Mitasova, I. Kosinovsky, D. Gerdes Fall 1993
 * US Army Construction Engineering Research Lab
 * Copyright 1993, H. Mitasova (University of Illinois),
 * I. Kosinovsky, (USA-CERL), and D.Gerdes (USA-CERL)
 *
 * modified by McCauley in August 1995
 * modified by Mitasova in August 1995
 *
 */

#include <stdio.h>
#include <math.h>
#include <unistd.h>
#include <grass/gis.h>
#include <grass/interpf.h>

double IL_dist_square(double *pt1, double *pt2, int dim)
{
    int i;
    double sum = 0, s;

    for (i = 0; i < dim; i++) {
        s = pt1[i] - pt2[i];
        sum += s * s;
    }
    return sum;
}