mirror of
https://github.com/HDFGroup/hdf5.git
synced 2024-11-27 02:10:55 +08:00
ce18afe663
Purpose:
added more tests to the H5set_extent function
Description:
there was a bug in the H5S_select_fill call when the fill value was not defined
added 2 more set of tests that call H5set_extent without the fill value being defined
there are now 4 sets of tests:
Test H5Dset_extent with chunks on the raw data cache, with and without fill value
Test H5Dset_extent with chunks written to file (b-tree routines), with and without fill value
Platforms tested:
w2000
linux
548 lines
14 KiB
C
548 lines
14 KiB
C
/*
|
|
* Copyright (C) 2002 NCSA
|
|
* All rights reserved.
|
|
*
|
|
* Programmer: Pedro Vicente <pvn@ncsa.uiuc.edu>
|
|
* April 12, 2002
|
|
*
|
|
* Purpose: Tests the H5Dset_extent call
|
|
*/
|
|
|
|
|
|
#include "hdf5.h"
|
|
#include "h5test.h"
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
*
|
|
* Tests the function H5Dset_extent. In the current version of the library
|
|
* the dataset MUST be chunked.
|
|
*
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
#define RANK 2
|
|
|
|
|
|
int main( void )
|
|
{
|
|
|
|
hid_t file_id;
|
|
hid_t dataset_id;
|
|
hid_t space_id;
|
|
hid_t plist_id;
|
|
hsize_t dims[RANK] = { 90, 90 };
|
|
hsize_t dims_new[RANK] = { 70, 70 };
|
|
hsize_t dims_chunk[RANK] = { 20, 20 };
|
|
hsize_t dims_out[RANK];
|
|
hsize_t maxdims[RANK] = { H5S_UNLIMITED, H5S_UNLIMITED };
|
|
int data[ 90 ][ 90 ];
|
|
int buf1[ 70 ][ 70 ];
|
|
int buf2[ 90 ][ 90 ];
|
|
int i, j, n = 0;
|
|
int fillvalue = 1; /* Fill value for the dataset */
|
|
|
|
|
|
for( i = 0; i < 90; i++ )
|
|
{
|
|
for( j = 0; j < 90; j++ )
|
|
{
|
|
data[i][j] = n++;
|
|
}
|
|
}
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Test H5Dset_extent with chunks on the raw data cache
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
|
|
/* Create a new file using default properties. */
|
|
if ((file_id = H5Fcreate( "set_extent_create.h5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT ))<0) goto out;
|
|
|
|
TESTING("extend dataset create with fill value");
|
|
|
|
/* Create the data space with unlimited dimensions. */
|
|
if ((space_id = H5Screate_simple( RANK, dims, maxdims ))<0) goto out;
|
|
|
|
/* Modify dataset creation properties, i.e. enable chunking. */
|
|
if ((plist_id = H5Pcreate (H5P_DATASET_CREATE ))<0) goto out;
|
|
if (H5Pset_chunk( plist_id, RANK, dims_chunk )<0) goto out;
|
|
if (H5Pset_fill_value( plist_id, H5T_NATIVE_INT, &fillvalue )<0) goto out;
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Create and write one dataset
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
/* Create a new dataset */
|
|
if ((dataset_id = H5Dcreate( file_id , "Dataset1", H5T_NATIVE_INT, space_id, plist_id ))<0) goto out;
|
|
|
|
/* Write the data. */
|
|
if (H5Dwrite( dataset_id , H5T_NATIVE_INT, space_id, H5S_ALL, H5P_DEFAULT, data )<0) goto out;
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Set new dimensions for the array; shrink it
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
/* Set new dimensions for the array. */
|
|
if (H5Dset_extent( dataset_id , dims_new )<0) goto out;
|
|
|
|
/* Get the space. */
|
|
if ((space_id = H5Dget_space( dataset_id ))<0) goto out;
|
|
|
|
/* Get dimensions. */
|
|
if (H5Sget_simple_extent_dims( space_id, dims_out, NULL )<0) goto out;
|
|
|
|
if ( dims_out[0] != dims_new[0] )
|
|
goto out;
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Read
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
/* Read the new dataset. */
|
|
if (H5Dread( dataset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf1 )<0) goto out;
|
|
|
|
|
|
/* Compare the read array with the original array */
|
|
for( i = 0; i < (int)dims_out[0]; i++ )
|
|
{
|
|
for( j = 0; j < (int)dims_out[1]; j++ )
|
|
{
|
|
if ( buf1[i][j] != data[i][j] ) {
|
|
goto out;
|
|
}
|
|
}
|
|
}
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Set new dimensions for the array; expand it again
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
/* Set new dimensions for the array. */
|
|
if (H5Dset_extent( dataset_id , dims )<0) goto out;
|
|
|
|
/* Get the space. */
|
|
if ((space_id = H5Dget_space( dataset_id ))<0) goto out;
|
|
|
|
/* Get dimensions. */
|
|
if (H5Sget_simple_extent_dims( space_id, dims_out, NULL )<0) goto out;
|
|
|
|
if ( dims_out[0] != dims[0] )
|
|
goto out;
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Read
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
/* Read the new dataset. */
|
|
if (H5Dread( dataset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf2 )<0) goto out;
|
|
|
|
/* Compare the read array with the original array */
|
|
for( i = 0; i < (int)dims_out[0]; i++ )
|
|
{
|
|
for( j = 0; j < (int)dims_out[1]; j++ )
|
|
{
|
|
if ( i >= 70 || j >= 70 )
|
|
{
|
|
if ( buf2[i][j] != fillvalue ) {
|
|
goto out;
|
|
}
|
|
}
|
|
else
|
|
{
|
|
if ( buf2[i][j] != data[i][j] ) {
|
|
goto out;
|
|
}
|
|
}
|
|
|
|
}
|
|
}
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Close/release resources
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
H5Dclose( dataset_id );
|
|
H5Sclose( space_id );
|
|
H5Pclose( plist_id );
|
|
|
|
PASSED();
|
|
TESTING("extend dataset create without fill value");
|
|
|
|
/* Create the data space with unlimited dimensions. */
|
|
if ((space_id = H5Screate_simple( RANK, dims, maxdims ))<0) goto out;
|
|
|
|
/* Modify dataset creation properties, i.e. enable chunking. */
|
|
if ((plist_id = H5Pcreate (H5P_DATASET_CREATE ))<0) goto out;
|
|
if (H5Pset_chunk( plist_id, RANK, dims_chunk )<0) goto out;
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Create and write one dataset
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
/* Create a new dataset */
|
|
if ((dataset_id = H5Dcreate( file_id , "Dataset2", H5T_NATIVE_INT, space_id, plist_id ))<0) goto out;
|
|
|
|
/* Write the data. */
|
|
if (H5Dwrite( dataset_id , H5T_NATIVE_INT, space_id, H5S_ALL, H5P_DEFAULT, data )<0) goto out;
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Set new dimensions for the array; shrink it
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
/* Set new dimensions for the array. */
|
|
if (H5Dset_extent( dataset_id , dims_new )<0) goto out;
|
|
|
|
/* Get the space. */
|
|
if ((space_id = H5Dget_space( dataset_id ))<0) goto out;
|
|
|
|
/* Get dimensions. */
|
|
if (H5Sget_simple_extent_dims( space_id, dims_out, NULL )<0) goto out;
|
|
|
|
if ( dims_out[0] != dims_new[0] )
|
|
goto out;
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Read
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
/* Read the new dataset. */
|
|
if (H5Dread( dataset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf1 )<0) goto out;
|
|
|
|
|
|
/* Compare the read array with the original array */
|
|
for( i = 0; i < (int)dims_out[0]; i++ )
|
|
{
|
|
for( j = 0; j < (int)dims_out[1]; j++ )
|
|
{
|
|
if ( buf1[i][j] != data[i][j] ) {
|
|
goto out;
|
|
}
|
|
}
|
|
}
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Set new dimensions for the array; expand it again
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
/* Set new dimensions for the array. */
|
|
if (H5Dset_extent( dataset_id , dims )<0) goto out;
|
|
|
|
/* Get the space. */
|
|
if ((space_id = H5Dget_space( dataset_id ))<0) goto out;
|
|
|
|
/* Get dimensions. */
|
|
if (H5Sget_simple_extent_dims( space_id, dims_out, NULL )<0) goto out;
|
|
|
|
if ( dims_out[0] != dims[0] )
|
|
goto out;
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Read
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
/* Read the new dataset. */
|
|
if (H5Dread( dataset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf2 )<0) goto out;
|
|
|
|
/* Compare the read array with the original array */
|
|
for( i = 0; i < (int)dims_out[0]; i++ )
|
|
{
|
|
for( j = 0; j < (int)dims_out[1]; j++ )
|
|
{
|
|
if ( i >= 70 || j >= 70 )
|
|
{
|
|
if ( buf2[i][j] != 0 ) {
|
|
goto out;
|
|
}
|
|
}
|
|
else
|
|
{
|
|
if ( buf2[i][j] != data[i][j] ) {
|
|
goto out;
|
|
}
|
|
}
|
|
|
|
}
|
|
}
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Close/release resources
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
H5Dclose( dataset_id );
|
|
H5Sclose( space_id );
|
|
H5Pclose( plist_id );
|
|
|
|
H5Fclose( file_id );
|
|
|
|
PASSED();
|
|
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Test H5Dset_extent with chunks written to file
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
|
|
/* Create a new file using default properties. */
|
|
if ((file_id = H5Fcreate( "set_extent_read.h5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT ))<0) goto out;
|
|
|
|
|
|
TESTING("extend dataset read with fill value");
|
|
|
|
/* Create the data space with unlimited dimensions. */
|
|
if ((space_id = H5Screate_simple( RANK, dims, maxdims ))<0) goto out;
|
|
|
|
/* Modify dataset creation properties, i.e. enable chunking. */
|
|
if ((plist_id = H5Pcreate (H5P_DATASET_CREATE ))<0) goto out;
|
|
if (H5Pset_chunk( plist_id, RANK, dims_chunk )<0) goto out;
|
|
if (H5Pset_fill_value( plist_id, H5T_NATIVE_INT, &fillvalue )<0) goto out;
|
|
|
|
/* Create a new dataset within the file using cparms creation properties. */
|
|
if ((dataset_id = H5Dcreate( file_id , "Dataset1", H5T_NATIVE_INT, space_id, plist_id ))<0) goto out;
|
|
|
|
/* Write the data. */
|
|
if (H5Dwrite( dataset_id , H5T_NATIVE_INT, space_id, H5S_ALL, H5P_DEFAULT, data )<0) goto out;
|
|
|
|
/* Close/release resources. */
|
|
H5Dclose( dataset_id );
|
|
H5Sclose( space_id );
|
|
H5Pclose( plist_id );
|
|
H5Fclose( file_id );
|
|
|
|
|
|
/* Open the file */
|
|
if ((file_id = H5Fopen( "set_extent_read.h5", H5F_ACC_RDWR, H5P_DEFAULT ))<0) goto out;
|
|
|
|
/* Open the dataset */
|
|
if ((dataset_id = H5Dopen( file_id , "Dataset1" ))<0) goto out;
|
|
|
|
/* Set new dimensions for the array. */
|
|
if (H5Dset_extent( dataset_id, dims_new )<0) goto out;
|
|
|
|
/* Get the space. */
|
|
if ((space_id = H5Dget_space( dataset_id ))<0) goto out;
|
|
|
|
/* Get dimensions. */
|
|
if (H5Sget_simple_extent_dims( space_id, dims_out, NULL )<0) goto out;
|
|
|
|
if ( dims_out[0] != dims_new[0] )
|
|
goto out;
|
|
|
|
/* Read the new dataset. */
|
|
if (H5Dread( dataset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf1 )<0) goto out;
|
|
|
|
/* Compare the read array with the original array */
|
|
for( i = 0; i < (int)dims_out[0]; i++ )
|
|
{
|
|
for( j = 0; j < (int)dims_out[1]; j++ )
|
|
{
|
|
if ( buf1[i][j] != data[i][j] ) {
|
|
goto out;
|
|
}
|
|
}
|
|
}
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Set new dimensions for the array; expand it again
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
/* Set new dimensions for the array. */
|
|
if (H5Dset_extent( dataset_id , dims )<0) goto out;
|
|
|
|
/* Get the space. */
|
|
if ((space_id = H5Dget_space( dataset_id ))<0) goto out;
|
|
|
|
/* Get dimensions. */
|
|
if (H5Sget_simple_extent_dims( space_id, dims_out, NULL )<0) goto out;
|
|
|
|
if ( dims_out[0] != dims[0] )
|
|
goto out;
|
|
|
|
/* Read the new dataset. */
|
|
if (H5Dread( dataset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf2 )<0) goto out;
|
|
|
|
/* Compare the read array with the original array */
|
|
for( i = 0; i < (int)dims_out[0]; i++ )
|
|
{
|
|
for( j = 0; j < (int)dims_out[1]; j++ )
|
|
{
|
|
if ( i >= 70 || j >= 70 )
|
|
{
|
|
if ( buf2[i][j] != fillvalue ) {
|
|
goto out;
|
|
}
|
|
}
|
|
else
|
|
{
|
|
if ( buf2[i][j] != data[i][j] ) {
|
|
goto out;
|
|
}
|
|
}
|
|
|
|
}
|
|
}
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Close/release resources
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
H5Dclose( dataset_id );
|
|
H5Sclose( space_id );
|
|
|
|
PASSED();
|
|
|
|
|
|
TESTING("extend dataset read without fill value");
|
|
|
|
/* Create the data space with unlimited dimensions. */
|
|
if ((space_id = H5Screate_simple( RANK, dims, maxdims ))<0) goto out;
|
|
|
|
/* Modify dataset creation properties, i.e. enable chunking. */
|
|
if ((plist_id = H5Pcreate (H5P_DATASET_CREATE ))<0) goto out;
|
|
if (H5Pset_chunk( plist_id, RANK, dims_chunk )<0) goto out;
|
|
|
|
/* Create a new dataset within the file using cparms creation properties. */
|
|
if ((dataset_id = H5Dcreate( file_id , "Dataset2", H5T_NATIVE_INT, space_id, plist_id ))<0) goto out;
|
|
|
|
/* Write the data. */
|
|
if (H5Dwrite( dataset_id , H5T_NATIVE_INT, space_id, H5S_ALL, H5P_DEFAULT, data )<0) goto out;
|
|
|
|
/* Close/release resources. */
|
|
H5Dclose( dataset_id );
|
|
H5Sclose( space_id );
|
|
H5Pclose( plist_id );
|
|
H5Fclose( file_id );
|
|
|
|
|
|
/* Open the file */
|
|
if ((file_id = H5Fopen( "set_extent_read.h5", H5F_ACC_RDWR, H5P_DEFAULT ))<0) goto out;
|
|
|
|
/* Open the dataset */
|
|
if ((dataset_id = H5Dopen( file_id , "Dataset2" ))<0) goto out;
|
|
|
|
/* Set new dimensions for the array. */
|
|
if (H5Dset_extent( dataset_id, dims_new )<0) goto out;
|
|
|
|
/* Get the space. */
|
|
if ((space_id = H5Dget_space( dataset_id ))<0) goto out;
|
|
|
|
/* Get dimensions. */
|
|
if (H5Sget_simple_extent_dims( space_id, dims_out, NULL )<0) goto out;
|
|
|
|
if ( dims_out[0] != dims_new[0] )
|
|
goto out;
|
|
|
|
/* Read the new dataset. */
|
|
if (H5Dread( dataset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf1 )<0) goto out;
|
|
|
|
/* Compare the read array with the original array */
|
|
for( i = 0; i < (int)dims_out[0]; i++ )
|
|
{
|
|
for( j = 0; j < (int)dims_out[1]; j++ )
|
|
{
|
|
if ( buf1[i][j] != data[i][j] ) {
|
|
goto out;
|
|
}
|
|
}
|
|
}
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Set new dimensions for the array; expand it again
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
/* Set new dimensions for the array. */
|
|
if (H5Dset_extent( dataset_id , dims )<0) goto out;
|
|
|
|
/* Get the space. */
|
|
if ((space_id = H5Dget_space( dataset_id ))<0) goto out;
|
|
|
|
/* Get dimensions. */
|
|
if (H5Sget_simple_extent_dims( space_id, dims_out, NULL )<0) goto out;
|
|
|
|
if ( dims_out[0] != dims[0] )
|
|
goto out;
|
|
|
|
/* Read the new dataset. */
|
|
if (H5Dread( dataset_id, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf2 )<0) goto out;
|
|
|
|
/* Compare the read array with the original array */
|
|
for( i = 0; i < (int)dims_out[0]; i++ )
|
|
{
|
|
for( j = 0; j < (int)dims_out[1]; j++ )
|
|
{
|
|
if ( i >= 70 || j >= 70 )
|
|
{
|
|
if ( buf2[i][j] != 0 ) {
|
|
goto out;
|
|
}
|
|
}
|
|
else
|
|
{
|
|
if ( buf2[i][j] != data[i][j] ) {
|
|
goto out;
|
|
}
|
|
}
|
|
|
|
}
|
|
}
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Close/release resources
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
H5Dclose( dataset_id );
|
|
H5Sclose( space_id );
|
|
|
|
|
|
H5Fclose( file_id );
|
|
|
|
PASSED();
|
|
|
|
|
|
return 0;
|
|
|
|
|
|
out:
|
|
H5Dclose( dataset_id );
|
|
H5Sclose( space_id );
|
|
H5Pclose( plist_id );
|
|
H5Fclose( file_id );
|
|
H5_FAILED();
|
|
return 1;
|
|
|
|
}
|
|
|
|
|
|
|