F5T.c File Reference

#include "F5T.h"
#include "F5A.h"
#include "F5private.h"

Go to the source code of this file.

Functions
hid_t	F5Tcommit (hid_t save_location, const char *name, hid_t TransientTypeID)
hid_t	F5Tappend (hid_t loc_id, const char *name, hid_t type_id)
void	F5Tclose (hid_t type_id)
const char *	F5Tget_precision_name (int TypeSize)
hid_t	F5Tsave_tensor (hid_t Domain_hid, const char *domain_name, const char *type_name, hid_t transient_type, int rank, int indices[], int grade)
hid_t	F5Tmake_precise_tensor (F5Path *fpath, const char *type_name, int TypeBytePrecision, hid_t transient_type, int rank, int indices[], int grade)

Function Documentation

F5Tappend()

hid_t F5Tappend	(	hid_t	loc_id,
		const char *	name,
		hid_t	type_id )

Open some type if it exists, otherwise create one.

Definition at line 25 of file F5T.c.

{
hid_t   FileType_id;
        H5E_BEGIN_TRY
        FileType_id = H5Topen2(loc_id, name, H5P_DEFAULT);
        H5E_END_TRY
 
 
        if (FileType_id<0)
        {
                FileType_id = F5Tcommit(loc_id, name, type_id);
        }
        else
        {
                F5printf(1,"F5Tappend(,name=%s,)", name);
 
                if (H5Tequal(FileType_id, type_id) != 1)
                        return 0;
        }
        return FileType_id;
}

References F5printf, F5Tcommit(), and name.

Referenced by F5file_type().

F5Tclose()

void F5Tclose

(

hid_t

type_id

)

Close a HDF5 Type ID if it's a valid identifier, otherwise silently ignore.

Definition at line 49 of file F5T.c.

{
        if (type_id>0)
                H5Tclose(type_id);
}

References H5Tclose().

Referenced by F5B_delete_tangential_types().

F5Tcommit()

hid_t F5Tcommit	(	hid_t	save_location,
		const char *	name,
		hid_t	TransientTypeID )

Save a type to a file and return a copy of the type ID.

Note

The returned HDF5 id must be released via H5Tclose().

Parameters

save_location	Where to save the type, must be some group ID
name	The name of the type
TransientTypeID	The intrinsic type ID

Returns

A type Id which is bound to the file to which save_location belongs

Definition at line 6 of file F5T.c.

{
hid_t   hid = H5Tcopy( TransientTypeID );
herr_t  herr;
 
/*      printf("F5Tcommit(): %s\n", name); */
        H5E_BEGIN_TRY
        herr = H5Tcommit2 (save_location, name, hid, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
        H5E_END_TRY
 
        if (herr<0)
        {
/*              printf("F5Tcommit(): %s FAILED\n", name); */
                H5Tclose( hid );
                hid = H5Topen2(save_location, name, H5P_DEFAULT);
        }
        return hid;
}

References H5Tclose(), and name.

Referenced by F5Tappend(), and F5Tsave_tensor().

F5Tget_precision_name()

const char * F5Tget_precision_name

(

int

TypeSize

)

Retrieve a string corresponding to the type size, used for coordinate sub-precision.

Definition at line 56 of file F5T.c.

{
        switch(TypeSize)
        {
        case 1:         return FIBER_HDF5_CHART_ONEBYTE_PRECISION;
        case 2:         return FIBER_HDF5_CHART_HALF_PRECISION  ;
        case 4:         return FIBER_HDF5_CHART_SINGLE_PRECISION;
        case 8:         return FIBER_HDF5_CHART_DOUBLE_PRECISION;
        case 16:        return FIBER_HDF5_CHART_EXTENDED_PRECISION;
 
        default:        return 0;
        }
}

References FIBER_HDF5_CHART_DOUBLE_PRECISION, FIBER_HDF5_CHART_EXTENDED_PRECISION, FIBER_HDF5_CHART_HALF_PRECISION, FIBER_HDF5_CHART_ONEBYTE_PRECISION, and FIBER_HDF5_CHART_SINGLE_PRECISION.

Referenced by F5Tmake_precise_tensor().

F5Tmake_precise_tensor()

hid_t F5Tmake_precise_tensor	(	F5Path *	fpath,
		const char *	type_name,
		int	TypeBytePrecision,
		hid_t	transient_type,
		int	rank,
		int	indices[],
		int	grade )

Definition at line 115 of file F5T.c.

{
const char*precision_name = F5Tget_precision_name( TypeBytePrecision );
        if (!precision_name)
        {
                puts("F5Tmake_precise_tensor(): Invalid Precision byte size specified!");
                return -1;
        }
 
        if (transient_type <= 0)
        {
                puts("F5Tmake_precise_tensor(): Invalid transient type given.");
                return -1;
        }
 
        {
        hid_t   type_id;
        hid_t   Domain_hid = H5Gopen2(fpath->GlobalChart_hid, fpath->myChart->domain_name, H5P_DEFAULT);
        hid_t   Precision_hid;
                if (Domain_hid<0)
                {
                        puts("F5Tmake_precise_tensor(): Opening DOMAIN NAME failed!");
                        return -1;
                }
 
                Precision_hid = H5Gopen2(Domain_hid,precision_name, H5P_DEFAULT);
                if (Precision_hid<0)
                {
                        puts("F5Tmake_tensor(): Opening PRECISION DOMAIN failed!");
                        return -1;
                }
 
                F5printf(20, "   F5Tmake_tensor(,%s, %x, rank=%d,, grade=%d)", type_name, transient_type, rank, grade);
 
                type_id = F5Tsave_tensor(Precision_hid, fpath->myChart->domain_name, type_name, transient_type, rank, indices, grade);
 
                H5Gclose(Precision_hid);
                H5Gclose(Domain_hid);
 
                return type_id;
        }
}

References ChartDomain_IDs::domain_name, F5printf, F5Tget_precision_name(), F5Tsave_tensor(), F5Path::GlobalChart_hid, H5Gclose, and F5Path::myChart.

Referenced by F5file_type().

F5Tsave_tensor()

hid_t F5Tsave_tensor	(	hid_t	Domain_hid,
		const char *	domain_name,
		const char *	type_name,
		hid_t	transient_type,
		int	rank,
		int	indices[],
		int	grade )

@TODO Specify symmetricity of the indices. @TODO Specify non-tensorial transformation rules, e.g. Christoffel symbols

Definition at line 74 of file F5T.c.

{
hid_t   type_id = F5Tcommit(Domain_hid, type_name, transient_type);
 
        /*
          THIS would be nice!
        H5Glink(type_id, H5G_LINK_SOFT,  FIBER_HDF5_CHART_DOMAIN, "/Charts/" domain_name);
        */
 
        F5Asave_string(type_id, FIBER_HDF5_CHART_DOMAIN, domain_name);
 
        if (rank>0 && indices)
        {
                if (H5Aexists(type_id, FIBER_HDF5_CHART_TENSOR_TYPE) == 0)
                {
                hsize_t Tdims = rank;
                hid_t   index_ds = H5Screate_simple(1, &Tdims, &Tdims ) ;
                hid_t   Rank = H5Acreate2(type_id, FIBER_HDF5_CHART_TENSOR_TYPE, H5T_NATIVE_INT, index_ds, H5P_DEFAULT, H5P_DEFAULT);
                        H5Awrite(Rank, H5T_NATIVE_INT, indices );
                        H5Aclose(Rank);
                        H5Sclose(index_ds);
                }
                else
                {
                        F5printf(20, "   F5Tsave_tensor(domain_name=%s, type_name=%s, %x, rank=%d,, grade=%d): Attribute \"%s\" already exists!", domain_name, type_name, transient_type, rank, grade, FIBER_HDF5_CHART_TENSOR_TYPE);
                }
                
        }
 
        if (rank>0 && grade>0)
        {
                F5Asave_ints( type_id, FIBER_HDF5_GRADE_ATTRIBUTE, &grade, 1);
        }
 
        return type_id;
}

References F5Asave_ints(), F5Asave_string(), F5printf, F5Tcommit(), FIBER_HDF5_CHART_DOMAIN, FIBER_HDF5_CHART_TENSOR_TYPE, and FIBER_HDF5_GRADE_ATTRIBUTE.

Referenced by F5B_save_global_chart(), and F5Tmake_precise_tensor().