F5AMR.c

Go to the documentation of this file.

#include "F5AMR.h"
#include "F5B.h"
#include "F5X.h"
#include "F5defs.h"
#include "F5L.h"
#include "F5F.h"
#include "F5R.h"
#include "F5private.h"
#include "F5coordinates.h"
 
#include <assert.h>
 
/*TODO: Implement these functions.
 
AMRInfo*F5AMRopen(hid_t fileID,
                         const int rank,
                         double rootOrigin[3],
                         double rootDelta[3],
                         const char*description);
 
void    F5AMRclose(AMRInfo*);
*/
 
 
/*TODO: rewrite these function to use the F5Path/F5F interfaces */
 
#if 0
herr_t F5AMRwrite_block( AMRInfo*AI,
                         const int    timestep,
                         const double physTime,
                         const char*  fieldname,
                         const int    level,
                         int          dimension,
                         const int    *location,
                         const int    *size,
                         hid_t fieldtype,
                         const void* dataPtr,
                         int   Centering,
                         const char*blockname,
                         const char*coordinate_system,
                         hid_t property_id)
{
hid_t   S_id, G_id, T_id, R_id, F_id, B_id;
char    Tname[1024];
int     i;
hsize_t dims[FIBER_MAX_RANK], start[FIBER_MAX_RANK];
hsize_t levels[3];
 
        if (!coordinate_system)
                coordinate_system = FIBER_HDF5_DEFAULT_CHART;
 
        assert(AI);
 
        /*  /T=0.345/                           */
        S_id = F5Bappend_slice(AI->fileID, &physTime);
 
        /*  /T=0.345/"mydata"                   */
        G_id = F5Gappend(S_id, AI->description, 0);
 
        levels[0] = levels[1] = levels[2] = level;
        assert(levels);
        TopologyName(Tname, sizeof(Tname), levels, Centering, dimension);
 
        /*  /T=0.345/"mydata"/VertexLevel 23/   */
        T_id = F5Gappend(G_id, Tname, 0);
        R_id = F5Gappend(T_id,  coordinate_system, 0) ;
 
        assert(dimension<=32);
        for(i=0; i<dimension; i++)
        {
                start[i] = location[i];
                dims[i]  = size[i];
        }
 
 
        B_id = F5Lwrite_fraction(R_id, fieldname, dimension,
                                 NULL, /* full grid dimension !? */
                                 dims,
                                 fieldtype, 0, /* memory field types not supported here yet. */
                                 dataPtr, start,
                                 0, 0, /* Ghostzones/ border information */
                                 blockname,
                                 F5Tsave_F5field_enum(G_id),
                                 property_id);
 
        /* TODO: Add attributes and relative positional information (parent level/child level information) */
 
        H5Dclose(B_id);
/*      H5Gclose(F_id);  */
        H5Gclose(R_id);
        H5Gclose(T_id);
        H5Gclose(G_id);
        H5Gclose(S_id);
        return 0;
}
 
/**
@todo   Need to adjust this function to the modernized F5Path
        interface, or remove it.
 */
herr_t F5CAMRwrite_block_v(AMRInfo*AI,
                           const int    timestep,
                           const double physTime,
                           const char*  fieldname,
                           const int    level,
                           int          dimension,
                           const int    *location,
                           const int    *size,
                           hid_t fieldtype,
                           const void**dataPtr,
                           int   Centering,
                           const char*blockname,
                           const char*coordinate_system,
                           hid_t property_id)
{
hid_t   S_id, G_id, T_id, R_id, F_id, B_id;
char    Tname[1024];
int     i;
hsize_t dims[32];
 
        assert(AI);
 
        /*  /T=0.345/                           */
        S_id = F5Bappend_slice(AI->fileID, &physTime);
 
        /*  /T=0.345/"mydata"                   */
        G_id = F5Gappend(S_id, AI->description, 0);
 
        {
        hsize_t levels[3];
        levels[0] = levels[1] = levels[2] = level;
        TopologyName(Tname, sizeof(Tname), levels, Centering, dimension);
        }
 
        /*  /T=0.345/"mydata"/VertexLevel 23/   */
        T_id = F5Gappend(G_id, Tname, 0);
        /*  /T=0.345/"mydata"/VertexLevel 23/Cartesian Chart 3D/        */
        R_id = F5Gappend(T_id,  coordinate_system, 0) ;
 
        assert(dimension<=32);
        for(i=0; i<dimension; i++)
        {
                dims[i] = size[i];
        }
 
        /*  /T=0.345/"mydata"/VertexLevel 23/Cartesian Chart 3D/"density"/      */
        F_id = F5Gappend(R_id,  fieldname, 0) ;
 
        /* TODO: Add reference to field type here */
 
        /* write field data block-wise */
 
        /*  /T=0.345/"mydata"/VertexLevel 23/Cartesian Chart 3D/"density"/block-000/    */
        B_id = F5Lwrites(F_id, blockname, dimension, dims, fieldtype, 0, dataPtr,
                         F5Tsave_F5field_enum(G_id),
                         property_id);
 
 
#if 0
        B_id = H5Gcreate(F_id, blockname, 0);
 
        /* write the data into the dataset */
        {
        hid_t   ds_id = H5Screate_simple(dimension, dims, NULL);
        hid_t   creator_id = H5Pcreate(H5P_DATASET_CREATE);
 
        int     Ncomponents = H5Tget_nmembers(fieldtype);
        int     comp;
                for(comp=0; comp<Ncomponents; comp++)
                {
                hid_t   id, ComponentType;
                char       * ComponentName;
                const void * data = dataPtr[comp];
 
                        if (!data) continue;
                        ComponentName = H5Tget_member_name(fieldtype, comp);
                        ComponentType = H5Tget_member_type(fieldtype, comp);
 
                        id = H5Dcreate(B_id, ComponentName, fieldtype, ds_id, creator_id);
                        H5Dwrite(id, fieldtype, H5S_ALL, H5S_ALL, H5P_DEFAULT, data);
                        H5Dclose(id);
 
                        free(ComponentName);
                        H5Tclose(ComponentType);
                }
 
                H5Sclose(ds_id);
                H5Pclose(creator_id);
        }
#endif
        H5Gclose(B_id);
        /* TODO: Add attributes and relative positional information (parent level/child level information) */
 
        H5Gclose(F_id);
        H5Gclose(R_id);
        H5Gclose(T_id);
        H5Gclose(G_id);
        H5Gclose(S_id);
        return 0;
}
#endif