diff options
Diffstat (limited to 'test/gass_read.c')
| -rw-r--r-- | test/gass_read.c | 212 |
1 files changed, 0 insertions, 212 deletions
diff --git a/test/gass_read.c b/test/gass_read.c deleted file mode 100644 index 4e6e1e1..0000000 --- a/test/gass_read.c +++ /dev/null @@ -1,212 +0,0 @@ -/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * - * Copyright by the Board of Trustees of the University of Illinois. * - * All rights reserved. * - * * - * This file is part of HDF5. The full HDF5 copyright notice, including * - * terms governing use, modification, and redistribution, is contained in * - * the files COPYING and Copyright.html. COPYING can be found at the root * - * of the source code distribution tree; Copyright.html can be found at the * - * root level of an installed copy of the electronic HDF5 document set and * - * is linked from the top-level documents page. It can also be found at * - * http://hdf.ncsa.uiuc.edu/HDF5/doc/Copyright.html. If you do not have * - * access to either file, you may request a copy from hdfhelp@ncsa.uiuc.edu. * - * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */ - -/* - * Programmer: Saurabh Bagchi (bagchi@uiuc.edu) - * Wednesday, August 11, 1999. - * - * Modifications: Saurabh Bagchi (Aug 17, 1999) - * Modified to work with VFL (HDF51.3). - */ - -/* Test the following functionality of the GASS driver. - 1. Open a remote file for read (the dataset was written using gass_write.c). - 2. Create a memory buffer to hold the dataset. - 3. Read the dataset into the memory buffer. - 4. Get some information about the dataset from the file. -*/ -#include "h5test.h" -#include <string.h> - -#ifndef H5_HAVE_GASS -int main(void) -{ - printf("Test skipped because GASS driver not available\n"); - return 0; -} -#else - -/*#define URL "ftp://gass:gass12@which/tmp/gass/junk.w"*/ -/*#define URL "http://paz.ncsa.uiuc.edu:8080/test/a.h5"*/ -#define URL "http://hdf.ncsa.uiuc.edu/GLOBUS/a.h5" -/* #define DATASETNAME "Int1Array" */ -#define NX_SUB 3 /* hyperslab dimensions */ -#define NY_SUB 4 -#define NX 7 /* output buffer dimensions */ -#define NY 7 -#define NZ 3 -#define RANK 2 -#define RANK_OUT 3 - -int -main (int argc, char **argv) -{ - hid_t fapl = -1, file, dataset; /* handles */ - char DATASETNAME[32]; - hid_t datatype, dataspace; - hid_t memspace; - H5T_class_t class; /* data type class */ - H5T_order_t order; /* data order */ - size_t size; /* - * size of the data element - * stored in file - */ - hsize_t dimsm[3]; /* memory space dimensions */ - hsize_t dims_out[2]; /* dataset dimensions */ - herr_t status; - - int data_out[NX][NY][NZ ]; /* output buffer */ - - hsize_t count[2]; /* size of the hyperslab in the file */ - hssize_t offset[2]; /* hyperslab offset in the file */ - hsize_t count_out[3]; /* size of the hyperslab in memory */ - hssize_t offset_out[3]; /* hyperslab offset in memory */ - int i, j, k, status_n, rank; - GASS_Info ginf; - - if (argc > 2) { - printf ("Incorrect command line. \n"); - printf ("Correct command line: %s [DATASET NAME] \n", argv[0]); - printf ("The two datasets present in the file are: \"IntArray\" and " - "\"Int1Array\"\n"); - exit(1); - } - - printf ("Correct command line: %s [DATASET NAME] \n", argv[0]); - printf ("The two datasets present in the file are: \"IntArray\" and " - "\"Int1Array\"\n"); - printf ("Default dataset is IntArray. \n"); - - if (argc == 1) - strcpy (DATASETNAME, "IntArray"); - else if (!strcmp(argv[1],"IntArray") || !strcmp(argv[1],"Int1Array")) - strcpy (DATASETNAME, argv[1]); - else - strcpy (DATASETNAME, "IntArray"); - printf ("\n Reading dataset %s \n\n", DATASETNAME); - - for (j = 0; j < NX; j++) { - for (i = 0; i < NY; i++) { - for (k = 0; k < NZ ; k++) - data_out[j][i][k] = 0; - } - } - - /* Create access property list and set the driver to GASS */ - fapl = H5Pcreate (H5P_FILE_ACCESS); - if (fapl < 0) { - printf (" H5Pcreate failed. \n"); - return -1; - } - - - ginf.block_size = 0; - ginf.max_length =0; - - /* ginf = GASS_INFO_NULL; */ - - status = H5Pset_fapl_gass (fapl, ginf); - if (status < 0) { - printf ("H5Pset_fapl_gass failed. \n"); - return -1; - } - - /* - * Open the file and the dataset. - */ - file = H5Fopen(URL, H5F_ACC_RDONLY, fapl); - dataset = H5Dopen(file, DATASETNAME); - - /* - * Get datatype and dataspace handles and then query - * dataset class, order, size, rank and dimensions. - */ - datatype = H5Dget_type(dataset); /* datatype handle */ - class = H5Tget_class(datatype); - if (class == H5T_INTEGER) printf("Data set has INTEGER type \n"); - order = H5Tget_order(datatype); - if (order == H5T_ORDER_LE) printf("Little endian order \n"); - - size = H5Tget_size(datatype); - printf(" Data size is %d \n", size); - - dataspace = H5Dget_space(dataset); /* dataspace handle */ - rank = H5Sget_simple_extent_ndims(dataspace); - status_n = H5Sget_simple_extent_dims(dataspace, dims_out, NULL); - printf("rank %d, dimensions %lu x %lu \n", rank, - (unsigned long)(dims_out[0]), (unsigned long)(dims_out[1])); - - /* - * Define hyperslab in the dataset. - */ - offset[0] = 1; - offset[1] = 2; - count[0] = NX_SUB; - count[1] = NY_SUB; - status = H5Sselect_hyperslab(dataspace, H5S_SELECT_SET, offset, NULL, - count, NULL); - - /* - * Define the memory dataspace. - */ - dimsm[0] = NX; - dimsm[1] = NY; - dimsm[2] = NZ ; - memspace = H5Screate_simple(RANK_OUT,dimsm,NULL); - - /* - * Define memory hyperslab. - */ - offset_out[0] = 3; - offset_out[1] = 0; - offset_out[2] = 0; - count_out[0] = NX_SUB; - count_out[1] = NY_SUB; - count_out[2] = 1; - status = H5Sselect_hyperslab(memspace, H5S_SELECT_SET, offset_out, NULL, - count_out, NULL); - - /* - * Read data from hyperslab in the file into the hyperslab in - * memory and display. - */ - status = H5Dread(dataset, H5T_NATIVE_INT, memspace, dataspace, - H5P_DEFAULT, data_out); - for (j = 0; j < NX; j++) { - for (i = 0; i < NY; i++) printf("%d ", data_out[j][i][0]); - printf("\n"); - } - /* - * 0 0 0 0 0 0 0 - * 0 0 0 0 0 0 0 - * 0 0 0 0 0 0 0 - * 3 4 5 6 0 0 0 - * 4 5 6 7 0 0 0 - * 5 6 7 8 0 0 0 - * 0 0 0 0 0 0 0 - */ - - /* - * Close/release resources. - */ - H5Tclose(datatype); - H5Dclose(dataset); - H5Sclose(dataspace); - H5Sclose(memspace); - H5Fclose(file); - H5Pclose(fapl); - - return 0; -} -#endif |