hdf5_wall.c

Go to the documentation of this file.

 
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <hdf5.h>
#include <hdf5_hl.h>
#include "../wall.h"
#include "../wall/wall_2d.h"
#include "../wall/wall_3d.h"
#include "hdf5_wall.h"
#include "hdf5_helpers.h"
 
#define WPATH 
int hdf5_wall_read_2D(hid_t f, wall_2d_data* data, char* qid);
int hdf5_wall_read_3D(hid_t f, wall_3d_data* data, char* qid);
 
int hdf5_wall_init(hid_t f, wall_data* data, char* qid) {
 
    char path[256];
    int err = 1;
 
    /* Read data the QID corresponds to */
    hdf5_gen_path("/wall/wall_2D_XXXXXXXXXX", qid, path);
    if(hdf5_find_group(f, path) == 0) {
        data->type = wall_type_2D;
        err = hdf5_wall_read_2D(f, &(data->w2d), qid);
    }
    hdf5_gen_path("/wall/wall_3D_XXXXXXXXXX", qid, path);
    if(hdf5_find_group(f, path) == 0) {
        data->type = wall_type_3D;
        err = hdf5_wall_read_3D(f, &(data->w3d), qid);
    }
    return err;
}
 
int hdf5_wall_read_2D(hid_t f, wall_2d_data* data, char* qid) {
    #undef WPATH
    #define WPATH "/wall/wall_2D_XXXXXXXXXX/"
 
    int nelements;
    if( hdf5_read_int(WPATH "nelements", &nelements,
                      f, qid, __FILE__, __LINE__) ) {return 1;}
 
    real* r = (real*) malloc(nelements * sizeof(real));
    real* z = (real*) malloc(nelements * sizeof(real));
    if( hdf5_read_double(WPATH "r", r, f, qid, __FILE__, __LINE__) ) {return 1;}
    if( hdf5_read_double(WPATH "z", z, f, qid, __FILE__, __LINE__) ) {return 1;}
 
    int err = wall_2d_init(data, nelements, r, z);
    free(r);
    free(z);
    return err;
}
 
int hdf5_wall_read_3D(hid_t f, wall_3d_data* data, char* qid) {
    #undef WPATH
    #define WPATH "/wall/wall_3D_XXXXXXXXXX/"
 
    int nelements;
    if( hdf5_read_int(WPATH "nelements", &nelements,
                      f, qid, __FILE__, __LINE__) ) {return 1;}
 
    /* Allocate temporary arrays for x1x2x3, y1y2y3, z1z2z3 for each triangle */
    real* x1x2x3 = (real*)malloc(3 * nelements * sizeof(real));
    real* y1y2y3 = (real*)malloc(3 * nelements * sizeof(real));
    real* z1z2z3 = (real*)malloc(3 * nelements * sizeof(real));
 
    if( hdf5_read_double(WPATH "x1x2x3", x1x2x3,
                         f, qid, __FILE__, __LINE__) ) {return 1;}
    if( hdf5_read_double(WPATH "y1y2y3", y1y2y3,
                         f, qid, __FILE__, __LINE__) ) {return 1;}
    if( hdf5_read_double(WPATH "z1z2z3", z1z2z3,
                         f, qid, __FILE__, __LINE__) ) {return 1;}
 
    int err = wall_3d_init(data, nelements, x1x2x3, y1y2y3, z1z2z3);
    free(x1x2x3);
    free(y1y2y3);
    free(z1z2z3);
    return err;
}