/*!\file:  DomainOutlineRead.cpp
 * \brief DomainOutlineRead.c: read the vertex coordinates defined in a domain 
 * outline from Argus (.exp file). The first contour in the file is for 
 * the outside domain outline. The following contours represent holes in
 * the domain.
 */

#include <stdio.h>
#include "../Alloc/alloc.h"
#include "../../include/include.h"
#include "../../objects/objects.h"
#include "../Exceptions/exceptions.h"
#include "../../Container/DataSet.h"

int DomainOutlineRead(int* pnprof,int** pprofnvertices,double*** ppprofx,double*** ppprofy,bool** pclosed,char* domainname,bool whole=true){

	
	/*indexing: */
	int i,counter;

	/*I/O: */
	FILE* fid=NULL;
	char chardummy[256];
	double  ddummy;

	/*output: */
	int nprof; //number of profiles in the domainname file
	int* profnvertices=NULL; //array holding the number of vertices for the nprof profiles
	double** pprofx=NULL; //array of profiles x coordinates
	double** pprofy=NULL; //array of profiles y coordinates
	bool* closed=NULL; //array holding closed flags for the nprof profiles

	/*For each profile: */
	int n;
	double* x=NULL;
	double* y=NULL;
	bool cl;

	/*open domain outline file for reading: */
	if ((fid=fopen(domainname,"r"))==NULL){
		_error_("%s%s","could not find domain file ",domainname); 
	}

	/*Do a first pass through the domainname file, to figure out how many profiles 
	 *we need to read: */
	nprof=1;
	for(;;){
		fscanf(fid,"%256s %256s\n",chardummy,chardummy);
		fscanf(fid,"%256s %256s\n",chardummy,chardummy);
		fscanf(fid,"%256s %256s %256s %256s\n",chardummy,chardummy,chardummy,chardummy);
		fscanf(fid,"%20u %256s\n",&n,chardummy);
		fscanf(fid,"%256s %256s %256s %256s %256s\n",chardummy,chardummy,chardummy,chardummy,chardummy);
		for (i=0;i<n;i++){
			fscanf(fid,"%20lf %20lf\n",&ddummy,&ddummy);
		}
		/*Ok, we have faked one profile reading, check whether we are at the end of the file, otherwise, keep fake reading next profile:*/
		if (feof(fid)){
			break;
		}
		nprof++;
	}
	
	/*Allocate and initialize all the profiles: */
	profnvertices=(int*)xmalloc(nprof*sizeof(int));
	pprofx=(double**)xmalloc(nprof*sizeof(double*));
	pprofy=(double**)xmalloc(nprof*sizeof(double*));
	for (i=0;i<nprof;i++){
		pprofx[i]=NULL;
		pprofy[i]=NULL;
	}
	closed=(bool*)xmalloc(nprof*sizeof(bool));

	/*Reaset file pointer to beginning of file: */
	fseek(fid,0,SEEK_SET);

	/*Start reading profiles: */
	for(counter=0;counter<nprof;counter++){

		/*Skip header: */
		fscanf(fid,"%256s %256s\n",chardummy,chardummy);
		fscanf(fid,"%256s %256s\n",chardummy,chardummy);
		fscanf(fid,"%256s %256s %256s %256s\n",chardummy,chardummy,chardummy,chardummy);
		
		/*Get number of profile vertices: */
		fscanf(fid,"%20u %256s\n",&n,chardummy);
	
		/*Skip next line: */
		fscanf(fid,"%256s %256s %256s %256s %256s\n",chardummy,chardummy,chardummy,chardummy,chardummy);

		/*Allocate vertices: */
		printf("number of n: %i\n",n);
		x=(double*)xmalloc(n*sizeof(double));
		y=(double*)xmalloc(n*sizeof(double));
		

		/*Read vertices: */
		for (i=0;i<n;i++){
			fscanf(fid,"%20lf %20lf\n",&x[i],&y[i]);
		}

		/*Now check that we are dealing with open contours: */
		cl=false;
		if((x[0]==x[n-1]) && (y[0]==y[n-1])){
			cl=true;
			if (!whole) {
				n=n-1;
			}
		}

		/*Assign pointers: */
		profnvertices[counter]=n;
		pprofx[counter]=x;
		pprofy[counter]=y;
		closed[counter]=cl;
	}

	/*close domain outline file: */
	fclose(fid);

	/*Assign output pointers: */
	*pnprof=nprof;
	*pprofnvertices=profnvertices;
	*ppprofx=pprofx;
	*ppprofy=pprofy;
	if(pclosed)*pclosed=closed;
	else       xfree((void**)&closed);
}

DataSet* DomainOutlineRead(char* domainname,bool whole=true){

	/*indexing: */
	int i;

	/*intermediary: */
	int nprof;
	int* profnvertices=NULL;
	double** pprofx=NULL;
	double** pprofy=NULL;

	Contour* contour=NULL;

	/*output: */
	DataSet* domain=NULL;

	/*get domain outline from intermediary function:*/
	DomainOutlineRead(&nprof,&profnvertices,&pprofx, &pprofy, NULL,domainname,whole);

	double* a=pprofx[0]; for(i=0;i<4;i++)printf("%g\n",a[i]);
	a=pprofy[0]; for(i=0;i<4;i++)printf("%g\n",a[i]);
	printf("%i %i\n",profnvertices[0],nprof);

	/*now create dataset of contours: */
	domain=new DataSet(0);
	printf("ok1\n");

	for(i=0;i<nprof;i++){
		domain->AddObject(new Contour(i,profnvertices[i],pprofx[i],pprofy[i],1));
	}
	printf("ok2\n");

	return domain;
}