F5LT.c File Reference

#include "F5L.h"
#include "F5B.h"
#include "F5X.h"
#include "F5T.h"
#include "F5A.h"
#include "F5defs.h"
#include "F5private.h"
#include "F5coordinates.h"
#include <stdlib.h>
#include <assert.h>
#include <string.h>

Go to the source code of this file.

Functions
ChartDomain_IDs *	F5FB_read_global_chart (F5Path *f)
int	F5LTget_index_depth (hid_t Top_hid)
int	F5LTget_skeleton_dimensionality (hid_t Top_hid)
int	F5LTset_dataspace (F5Path *f, int rank, hsize_t *dims)
int	F5LTexpand_dataspace (F5Path *f, int rank, hsize_t *dims)
hid_t	F5LTmake_enum_type (F5Path *fpath)
int	F5LTget_maximal_grid_refinement (hid_t ContentsGroup_id, const char *gridname, hsize_t *refinement)
int	F5LTget_maximal_skeleton_refinement (hid_t ContentsGroup_id, const char *gridname, hsize_t *refinement, int IndexDepth, int Dimensionality)
int	F5LTset_maximal_grid_refinement (hid_t ContentsGroup_id, const char *gridname, int IndexDepth, int Dimensionality, const hsize_t *refinement)
hid_t	F5LTcreate_topology (hid_t Grid_hid, const char *TopologyName, int IndexDepth, int SkeletonDimensionality, int DataspaceDimensionality, const hsize_t *refinement)
F5Path *	F5LTcreateV (hid_t File_id, const double *time, const char *gridname, ChartDomain_IDs *ChartDomain, ChartDomain_IDs *coord_creator(void *udata), void *udata, const char *coordinate_system, const char *TopologyName, int IndexDepth, int SkeletonDimensionality, int Dimensionality, const hsize_t *refinement)
	if (return f)
F5Path *	F5LTcreate (hid_t File_id, const double *time, const char *gridname, ChartDomain_IDs *ChartDomain, ChartDomain_IDs *coord_creator(), const char *coordinate_system, const char *TopologyName, int IndexDepth, int SkeletonDimensionality, int Dimensionality, const hsize_t *refinement)

Variables
	else

Function Documentation

F5FB_read_global_chart()

ChartDomain_IDs * F5FB_read_global_chart ( F5Path * f )

extern

Definition at line 14 of file F5B.c.

{
        return  F5B_read_global_chart(f->GlobalChart_hid, f->myChart->domain_name);
}

References ChartDomain_IDs::domain_name, F5B_read_global_chart(), F5Path::GlobalChart_hid, and F5Path::myChart.

Referenced by F5Bopen_coordinate_domain(), and F5LTcreateV().

F5LTexpand_dataspace()

int F5LTexpand_dataspace	(	F5Path *	f,
		int	rank,
		hsize_t *	dims )

Definition at line 120 of file F5LT.c.

{
hsize_t current_dims[FIBER_MAX_RANK];
int     current_rank, Mrank;
int     i;
 
        F5printf(30, "F5LTexpand_dataspace(,%d,)", rank);
        if (!f)
                return 1;
 
        F5Aget_ints( f->Topology_hid, FIBER_HDF5_RANK_ATTRIBUTE, &Mrank, 1);
        if (Mrank<0)
        {
                F5printf(30, "F5LTexpand_dataspace(,%d,): INVALID FILE: RANK attribute is not integer - cannot set dataspace!", rank);
                return 0;
        }
        if (Mrank>=0)
        {
                if (Mrank != rank)
                {
                        F5printf(30, "F5LTexpand_dataspace(,%d,): Dimensions incompatible with RANK attribute %d - inconsistent dataspace!", rank, Mrank);
                        puts("Set DATAPACE: Inconsistent rank!");
                        return 0;
                }
        }
 
        current_rank = F5LAget_dimensions(f->Topology_hid,
                                          FIBER_FIELD_DATASPACE_DIMENSIONS_ATTRIBUTE,
                                          current_dims);
 
        if (current_rank>0)
        {
                if (rank != current_rank)
                {
                        puts("Set DATAPACE: Inconsistent rank!");
                        return 0;
                }
 
                for(i=0; i<rank; i++)
                {
                        if (current_dims[i] < dims[i])
                                current_dims[i] = dims[i];
                }
        }
        else
                for(i=0; i<rank; i++)
                        current_dims[i] = dims[i];
 
        F5LAsave_dimensions(f->Topology_hid, FIBER_FIELD_DATASPACE_DIMENSIONS_ATTRIBUTE, rank, current_dims);
        F5printf(30, "F5LTexpand_dataspace(,rank=%d,): expand dimensions of topology.", rank);
 
        for(i=0; i<rank; i++)
                dims[i] = current_dims[i];
 
        return 1;
}

References F5Aget_ints(), F5LAget_dimensions(), F5LAsave_dimensions(), F5printf, FIBER_FIELD_DATASPACE_DIMENSIONS_ATTRIBUTE, FIBER_HDF5_RANK_ATTRIBUTE, FIBER_MAX_RANK, and F5Path::Topology_hid.

Referenced by F5Fexpand_dataspace().

F5LTget_skeleton_dimensionality()

int F5LTget_skeleton_dimensionality

(

hid_t

Top_hid

)

Definition at line 34 of file F5LT.c.

{
int     result = FIBER_HDF5_INVALID_SKELETON_DIMENSIONALITY;
        if (F5Aget_ints(Top_hid, FIBER_HDF5_SKELETON_DIMENSIONALITY_ATTRIB, &result, 1) <= 0)
                F5printf(5," *Warning* F5LTget_F5LTget_skeleton_dimensionality() could not find skeleton dimensionality,setting unknown -1");
 
        return result;
}

References F5Aget_ints(), F5printf, FIBER_HDF5_INVALID_SKELETON_DIMENSIONALITY, and FIBER_HDF5_SKELETON_DIMENSIONALITY_ATTRIB.

Referenced by F5iterate_topologies().

F5LTmake_enum_type()

hid_t F5LTmake_enum_type

(

F5Path *

fpath

)

Get the field enum type associated with this F5Path. May be used to check versioning information of field descriptions. \TODO Retrieve versioning from existing file, current implementation just writes.

Definition at line 179 of file F5LT.c.

{
        if (!fpath) return 0;
        if (fpath->field_enum_type_hid<=0)
        {
                assert(fpath->ContentsGroup_hid);
                F5printf(30, " F5LTmake_enum_type() - save type.");
                fpath->field_enum_type_hid = F5Tsave_F5field_enum(fpath->ContentsGroup_hid);
                F5printf(30, " F5LTmake_enum_type() - DONE saving type.");
        }
        return fpath->field_enum_type_hid;
}

References F5Path::ContentsGroup_hid, F5printf, F5Tsave_F5field_enum(), and F5Path::field_enum_type_hid.

Referenced by F5Fcreate(), F5FSwrite_fraction(), F5Fwrite_1D(), F5Fwrite_flexible(), F5Fwrite_fraction(), F5Fwrite_fraction_external(), F5Fwrite_fractionS(), F5Fwrite_linear(), F5Fwrite_linear_fraction(), F5FwriteIMAGE(), F5Fwrites(), F5FwriteX(), F5Rcreate_coordinate_topology(), F5Rcreate_edges_as_vertices_topology(), F5Rcreate_hexaedrons_as_vertices_topology(), F5Rcreate_rectilinear(), F5Rcreate_tetrahedrons_as_vertices_topology(), F5Rcreate_triangles_as_vertices_topology(), F5Rcreate_uniform_sparse(), F5write_particle_cartesian3Dv(), F5write_particle_cartesian3DvU(), and F5write_particle_double_cartesian3Dv().

F5LTset_dataspace()

int F5LTset_dataspace	(	F5Path *	f,
		int	rank,
		hsize_t *	dims )

Definition at line 43 of file F5LT.c.

{
hsize_t current_dims[FIBER_MAX_RANK];
int     current_rank, Mrank = -1;
 
        if (!f)
        {
                return 0;
        }
 
        F5printf(30, "F5LTset_dataspace(,rank=%d,)", rank);
 
        F5Aget_ints( f->Topology_hid, FIBER_HDF5_RANK_ATTRIBUTE, &Mrank, 1);
        if (Mrank<0)
        {
                F5printf(30, "F5LTset_dataspace(,%d,): INVALID FILE: RANK attribute is not integer - cannot set dataspace!", rank);
                return 0;
        }
 
        if (Mrank>=1)
        {
                if (Mrank != rank)
                {
                        if (rank==4)
                                F5printf(0, "F5LTset_dataspace(,%d,): Dimensions %dx%dx%dx%d incompatible with RANK attribute %d - inconsistent dataspace!",
                                         rank, (int)dims[0],(int)dims[1],(int)dims[2],(int)dims[3],
                                         Mrank);
                        else
                                F5printf(0, "F5LTset_dataspace(,%d,): Dimensions incompatible with RANK attribute %d - inconsistent dataspace!", rank, Mrank);
 
                        return 0;
                }
        }
 
        current_rank = F5LAget_dimensions(f->Topology_hid,
                                          FIBER_FIELD_DATASPACE_DIMENSIONS_ATTRIBUTE,
                                          current_dims);
 
        if (current_rank>0)
        {
        int     i;
                if (rank != current_rank)
                {
                        puts("F5LTset_dataspace(): Set DATAPACE: Inconsistent rank with DATASPACE_DIMENSIONS!");
                        return 0;
                }
 
                for(i=0; i<rank; i++)
                {
                        if (current_dims[i] != dims[i])
                        {
                                F5printf(5, "F5LTset_dataspace(,rank=%d,): new dimensions[%d] is %d, but already set to %d in the file as  %dx%dx%dx%d - aborting.",
                                         rank,              i, (int)(current_dims[i]), (int)(dims[i]),
                                         rank>0?(int)dims[0]:1, rank>1?(int)dims[1]:1, rank>2?(int)dims[2]:1, rank>3?(int)dims[3]:1);
                                return 0;
                        }
                }
                F5printf(5, "F5LTset_dataspace(,rank=%d,): dimensions are ok.", rank);
                return 1;
        }
        {int i;
        for(i=0; i<rank; i++)
        {
                if (dims[i] < 1)
                {
                        F5printf(10, "F5LTset_dataspace(,rank=%d,): do NOT set new topology dimensions, because dims[%d]=%d\n", rank, i, (int)dims[i] );
                        return 0;
                }
        }
        }
 
        F5LAsave_dimensions(f->Topology_hid, FIBER_FIELD_DATASPACE_DIMENSIONS_ATTRIBUTE, rank, dims);
        F5printf(30, "F5LTset_dataspace(,rank=%d,): set new dimensions to skeleton as %dx%dx%dx%d.", rank,
                 rank>0?(int)dims[0]:1, rank>1?(int)dims[1]:1, rank>2?(int)dims[2]:1, rank>3?(int)dims[3]:1);
        return 1;
}

References F5Aget_ints(), F5LAget_dimensions(), F5LAsave_dimensions(), F5printf, FIBER_FIELD_DATASPACE_DIMENSIONS_ATTRIBUTE, FIBER_HDF5_RANK_ATTRIBUTE, FIBER_MAX_RANK, and F5Path::Topology_hid.

Referenced by F5Fcreate(), F5FSwrite_fraction(), F5Fwrite_flexible(), F5Fwrite_fraction(), F5Fwrite_fraction_external(), F5Fwrite_fractionS(), F5Fwrite_linear(), F5FwriteIMAGE(), F5Fwrites(), and F5FwriteX().

if()

if

(

return

f

)

Definition at line 681 of file F5LT.c.

{
ChartDomain_IDs* (*coord_creator)() = (ChartDomain_IDs*(*)() )udata;
        assert( coord_creator );
        return coord_creator();
}

References F5printf_indent, and TopologyName().

Variable Documentation

else

else

Initial value:

{
                F5printf(1, "F5 Creating a Topological %d-dimensional Skeleton of Index Depth=%d, "
                         "Info - No chart used for %s, seems to be a relative representation in [%s]"
                         , Dimensionality, IndexDepth
                         , TopologyName, coordinate_system)

Definition at line 673 of file F5LT.c.