#include "diabolo.h"
#include "fichier.h"
#include "bolo.h"
#include "arcunit.h"






void	lit_fichier_DIA2(int fen)
{
//int i,f;
F_entete			*E;
F_block_desc 		*BD;
F_block_param 		*BP;
F_donne		*	T;							// tableau courant de donnee
long tt=1397151033;
long total=0;
long val_controleur;
long blockold=0;
double total_controleur=-1;
int bb=0;
int i;
def_gains;

T=malloc(max_nb_donne*sizeof(F_donne));
E=malloc(sizeof(F_entete));
BD=malloc(sizeof(F_block_desc));
BP=malloc(sizeof(F_block_param));
printf("lecture du fichier  \n");
position(fen,0L);
select(printf_ref);
litF(fen,sizeof(F_entete),E);
printf(" version :  =%d     ",E->version);
printf("temps debut, fin =%ld , %ld  \n",tt+E->tempsdebut,tt+E->tempsfin);
E->commentaires[255]=0;
printf("commentaires à la main : %s \n",E->commentaires);

while ( litF(fen,sizeof(F_block_desc),BD) != -1)
	{
//	printf(" BD.change = %d    nb ligne= %d \n",BD->change,BD->n_ligne_don);
	if(BD->change>10)
		{
		char buf[20000];
		printf("block  type  %d,  longueur %d ",BD->change,BD->n_ligne_don);
		litF(fen,BD->n_ligne_don,buf);
		}
	else
	
	if(BD->change==10)
		{
		char comment[10000];
		printf(" commentaire longueur %d ",BD->n_ligne_don);
		litF(fen,BD->n_ligne_don,comment);
//		printf("commentaire : \n");
//		printf(comment);
		printf("\n");
		}
	else
		{
			
		bb++;
		litF(fen,sizeof(F_block_param),BP);
		if(total_controleur==-1)	{val_controleur=BP->temps_controleur_debut;total_controleur=0;}
		total_controleur+=(double)( (BP->temps_controleur_fin - val_controleur) & 0x7fffffff ) ;
		val_controleur=BP->temps_controleur_fin;
		printf(" block %d :chg=%d , %d lignes , horloge %x ,temps debut %ld , fin =%ld ,  ereur%ld ",		
			BD->n_blok,
			BD->change,
			BD->n_ligne_don,
			BP->reglage.horloge,
			BP->temps_controleur_debut,BP->temps_controleur_fin,
			BP->temps_controleur_fin-(blockold+(long)BD->n_ligne_don));
		blockold=BP->temps_controleur_fin;
//		printf("bololu =%d \n",BD->bololut);
		if(BD->n_ligne_don>max_nb_donne) {printf("\n trop de ligne a lire : erreur \n\n");break;}
		litF(fen,sizeof(F_donne)*(long)BD->n_ligne_don,T);
		total+=BD->n_ligne_don;
		if(bb==5) // pour le 5 eme block uniquement
				{
				for(i=0;i<nb_max_bolo;i++)
					{
					printf("\nbolo%d : gain=%d -> %g   V=%d  I=%d  T=%d  L=%d  phase=%d "
						,i+1,gainbrut(BP->reglage.bolo[i]),gain_ampli(BP->reglage.bolo[i])
						,dac_V(BP->reglage.bolo[i]),dac_I(BP->reglage.bolo[i])
						,dac_T(BP->reglage.bolo[i]),dac_L(BP->reglage.bolo[i])
						,phase(BP->reglage.bolo[i])   );
					}
				printf("\n periode horloge=%g ,   nb_mesures=%d  ,  frequence=%g  , temps mort=%d  \n"
					,2./(double)(BP->reglage.horloge.nb_mesures*24),BP->reglage.horloge.nb_mesures/2 ,500000.,BP->reglage.horloge.temp_mort);
				}	  
		}
	
	
	
	printf("\n");
	}

printf(" total fichier %ld lignes  \n",total);
printf("temps  fichier %ld sec     \n",E->tempsfin-E->tempsdebut);
printf("temps  controleur %8.3f sec     \n",1e-6*total_controleur);
printf("frequence de modulation calculée a partir du temps controleur = %g Hz  \n",((double)(total-1))/(2*1e-6*total_controleur));
printf("periode calculee  = %g  ms \n",(1e-3*total_controleur/((double)(total-1))));
printf("frequence moyenne reelle 	 = %g Hz  \n",((double)total)/(double)(2*(E->tempsfin-E->tempsdebut)));
printf("frequence generateur   	 = %g MHz  \n",
		15.0*(((double)(total-1))/(double)(E->tempsfin-E->tempsdebut))/(((double)(total-1))/(1e-6*total_controleur)));
/*
f=ouvreD(0,0,"",0);
for(i=0;i<1000;i++)
	{
	double x,y;
	x=i;
	y=T[i].bolo[0];
//	y=BP->bol_per[0][i];
	trace(f,x,y);
	}
*/	
}