This commit is contained in:
Rodolfo Cuevas 2018-12-12 12:33:06 -03:00
parent 00444c46f2
commit b1c3c1aa31
1 changed files with 17 additions and 58 deletions

View File

@ -39,6 +39,20 @@ float min(float x, float y){ // Función para encontrar el mayor entre dos valo
return (x < y)? x : y; return (x < y)? x : y;
} }
float stripClosest(point_t strip[], int size, float d) // Función para encontrar la distancia entre los puntos más cerca del arreglo dado
{
float min = d; // inicializa en la distancia minima d
qsort(strip, size, sizeof(point_t), compareY);
for (int i = 0; i < size; ++i)
for (int j = i+1; j < size && (strip[j].y - strip[i].y) < min; ++j)
if (distance(strip[i],strip[j]) < min)
min = dist(strip[i], strip[j]);
return min;
}
float closestUtil(point_t P[], int n){ // Función para encontrar la distancia más corta entre puntos float closestUtil(point_t P[], int n){ // Función para encontrar la distancia más corta entre puntos
int mid = n/2; int mid = n/2;
point_t midpoint_t = P[mid]; point_t midpoint_t = P[mid];
@ -62,66 +76,11 @@ float closestUtil(point_t P[], int n){ // Función para encontrar la distancia m
distance es strip[]*/ distance es strip[]*/
return min(d, stripClosest(strip, j, d) ); return min(d, stripClosest(strip, j, d) );
} }
float stripClosest(point_t strip[], int size, float d) // Función para encontrar la distancia entre los puntos más cerca del arreglo dado
{
float min = d; // inicializa en la distancia minima d
qsort(strip, size, sizeof(point_t), compareY);
for (int i = 0; i < size; ++i)
for (int j = i+1; j < size && (strip[j].y - strip[i].y) < min; ++j)
if (distance(strip[i],strip[j]) < min)
min = dist(strip[i], strip[j]);
return min;
}
point_t * divide_and_conquer(point_t *points, unsigned int n, double *minimum_dist) { point_t * divide_and_conquer(point_t *points, unsigned int n, double *minimum_dist) {
point_t *closest_pair = malloc(sizeof(point_t) * 2); point_t *closest_pair = malloc(sizeof(point_t) * 2);
qsort(P, n, sizeof(point_t), compareX);
return closestUtil(P, n); //Usa la funcion recursiva closestUtil para encontrar la distancia minima
}
float closest(point_t P[], int n)
{
qsort(P, n, sizeof(point_t), compareX);
return closestUtil(P, n); /*Uso recursivo de la funcion closestUtil() para encontrar la distancia más pequeña*/
}
return closest_pair;
}
float closestUtil(point_t P[], int n) /*Función recursiva para encontrar la distancia más pequeña. El arreglo
P contiene todos los puntos ordenados respecto a la cordenada X */
{
int mid = n/2;
point_t midpoint_t = P[mid];
/*Considerando la linea vertical que pasa a travez del punto medio, calcula el distanca mas corta en dl
de la izquierda de la mitad y dr en el lado derecho*/
float dl = closestUtil(P, mid);
float dr = closestUtil(P + mid, n-mid);
// Find the smaller of two distances
float d = min(dl, dr);
// Build an array strip[] that contains point_ts close (closer than d)
// to the line passing through the middle point_t
point_t strip[n];
int j = 0;
for (int i = 0; i < n; i++)
if (abs(P[i].x - midpoint_t.x) < d)
strip[j] = P[i], j++;
// Find the closest point_ts in strip. Return the minimum of d and closest
// distance is strip[]
return min(d, stripClosest(strip, j, d) );
}
// The main functin that finds the smallest distance
// This method mainly uses closestUtil()
float closest(point_t P[], int n)
{
qsort(P, n, sizeof(point_t), compareX);
// Use recursive function closestUtil() to find the smallest distance
return closestUtil(P, n);
}