void get_3d(double ***  );
void get_1d(double *  );

int test(){

  //declaring a 1d array

  double *nu;
  nu = malloc(6*sizeof(double));

  for(int j=0;j<6;j++){
    nu[j]=j;
  }
  

  //declaring a 3d array

 double ***nu_iso;

  nu_iso = malloc(6*sizeof(double**));
  for(int i=0;i<6;i++)
    {
      nu_iso[i] = malloc(13*sizeof(double *));

      for(int j=0;j<13;j++)
	{
	  nu_iso[i][j] = malloc(8*sizeof(double));
	}
    }


  for(int j=0;j<6;j++){
    for(int k=0;k<13;k++){
      for(int l=0;l<8;l++){

	nu_iso[j][k][l] = j+k+l;
	//printf("%d %d %d %lf\n",j,k,l,nu_iso[j][k][l]);

      }
    }
  }


  get_1d(nu);
  get_3d(nu_iso);


  return 0;
}



void get_1d(double *bin_corr ){

  cout<<"inside get_1d"<<endl;

  for(int j=0;j<6;j++){

    printf("%d %lf\n",j,bin_corr[j]);

  }

void get_3d(double ***bin_corr ){

  cout<<"inside get_3d"<<endl;

  for(int j=0;j<6;j++){
    for(int k=0;k<13;k++){
      for(int l=0;l<8;l++){

	printf("%d %d %d %lf\n",j,k,l,bin_corr[j][k][l]);

      }
    }
  }

}




}