prova

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <float.h>

#define FILE_PUNTI  "hits.bin"
#define FILE_PUNTI2 "nostro.bin"
#define FILE_RETTE  "rette.bin"
#define CONST       2.829

typedef struct _dot
  {
    int x,y;
  }dot;

typedef struct _retta
  {
    float a,b,c;
  }retta;

typedef float* pfloat;
typedef dot*   pdot;
typedef retta* pretta;
typedef FILE*  pfile;

int main (void)
  {
    int    i, j, n_r,n_d;
    dot    d;
    float  a,b,c,min1=FLT_MAX,min2=FLT_MAX;
    pfloat p_dist;
    pretta p_r;
    pfile  pf,pf2;

   /*Immagino che i punti siano troppi per essere acquisiti con un fread tutti in una volta*/
    if((pf=fopen(FILE_RETTE,"rb"))==NULL)                      goto fopen_err;
    if((fread(&n_r,sizeof(int),1,pf))!=1)                      goto fread_err;
    if((p_r=(pretta)malloc(sizeof(retta)*n_r))==NULL)          goto mem_err;
    if((fread(p_r,sizeof(retta),n_r,pf))!=n_r)                 goto fread_err;
    fclose(pf);
/*
    for(i=0;i<n_r;i++)
      {
        printf("La retta %d e': (%f)x+(%f)y+(%f)=0\n",(i+1),p_r[i].a,p_r[i].b,p_r[i].c);
      }
*/
    if((pf=fopen(FILE_PUNTI,"rb"))==NULL)                      goto fopen_err;
    if((fread(&n_d,sizeof(int),1,pf))!=1)                      goto fread_err;
    if((pf2=fopen(FILE_PUNTI2,"wb"))==NULL)                    goto fopen_err;
    for (i=0;i<n_d;i++)
      {
        if((fread(&d,sizeof(dot),1,pf))!=1)                    goto fread_err;
        fwrite(&d,sizeof(dot),1,pf2);
      }

    fclose(pf);
    fclose(pf2);

    if((p_dist=(pfloat)malloc(sizeof(float)*n_r))==NULL)       goto mem_err;

    for(i=0;i<n_r;i++)
      p_dist[i]=0;

    printf("Ti stampero' le rette la cui distanza media dai punti sia < %f\n",CONST);

    for (i=0;i<n_r;i++)
      {
        if((pf=fopen(FILE_PUNTI2,"rb"))==NULL)                 goto fopen_err;
        a=p_r[i].a; /*vengono richiamate tante volte*/
        b=p_r[i].b;
        c=p_r[i].c;
        for (j=0;j<n_d;j++)
          {
            if((fread(&d,sizeof(dot),1,pf))!=1)                goto fread_err;
            p_dist[i]+=fabs((a*d.x)+(b*d.y)+c)/(sqrt((a*a)+(b*b)));
          }
        p_dist[i]/=n_d;  /*media distanze*/
        if(p_dist[i]<CONST)
          printf ("La retta %d°: (%f)x+(%f)y+(%f)=0  con distanza %f\n",(i+1),a,b,c,p_dist[i]);
        fclose(pf);
      }

    for (i=0;i<n_r; i++)
      {
        if(p_dist[i]<min1)
          {
            min2=min1;
            min1=p_dist[i];
          }
      }

    printf ("Il primo minimo e' %f\nIl secondo minimo e'%f\n",min1,min2);
    free(p_r);
    free(p_dist);
    return EXIT_SUCCESS;

    fopen_err: printf("C'e' un problema nell'apertura del file.\n");              return EXIT_FAILURE;
    fread_err: printf("Il numero acquisito non coincide con quello inserito.\n"); return EXIT_FAILURE;
    mem_err:   printf("Non c'e' abbastanza spazio in memoria.\n");                return EXIT_FAILURE;

  }