HDF5Float64.cc

Go to the documentation of this file.

// This file is part of hdf5_handler a HDF5 file handler for the OPeNDAP
// data server.

// Author: Hyo-Kyung Lee <hyoklee@hdfgroup.org> and Muqun Yang
// <myang6@hdfgroup.org> 

// Copyright (c) 2009-2013 The HDF Group, Inc. and OPeNDAP, Inc.
//
// This is free software; you can redistribute it and/or modify it under the
// terms of the GNU Lesser General Public License as published by the Free
// Software Foundation; either version 2.1 of the License, or (at your
// option) any later version.
//
// This software is distributed in the hope that it will be useful, but
// WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
// or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
// License for more details.
//
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
//
// You can contact OPeNDAP, Inc. at PO Box 112, Saunderstown, RI. 02874-0112.
// You can contact The HDF Group, Inc. at 1901 South First Street,
// Suite C-2, Champaign, IL 61820  

// #define DODS_DEBUG




#include <string>
#include <ctype.h>
#include "config_hdf5.h"
#include "InternalErr.h"
#include "h5dds.h"
#include "HDF5Float64.h"
#include "HDF5Structure.h"
#include "debug.h"

typedef struct s2_float64_t {
    dods_float64 a;
} s2_float64_t;



HDF5Float64::HDF5Float64(const string & n, const string &d) : Float64(n, d)
{
    ty_id = -1;
    dset_id = -1;
}

BaseType *HDF5Float64::ptr_duplicate()
{
    return new HDF5Float64(*this);
}

bool HDF5Float64::read()
{
    if (read_p())
        return true;

    if (get_dap_type(ty_id) == "Float64") {
        dods_float64 buf;
        get_data(dset_id, (void *) &buf);
        set_read_p(true);
        set_value(buf);

        // Release the handles.
        if (H5Tclose(ty_id) < 0) {
            throw InternalErr(__FILE__, __LINE__, "Unable to close the datatype.");
        }
        if (H5Dclose(dset_id) < 0) {
            throw InternalErr(__FILE__, __LINE__, "Unable to close the dset.");
        }

    }

    if (get_dap_type(ty_id) == "Structure") {
        DBG(cerr << "=read(): Structure" << endl);
        BaseType *q = get_parent();
        if (!q)
            throw InternalErr(__FILE__, __LINE__, "null pointer");
        HDF5Structure &p = static_cast<HDF5Structure &> (*q);
        DBG(cerr << "=read(): Size = " << p.get_entire_array_size() << endl);

#ifdef DODS_DEBUG
        int i = H5Tget_nmembers(ty_id);
        if(i < 0) {
            throw InternalErr(__FILE__, __LINE__, "H5Tget_nmembers() failed.");
        }
#endif
        int j = 0;
        int k = 0;

        vector<s2_float64_t> buf(p.get_entire_array_size());

        string myname = name();
        string parent_name;

        hid_t s2_float64_tid = H5Tcreate(H5T_COMPOUND, sizeof(s2_float64_t));
        hid_t stemp_tid;

        if (s2_float64_tid < 0) {
            throw InternalErr(__FILE__, __LINE__, "cannot create a new datatype");
        }

        DBG(cerr << "=read() ty_id=" << ty_id << " name=" << myname <<
                " size=" << i << endl);
        while (q != NULL) {
            if (q->is_constructor_type()) { // Grid, structure or sequence
                if (k == 0) {
                    // Bottom level structure
                    if (H5Tinsert(s2_float64_tid, myname.c_str(), HOFFSET(s2_float64_t, a), H5T_NATIVE_DOUBLE) < 0) {
                        throw InternalErr(__FILE__, __LINE__, "Unable to add datatype.");
                    }
                }
                else {
                    stemp_tid = H5Tcreate(H5T_COMPOUND, sizeof(s2_float64_t));
                    if (stemp_tid < 0) {
                        throw InternalErr(__FILE__, __LINE__, "cannot create a new datatype");
                    }
                    if (H5Tinsert(stemp_tid, parent_name.c_str(), 0, s2_float64_tid) < 0) {
                        throw InternalErr(__FILE__, __LINE__, "Unable to add datatype.");
                    }
                    s2_float64_tid = stemp_tid;
                }
                parent_name = q->name();
                p = static_cast<HDF5Structure &> (*q);
                // Remember the index of array from the last parent.
                j = p.get_array_index();
                q = q->get_parent();
            }
            else {
                q = NULL;
            }
            k++;
        }

        if (H5Dread(dset_id, s2_float64_tid, H5S_ALL, H5S_ALL, H5P_DEFAULT, &buf[0]) < 0) {
            // buf is deleted in the catch ... block below and
            // should not be deleted here. pwest Mar 18, 2009
            //delete[] buf;
            throw InternalErr(__FILE__, __LINE__, "hdf5_dods server failed when getting int32 data for structure");
            //string
            //()
            //+ Msgi);
        }
        set_read_p(true);
        DBG(cerr << "index " << j << endl);

        set_value(buf[j].a);
    }

    return true;
}

void HDF5Float64::set_did(hid_t dset)
{
    dset_id = dset;
}

void HDF5Float64::set_tid(hid_t type)
{
    ty_id = type;
}

hid_t HDF5Float64::get_did()
{
    return dset_id;
}

hid_t HDF5Float64::get_tid()
{
    return ty_id;
}