/* This is part of the netCDF package. Copyright 2005 University
   Corporation for Atmospheric Research/Unidata See COPYRIGHT file for
   conditions of use. See www.unidata.ucar.edu for more info.

   Test netcdf-4 variable length code. 

   $Id: tst_vl.c,v 1.24 2008/06/10 19:38:03 russ Exp $
*/

#include <nc_tests.h>
#include "netcdf.h"

#define FILE_NAME "tst_vl.nc"
#define DIM_LEN 3
#define VLEN_NAME "vlen_name"
#define ATT_NAME "att_name"
#define DIM_NAME "dim"
#define VAR_NAME "var"
#define PHONY_VAL -99

/* Use this struct to hold vlen data: */
/* typedef struct  { */
/* size_t len;  Length of VL data (in base type units)        */
/* void *p;     Pointer to VL data          */
/* } hvl_t; */
int
main(int argc, char **argv)
{

#ifdef USE_PARALLEL
   MPI_Init(&argc, &argv);
#endif

   printf("\n*** Testing netcdf-4 variable length type functions.\n");
   printf("*** testing vlen (of ints) attribute...");
   {
      int ncid, typeid;
      nc_vlen_t data[DIM_LEN];
      size_t size_in;
      nc_type base_nc_type_in;
      char name_in[NC_MAX_NAME + 1];
      int *phony, class_in;
      int i, j;
      
      /* Create phony data. */
      for (i=0; i<DIM_LEN; i++)
      {
	 if (!(phony = malloc(sizeof(int) * (i+1))))
	    return NC_ENOMEM;
	 for (j = 0; j < i + 1; j++)
	    phony[j] = PHONY_VAL;
	 data[i].p = phony;
	 data[i].len = i+1;
      }
      
      /*nc_set_log_level(5);*/

      if (nc_create(FILE_NAME, NC_NETCDF4, &ncid)) ERR;
      if (nc_def_vlen(ncid, VLEN_NAME, NC_INT, &typeid)) ERR;
      if (nc_inq_vlen(ncid, typeid, name_in, &size_in, &base_nc_type_in)) ERR;
      if (base_nc_type_in != NC_INT || (size_in != sizeof(int) || strcmp(name_in, VLEN_NAME))) ERR;
      if (nc_inq_user_type(ncid, typeid, name_in, &size_in, &base_nc_type_in, NULL, &class_in)) ERR;
      if (base_nc_type_in != NC_INT || (size_in != sizeof(int) || strcmp(name_in, VLEN_NAME))) ERR;
      if (nc_inq_compound(ncid, typeid, name_in, &size_in, NULL) != NC_EBADTYPE) ERR;
      if (nc_put_att(ncid, NC_GLOBAL, ATT_NAME, typeid, DIM_LEN, data)) ERR;
      if (nc_close(ncid)) ERR;
      
      /* Reset data to make sure it was really copied. */
      for (i=0; i<DIM_LEN; i++)
	 for (j = 0; j < i + 1; j++)
	    ((int *)(data[i].p))[j] = 0;

      /* Check it out. */
      if (nc_open(FILE_NAME, NC_NOWRITE, &ncid)) ERR;
/*      if (nc_get_att(ncid, NC_GLOBAL, ATT_NAME, data_in)) ERR;*/
/*      for (i = 0; i < DIM_LEN; i++)
      {
	 if (data_in[i].len != i + 1) ERR;
	 for (j = 0; j < data_in[i].len; j++)
	    if (*((int *)data_in[i].p) != PHONY_VAL)
	    {
	       printf("*((int *)data_in[%d].p = %d (0x%x)\n", i, *((int *)data_in[i].p),
		      *((int *)data_in[i].p));
	       ERR;
	       return 2;
	    }
      }
      if (nc_inq_att(ncid, NC_GLOBAL, ATT_NAME, &typeid, &len_in)) ERR;
      if (len_in != DIM_LEN) ERR;
      if (nc_inq_vlen(ncid, typeid, name_in, &size_in, &base_nc_type_in)) ERR;
      if (base_nc_type_in != NC_INT || (size_in != 4 && strcmp(name_in, VLEN_NAME)))
      {
	 printf("base_nc_type_in = %d size_in = %d name_in = %s\n", base_nc_type_in,
		(int)size_in, name_in);
	 ERR;
      }
      if (nc_inq_user_type(ncid, typeid, name_in, &size_in, &base_nc_type_in, NULL, &class_in)) ERR;
      if (base_nc_type_in != NC_INT || (size_in != 4 || strcmp(name_in, VLEN_NAME))) ERR;*/
      if (nc_close(ncid)) ERR;

      /* Free the memory used in our phony data. */
      for (i=0; i<DIM_LEN; i++)
      {
/*	 if (nc_free_vlen(&data[i]) || nc_free_vlen(&data_in[i])) ERR;*/
	 free(data[i].p);
/*	 free(data_in[i].p);*/
      }
   }

   SUMMARIZE_ERR;

/*    printf("*** testing vlen (of ints) variable..."); */
/*    { */
/*       int dimid, varid, dimids[] = {0}; */
/*       char var_name[NC_MAX_NAME+1]; */
/*       int nvars, natts, ndims, unlimdimid, dimids_var[1]; */
/*       nc_type var_type; */
/*       int ncid, typeid; */
/*       nc_vlen_t data[DIM_LEN], data_in[DIM_LEN]; */
/*       int *phony; */
/*       int i, j; */
      
/*       /\* Create phony data. *\/ */
/*       for (i=0; i<DIM_LEN; i++) */
/*       { */
/* 	 if (!(phony = malloc(sizeof(int) * (i+1)))) */
/* 	    return NC_ENOMEM; */
/* 	 for (j=0; j<i+1; j++) */
/* 	    phony[j] = PHONY_VAL; */
/* 	 data[i].p = phony; */
/* 	 data[i].len = i+1; */
/*       } */

/*       if (nc_create(FILE_NAME, NC_NETCDF4, &ncid)) ERR; */
/*       if (nc_def_vlen(ncid, "name1", NC_INT, &typeid)) ERR; */
/*       if (nc_def_dim(ncid, DIM_NAME, DIM_LEN, &dimid)) ERR; */
/*       if (nc_def_var(ncid, VAR_NAME, typeid, 1, dimids, &varid)) ERR; */
/*       if (nc_put_var(ncid, varid, data)) ERR; */
/*       if (nc_close(ncid)) ERR; */
      
/*       /\* Check it out. *\/ */
/*       if (nc_open(FILE_NAME, NC_NOWRITE, &ncid)) ERR; */
/*       if (nc_inq(ncid, &ndims, &nvars, &natts, &unlimdimid)) ERR; */
/*       if (ndims != 1 || nvars != 1 || natts != 0 || unlimdimid != -1) ERR; */
/*       if (nc_inq_var(ncid, 0, var_name, &var_type, &ndims, dimids_var, &natts)) ERR; */
/*       if (ndims != 1 || strcmp(var_name, VAR_NAME) ||  */
/* 	  dimids_var[0] != dimids[0] || natts != 0) ERR; */
/*       if (nc_get_var(ncid, 0, data_in)) ERR; */
/*       for (i=0; i<DIM_LEN; i++) */
/* 	 for (j=0; j<data_in[i].len; j++) */
/* 	    if (*((int *)data_in[i].p) != *((int *)data[i].p)) */
/* 	    { */
/* 	       printf("*((int *)data_in[%d].p = %d (0x%x)\n", i, *((int *)data_in[i].p), */
/* 		      *((int *)data_in[i].p)); */
/* 	       ERR; */
/* 	       return 2; */
/* 	    } */
/*       if (nc_close(ncid)) ERR; */

/*       /\* Free the memory used in our phony data. *\/ */
/*       for (i=0; i<DIM_LEN; i++) */
/* 	 if (nc_free_vlen(&data[i]) || nc_free_vlen(&data_in[i])) ERR; */
/*    } */

/*    SUMMARIZE_ERR; */

   FINAL_RESULTS;
#ifdef USE_PARALLEL
   MPI_Finalize();
#endif   
}