F5T.c

Go to the documentation of this file.

#include "F5T.h"
#include "F5A.h"
#include "F5private.h"
 
 
hid_t   F5Tcommit(hid_t save_location, const char*name, hid_t TransientTypeID)
{
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;
}
 
hid_t   F5Tappend(hid_t loc_id, const char *name, hid_t type_id)
{
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;
}
 
 
 
void    F5Tclose(hid_t type_id)
{
        if (type_id>0)
                H5Tclose(type_id);
}
 
 
const char*F5Tget_precision_name(int TypeSize)
{
        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;
        }
}
 
/**
 @TODO Specify symmetricity of the indices.
 @TODO Specify non-tensorial transformation rules, e.g. Christoffel symbols
*/
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   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;
}
 
 
 
hid_t   F5Tmake_precise_tensor(F5Path*fpath, const char*type_name, int TypeBytePrecision,
                               hid_t transient_type, int rank, int indices[], int grade)
{
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;
        }
}