F5coordinates.c

Go to the documentation of this file.

#include "F5Bchart.h"
#include "F5defs.h"
#include "F5coordinates.h"
 
hid_t   F5T_IntegerFraction32()
{
static  hid_t IF_id = -1;
 
        if (IF_id<0)
        {
                IF_id = H5Tcreate(H5T_COMPOUND, sizeof(IntegerFraction32) );
                H5Tinsert(IF_id, "num"  ,                   0, H5T_NATIVE_INT32);
                H5Tinsert(IF_id, "denom", sizeof(F5_uint32_t), H5T_NATIVE_INT32);
        }
        return IF_id;
}
 
hid_t   F5T_IntegerFraction64()
{
static  hid_t IF_id = -1;
 
        if (IF_id<0)
        {
                IF_id = H5Tcreate(H5T_COMPOUND, sizeof(IntegerFraction64) );
                H5Tinsert(IF_id, "num"  ,                   0, H5T_NATIVE_INT64);
                H5Tinsert(IF_id, "denom", sizeof(F5_uint64_t), H5T_NATIVE_INT64);
        }
        return IF_id;
}
 
#define STANDARD_CHART(what)                            \
ChartDomain_IDs*F5B_standard_##what()                   \
{                                                       \
static ChartDomain_IDs*myChart = 0;                     \
                                                        \
        if (!myChart)                                   \
                myChart = F5B_new_global_##what();      \
                                                        \
        return myChart;                                 \
}
 
 
/*
Note: two main functions co-exist for the creation of coordinate systems for Representations:
 
        F5B_new_global_domain:  Only positional type
        F5B_new_global_chart:   Also tangential/tensor types
*/
 
 
 
/*
TODO: Add standard chart and cleanup function for each coordinate system
 */
 
/************* 4D charts ************************/
 
ChartDomain_IDs*F5B_new_global_cartesian_chart4D()
{
const char*coordnames[] = { "t", "x", "y", "z" };
        return F5B_new_global_float_chart(coordnames, 4, "Cartesian 4D", F5_FORTRAN_ORDER);
}
 
ChartDomain_IDs*F5B_new_global_polar_chart4D()
{
const char*coordnames[] = { "t", "r", "theta", "phi" };
        return F5B_new_global_float_chart(coordnames, 4, "Polar 4D", F5_FORTRAN_ORDER);
}
 
/************* 3D charts ************************/
 
ChartDomain_IDs*F5B_new_integer_regular_domain3D()
{
const char*coordnames[] = {"x", "y", "z" };
        return F5B_new_global_int_domain(coordnames, 3, "Integer3D", F5_FORTRAN_ORDER, 3);
}
 
ChartDomain_IDs*F5B_new_rational_regular_domain3D()
{
const char*coordnames[] = {"x", "y", "z" };
        return F5B_new_global_fractional_domain(coordnames, 3, "Rational3D", F5_FORTRAN_ORDER, 3);
}
 
ChartDomain_IDs*F5B_new_global_cartesian_chart3D()
{
const char*coordnames[] = {"x", "y", "z" };
        return F5B_new_global_float_chart(coordnames, 3, FIBER_HDF5_CARTESIAN_CHART_DOMAIN, F5_FORTRAN_ORDER);
}
 
/********************************************************/
 
ChartDomain_IDs*F5B_new_global_triangle_cells_ijk()
{
const char*coordnames[] = {"i", "j", "k" };
        return F5B_new_global_int_domain(coordnames, 3, F5_TRIANGLE_CELL_DOMAIN, F5_FORTRAN_ORDER, 2);
}
STANDARD_CHART(triangle_cells_ijk)
 
ChartDomain_IDs*F5B_new_global_triangle_cells()
{
const char*coordnames[] = {"ii","ij","jj"};
        return F5B_new_global_int_domain(coordnames, 3, F5_TRIANGLE_CELL_DOMAIN, F5_FORTRAN_ORDER, 2);
}
 
STANDARD_CHART(triangle_cells)
 
ChartDomain_IDs*F5B_new_global_edge_cells()
{
const char*coordnames[] = {"i", "j" };
        return F5B_new_global_int_domain(coordnames, 2, F5_EDGE_CELL_DOMAIN, F5_FORTRAN_ORDER, 1);
}
 
STANDARD_CHART(edge_cells)
 
 
ChartDomain_IDs*F5B_new_global_quad_cells()
{
/*const char*coordnames[4] = {"i00", "i01", "i11", "i10" };*/
const char*coordnames[4] = { "ii","ij","ji","jj"};
        return F5B_new_global_int_domain(coordnames, 4, F5_QUAD_CELL_DOMAIN, F5_FORTRAN_ORDER, 2);
}
STANDARD_CHART(quad_cells)
 
ChartDomain_IDs*old_F5B_new_global_tetrahedral_cells()
{
const char*coordnames[4] = {"i0", "i01", "i10", "i11" };
        return F5B_new_global_int_domain(coordnames, 4, F5_TETRAHEDRON_CELL_DOMAIN, F5_FORTRAN_ORDER, 3);
}
 
 
ChartDomain_IDs*F5B_new_global_tetrahedral_cells()
{
const char*coordnames[4] = { "iii", "iij", "ijj", "jjj" };
        return F5B_new_global_int_domain(coordnames, 4, F5_TETRAHEDRON_CELL_DOMAIN, F5_FORTRAN_ORDER, 3);
}
STANDARD_CHART(tetrahedral_cells)
 
 
ChartDomain_IDs*F5B_new_global_hexahedral_cells()
{
static  const char*coordnames[8] = { "iii","iij","iji","ijj","jii","jij","jji","jjj" };
        return F5B_new_global_int_domain(coordnames, 8, F5_HEXAHEDRON_CELL_DOMAIN, F5_FORTRAN_ORDER, 3);
}
STANDARD_CHART(hexahedral_cells)
 
extern
ChartDomain_IDs*F5B_new_global_fragment_neighbour_domain(const char*coordnames[],
                                                         int dimension,
                                                         const char*name,
                                                         const int perm_vector[],
                                                         int cell_dimensionality);
ChartDomain_IDs*F5B_new_global_fragname()
{
static const char*coordnames[] = { "a", "b", "c", "d", "e", "f", "g", "h" };
        return F5B_new_global_fragment_neighbour_domain(coordnames, 8, F5_FRAGMENT_NEIGHBOUR_CELL_DOMAIN, F5_FORTRAN_ORDER, 3);
}
 
/********************************************************/
 
static ChartDomain_IDs*standard_cartesian_chart3D = 0;
 
ChartDomain_IDs*F5B_standard_cartesian_chart3D()
{
        if (!standard_cartesian_chart3D)
                standard_cartesian_chart3D = F5B_new_global_cartesian_chart3D();
 
        return standard_cartesian_chart3D;
}
 
 
ChartDomain_IDs*F5B_new_global_polar_chart3D()
{
const char*coordnames[] = {"r", "theta", "phi" };
        return F5B_new_global_float_chart(coordnames, 3, "Polar 3D", F5_FORTRAN_ORDER);
}
 
ChartDomain_IDs*F5B_new_global_cylindrical_chart3D()
{
const char*coordnames[] = {"r", "phi", "z"};
        return F5B_new_global_float_chart(coordnames, 3, "Cylindrical 3D", F5_FORTRAN_ORDER);
}
 
 
ChartDomain_IDs*F5B_new_global_texture_chart3D()
{
const char*coordnames[] = {"u", "v", "w" };
        return F5B_new_global_float_chart(coordnames, 3, "Texture 3D", F5_FORTRAN_ORDER);
}
 
ChartDomain_IDs*F5B_new_global_t_polar2_chart3D()
{
const char*coordnames[] = {"t", "r", "phi" };
        return F5B_new_global_float_chart(coordnames, 3, "Time Polar 2D", F5_FORTRAN_ORDER);
}
 
/************* 2D charts ************************/
 
ChartDomain_IDs*F5B_new_global_cartesian_chart2D()
{
const char*coordnames[] = {"x", "y" };
        return F5B_new_global_float_chart(coordnames, 2, "Cartesian 2D", F5_FORTRAN_ORDER);
}
 
ChartDomain_IDs*F5B_new_global_polar_chart2D()
{
const char*coordnames[] = {"r", "phi" };
        return F5B_new_global_float_chart(coordnames, 2, "Polar 2D", F5_FORTRAN_ORDER);
}
 
ChartDomain_IDs*F5B_new_global_spherical_chart2D()
{
const char*coordnames[] = {"phi", "theta" };
        return F5B_new_global_float_chart(coordnames, 2, "Spherical 2D", F5_FORTRAN_ORDER);
}
 
ChartDomain_IDs*F5B_new_global_axial_chart2D()
{
const char*coordnames[] = {"r", "theta" };
        return F5B_new_global_float_chart(coordnames, 2, "Axial 2D", F5_FORTRAN_ORDER);
}
 
/************* 1D charts ************************/
ChartDomain_IDs* F5B_new_global_cartesian_chart1D()
{
const char*coordnames[] = {"x"};
        return F5B_new_global_float_chart(coordnames, 1, "Cartesian 1D", F5_FORTRAN_ORDER);
}
 
ChartDomain_IDs* F5B_new_global_spherical_chart1D()
{
const char*coordnames[] = {"phi"};
        return F5B_new_global_float_chart(coordnames, 1, "Spherical 1D", F5_FORTRAN_ORDER);
}
 
ChartDomain_IDs* F5B_new_global_time_chart1D()
{
const char*coordnames[] = {"t"};
        return F5B_new_global_float_chart(coordnames, 1, "Time 1D", F5_FORTRAN_ORDER);
}
 
 
ChartDomain_IDs*F5B_new_global_texture_chart2D()
{
const char*coordnames[] = {"u", "v" };
        return F5B_new_global_float_chart(coordnames, 2, "Texture 2D",  F5_FORTRAN_ORDER);
}
 
static ChartDomain_IDs*standard_texture_chart2D = 0;
 
ChartDomain_IDs*F5B_standard_texture_chart2D()
{
        if (!standard_texture_chart2D)
                standard_texture_chart2D = F5B_new_global_texture_chart2D();
 
        return standard_texture_chart2D;
}
 
/******************** Color Charts **************/
ChartDomain_IDs*F5B_new_global_color_chart_rgba_real()
{
const char*coordnames[] = {"r", "g", "b", "a" };
        return F5B_new_global_float_chart(coordnames, 4, "RGBA", F5_FORTRAN_ORDER);
}
 
ChartDomain_IDs*F5B_new_global_color_chart_rgba_int()
{
const char*coordnames[] = {"r", "g", "b", "a" };
        return F5B_new_global_int_chart(coordnames, 4, "RGBA", F5_FORTRAN_ORDER);
}
 
 
/********************  CLEANUP ******************/
 
void F5Bcoordinate_cleanup()
{
        F5B_delete_global_chart(standard_cartesian_chart3D);
        F5B_delete_global_chart(standard_texture_chart2D);
}