F5Bchart.h

/*
 $Id: F5Bchart.h,v 1.23 2008/03/19 21:21:43 werner Exp $
*/

#ifndef __F5_F5Bchart_H
#define __F5_F5Bchart_H

#include <F5/hdf5inc.h>
#include <F5/F5defs.h>
#include <F5/F5types.h>

#include "F5WinDLLApi.h"

#ifdef  __cplusplus
extern "C"
{
#endif

#if defined(_WINDOWS) || defined(_WIN32) || defined(__MSC_VER__)
#define H5Tname_free(x)
#else
#define H5Tname_free(x) free(x)
#endif

struct  _F5_ChartPrecisionTypes
{
        hid_t   Point_hid_t,
                
                Vector_hid_t,

                Covector_hid_t,

                Bivector_hid_t,

                Bicovector_hid_t,

                Pseudoscalar_hid_t,
                
                Pseudocoscalar_hid_t,
                
                Metric_hid_t,

                Cometric_hid_t;
};

typedef struct _F5_ChartPrecisionTypes F5_ChartPrecisionTypes;

F5_API void F5Bchart_init_tangential_types(F5_ChartPrecisionTypes*H5Types, 
                                           const char*coordnames[],
                                           int dimension,
                                           hid_t coord_type);

        
F5_API void F5Bchart_add_tangential_types(F5_ChartPrecisionTypes*H5Types, 
                                          const char*coordnames[],
                                          int dimension,
                                          hid_t coord_type);

F5_API  void F5B_delete_tangential_types(F5_ChartPrecisionTypes*IDs);

F5_API void F5Bchart_read_tangential_types(F5_ChartPrecisionTypes*H5Types,
                                           hid_t        location);

        
typedef struct 
{
        const char*domain_name;

        unsigned  refs;                 
        int     perm_vector[FIBER_MAX_RANK];
        int     cell_dimensionality;

        F5_ChartPrecisionTypes
/*              HalfPrecision, */
                SinglePrecision,
                DoublePrecision;
/*              QuadPrecision;  */

        F5Ttensor_t**TensorTypes;

}
        ChartDomain_IDs;


F5_API void F5Binsert_tensor_type(ChartDomain_IDs*CD, F5Ttensor_t*);
F5_API F5Ttensor_t*F5Bget_tensor_type_byname(ChartDomain_IDs*CD, const char*name);


        
typedef ChartDomain_IDs*(*new_global_chart_f)();

F5_API ChartDomain_IDs*F5B_new_global_domain(const char*coordnames[],
                                             int dimension,
                                             const char*name,
                                             const int perm_vector[],
                                             int cell_dimensionality,
                                             hid_t SinglePrecisionType,
                                             hid_t DoublePrecisionType
        );


F5_API ChartDomain_IDs*F5B_new_global_float_chart(const char*coordnames[],
                                                  int dimension,
                                                  const char*name,
                                                  const int perm_vector[FIBER_MAX_RANK] );

F5_API ChartDomain_IDs*F5B_new_global_int_domain(const char*coordnames[],
                                                 int dimension,
                                                 const char*name,
                                                 const int perm_vector[FIBER_MAX_RANK],
                                                 int cell_dimensionality);

F5_API ChartDomain_IDs*F5B_new_global_fractional_domain(const char*coordnames[],
                                                        int dimension,
                                                        const char*name,
                                                        const int perm_vector[],
                                                        int cell_dimensionality);

        
/*F5_API ChartDomain_IDs*F5Bcopy_chart(ChartDomain_IDs*oldchart, const char*newname); */

F5_API ChartDomain_IDs*F5Bcopy_chart_domain(const ChartDomain_IDs*IDs);


F5_API int      F5Bget_chart_dims(const ChartDomain_IDs*IDs);
        
#define F5_FORTRAN_ORDER        ((int*)(0))
#define F5_C_ORDER              ((int*)(-1))

F5_API hsize_t*F5B_permute_dimensions(hsize_t*target_dims,
                                      int rank,
                                      const hsize_t*source_dims,
                                      int perm_vector[FIBER_MAX_RANK],
                                      int perm_vector_size);

F5_API const int*F5B_fortran_permute_vector();
        
F5_API void F5B_delete_global_chart(ChartDomain_IDs*ID);

F5_API void F5Bincrement_domain_refs(ChartDomain_IDs*ChartDomain);
        
F5_API ChartDomain_IDs*F5B_save_global_chart(const ChartDomain_IDs*IDs, hid_t Domain_hid);
F5_API ChartDomain_IDs*F5B_read_global_chart(hid_t ChartsID, const char*name);

F5_API ChartDomain_IDs*F5B_new_global_cartesian_chart3D();
        
#ifdef  __cplusplus
}    /* extern "C" */
#endif

#endif  /* __F5_F5Bchart_H */