F5B.c

Go to the documentation of this file.

#include "F5B.h"
#include "F5X.h"
#include "F5F.h"
#include "F5defs.h"
#include "F5coordinates.h"
#include "F5private.h"
 
 
#include <stdlib.h>
 
#include <string.h>
#include <assert.h>
 
ChartDomain_IDs*F5FB_read_global_chart(F5Path*f)
{
        return  F5B_read_global_chart(f->GlobalChart_hid, f->myChart->domain_name);
}
 
 
struct  mkDomain
{
        const char*name;
        hid_t domain_type;
};
 
 
#if 0
/**\ingroup F5L
  Create a Topology of the given name with a representation in
  another Topology.
  The types in the chart domain (F5Path->myChart) will only contain a copy
  of the domain_type in the Point_hid_t, all other types like the
  Vector, Covector etc. types are invalid (-1).
 
  @note A new ChartDomain object is created, even if there might
        be one already somewhere in memory elsewhere.
  @param domainname          A description for the refinement type.
  @param domain_type         An hdf5 type that describes the pointwise refinement
  @param representation_name The name of the representation reference, which
                             will be of the type domainname.
 */
F5Path* F5LTcreateT(hid_t File_id, const double*time,
                    const char*gridname,
                    const char*domainname,
                    hid_t      domain_type,
                    const char*representation_name,
                    const char*TopologyName,
                    int IndexDepth,
                    int Dimensionality,
                    const hsize_t*refinement)
{
struct  mkDomain D;
        D.name = domainname;
        D.domain_type = domain_type;
 
        return F5LTcreateV(File_id, time, gridname,
                           0,
                           F5B_new_global_domain,
                           &D,
                           representation_name,
                           TopologyName,
                           IndexDepth,
                           Dimensionality,
                           refinement);
}
#endif
 
void F5Bincrement_domain_refs(ChartDomain_IDs*ChartDomain)
{
        if (ChartDomain)
                ChartDomain->refs++;
}
 
F5Path*F5Bopen_coordinate_domain(hid_t File_id,
                                 ChartDomain_IDs*ChartDomain)
{
F5Path  *f;
        if (File_id<1)
        {
                F5printf(10, "*F5Bopen_coordinate_domain(): Invalid file ID given!");
                return NULL;
        }
 
        f = newF5Path(File_id);
 
        F5printf(30, "F5Bopen_coordinate_domain()");
 
        H5E_BEGIN_TRY
        f->ContentsGroup_hid = H5Gopen2(File_id, FIBER_CONTENTS,H5P_DEFAULT);
 
        if (f->ContentsGroup_hid<=0)
                F5printf(10, "File contains no "FIBER_CONTENTS" contents group (old F5 file?).");
 
        /* NOTE: This should actually retrieve the enum types plus
           field versioning information from the file!
         */
        if (f->ContentsGroup_hid>=0)
        {
                f->field_enum_type_hid = H5Topen2(f->ContentsGroup_hid, FIBER_ARRAY_TYPE_ENUM, H5P_DEFAULT);
                if (f->field_enum_type_hid < 0)
                        F5printf(10, "File Contents group contains no field enumeration type (old F5 file?)");
        }
        H5E_END_TRY
 
        f->Slice_hid            = 0;
        f->Grid_hid             = 0;
        f->Charts_hid           = 0;
        f->Chart_hid            = 0;
        f->Topology_hid         = 0;
        f->Representation_hid   = 0;
 
        /* preparing chart */
 
        f->GlobalChart_hid = F5Gtry_to_open(File_id, FIBER_HDF5_GLOBAL_CHARTS) ;
        if (f->GlobalChart_hid<0)
        {
                F5printf(10, "File contains no "FIBER_HDF5_GLOBAL_CHARTS" entry (invalid F5 file?)");
                deleteF5Path(f);
                F5printf(30, "~F5Bopen_coordinate_domain() failed");
                return NULL;
        }
 
        assert(ChartDomain);
        {
                f->myChart = ChartDomain;
                f->myChart->refs++;
        }
 
        /* Creates an entry in the "global charts" group for all cartesian charts,
           this is the place to store the named types of cartesian charts.
 
        NOTE: This coordinate/chart reference stuff is planned to be replaced by shared
              types. This revision will work without links.
        */
        if (F5Gexist(f->GlobalChart_hid, f->myChart->domain_name) )
        {
                f->FileIDs = F5FB_read_global_chart(f);
                F5printf(30, "~F5Bopen_coordinate_domain() is ok");
                return f;
        }
#if     1
        F5printf(10, "F5Bopen_coordinate_domain() failed, but still trying to make sense out of this file.");
        f->FileIDs = 0;
        return f;
#else
        H5Gclose( f->GlobalChart_hid );
        deleteF5Path(f);
 
        F5printf(30, "~F5Bopen_coordinate_domain() failed");
        return NULL;
#endif
}
 
 
 
F5Path*F5Binitialize_path(hid_t File_id)
{
/*const char*cs = coordinate_system ? coordinate_system : FIBER_HDF5_DEFAULT_CHART; */
const char*cs = FIBER_HDF5_DEFAULT_CHART;
ChartDomain_IDs*ChartDomain = (strcmp(cs, FIBER_HDF5_DEFAULT_CHART)==0) ? F5B_standard_cartesian_chart3D() : 0;
 
F5Path*f = F5Bopen_coordinate_domain( File_id, ChartDomain);
        if (!f) /* hack for backward compatibility. */
        {
        ChartDomain_IDs*oldChartDomain = F5B_standard_cartesian_chart3D();
                oldChartDomain->domain_name = FIBER_HDF5_DEPRECATED_OLD_CARTESIAN_DOMAIN;
 
                F5printf(10, "*F5Binitialize_path(): Could not initialize with "
                         FIBER_HDF5_DEFAULT_CHART ", now trying " FIBER_HDF5_DEPRECATED_OLD_CARTESIAN_DOMAIN "...");
 
                f = F5Bopen_coordinate_domain(File_id, ChartDomain);
                oldChartDomain->domain_name = FIBER_HDF5_CARTESIAN_CHART_DOMAIN;
 
                if (!f)
                {
                        F5printf(10, "*F5Binitialize_path(): Could not initialize with "
                                 FIBER_HDF5_DEPRECATED_OLD_CARTESIAN_DOMAIN " either... not good.");
 
                        f = F5Bopen_coordinate_domain( File_id, F5B_standard_cartesian_chart3D() );
                }
        }
        return f;
}
 
 
void    F5close(F5Path*f)
{
        if (!f) return;
 
        F5B_delete_global_chart(f->myChart);
        F5B_delete_global_chart(f->FileIDs);
        if (f->field_enum_type_hid>0)   H5Tclose(f->field_enum_type_hid);
        F5printf(100, "F5close(): Closing representation... %d", f->Representation_hid );
        if (f->Representation_hid)      H5Gclose(f->Representation_hid);
        F5printf(100, "F5close(): Closing global chart...");
        if (f->GlobalChart_hid)         H5Gclose(f->GlobalChart_hid);
        F5printf(100, "F5close(): Closing chart group...");
        if (f->Charts_hid)              H5Gclose(f->Charts_hid);
        F5printf(100, "F5close(): Closing chart...");
        if (f->Chart_hid)               H5Gclose(f->Chart_hid);
        if (f->Topology_hid)            H5Gclose(f->Topology_hid);
        F5printf(100, "F5close(): Closing grid...");
        if (f->Grid_hid)                H5Gclose(f->Grid_hid);
        F5printf(100, "F5close(): Closing slice...");
        if (f->Slice_hid)               H5Gclose(f->Slice_hid);
        F5printf(100, "F5close(): Closing ContentsGroup...");
        if (f->ContentsGroup_hid>0)     H5Gclose(f->ContentsGroup_hid);
 
        if (f->Field_hid>0)
        {
                F5printf(40, "F5close(): field is %s.", F5Fis_group(f)?"a group":"a dataset");
 
                if (F5Fis_group(f) )
                        H5Gclose( f->Field_hid );
                else
                        H5Dclose( f->Field_hid );
        }
 
        deleteF5Path(f);
}
 
void F5closep(F5Path**f)
{
        if (!f)  return;
        if (!*f) return;
        F5close(*f);
        *f = NULL;
}
 
void F5Bswap_dims (int dimensions, const hsize_t* in, hsize_t* out)
{
int i;
 
    for (i = 0; i < dimensions; i++) {
        out[i] = in[dimensions - 1 - i];
    }
}