F5particles.c

Go to the documentation of this file.

#include "F5particles.h"
#include "F5B.h"
#include "F5R.h"
#include "F5F.h"
#include "F5X.h"
#include "F5Fmeta.h"
 
#include "F5defs.h"
#include "F5Bchart.h"
#include "F5L.h"
#include "F5private.h"
 
#include <stdarg.h>
#include <assert.h>
 
F5ErrorCode F5write_particle_positions(hid_t File_id, double time,
                                       const char*gridname,
                                       const F5_vec3_point_t*Coords, int nCoords,
                                       const char*coordinate_system)
{
F5Path*f = F5Rcreate_cartesian_nD(File_id, time, gridname, 0, coordinate_system);
        assert(f);
 
F5ErrorCode rv = F5Fwrite_1D( f, FIBER_HDF5_POSITIONS_STRING, nCoords,  F5T_COORD3_FLOAT, F5T_COORD3_FLOAT, Coords, H5P_DEFAULT);
 
        F5close(f);
 
        return rv;
}
 
 
/** Creates a grid object of a point cloud (vertices) and additional fields.
    Additional fields can be added in the ellipses ( ... ) by specifying
    the field name, the data type and the pointer to the data array:
    e.g.: ( ... ,"Velocity", F5T_VEC3_FLOAT, V,
                ,"Acceleration", F5T_VEC3_FLOAT, A, ... )
*/
int F5write_particle_cartesian3Dv(hid_t File_id, double time, const char*gridname,
                                  const F5_vec3_point_t*Coords, int nCoords,
                                  const char*coordinate_system,
                                  ...)
{
const char* fieldname;
hid_t fieldtype;
const void* dataPtr;
va_list vl;
int NumOfWrittenFields = 0;
        //const char* unittype;
 
        /* Create the Field */
F5Path*f = F5Rcreate_cartesian_nD(File_id, time, gridname, 0, coordinate_system);
 
        assert(f);
 
hid_t   creator_hid = H5Pcreate(H5P_DATASET_CREATE);
F5ErrorCode rv;
 
        F5printf(20, "F5write_particle_cartesian3Dv(): write positional information");
 
        rv = F5Fwrite_1D( f, FIBER_HDF5_POSITIONS_STRING, nCoords,  F5T_COORD3_FLOAT, F5T_COORD3_FLOAT, Coords, H5P_DEFAULT);
 
        /* Init the step through the argument list */
        va_start(vl, coordinate_system);
 
        F5printf(20, "F5write_particle_cartesian3Dv(): Step through the argument list");
        /* Step through the argument list */
        for(;;)
        {
        hid_t   file_type;
 
                /* Get the Name of the field */
                fieldname = va_arg(vl, const char*);
                if (fieldname == NULL)  break;
 
                /* Get the type of the field */
                fieldtype = va_arg(vl, hid_t);
                /* Get the pointer to the field */
                dataPtr = va_arg(vl, const void*);
                /* Get the unit of the field*/
                //unittype = va_arg(vl, const void*);
                if (dataPtr == NULL) break; 
 
                F5printf(20, "F5write_particle_cartesian3Dv(): write data for field %s", fieldname);
 
                file_type = F5file_type(f, fieldtype);
 
                /* Write the stuff */
 
                /* WHY NOT USE F5Fwrite here?? */
                f->Field_hid = F5Lwrite1D(f->Representation_hid,
                                          fieldname, nCoords, file_type, 0, dataPtr,
                                          F5LTmake_enum_type(f), creator_hid); 
 
                F5I_add_field(f, fieldname);
                H5Dclose(f->Field_hid);
                f->Field_hid = 0;
                NumOfWrittenFields++;
        }
 
        /* Close the Step through */
        va_end(vl);
 
        H5Pclose(creator_hid);
 
        /* Close the field */
        F5close(f);
 
        return NumOfWrittenFields;
}
 
int F5write_particle_double_cartesian3Dv(hid_t File_id, double time, const char*gridname,
                                  const F5_point3_double_t*Coords, int nCoords,
                                  const char*coordinate_system,
                                  ...)
{
const char* fieldname;
hid_t fieldtype;
const void* dataPtr;
va_list vl;
int NumOfWrittenFields = 0;
        //const char* unittype;
 
        /* Create the Field */
F5Path*f = F5Rcreate_cartesian_nD(File_id, time, gridname, 0, coordinate_system);
 
        assert(f);
 
hid_t   creator_hid = H5Pcreate(H5P_DATASET_CREATE);
F5ErrorCode rv;
 
        F5printf(20, "F5write_particle_cartesian3Dv(): write positional information");
 
        rv = F5Fwrite_1D( f, FIBER_HDF5_POSITIONS_STRING, nCoords,  F5T_COORD3_DOUBLE, F5T_COORD3_DOUBLE, Coords, H5P_DEFAULT);
 
        /* Init the step through the argument list */
        va_start(vl, coordinate_system);
 
        F5printf(20, "F5write_particle_cartesian3Dv(): Step through the argument list");
        /* Step through the argument list */
        for(;;)
        {
        hid_t   file_type;
 
                /* Get the Name of the field */
                fieldname = va_arg(vl, const char*);
                /* Get the type of the field */
                fieldtype = va_arg(vl, hid_t);
                /* Get the pointer to the field */
                dataPtr = va_arg(vl, const void*);
                /* Get the unit of the field*/
                //unittype = va_arg(vl, const void*);
 
                /* Break? */
                if ( (fieldname == NULL) || (dataPtr == NULL) )
                        break;
 
                F5printf(20, "F5write_particle_cartesian3Dv(): write data");
 
                file_type = F5file_type(f, fieldtype);
 
                /* Write the stuff */
 
                /* WHY NOT USE F5Fwrite here?? */
                f->Field_hid = F5Lwrite1D(f->Representation_hid,
                                          fieldname, nCoords, file_type, 0, dataPtr,
                                          F5LTmake_enum_type(f), creator_hid); 
 
                F5I_add_field(f, fieldname);
                H5Dclose(f->Field_hid);
                f->Field_hid = 0;
                NumOfWrittenFields++;
        }
 
        /* Close the Step through */
        va_end(vl);
 
        H5Pclose(creator_hid);
 
        /* Close the field */
        F5close(f);
 
        return NumOfWrittenFields;
}
 
 
 
/** added unit support in convenience function */
int F5write_particle_cartesian3DvU(hid_t File_id, double time, const char*gridname,
                                   const F5_vec3_point_t*Coords, const char*CoordUnit, int nCoords,
                                   const char*coordinate_system,
                                   ...)
{
        const char* fieldname;
        hid_t fieldtype;
        const void* dataPtr;
        va_list vl;
        int NumOfWrittenFields = 0;
        const char* unittype;
 
 
        /* Create the Field */
F5Path*f = F5Rcreate_cartesian_nD(File_id, time, gridname, 0, coordinate_system);
 
hid_t   creator_hid = H5Pcreate(H5P_DATASET_CREATE);
 
        F5printf(20, "F5write_particle_cartesian3DvU(): write positional information");
        F5Fwrite_1D( f, FIBER_HDF5_POSITIONS_STRING, nCoords,  F5T_COORD3_FLOAT, F5T_COORD3_FLOAT, Coords, H5P_DEFAULT );
 
        F5Fset_metric_unit_description(f, CoordUnit );
 
#if     0
        /* write the positional information.
           yack. should use F5F functions here...!
         */
{
hsize_t dims = nCoords;
hid_t   coord_ds = H5Screate_simple(1, &dims, NULL);
hid_t   id = F5Dappend(f->Representation_hid, FIBER_HDF5_POSITIONS_STRING, f->FileIDs->Point_hid_t, coord_ds, creator_hid);
        H5Dwrite(id, f->FileIDs->Point_hid_t, H5S_ALL, H5S_ALL, H5P_DEFAULT, Coords);
 
        f->Field_hid = id;
        F5I_add_field(f, FIBER_HDF5_POSITIONS_STRING, FIBER_HDF5_SIMPLE_FIELD);
        f->Field_hid = 0;
        H5Dclose(id);
        H5Sclose(coord_ds);
}
#endif
 
        /* Init the step through the argument list */
        va_start(vl, coordinate_system);
 
        F5printf(20, "F5write_particle_cartesian3DvU(): Step through the argument list");
        /* Step through the argument list */
 
        for(;;)
        {
        hid_t   file_type; // MSVC wants it to be at the beginning of the scope
 
                F5printf(20, "F5write_particle_cartesian3DvU(): Step %d", NumOfWrittenFields);
 
 
                /* Get the Name of the field */
                fieldname = va_arg(vl, const char*);
                /* Get the type of the field */
                fieldtype = va_arg(vl, hid_t);
                /* Get the pointer to the field */
                dataPtr = va_arg(vl, const void*);
                /* Get the unit of the field*/ 
                unittype = (const char*)va_arg(vl, const void*);
 
                /* Break? */
                if ( (fieldname == NULL) || (dataPtr == NULL) || unittype == NULL )
                        break;
 
                F5printf(20, "F5write_particle_cartesian3DvU(): write data");
 
                file_type = F5file_type(f, fieldtype);
 
                /* Write the stuff */
 
                /* WHY NOT USE F5Fwrite here?? */
                f->Field_hid = F5Lwrite1D(f->Representation_hid,
                                          fieldname, nCoords, file_type, 0, dataPtr,
                                          F5LTmake_enum_type(f), creator_hid); 
 
                F5I_add_field(f, fieldname);
 
                F5Fset_metric_unit_description(f, unittype );
 
                H5Dclose(f->Field_hid);
 
 
 
                f->Field_hid = 0;
 
                F5printf(20, "F5write_particle_cartesian3DvU():4 fieldid=%d", (int)f->Field_hid );
                NumOfWrittenFields++;
        }
 
        /* Close the Step through */
        va_end(vl);
 
        H5Pclose(creator_hid);
 
        /* Close the field */
        F5close(f);
 
        return NumOfWrittenFields;
}