diff options
Diffstat (limited to 'src/H5Smpio.c')
| -rw-r--r-- | src/H5Smpio.c | 303 |
1 files changed, 0 insertions, 303 deletions
diff --git a/src/H5Smpio.c b/src/H5Smpio.c index b46d347..0a9a7bf 100644 --- a/src/H5Smpio.c +++ b/src/H5Smpio.c @@ -191,7 +191,6 @@ H5S_mpio_none_type( const H5S_t UNUSED *space, size_t UNUSED elmt_size, * *------------------------------------------------------------------------- */ -#ifndef AKC_OLD static herr_t H5S_mpio_hyper_type( const H5S_t *space, size_t elmt_size, /* out: */ @@ -453,308 +452,6 @@ done: #endif FUNC_LEAVE_NOAPI(ret_value); } -#else -/* keep this old code for now. */ -static herr_t -H5S_mpio_hyper_type( const H5S_t *space, size_t elmt_size, - /* out: */ - MPI_Datatype *new_type, - size_t *count, - hsize_t *extra_offset, - hbool_t *is_derived_type ) -{ - H5S_sel_iter_t sel_iter; /* Selection iteration info */ - hbool_t sel_iter_init=0; /* Selection iteration info has been initialized */ - - struct dim { /* less hassle than malloc/free & ilk */ - hssize_t start; - hsize_t strid; - hsize_t block; - hsize_t xtent; - hsize_t count; - } d[H5S_MAX_RANK]; - - int i; - int offset[H5S_MAX_RANK]; - int max_xtent[H5S_MAX_RANK]; - H5S_hyper_dim_t *diminfo; /* [rank] */ - int rank; - int block_length[2]; - MPI_Datatype inner_type, outer_type, old_type[2]; - MPI_Aint extent_len, displacement[2]; - int mpi_code; /* MPI return code */ - herr_t ret_value = SUCCEED; - - FUNC_ENTER_NOAPI_NOINIT(H5S_mpio_hyper_type); - - /* Check args */ - assert (space); - assert(sizeof(MPI_Aint) >= sizeof(elmt_size)); - if (0==elmt_size) - goto empty; - - /* Initialize selection iterator */ - if (H5S_select_iter_init(&sel_iter, space, elmt_size)<0) - HGOTO_ERROR (H5E_DATASPACE, H5E_CANTINIT, FAIL, "unable to initialize selection iterator"); - sel_iter_init=1; /* Selection iteration info has been initialized */ - - /* Abbreviate args */ - diminfo=sel_iter.u.hyp.diminfo; - assert (diminfo); - - /* make a local copy of the dimension info so we can operate with them */ - - /* Check if this is a "flattened" regular hyperslab selection */ - if(sel_iter.u.hyp.iter_rank!=0 && sel_iter.u.hyp.iter_rank<space->extent.rank) { - /* Flattened selection */ - rank=sel_iter.u.hyp.iter_rank; - assert (rank >= 0 && rank<=H5S_MAX_RANK); /* within array bounds */ - if (0==rank) - goto empty; - -#ifdef H5Smpi_DEBUG - HDfprintf(stderr, "%s: Flattened selection\n",FUNC); -#endif - for ( i=0; i<rank; ++i) { - d[i].start = diminfo[i].start+sel_iter.u.hyp.sel_off[i]; - d[i].strid = diminfo[i].stride; - d[i].block = diminfo[i].block; - d[i].count = diminfo[i].count; - d[i].xtent = sel_iter.u.hyp.size[i]; -#ifdef H5Smpi_DEBUG - HDfprintf(stderr, "%s: start=%Hd stride=%Hu count=%Hu block=%Hu xtent=%Hu", - FUNC, d[i].start, d[i].strid, d[i].count, d[i].block, d[i].xtent ); - if (i==0) - HDfprintf(stderr, " rank=%d\n", rank ); - else - HDfprintf(stderr, "\n" ); -#endif - if (0==d[i].block) - goto empty; - if (0==d[i].count) - goto empty; - if (0==d[i].xtent) - goto empty; - } - } /* end if */ - else { - /* Non-flattened selection */ - rank = space->extent.rank; - assert (rank >= 0 && rank<=H5S_MAX_RANK); /* within array bounds */ - if (0==rank) - goto empty; - -#ifdef H5Smpi_DEBUG - HDfprintf(stderr, "%s: Non-flattened selection\n",FUNC); -#endif - for ( i=0; i<rank; ++i) { - d[i].start = diminfo[i].start+space->select.offset[i]; - d[i].strid = diminfo[i].stride; - d[i].block = diminfo[i].block; - d[i].count = diminfo[i].count; - d[i].xtent = space->extent.size[i]; -#ifdef H5Smpi_DEBUG - HDfprintf(stderr, "%s: start=%Hd stride=%Hu count=%Hu block=%Hu xtent=%Hu", - FUNC, d[i].start, d[i].strid, d[i].count, d[i].block, d[i].xtent ); - if (i==0) - HDfprintf(stderr, " rank=%d\n", rank ); - else - HDfprintf(stderr, "\n" ); -#endif - if (0==d[i].block) - goto empty; - if (0==d[i].count) - goto empty; - if (0==d[i].xtent) - goto empty; - } - } /* end else */ - -/********************************************************************** - Compute array "offset[rank]" which gives the offsets for a multi- - dimensional array with dimensions "d[i].xtent" (i=0,1,...,rank-1). -**********************************************************************/ - offset[rank-1] = 1; - max_xtent[rank-1] = d[rank-1].xtent; -#ifdef H5Smpi_DEBUG - i=rank-1; - HDfprintf(stderr, " offset[%2d]=%d; max_xtent[%2d]=%d\n", - i, offset[i], i, max_xtent[i]); -#endif - for (i=rank-2; i>=0; --i) { - offset[i] = offset[i+1]*d[i+1].xtent; - max_xtent[i] = max_xtent[i+1]*d[i].xtent; -#ifdef H5Smpi_DEBUG - HDfprintf(stderr, " offset[%2d]=%d; max_xtent[%2d]=%d\n", - i, offset[i], i, max_xtent[i]); -#endif - } - - /* Create a type covering the selected hyperslab. - * Multidimensional dataspaces are stored in row-major order. - * The type is built from the inside out, going from the - * fastest-changing (i.e., inner) dimension * to the slowest (outer). */ - -/******************************************************* -* Construct contig type for inner contig dims: -*******************************************************/ -#ifdef H5Smpi_DEBUG - HDfprintf(stderr, "%s: Making contig type %d MPI_BYTEs\n", FUNC,elmt_size ); - for (i=rank-1; i>=0; --i) - HDfprintf(stderr, "d[%d].xtent=%Hu \n", i, d[i].xtent); -#endif - if (MPI_SUCCESS != (mpi_code= MPI_Type_contiguous( (int)elmt_size, MPI_BYTE, &inner_type ))) - HMPI_GOTO_ERROR(FAIL, "MPI_Type_contiguous failed", mpi_code); - -/******************************************************* -* Construct the type by walking the hyperslab dims -* from the inside out: -*******************************************************/ - for ( i=rank-1; i>=0; --i) { -#ifdef H5Smpi_DEBUG - HDfprintf(stderr, "%s: Dimension i=%d \n" - "count=%Hu block=%Hu stride=%Hu\n", - FUNC, i, d[i].count, d[i].block, d[i].strid ); -#endif - -#ifdef H5Smpi_DEBUG - HDfprintf(stderr, "%s: i=%d Making vector-type \n", FUNC,i); -#endif - /**************************************** - * Build vector in current dimension: - ****************************************/ - mpi_code =MPI_Type_vector((int)(d[i].count), /* count */ - (int)(d[i].block), /* blocklength */ - (int)(d[i].strid), /* stride */ - inner_type, /* old type */ - &outer_type ); /* new type */ - - MPI_Type_free( &inner_type ); - if (mpi_code!=MPI_SUCCESS) - HMPI_GOTO_ERROR(FAIL, "couldn't create MPI vector type", mpi_code); - - displacement[1] = (MPI_Aint)elmt_size * max_xtent[i]; - if(MPI_SUCCESS != (mpi_code = MPI_Type_extent(outer_type, &extent_len))) - HMPI_GOTO_ERROR(FAIL, "MPI_Type_extent failed", mpi_code); - - /************************************************* - * Restructure this datatype ("outer_type") - * so that it still starts at 0, but its extent - * is the full extent in this dimension. - *************************************************/ - if ((int)extent_len < displacement[1]) { - -#ifdef H5Smpi_DEBUG - HDfprintf(stderr, "%s: i=%d Extending struct type\n" - "***displacements: 0, %d\n", FUNC, i, displacement[1]); -#endif - -#ifdef H5_HAVE_MPI2 /* have MPI-2 (this function is not included in MPICH) */ - mpi_code = MPI_Type_create_resized - ( outer_type, /* old type */ - 0, /* blocklengths */ - displacement[1], /* displacements */ - &inner_type); /* new type */ -#else /* do not have MPI-2 */ - block_length[0] = 1; - block_length[1] = 1; - - displacement[0] = 0; - - old_type[0] = outer_type; - old_type[1] = MPI_UB; - mpi_code = MPI_Type_struct ( 2, /* count */ - block_length, /* blocklengths */ - displacement, /* displacements */ - old_type, /* old types */ - &inner_type); /* new type */ -#endif - - MPI_Type_free (&outer_type); - if (mpi_code!=MPI_SUCCESS) - HMPI_GOTO_ERROR(FAIL, "couldn't resize MPI vector type", mpi_code); - } - else { - inner_type = outer_type; - } - } /* end for */ -/*************************** -* End of loop, walking -* thru dimensions. -***************************/ - - - /* At this point inner_type is actually the outermost type, even for 0-trip loop */ - -/*************************************************************** -* Final task: create a struct which is a "clone" of the -* current struct, but displaced according to the d[i].start -* values given in the hyperslab description: -***************************************************************/ - displacement[0] = 0; - for (i=rank-1; i>=0; i--) - displacement[0] += d[i].start * offset[i]; - -printf("dumping MPI_BYTE\n"); -printdatatype(MPI_INT); -printdatatype(MPI_BYTE); - if (displacement[0] > 0) { - displacement[0] *= elmt_size; - block_length[0] = 1; - old_type[0] = inner_type; - -#ifdef H5Smpi_DEBUG - HDfprintf(stderr, "%s: Making final struct\n***count=1:\n", FUNC); - HDfprintf(stderr, "\tblocklength[0]=%d; displacement[0]=%d\n", - block_length[0], displacement[0]); -#endif - - - if (MPI_SUCCESS != (mpi_code= MPI_Type_struct( 1, /* count */ - block_length, /* blocklengths */ - displacement, /* displacements */ - old_type, /* old type */ - new_type )) /* new type */ - ) - HMPI_GOTO_ERROR(FAIL, "couldn't create MPI struct type", mpi_code); - - if (MPI_SUCCESS != (mpi_code= MPI_Type_free (&old_type[0]))) - HMPI_GOTO_ERROR(FAIL, "MPI_Type_free failed", mpi_code); - } - else { - *new_type = inner_type; - } - - if (MPI_SUCCESS != (mpi_code= MPI_Type_commit( new_type ))) - HMPI_GOTO_ERROR(FAIL, "MPI_Type_commit failed", mpi_code); - - /* fill in the remaining return values */ - *count = 1; /* only have to move one of these suckers! */ - *extra_offset = 0; - *is_derived_type = 1; - HGOTO_DONE(SUCCEED); - -empty: - /* special case: empty hyperslab */ - *new_type = MPI_BYTE; - *count = 0; - *extra_offset = 0; - *is_derived_type = 0; - -done: - /* Release selection iterator */ - if(sel_iter_init) { - if (H5S_SELECT_ITER_RELEASE(&sel_iter)<0) - HDONE_ERROR (H5E_DATASPACE, H5E_CANTRELEASE, FAIL, "unable to release selection iterator"); - } /* end if */ - -#ifdef H5Smpi_DEBUG - HDfprintf(stderr, "Leave %s, count=%ld is_derived_type=%d\n", - FUNC, *count, *is_derived_type ); -#endif - FUNC_LEAVE_NOAPI(ret_value); -} -#endif /*------------------------------------------------------------------------- |