c_sf.c

Go to the documentation of this file.

/*
 * c_sf.c - wrappers for the table reading in the chimaera code
 *
 *  W. Robinson, SAIC, 12 Dec 2018
 */
#include <string.h>
#include "mfhdf.h"
 
/*
 * first, the sfstart replacement
 */
int32_t c_sfstart_( char *name, int32_t *access, int nam_len )
 /*  return the file ID
     name  I file name
     access  I access mode read is DFACC_READ (also_WRITE, CREATE)
     nam_len  I  length of ftn string
  */
     
  {
  char *fn, *bptr;
  int32_t ret;
 
  fn = (char *) malloc( nam_len + 1 * sizeof( char ) );
  fn[ nam_len ] = (char) 0;
  strncpy( fn, name, nam_len );
  bptr = strchr( fn, ' ' );
  if( bptr != NULL ) { 
    *bptr = 0;
  }
 
  ret = (int32_t) SDstart( fn, *access);
  free(fn);
 
  return(ret);
  }
/*
 * next the sfselect
 */
int32_t c_sfselect_( int32_t *file_id, int32_t *index )
 /*  return the SDS ID
     file_id  I  file ID
     index  I  index from c_sfn2index_
  */
  {
  int32_t fid;
 
  fid = SDselect( *file_id, *index );
  return( fid );
  }
/*
 * next the sfn2index
 */
int32_t c_sfn2index_( int32_t *fid, char *name, int32_t nam_len )
 /*  return the index of the SDS name
     name  I  SDS name
     nam_len  I  length of SDS name
  */
  {
  char *fn, *bptr;
  int32_t ret;
 
  fn = (char *) malloc( nam_len + 1 * sizeof( char ) );
  fn[ nam_len ] = (char) 0;
  strncpy( fn, name, nam_len );
  bptr = strchr( fn, ' ' );
  if( bptr != NULL ) { 
    *bptr = 0;
  }
 
  ret = SDnametoindex( *fid, fn );
  free(fn);
 
  return(ret);
  }
/*
 * next the sfginfo
 */
int32_t c_sfginfo_( int32_t *vid, char *name, int32_t *rank, int32_t *dims, 
  int32_t *nt, int32_t *nattr, int32_t nam_len )
 /*  get info about a SDS
     return 0 if good 
     vid  I SDS ID
     name  O  SDS name
     rank  O  # dims
     dims  O  length of each dim
     nt   O  data type
     nattr  O  # attributes
     nam_len  I  len of name string  
  */
  {
  char *fn, *bptr;
  int32_t ret;
  int32 cdims[100], i;
  int32 rank32, nt32, nattr32;
 
  fn = NULL;
  if( nam_len != 0 ) {
 
    fn = (char *) malloc( nam_len + 1 * sizeof( char ) );
    fn[ nam_len ] = (char) 0;
    strncpy( fn, name, nam_len );
    bptr = strchr( fn, ' ' );
    if( bptr != NULL ) {
      *bptr = 0;
    }
  }
 
  ret = SDgetinfo((int32) *vid, fn, &rank32, cdims, &nt32, &nattr32);
 
    for(i = 0; i < rank32; i++)
        dims[i] = cdims[rank32 - i - 1];
  if( nam_len != 0 ) free(fn);
  *rank  = (int32_t) rank32;
  *nt    = (int32_t) nt32;
  *nattr = (int32_t) nattr32;
 
  return( ret );
  } 
/*
 * next the sfrdata
 */
int32_t c_sfrdata_( int32_t *vid, int32_t *start, int32_t *stride, 
  int32_t *end, VOIDP values )
 /*  read the data slice
     return 0 if good
     vid  I SDS ID
     start, stride, end  I  the start, stride and # too read in each dim
     values  O  space for read array
  */
  {
  int32_t ret;
  int32_t   i, rank, dims[100], nt, nattrs, status;
  int32_t   cstart[100], cstride[100], cend[100];
  int     nostride;
 
  status = SDgetinfo( *vid, NULL, &rank, dims, &nt, &nattrs);
  if(status == FAIL) return FAIL;
 
  nostride = TRUE;
 
  for(i = 0; i < rank; i++) {
    cstart[i] = start[rank - i - 1];
    cend[i]   = end[rank - i - 1];
    if((cstride[i] = stride[rank - i - 1]) != 1) nostride = FALSE;
  }
 
  ret = (int32_t) SDreaddata(*vid, cstart, (nostride? NULL : cstride), 
    cend, values);
 
  return( ret );
  }
/*
 * next the sfendacc
 */
int32_t c_sfendacc_( int32_t *vid )
 /*  close SDS id
     return 0 if good
     vid  I SDS ID
  */
  {
  int32_t ret;
 
  ret = SDendaccess( *vid );
 
  return( ret );
  }
/*
 * next the sfend
 */
int32_t c_sfend_( int32_t *id )
 /*  close file
     return 0 if good
     id  I  file id
  */
  {
  int32_t ret;
 
  ret = SDend( *id );
 
  return( ret );
  }