hdf5/test/srb_read.c
Raymond Lu 86ba18c0a7 [svn-r2800]
Purpose:
    Fix a typo
Description:
    Just a typo in print output.
Solution:
    [details about the changes, algorithm, etc...]
    [Please as detail as you can since your own explanation is
    better than others guessing it from the code.]
Platforms tested:
    [machines you have tested the changed version.  This is absolute
    important.  Test it out on at least two or three different platforms
    such as Big-endian-32bit (SUN/IRIX), little-endian-32(LINUX) and
    64-bit (IRIX64/UNICOS/DEC-ALPHA) would be good.]
2000-11-03 15:41:13 -05:00

215 lines
7.9 KiB
C

/*
* Copyright (C) 2000 NCSA
* All rights reserved.
*
* Programmer: Raymond Lu <slu@ncsa.uiuc.edu>
* April 14, 2000
*
* Purpose: Test HDF-SRB if it is configured and compiled. Read a chunk of
* dataset of an HDF5 file with an integer array on SRB server,
* using hyperslab.
* Usage: The HDF5 is built on the top of SRB. So you just need to
* activate a connection to SRB server by calling function
* H5Pset_fapl_srb()(after creating file property list by calling
* H5Pcreate()). All server information is passed in through its
* parameter SRB_Info(a structure). Its fields are
* char *srbHost: SRB host address of server. If the input is
* NULL, it will use the environment varible "srbHost" if
* it is defined. If it is not defined, the
* ~/.srb/.MdasEnd file will be checked next. If not, it
* will use the hostname of the client machine.
* char *srbPort: SRB host port number of server. If the
* input value is NULL, it will use the env variable
* "srbPort" if it is defined. If not, it defaults to
* 5558.
* char *srbAuth: SRB Authentication-password. It is used to
* define password for MDAS or SEA authentication.
* For SEA authentication, this is the password used by
* the SEA library to decrypt the encrypted private key
* stored in the file ~/.SEAuuuuu@ddddd(where uuuuu is
* the user ID and ddddd is the user domain name). This
* input is not needed if an unencrypted private key is
* available in the /tmp directory(generated using the
* 'seaauth auto' command). To provide additional
* flexibility, a client may also use the environment
* variable "srbAuth" to specify the password. A client
* may also supply the password in the ~/.srb/.MdasAuth
* file. If a client uses more than one method to specfy
* the password, the value given in this function call
* will take precedent, then the environment variable
* "srbAuth", and lastly, the ~/.srb/.MdasAuth file.
* int storSysType: Storage system type on SRB server.
* 0=Unix, 1=UniTree, 2=HPSS, 3=FTP, 4=HTTP
* int mode: File access mode, same definition with Unix.
* int size: File Size, only valid for HPSS, -1 is default
*
* Then you can open an HDF5 file by calling H5Fopen(). When
* you pass in the file name, it has to be file name with an
* absolute path. It you use SDSC server(ghidorah.sdsc.edu), your
* home directory is possibly in /projects/mdas/srb/SRBVault/.
*/
#include <h5test.h>
#ifndef H5_HAVE_SRB
int main(void)
{
printf("Test skipped because SRB driver not available\n");
return 0;
}
#else
#define fileName "/projects/mdas/srb/SRBVault/slu.ncsa/a.h5"
#define DATASETNAME "IntArray"
#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(void)
{
hid_t fapl=-1, fid = -1, dataset;
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;
SRB_Info srb_info={ NULL, /* Use host name in ~/.srb/.MdasEnv */
NULL, /* Use environment variable "srbPort" */
NULL, /* Use password in ~/.srb/.MdasEnv */
0, /* Use Unix storage system. */
0600, /* File was created for read and write for
* owner */
-1 /* default */
};
for (j = 0; j < NX; j++) {
for (i = 0; i < NY; i++) {
for (k = 0; k < NZ ; k++)
data_out[j][i][k] = 0;
}
}
fapl = H5Pcreate(H5P_FILE_ACCESS);
if (fapl < 0) {
printf (" H5Pcreate failed. \n");
return -1;
}
status = H5Pset_fapl_srb(fapl, srb_info);
if (status < 0) {
printf ("H5Pset_fapl_srb failed. \n");
return -1;
}
fid = H5Fopen(fileName, H5F_ACC_RDONLY, fapl);
if (fid < 0) {
printf ("H5Fopen failed. \n");
return -1;
}
dataset = H5Dopen(fid, DATASETNAME);
if(dataset<0) {
printf ("H5Dopen failed. \n");
return -1;
}
/*
* 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(fid);
H5Pclose(fapl);
printf("Test finished!\n");
return 0;
}
#endif