F5L.h

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

#ifndef __F5Lield_H
#define __F5Lield_H

#include <F5/hdf5inc.h>

#include "F5WinDLLApi.h"
#include <F5/F5defs.h>

#ifdef  __cplusplus
#ifdef  F5_CPLUSPLUS
namespace F5
{

template <int N, typename T, typename Trait=void> 
  struct Extent
{
        T[N]    data;
};

} /* namespace F5 */
#endif /* F5_CPLUSPLUS */

extern "C"
{
#endif

#define F5P_DEFAULT     (-1)


F5_API hid_t    F5Lcreate(hid_t R_id, const char*fieldname, int dimension, const hsize_t*dims,
                          hid_t fieldtype, hid_t property_id);


F5_API hid_t    F5LTcreate_topology(hid_t Grid_hid,
                                    const char*TopologyName,
                                    int IndexDepth,
                                    int SkeletonDimensionality,
                                    int DataspaceDimensionality,
                                    const hsize_t*refinement);

F5_API  int     F5LTset_maximal_grid_refinement(hid_t ContentsGroup_id,
                                                const char*gridname,
                                                int     Dimensionality,
                                                const hsize_t*refinement);

F5_API  int     F5LTget_maximal_grid_refinement(hid_t ContentsGroup_id,
                                                const char*gridname,
                                                hsize_t*refinement);
        

F5_API int   F5LTget_index_depth(hid_t Top_hid);


F5_API hid_t    F5Lwrite(hid_t R_id, const char*fieldname,
                         int dimension, const hsize_t*dims,
                         hid_t fieldtype, hid_t memtype, const void * dataPtr,
                         hid_t enum_type, hid_t property_id);


F5_API hid_t    F5Lwrite_fraction(hid_t R_id, const char*fieldname,
                                  int dimension, const hsize_t*full_dims,
                                  const hsize_t*datasize,
                                  hid_t fieldtype, hid_t memtype, const void * dataPtr,
                                  const hsize_t*datastart,
                                  const hsize_t*start_border, const hsize_t*end_border,
                                  const char*fraction_name, hid_t enum_type, hid_t property_id);        

F5_API hid_t    F5LSwrite_fraction(hid_t R_id, const char*fieldname,
                                   int dimension, const hsize_t*full_dims,
                                   const hsize_t*datasize,
                                   hid_t fieldtype, hid_t memtype, const void *const* dataPtr,
                                   const hsize_t*datastart,
                                   const hsize_t*start_border, const hsize_t*end_border,
                                   const char*fraction_name,
                                   hid_t enum_type, hid_t property_id); 



F5_API hid_t    F5Lwrite1D(hid_t R_id, const char*fieldname, hsize_t nElements,
                           hid_t fieldtype, hid_t memtype, const void * dataPtr,
                           hid_t enum_type, hid_t property_id);


F5_API hid_t    F5Lwrite_linear(hid_t R_id, const char*fieldname,
                                int dimension, const hsize_t*dataspace_dims,
                                hid_t fieldtype, const void * base, const void * delta);


F5_API int F5Lis_linear(hid_t Rep_id, const char*fieldname);

F5_API hid_t    F5Lread_linear(hid_t F_id,
                               hsize_t*dims,
                               hid_t fieldtype,
                               void * base, void * delta);


F5_API hid_t    F5Lwrite_entity(hid_t R_id, const char*fieldname,
                                hid_t fieldtype,
                                const void * value);



F5_API hid_t F5Lwrites(hid_t F_id, const char*fieldname,
                       int dimension, const hsize_t*dims,
                       hid_t fieldtype, hid_t memtype,
                       const void *const*dataPtr,
                       hid_t enum_type,
                       hid_t property_id);

F5_API hid_t F5LwriteX(hid_t R_id, const char*fieldname,
                       int Mrank, const hsize_t dims[],
                       hid_t fieldtype, hid_t memtype, 
                       const void *const dataPtr[], const int polynom_order[],
                       const int component_map[],
                       hid_t property_id, hid_t ProcIDLocation);


F5_API hid_t    F5Lwrite_permutation(hid_t R_id, const char*fieldname,
                                     hsize_t size,
                                     hsize_t*permutation_indices,
                                     hsize_t nvalues,
                                     hid_t fieldtype, const void * dataPtr);


F5_API hid_t F5Lget_type(hid_t Field_hid, int FieldIDisGroup);


F5_API int      F5Lget_field_dimension_and_type(hid_t Representation_hid,
                                             const char*fieldname,
                                             hsize_t dims[FIBER_MAX_RANK], hid_t*type_id);

F5_API int      F5LAsave_dimensions(hid_t Field_id, const char*aname, int rank, const hsize_t*dims);
F5_API int      F5LAget_dimensions(hid_t Field_id, const char*aname, hsize_t dims[FIBER_MAX_RANK]);
F5_API void     F5Lset_space(hid_t Field_id, int rank, const hsize_t*dims);

F5_API void     F5Dset_type_attribute(hid_t F_id, hid_t fieldtype);
        
/*
  already declared
  
F5_API hid_t    F5LSwrite_fraction(hid_t R_id, const char*fieldname,
                                   int rank, const hsize_t*full_dims,
                                   const hsize_t*datasize,
                                   hid_t fieldtype, hid_t memtype, void *const*dataPtr,
                                   const hsize_t*datastart,
                                   const hsize_t*start_border, const hsize_t*end_border,
                                   const char*fraction_name,
                                   hid_t enum_type, hid_t property_id);
*/

        
#ifdef  __cplusplus
}    /* extern "C" */

#ifdef  F5_CPLUSPLUS

namespace F5
{

template <int N, typename Trait>
inline hid_t LwriteX(hid_t R_id, const char*fieldname, 
                     const Extent<N, hsize_t, Trait>&dims,
                     hid_t fieldtype, hid_t memtype, 
                     const Extent<N, const void*, Trait>&dataPtr,
                     const Extent<N, int, Trait>&polynom_order,
                     hid_t property_id = F5P_DEFAULT)
{
        return F5LwriteX(R_id, fieldname, N, dims.data,
                        fieldtype, memtype, 
                        dataPtr.data, polynom_order.data,
                        property_id);

}

template <int N, typename Trait>
inline hid_t LwriteX(hid_t R_id, const char*fieldname, 
                     const Extent<N, hsize_t, Trait>&dims,
                     hid_t fieldtype,
                     const Extent<N, const void*, Trait>&dataPtr,
                     const Extent<N, int, Trait>&polynom_order,
                     hid_t property_id = F5P_DEFAULT)
{
        return F5LwriteX(R_id, fieldname, N, dims.data,
                        fieldtype, 0, 
                        dataPtr.data, polynom_order.data,
                        property_id = F5P_DEFAULT);

}

template <typename T, typename value>
inline T evalPolynom(const value&v, const T*data, int n)
{
        if (n==0) return 0;
        if (n==1) return T[0];
        if (n==1) return T[0] + v*T[1];
        return T[0] + v*evalPolynom(v, data+1, n-1);
}


} /* namespace F5 */

#endif  /* F5_CPLUSPLUS */

#endif

#endif  /* __F5Lield_H */