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

typedef struct results_by_phase{
  double k_arrival[2];
  double k_delta;
  double k_estimate;
  double u_arrival[2];
  double u_delta;
  double u_estimate;
} results_by_phase_t;

results_by_phase_t phase_array[MAX_n_replies]; /* Should be max_phase */

/*........................................................................*/
void
print_estimates(void){
  /* Prints the estimate fields of the phase_rrray */
  int i;
  printf("The End-to-End bottleneck bandwidth estimates\n");
  printf("according to the packet pair algorithm are:\n");
  for (i=0;i<max_phase;i++){
      printf("%d. %f Kb/s\t%f Kb/s\n",i,phase_array[i].k_estimate,
	     phase_array[i].u_estimate);
  }
}

/*........................................................................*/
void
fill_phase_array(int marker_0,int marker_1){
  /* Reads results_array (one record represents a probe packet response) *
   * and fills phase_array (one record represents a full probe response) *
   * marker_0 and marker_1 correspond to the packet_nos to fill the array*/
  int i,phase_no,packet_no,size;
  double k_arrival,u_arrival;

  /* Initialize phase_array */
  for (i=0;i<max_phase;i++){
    phase_array[i].k_arrival[0]=-1.0;
    phase_array[i].k_arrival[1]=-1.0;
    phase_array[i].k_delta=-1.0;
    phase_array[i].k_estimate=-1.0;
    phase_array[i].u_arrival[0]=-1.0;
    phase_array[i].u_arrival[1]=-1.0;
    phase_array[i].u_delta=-1.0;
    phase_array[i].u_estimate=-1.0;
  }

  /* Now fill phase_array */
  for (i=0;i<n_replies;i++){
    phase_no=results_array[i].phase_no;
    packet_no= results_array[i].packet_no;
    k_arrival=results_array[i].k_arrival;
    u_arrival=results_array[i].u_arrival;
    size=results_array[i].size;
    if (packet_no==marker_0) {
      phase_array[phase_no].k_arrival[0]=k_arrival;
      phase_array[phase_no].u_arrival[0]=u_arrival;
    } else if (packet_no==marker_1) {
      phase_array[phase_no].k_arrival[1]=k_arrival;
      phase_array[phase_no].u_arrival[1]=u_arrival;
    }
      
  }

  /* Now some processing of phase_array filling */
  for (i=0;i<max_phase;i++){
    if ((phase_array[i].k_arrival[0]>0)
	&&(phase_array[i].k_arrival[1]>0)){
      phase_array[i].k_delta= 
	phase_array[i].k_arrival[1]-
	phase_array[i].k_arrival[0];
      phase_array[i].u_delta= 
	phase_array[i].u_arrival[1]-
	phase_array[i].u_arrival[0];       
    } 
  }

}

/*........................................................................*
 *                                                                        *
 *                Bandwidth Inference Functions                           * 
 *                                                                        *
 *........................................................................*/

/*........................................................................*/
double 
BW_e2e(void){
  /* Use phase_array to estimate BW using the packet-pair technique       *
   * fill_phase_array should be called before calling this function       */
  int i;
  for (i=0;i<max_phase;i++){
    if (phase_array[i].k_delta>0){
      phase_array[i].k_estimate=markers_size[0]*8/phase_array[i].k_delta;
      phase_array[i].u_estimate=markers_size[0]*8/phase_array[i].u_delta;
    }
  }
}