Uniform Grids

Functions
F5_API F5Path *	F5Rcreate_uniform (hid_t File_id, double time, const char *gridname, const void *origin, const void *spacing, hsize_t *dims, const char *coordinate_system)
F5_API F5Path *	F5Rcreate_uniform_cartesian3D (hid_t file_id, double time, const char *gridname, const F5_vec3_point_t *origin, const F5_vec3_float_t *spacing, hsize_t dims[3], const char *coordinate_system)
F5_API F5Path *	F5Rcreate_uniform_cartesian3Dbbox (hid_t File_id, double time, const char *gridname, const F5_vec3_point_t *start, const F5_vec3_point_t *end, hsize_t dims[3], const char *coordinate_system)
F5_API F5Path *	F5Fwrite_uniform_cartesian3D (hid_t file_id, double time, const char *gridname, const F5_vec3_point_t *origin, const F5_vec3_float_t *spacing, hsize_t dims[3], const char *fieldname, hid_t fieldtype, const void *dataPtr, const char *coordinate_system, hid_t property_id)
F5_API int	F5write_uniform_cartesian3Dv (hid_t file_id, double time, const char *gridname, const F5_vec3_point_t *origin, const F5_vec3_float_t *spacing, hsize_t dims[3], const char *coordinate_system, hid_t property_id,...)
F5_API int	F5write_uniform_cartesian3Dvs (hid_t file_id, double time, const char *gridname, const F5_vec3_point_t *origin, const F5_vec3_float_t *spacing, hsize_t dims[3], const char *coordinate_system, hid_t property_id,...)
F5_API F5Path *	F5Fwrite_uniform_cartesian3Ds (hid_t file_id, double time, const char *gridname, const F5_vec3_point_t *origin, const F5_vec3_float_t *spacing, hsize_t dims[3], const char *fieldname, hid_t fieldtype, const void **dataPtr, const char *coordinate_system, hid_t property_id)
F5_API F5Path *	F5Rcreate_uniform_sparse (hid_t file_id, double time, const char *gridname, const void *origin, const void *spacing, hsize_t *dims, const char *coordinate_system, const hsize_t *sparse_idx_Ptr, const hsize_t sparse_size, hid_t property_id)
F5_API F5Path *	F5Rcreate_uniform_sparse2 (hid_t file_id, double time, const char *gridname, const F5_vec3_double_t *origin, const F5_vec3_double_t *spacing, hsize_t *dims, const char *coordinate_system, const hsize_t *sparse_idx_Ptr, hsize_t sparse_size, hid_t property_id)

Detailed Description

Function Documentation

F5Fwrite_uniform_cartesian3D()

F5_API F5Path * F5Fwrite_uniform_cartesian3D	(	hid_t	file_id,
		double	time,
		const char *	gridname,
		const F5_vec3_point_t *	origin,
		const F5_vec3_float_t *	spacing,
		hsize_t	dims[3],
		const char *	fieldname,
		hid_t	fieldtype,
		const void *	dataPtr,
		const char *	coordinate_system,
		hid_t	property_id )

Write a data field given on uniform cartesian grid into file.

@param coordinate_system The coordinate system.
            May be NULL to refer to a standard chart.

Returns

A new F5Path pointer. Call F5close() if it is no longer used.

Definition at line 271 of file F5uniform.c.

{
F5Path*f = F5Rcreate_uniform_cartesian3D(file_id, time,
                                         gridname,
                                         origin,
                                         spacing,
                                         dims,
                                         coordinate_system);
        /* write the data into the dataset */
 
         if ( F5Fwrite(f, fieldname, 3, dims, fieldtype, 0, dataPtr, prop_id) != F5_SUCCESS)
         {
                 F5close(f);
                 return NULL;
         }
         F5Cset (f, fieldname, FIBER_REGULAR_GRID_URL "StandardCartesian/Uniform/VertexCentered");
 
         return f;
}

References F5_SUCCESS, F5close(), F5Cset(), F5Fwrite(), F5Rcreate_uniform_cartesian3D(), and FIBER_REGULAR_GRID_URL.

F5Fwrite_uniform_cartesian3Ds()

F5_API F5Path * F5Fwrite_uniform_cartesian3Ds	(	hid_t	file_id,
		double	time,
		const char *	gridname,
		const F5_vec3_point_t *	origin,
		const F5_vec3_float_t *	spacing,
		hsize_t	dims[3],
		const char *	fieldname,
		hid_t	fieldtype,
		const void **	dataPtr,
		const char *	coordinate_system,
		hid_t	property_id )

Write a data field given on uniform cartesian grid for a given coordinate system into file. If the data is a compound data type, then it is given as distinct homogeneous arrays in memory.

Parameters

coordinate_system

The coordinate system. May be NULL to refer to a standard chart.

Returns

A new F5Path pointer. Call F5close() if it is no longer used.

Definition at line 301 of file F5uniform.c.

{
F5Path*f = F5Rcreate_uniform_cartesian3D(file_id, time,
                                         gridname,
                                         origin,
                                         spacing,
                                         dims,
                                         coordinate_system);
 
        /* write the data into the dataset */
 
        if (! F5Fwrites(f, fieldname, 3, dims, fieldtype, 0, dataPtr, property_id) )
        {
                F5close(f);
                return NULL;
        }
        return f;
}

References F5close(), F5Fwrites(), and F5Rcreate_uniform_cartesian3D().

F5Rcreate_uniform()

F5_API F5Path * F5Rcreate_uniform	(	hid_t	File_id,
		double	time,
		const char *	gridname,
		const void *	origin,
		const void *	spacing,
		hsize_t *	dims,
		const char *	coordinate_system )

Create a regular grid which is uniform in the specified coordinate system.

Parameters

origin	The origin of the grid, must be of the coordinate type that is associated with the coordinate system.
spacing	The distance among neighboring points of the grid, must be of the vector type of the coordinate system, i.e. the coordinate type that stores coordinate differences.
coordinate_system	The coordinate system. May be NULL to refer to the standard cartesian chart.

Returns

A new F5Path pointer. Call F5close() if it is no longer used.

References F5_API.

F5Rcreate_uniform_cartesian3D()

F5_API F5Path * F5Rcreate_uniform_cartesian3D	(	hid_t	file_id,
		double	time,
		const char *	gridname,
		const F5_vec3_point_t *	origin,
		const F5_vec3_float_t *	spacing,
		hsize_t	dims[3],
		const char *	coordinate_system )

Create a regular three-dimensional grid which is uniform in the specified coordinate system.

Parameters

coordinate_system

The coordinate system. May be NULL to refer to a standard chart.

Returns

A new F5Path pointer. Call F5close() if it is no longer used.

Definition at line 73 of file F5uniform.c.

{
F5Path*f = F5Rcreate_cartesian_3D(File_id, time, gridname, coordinate_system);
F5_vec3_point_t end, center;
 
        F5Cset (f, FIBER_HDF5_POSITIONS_STRING, FIBER_REGULAR_GRID_URL "StandardCartesian/Uniform/");
        F5Fwrite_linear(f, FIBER_HDF5_POSITIONS_STRING,
                                 3, dims, F5T_COORD3_FLOAT, origin, spacing);
 
        end.x = origin->x + spacing->x*(dims[0]-1);
        end.y = origin->y + spacing->y*(dims[1]-1);
        end.z = origin->z + spacing->z*(dims[2]-1);
 
        center.x = 0.5*(origin->x + end.x);
        center.y = 0.5*(origin->y + end.y);
        center.z = 0.5*(origin->z + end.z);
 
 
        {F5_vec3_point_t minmax[2];
                minmax[0].x = origin->x;        minmax[0].y = origin->y;        minmax[0].z = origin->z;
                minmax[1].x = end.x;    minmax[1].y = end.y;    minmax[1].z = end.z;
 
                F5Fset_range(f, NULL, minmax);
        }
 
        F5Fset_average(f, NULL, &center);
 
        F5Fclose(f);
 
        return f;
}

References F5Cset(), F5Fclose(), F5Fset_average(), F5Fset_range(), F5Fwrite_linear(), F5Rcreate_cartesian_3D(), F5T_COORD3_FLOAT, FIBER_HDF5_POSITIONS_STRING, FIBER_REGULAR_GRID_URL, F5_point3_float_t::x, F5_vec3_float_t::x, F5_point3_float_t::y, F5_vec3_float_t::y, F5_point3_float_t::z, and F5_vec3_float_t::z.

Referenced by F5Fwrite_uniform_cartesian3D(), F5Fwrite_uniform_cartesian3Ds(), F5Rcreate_uniform_sparse(), F5write_uniform_cartesian3Dv(), and F5write_uniform_cartesian3Dvs().

F5Rcreate_uniform_cartesian3Dbbox()

F5_API F5Path * F5Rcreate_uniform_cartesian3Dbbox	(	hid_t	File_id,
		double	time,
		const char *	gridname,
		const F5_vec3_point_t *	start,
		const F5_vec3_point_t *	end,
		hsize_t	dims[3],
		const char *	coordinate_system )

Definition at line 110 of file F5uniform.c.

{
F5Path*f = F5Rcreate_cartesian_3D(File_id, time, gridname, coordinate_system);
F5_vec3_float_t spacing;
F5_vec3_point_t center;
 
        F5Cset (f, FIBER_HDF5_POSITIONS_STRING, FIBER_REGULAR_GRID_URL "StandardCartesian/Uniform/");
 
        spacing.x = (end->x - start->x) / (dims[0]-1);
        spacing.y = (end->y - start->y) / (dims[1]-1);
        spacing.z = (end->z - start->z) / (dims[2]-1);
 
        center.x = 0.5*(start->x + end->x);
        center.y = 0.5*(start->y + end->y);
        center.z = 0.5*(start->z + end->z);
 
        F5Fwrite_linear(f, FIBER_HDF5_POSITIONS_STRING,
                                 3, dims, F5T_COORD3_FLOAT, start, &spacing);
 
        {
        F5_vec3_point_t minmax[2];
                minmax[0].x = start->x; minmax[0].y = start->y; minmax[0].z = start->z;
                minmax[1].x = end->x;   minmax[1].y = end->y;   minmax[1].z = end->z;
                F5Fset_range(f, NULL, minmax);
        }
        F5Fset_average(f, NULL, &center);
 
        F5Fclose(f);
 
        return f;
}

References F5Cset(), F5Fclose(), F5Fset_average(), F5Fset_range(), F5Fwrite_linear(), F5Rcreate_cartesian_3D(), F5T_COORD3_FLOAT, FIBER_HDF5_POSITIONS_STRING, FIBER_REGULAR_GRID_URL, F5_point3_float_t::x, F5_vec3_float_t::x, F5_point3_float_t::y, F5_vec3_float_t::y, F5_point3_float_t::z, and F5_vec3_float_t::z.

F5Rcreate_uniform_sparse()

F5_API F5Path * F5Rcreate_uniform_sparse	(	hid_t	file_id,
		double	time,
		const char *	gridname,
		const void *	origin,
		const void *	spacing,
		hsize_t *	dims,
		const char *	coordinate_system,
		const hsize_t *	sparse_idx_Ptr,
		const hsize_t	sparse_size,
		hid_t	property_id )

Definition at line 545 of file F5uniform.c.

{
F5Path*f = F5Rcreate_uniform_cartesian3D(file_id, time,
                                         gridname,
                                         (const F5_point3_float_t*)origin,
                                         (const F5_vec3_float_t*)spacing,
                                         dims,
                                         coordinate_system);
 
        /* write the data into the dataset, todo: code-review, how to handle rank? 3 vs 1 */
        f->Field_hid = F5Lwrite(f->Representation_hid, FIBER_HDF5_SPARSE,
                                1, &sparse_size,
                                F5file_type(f, H5T_NATIVE_HSIZE), H5T_NATIVE_HSIZE, sparse_idx_Ptr,
                                F5LTmake_enum_type(f),
                                property_id); 
 
        /*
        //      if (! F5Fwrite(f, FIBER_HDF5_SPARSE, 1, &sparse_size, H5T_NATIVE_HSIZE, 0, sparse_idx_Ptr, property_id) )
        //{
        //              F5close(f);
        //              return NULL;
                //}
                */
        return f;
}

References Dims::dims, F5_API, F5file_type(), F5LTmake_enum_type(), F5Lwrite(), F5Rcreate_uniform_cartesian3D(), FIBER_HDF5_SPARSE, F5Path::Field_hid, and F5Path::Representation_hid.

F5Rcreate_uniform_sparse2()

F5_API F5Path * F5Rcreate_uniform_sparse2	(	hid_t	file_id,
		double	time,
		const char *	gridname,
		const F5_vec3_double_t *	origin,
		const F5_vec3_double_t *	spacing,
		hsize_t *	dims,
		const char *	coordinate_system,
		const hsize_t *	sparse_idx_Ptr,
		hsize_t	sparse_size,
		hid_t	property_id )

Definition at line 579 of file F5uniform.c.

{
F5Path*f = F5Rcreate_cartesian_nD(file_id, time,
                                  gridname, 1, coordinate_system);
 
  if (! F5Fwrite(f, FIBER_HDF5_SPARSE, 1, &sparse_size, H5T_NATIVE_HSIZE, 0, sparse_idx_Ptr, property_id) )
    {
      //     F5close(f);
      // return NULL;
    }
 
  assert( F5Fset_attribute( f, FIBER_UNIFORMARRAY_BASE_ATTRIBUTE, origin, 3, H5T_NATIVE_DOUBLE, 0 ) );
  assert( F5Fset_attribute( f, FIBER_UNIFORMARRAY_DELTA_ATTRIBUTE, spacing, 3, H5T_NATIVE_DOUBLE, 0 ) );
  assert( F5Fset_attribute( f, FIBER_UNIFORMARRAY_ELEMENTS_ATTRIBUTE, dims, 3, H5T_NATIVE_HSIZE, 0 ) );
 
  return f;
}

References Dims::dims, F5_API, F5Fset_attribute(), F5Fwrite(), F5Rcreate_cartesian_nD(), FIBER_HDF5_SPARSE, FIBER_UNIFORMARRAY_BASE_ATTRIBUTE, FIBER_UNIFORMARRAY_DELTA_ATTRIBUTE, and FIBER_UNIFORMARRAY_ELEMENTS_ATTRIBUTE.

F5write_uniform_cartesian3Dv()

F5_API int F5write_uniform_cartesian3Dv	(	hid_t	file_id,
		double	time,
		const char *	gridname,
		const F5_vec3_point_t *	origin,
		const F5_vec3_float_t *	spacing,
		hsize_t	dims[3],
		const char *	coordinate_system,
		hid_t	property_id,
			... )

Write a bunch of data fields given on the same uniform cartesian grid. This variation of the uniform grid writer function uses a variable number of arguments list. Each field is defined by a triple of

  const char *fieldname
  hid_t fieldtype
  const void  *dataPtr

If either the fieldname or the dataPtr are NULL, then the traversal of the argument is terminated.

Parameters

coordinate_system

The coordinate system. May be NULL to refer to a standard chart.

Returns

The number of successfully written fields.

Definition at line 332 of file F5uniform.c.

{
        const char* fieldname;
        hid_t fieldtype;
        const void* dataPtr;
        va_list vl;
        int NumOfWrittenFields = 0;
 
        /* Create the Field */
        F5Path* f = F5Rcreate_uniform_cartesian3D(file_id, time,
                                                  gridname,
                                                  origin,
                                                  spacing,
                                                  dims,
                                                  coordinate_system);
 
        /* Init the step through the argument list */
        va_start(vl, property_id);
 
        /* Step through the argument list */
        for(;;)
        {
                /* Get the Name of the field */
                fieldname = va_arg(vl, const char*);
                /* Get the type of the field */
                fieldtype = va_arg(vl, hid_t);
                /* Get the pointer to the field */
                dataPtr = va_arg(vl, const void*);
 
                /* Break? */
                if ( (fieldname == NULL) || (dataPtr == NULL) )
                        break;
 
                /* Write the stuff */
                F5Fwrite(f, fieldname, 3, dims, fieldtype, 0, dataPtr, property_id);
                NumOfWrittenFields++;
        }
 
        /* Close the Step through */
        va_end(vl);
 
        /* Close the field */
        F5close(f);
 
        return NumOfWrittenFields;
}

References F5close(), F5Fwrite(), and F5Rcreate_uniform_cartesian3D().

F5write_uniform_cartesian3Dvs()

F5_API int F5write_uniform_cartesian3Dvs	(	hid_t	file_id,
		double	time,
		const char *	gridname,
		const F5_vec3_point_t *	origin,
		const F5_vec3_float_t *	spacing,
		hsize_t	dims[3],
		const char *	coordinate_system,
		hid_t	property_id,
			... )

Write a bunch of data fields given on the same uniform cartesian grid. This variation of the uniform grid writer function uses a variable number of arguments list. Each field is defined by a triple of

  const char *fieldname
  hid_t fieldtype
  const void **dataPtr

If either the fieldname or the dataPtr are NULL, then the traversal of the argument is terminated.

Parameters

coordinate_system

The coordinate system. May be NULL to refer to a standard chart.

Returns

The number of successfully written fields.

Definition at line 388 of file F5uniform.c.

{
        const char* fieldname;
        hid_t fieldtype;
        const void** dataPtr;
        va_list vl;
        int NumOfWrittenFields = 0;
 
        /* Create the Field */
        F5Path* f = F5Rcreate_uniform_cartesian3D(file_id, time,
                                                  gridname,
                                                  origin,
                                                  spacing,
                                                  dims,
                                                  coordinate_system);
 
        /* Init the step through the argument list */
        va_start(vl,  property_id);
 
        /* Step through the argument list */
        for(;;)
        {
                /* Get the Name of the field */
                fieldname = va_arg(vl, const char*);
                /* Get the type of the field */
                fieldtype = va_arg(vl, hid_t);
                /* Get the pointer to the field */
                dataPtr = va_arg(vl, const void**);
 
                /* Break? */
                if ( (fieldname == NULL) || (dataPtr == NULL) )
                        break;
 
                /* Write the stuff */
                F5Fwrites(f, fieldname, 3, dims, fieldtype, 0, dataPtr,  property_id);
                NumOfWrittenFields++;
        }
 
        /* Close the Step through */
        va_end(vl);
 
        /* Close the field */
        F5close(f);
 
        return NumOfWrittenFields;
}

References F5close(), F5Fwrites(), and F5Rcreate_uniform_cartesian3D().