// VC6.cpp*: d�finit le point d'entr�e pour l'application console.
//

/clr
#include "stdafx.h"
#include <fstream>
#include <iomanip>
#include <iostream>
#include <cstring>
#include "stdio.h"
#include "stdlib.h"
#include <vcclr.h>
#include <vector>
//#include <string>

using namespace std;
using namespace System;

FILE *fp1;
//string ^str;
std::string;

//System::String^ str;

//Variable
char * prt_info2;
int i,nb_step_Y,nb_step_X,j,k,cc,dd,n;
long recl,lenr,numf,numx,numy;
float xmin,xmax,ymin,ymax;
//MAX=100;
//Struct
//typedef struct {
//float abscisse;
//float ordonnee;
//} tPoint;

//float X,Y,Z,A,B,C,F;

struct info{
float *X;
float *Y;
float *Z;
float *A;
float *B;
float *C;
float *F;
};

info info2,bar,juju;

int main()
{
String ^str;
int num_rec = 0;

pin_ptr<const wchar_t> wch = PtrToStringChars(str);

// Convert to a char*
size_t origsize = wcslen(wch) + 1;
const size_t newsize = 1000;
size_t convertedChars = 0;
char nstring[newsize];
wcstombs_s(&convertedChars, nstring, origsize, wch, _TRUNCATE);
//strcat_s(nstring, " (char *)");



errno_t err;
//Read from data file
if((err=fopen_s(&fp1,"C:\\Sample\\project2DAT.dat","rb"))!=0)
fread(&recl,sizeof(recl),1,fp1);
lenr=recl/256/4;
fseek(fp1,28+24,0);
fread(&numf,sizeof(numf),1,fp1);
fseek(fp1,3*28+16,0);
fread(&ymin,4,1,fp1);
fread(&ymax,4,1,fp1);
fread(&numy,4,1,fp1);
fseek(fp1,4*28+16,0);
fread(&xmin,4,1,fp1);
fread(&xmax,4,1,fp1);
fread(&numx,4,1,fp1);

for (long i=0;i<2;i++)
{
fread(&info2.X[i],4,1,fp1);
fread(&info2.Y[i],4,1,fp1);
fread(&info2.Z[i],4,1,fp1);
fread(&info2.A[i],4,1,fp1);
fread(&info2.B[i],4,1,fp1);
fread(&info2.C[i],4,1,fp1);
fread(&info2.F[i],4,1,fp1);
}

fclose(fp1);

//build the step for smaller grid
nb_step_Y=(int)ymax-(int)ymin;
nb_step_X=(int)xmax-(int)xmin;

//build the array with the lines by smaller grid
cc=0;

for (long j=0;j<nb_step_Y;j++)
{
    for (long k=0;k<nb_step_X;k++)
        {
         cc=cc+1;
         dd=0;
      for (i=0;i<numx;i++)
          {
           if ((xmin-0.5*k < info2.X[i] < xmin+0.5*k)&&(ymin-0.5*j < info2.Y[i] < ymin+0.5*j))
           dd=dd+1;
           juju.X[dd]=info2.X[i];
           juju.Y[dd]=info2.Y[i];
           juju.Z[dd]=info2.Z[i];
           juju.A[dd]=info2.A[i];
           juju.B[dd]=info2.B[i];
           juju.C[dd]=info2.C[i];
           juju.F[dd]=info2.F[i];

                 for (cc=1;cc<dd;cc++)
                     {
                      n=k*j;
                      bar.X[n] = juju.X[cc];
                      bar.Y[n] = juju.Y[cc];
                      bar.Z[n] = juju.Z[cc];
                      bar.A[n] = juju.A[cc];
                      bar.B[n] = juju.B[cc];
                      bar.C[n] = juju.C[cc];
                      bar.F[n] = juju.F[cc];
                     };
                  bar.X[n] = bar.X[n]/dd;
                  bar.Y[n] = bar.Y[n]/dd;
                  bar.Z[n] = bar.Z[n]/dd;
          };
        };
};
// recording of result in file txt
FILE* fichier = NULL;
fichier = fopen("C:\\Sample\\donnees.txt", "a+");

if (fichier != NULL)
 {
  fputs("R�capitulatif :\n", fichier);
  for (cc=1;cc<dd;cc++)
   {
    fprintf(fichier, "*", bar.X[cc], "*", bar.Y[cc], "*", bar.Z[cc], "*", bar.A[cc], "*", bar.B[cc], "*", bar.C[cc], "*", bar.F[cc]);
   }
  fclose(fichier); // On ferme le fichier qui a �t� ouvert
  cout << "writing is ok\n";
 }
 else
 {
 // errror 
  printf("Impossible to open the file");
 }
system("pause");
return 0;
};