F5P.c

Go to the documentation of this file.

#include "F5P.h"
#include "F5X.h"
#include "F5uniform.h"
/* #if defined(__APPLE__) && defined(__MACH__) && defined(__POWERPC__) */
#include <stdlib.h>
/* #else */
/* #include <malloc.h> */
/* #endif */
 
int F5P_is_regular3D(F5Path*grid, const char*coordinate_system)
{
hsize_t dims[FIBER_MAX_RANK];
int     rank = F5get_extension(grid, dims);
 
        if (rank==3)
                return 1;
 
        return 0;
}
 
int F5P_is_uniform(F5Path*grid, const char*coordinate_system)
{
hid_t   Top_id, Rep_id;
int     ok;
 
        if (!grid->Grid_hid)
                return 0;
 
        if (grid->Topology_hid)
                Top_id = grid->Topology_hid;
        else
                Top_id = H5Gopen2(grid->Grid_hid, FIBER_HDF5_POINTS, H5P_DEFAULT);
 
        Rep_id = H5Gopen2( Top_id, coordinate_system, H5P_DEFAULT );
 
        ok = F5Lis_linear(Rep_id, FIBER_HDF5_POSITIONS_STRING);
 
        H5Gclose( Rep_id );
 
        if (!grid->Topology_hid)
                H5Gclose( Top_id );
 
        return ok;
}
 
 
/**
   A grid is `Z-stacked' if its coordinates in a certain chart
   are uniform in all coordinates but the last one, which is
   given as a one-dimensional dataset.
  */
int F5P_is_Zstacked(F5Path*grid, const char*coordinate_system)
{
hid_t   Top_id, Rep_id;
int     ok;
hsize_t dims[FIBER_MAX_RANK];
int     rank;
 
        if (!grid->Grid_hid)
                return 0;
 
        rank = F5get_extension(grid, dims);
 
        if (rank!=3) return 0;
 
        if (grid->Topology_hid)
                Top_id = grid->Topology_hid;
        else
                Top_id = H5Gopen2(grid->Grid_hid, FIBER_HDF5_POINTS, H5P_DEFAULT);
 
        Rep_id = H5Gopen2( Top_id, coordinate_system, H5P_DEFAULT );
 
        ok = F5Lis_linear(Rep_id, FIBER_HDF5_POSITIONS_STRING);
        if (ok)
        {
        char*name;
        hid_t   pos_id = H5Gopen2(Rep_id, FIBER_HDF5_POSITIONS_STRING, H5P_DEFAULT); 
        hid_t   type_id = F5Lget_type(pos_id, 1, 0);
        int     n = H5Tget_nmembers(type_id);
 
                /* Look for the last member */
                name = H5Tget_member_name(type_id, n-1);
                if (!F5Dexist(pos_id, name))
                        ok = 0;
                free(name);
 
                H5Tclose(type_id);
                H5Gclose(pos_id);
        }
        H5Gclose( Rep_id );
 
        if (!grid->Topology_hid)
                H5Gclose( Top_id );
 
        return ok;
}
 
 
 
/*
Rectilinear: The positions contains an base and delta attribute,
like in the uniform case, but also exactly one dataset for each
coordinate. The base and delta attributes are not used here; they
are just to distinguish this case from the curvilinear case, where
coordinates at each point are given via separate coordinate
component fields, e.g. three scalar fields {x,y,z}.
Another criteria is that here the three fields {x,y,z} are
one-dimensional, while for curvilinear separated coordinates they
might be n-dimensional datasets. However, we don't want to
insist of saving n-dimensional datasets as n-dimensional HDF5
datasets (maybe this should be changed?).
*/
int F5P_is_rectilinear(F5Path*grid, const char*coordinate_system)
{
hid_t   Top_id, Rep_id;
int     ok;
 
        if (!grid->Grid_hid)
                return 0;
 
        if (grid->Topology_hid)
                Top_id = grid->Topology_hid;
        else
                Top_id = H5Gopen2(grid->Grid_hid, FIBER_HDF5_POINTS, H5P_DEFAULT);
 
        Rep_id = H5Gopen2( Top_id, coordinate_system, H5P_DEFAULT );
 
        ok = F5Lis_linear(Rep_id, FIBER_HDF5_POSITIONS_STRING);
        if (ok)
        {
        int     i, n;
        char*name;
        hid_t   pos_id = H5Gopen2(Rep_id, FIBER_HDF5_POSITIONS_STRING, H5P_DEFAULT); 
        hid_t   type_id = F5Lget_type(pos_id, 1, 0);
                n = H5Tget_nmembers(type_id);
                for(i=0; i<n; i++)
                {
                        name = H5Tget_member_name(type_id, i);
                        if (!F5Dexist(pos_id, name))
                        {
                                ok = 0;
                                free(name);
                                break;
                        }
                        free(name);
                }
                H5Tclose(type_id);
                H5Gclose(pos_id);
        }
        H5Gclose( Rep_id );
 
        if (!grid->Topology_hid)
                H5Gclose( Top_id );
 
        return ok;
}
 
 
/** The grid is curvilinear if the positions is a dataset of rank n
    with n>1 or it is a group with n scalar field that are named like
    the coordinate entries.
  */
int F5P_is_curvilinear(F5Path*grid, const char*coordinate_system)
{
hid_t   Top_id, Rep_id, type_id;
hsize_t dims[FIBER_MAX_RANK];
int     rank;
int     ok = 0;
 
        if (!grid->Grid_hid)
                return 0;
 
        rank = F5get_extension(grid, dims);
 
        if (rank<2) return 0;
 
        if (grid->Topology_hid)
                Top_id = grid->Topology_hid;
        else
                Top_id = H5Gopen2(grid->Grid_hid, FIBER_HDF5_POINTS, H5P_DEFAULT);
 
        Rep_id = H5Gopen2( Top_id, coordinate_system, H5P_DEFAULT );
 
        ok = F5Dexist(Rep_id, FIBER_HDF5_POSITIONS_STRING);
 
        if (!ok)
        {
                ok = F5Gexist(Rep_id, FIBER_HDF5_POSITIONS_STRING);
                if (ok)
                {
                int     i, n;
                char*name;
                hid_t   pos_id = H5Gopen2(Rep_id, FIBER_HDF5_POSITIONS_STRING, H5P_DEFAULT); 
                        type_id = F5Lget_type(pos_id, 1, 0);
                        n = H5Tget_nmembers(type_id);
                        for(i=0; i<n; i++)
                        {
                                name = H5Tget_member_name(type_id, i);
                                if (!F5Dexist(pos_id, name))
                                {
                                        ok = 0;
                                        break;
                                }
                                free(name);
                        }
 
                        H5Tclose(type_id);
                        H5Gclose(pos_id);
                }
        }
 
        H5Gclose( Rep_id );
 
        if (!grid->Topology_hid)
                H5Gclose( Top_id );
 
        return ok;
}
 
 
int F5P_is_triangular_surface(F5Path*grid, const char*coordinate_system)
{
hid_t   Cell_Top_id, Cell_Vertices_id, Vertices_id, Vertices_Type_id;
int     ok;
 
        if (!F5Gexist(grid->Grid_hid, FIBER_HDF5_CELLS) )
                if( !F5Gexist(grid->Grid_hid, FIBER_HDF5_FACES) )
                        return 0;
                else
                        Cell_Top_id = H5Gopen2(grid->Grid_hid, FIBER_HDF5_FACES, H5P_DEFAULT);
        else
                Cell_Top_id = H5Gopen2(grid->Grid_hid, FIBER_HDF5_CELLS, H5P_DEFAULT);
 
        if (!F5Gexist(Cell_Top_id, FIBER_HDF5_POINTS) )
        {
                H5Gclose( Cell_Top_id );
                return 0;
        }
 
        Cell_Vertices_id = F5Gtry_to_open(Cell_Top_id, FIBER_HDF5_POINTS);
        if (Cell_Vertices_id<1)
        {
                H5Gclose( Cell_Top_id );
                return 0;
        }
 
        Vertices_id = F5Dtry_to_open(Cell_Vertices_id, FIBER_HDF5_POSITIONS_STRING);
        if (Vertices_id<1)
        {
                H5Gclose( Cell_Top_id );
                H5Gclose( Cell_Vertices_id );
                return 0;
        }
 
        Vertices_Type_id = H5Dget_type(Vertices_id);
 
        if (H5Tequal(Vertices_Type_id, F5T_TRIANGLE32)>0 )
        {
                ok = 1;
        }
        else if (H5Tequal(Vertices_Type_id, F5T_TRIANGLE64)>0 )
        {
                ok = 1;
        }
        else
                ok = 0;
 
        H5Tclose( Vertices_Type_id );
        H5Gclose( Cell_Top_id );
        H5Gclose( Cell_Vertices_id );
        H5Dclose( Vertices_id );
        return ok;
}
 
/** No Connectivity information is given, and the dimensionality of
    the Positions is 1.
 */
int F5P_is_particle_system(F5Path*grid, const char*coordinate_system)
{
hsize_t dims[FIBER_MAX_RANK];
int     rank = F5get_extension(grid, dims);
 
        if (rank==1)
        {
                if (!grid->Grid_hid)
                        return 0;
 
                if (!F5Gexist(grid->Grid_hid, FIBER_HDF5_CELLS) )
                        return 1;
 
                return 0;
        }
 
        return 0;
}
 
 
/** Check if the grid has positions in the specified coordinate system */
int F5P_has_vertices(F5Path*grid, const char*coordinate_system)
{
hid_t   Top_id, Rep_id;
int     ok;
 
        if (!grid->Grid_hid)
                return 0;
 
        if (grid->Topology_hid)
                Top_id = grid->Topology_hid;
        else
                Top_id = H5Gopen2(grid->Grid_hid, FIBER_HDF5_POINTS, H5P_DEFAULT);
 
        Rep_id = H5Gopen2( Top_id, coordinate_system, H5P_DEFAULT );
 
        ok = F5Dexist(Rep_id, FIBER_HDF5_POSITIONS_STRING) ||
             F5Gexist(Rep_id, FIBER_HDF5_POSITIONS_STRING);
 
        H5Gclose( Rep_id );
 
        if (!grid->Topology_hid)
                H5Gclose( Top_id );
 
        return ok;
}
 
 
/** Check if the grid has some fields in the specified coordinate system */
int F5P_has_vertex_fields(F5Path*grid, const char*coordinate_system);
/** Check if the grid has some data fields defined on its cells */
int F5P_has_cell_fields(F5Path*grid, const char*coordinate_system);