Blame Float64.cc

Packit a4aae4
Packit a4aae4
// -*- mode: c++; c-basic-offset:4 -*-
Packit a4aae4
Packit a4aae4
// This file is part of libdap, A C++ implementation of the OPeNDAP Data
Packit a4aae4
// Access Protocol.
Packit a4aae4
Packit a4aae4
// Copyright (c) 2002,2003 OPeNDAP, Inc.
Packit a4aae4
// Author: James Gallagher <jgallagher@opendap.org>
Packit a4aae4
//
Packit a4aae4
// This library is free software; you can redistribute it and/or
Packit a4aae4
// modify it under the terms of the GNU Lesser General Public
Packit a4aae4
// License as published by the Free Software Foundation; either
Packit a4aae4
// version 2.1 of the License, or (at your option) any later version.
Packit a4aae4
//
Packit a4aae4
// This library is distributed in the hope that it will be useful,
Packit a4aae4
// but WITHOUT ANY WARRANTY; without even the implied warranty of
Packit a4aae4
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
Packit a4aae4
// Lesser General Public License for more details.
Packit a4aae4
//
Packit a4aae4
// You should have received a copy of the GNU Lesser General Public
Packit a4aae4
// License along with this library; if not, write to the Free Software
Packit a4aae4
// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
Packit a4aae4
//
Packit a4aae4
// You can contact OPeNDAP, Inc. at PO Box 112, Saunderstown, RI. 02874-0112.
Packit a4aae4
Packit a4aae4
// (c) COPYRIGHT URI/MIT 1994-1999
Packit a4aae4
// Please read the full copyright statement in the file COPYRIGHT_URI.
Packit a4aae4
//
Packit a4aae4
// Authors:
Packit a4aae4
//      jhrg,jimg       James Gallagher <jgallagher@gso.uri.edu>
Packit a4aae4
Packit a4aae4
// Implementation for Float64.
Packit a4aae4
//
Packit a4aae4
// jhrg 9/7/94
Packit a4aae4
Packit a4aae4
#include "config.h"
Packit a4aae4
Packit a4aae4
//#define DODS_DEBUG
Packit a4aae4
Packit a4aae4
#include <sstream>
Packit a4aae4
#include <iomanip>
Packit a4aae4
Packit a4aae4
#include "Byte.h"           // synonymous with UInt8 and Char
Packit a4aae4
#include "Int8.h"
Packit a4aae4
#include "Int16.h"
Packit a4aae4
#include "UInt16.h"
Packit a4aae4
#include "Int32.h"
Packit a4aae4
#include "UInt32.h"
Packit a4aae4
#include "Int64.h"
Packit a4aae4
#include "UInt64.h"
Packit a4aae4
#include "Float32.h"
Packit a4aae4
#include "Float64.h"
Packit a4aae4
#include "Str.h"
Packit a4aae4
#include "Url.h"
Packit a4aae4
Packit a4aae4
#include "DDS.h"
Packit a4aae4
#include "Marshaller.h"
Packit a4aae4
#include "UnMarshaller.h"
Packit a4aae4
Packit a4aae4
#include "DMR.h"
Packit a4aae4
#include "D4StreamMarshaller.h"
Packit a4aae4
#include "D4StreamUnMarshaller.h"
Packit a4aae4
Packit a4aae4
#include "util.h"
Packit a4aae4
#include "parser.h"
Packit a4aae4
#include "Operators.h"
Packit a4aae4
#include "dods-limits.h"
Packit a4aae4
#include "InternalErr.h"
Packit a4aae4
Packit a4aae4
Packit a4aae4
using std::cerr;
Packit a4aae4
using std::endl;
Packit a4aae4
Packit a4aae4
namespace libdap {
Packit a4aae4
Packit a4aae4
/** The Float64 constructor requires only the name of the variable
Packit a4aae4
    to be created.  The name may be omitted, which will create a
Packit a4aae4
    nameless variable.  This may be adequate for some applications.
Packit a4aae4
Packit a4aae4
    @param n A string containing the name of the variable to be
Packit a4aae4
    created.
Packit a4aae4
Packit a4aae4
*/
Packit a4aae4
Float64::Float64(const string &n) : BaseType(n, dods_float64_c), d_buf(0)
Packit a4aae4
{}
Packit a4aae4
Packit a4aae4
/** The Float64 server-side constructor accepts the name of the variable and
Packit a4aae4
    the dataset name from which this instance is created.
Packit a4aae4
Packit a4aae4
    @param n A string containing the name of the variable to be created.
Packit a4aae4
    @param d A string containing the name of the dataset from which this
Packit a4aae4
    variable is created
Packit a4aae4
*/
Packit a4aae4
Float64::Float64(const string &n, const string &d) : BaseType(n, d, dods_float64_c), d_buf(0)
Packit a4aae4
{}
Packit a4aae4
Packit a4aae4
Float64::Float64(const Float64 &copy_from) : BaseType(copy_from)
Packit a4aae4
{
Packit a4aae4
    d_buf = copy_from.d_buf;
Packit a4aae4
}
Packit a4aae4
Packit a4aae4
BaseType *
Packit a4aae4
Float64::ptr_duplicate()
Packit a4aae4
{
Packit a4aae4
    return new Float64(*this);
Packit a4aae4
}
Packit a4aae4
Packit a4aae4
Float64 &
Packit a4aae4
Float64::operator=(const Float64 &rhs)
Packit a4aae4
{
Packit a4aae4
    if (this == &rhs)
Packit a4aae4
        return *this;
Packit a4aae4
Packit a4aae4
    dynamic_cast<BaseType &>(*this) = rhs;
Packit a4aae4
Packit a4aae4
    d_buf = rhs.d_buf;
Packit a4aae4
Packit a4aae4
    return *this;
Packit a4aae4
}
Packit a4aae4
Packit a4aae4
unsigned int
Packit a4aae4
Float64::width(bool) const
Packit a4aae4
{
Packit a4aae4
    return sizeof(dods_float64);
Packit a4aae4
}
Packit a4aae4
Packit a4aae4
bool
Packit a4aae4
Float64::serialize(ConstraintEvaluator &eval, DDS &dds, Marshaller &m, bool ce_eval)
Packit a4aae4
{
Packit a4aae4
#if USE_LOCAL_TIMEOUT_SCHEME
Packit a4aae4
    dds.timeout_on();
Packit a4aae4
#endif
Packit a4aae4
    if (!read_p())
Packit a4aae4
        read();  // read() throws Error and InternalErr
Packit a4aae4
Packit a4aae4
    if (ce_eval && !eval.eval_selection(dds, dataset()))
Packit a4aae4
        return true;
Packit a4aae4
#if USE_LOCAL_TIMEOUT_SCHEME
Packit a4aae4
    dds.timeout_off();
Packit a4aae4
#endif
Packit a4aae4
    m.put_float64( d_buf ) ;
Packit a4aae4
Packit a4aae4
    return true;
Packit a4aae4
}
Packit a4aae4
Packit a4aae4
bool
Packit a4aae4
Float64::deserialize(UnMarshaller &um, DDS *, bool)
Packit a4aae4
{
Packit a4aae4
    um.get_float64( d_buf ) ;
Packit a4aae4
Packit a4aae4
    return false;
Packit a4aae4
}
Packit a4aae4
Packit a4aae4
void
Packit a4aae4
Float64::compute_checksum(Crc32 &checksum)
Packit a4aae4
{
Packit a4aae4
	checksum.AddData(reinterpret_cast<uint8_t*>(&d_buf), sizeof(d_buf));
Packit a4aae4
}
Packit a4aae4
Packit a4aae4
/**
Packit a4aae4
 * @brief Serialize an Int8
Packit a4aae4
 * @param m
Packit a4aae4
 * @param dmr Unused
Packit a4aae4
 * @param eval Unused
Packit a4aae4
 * @param filter Unused
Packit a4aae4
 * @exception Error is thrown if the value needs to be read and that operation fails.
Packit a4aae4
 */
Packit a4aae4
void
Packit a4aae4
Float64::serialize(D4StreamMarshaller &m, DMR &, /*ConstraintEvaluator &,*/ bool)
Packit a4aae4
{
Packit a4aae4
    if (!read_p())
Packit a4aae4
        read();          // read() throws Error
Packit a4aae4
Packit a4aae4
    m.put_float64( d_buf ) ;
Packit a4aae4
}
Packit a4aae4
Packit a4aae4
void
Packit a4aae4
Float64::deserialize(D4StreamUnMarshaller &um, DMR &)
Packit a4aae4
{
Packit a4aae4
    um.get_float64( d_buf ) ;
Packit a4aae4
}
Packit a4aae4
Packit a4aae4
unsigned int
Packit a4aae4
Float64::val2buf(void *val, bool)
Packit a4aae4
{
Packit a4aae4
    // Jose Garcia
Packit a4aae4
    // This method is public therefore and I believe it has being designed
Packit a4aae4
    // to be use by read which must be implemented on the surrogated library,
Packit a4aae4
    // thus if the pointer val is NULL, is an Internal Error.
Packit a4aae4
    if (!val)
Packit a4aae4
        throw InternalErr(__FILE__, __LINE__,
Packit a4aae4
                          "The incoming pointer does not contain any data.");
Packit a4aae4
Packit a4aae4
    d_buf = *(dods_float64 *)val;
Packit a4aae4
Packit a4aae4
    return width();
Packit a4aae4
}
Packit a4aae4
Packit a4aae4
unsigned int
Packit a4aae4
Float64::buf2val(void **val)
Packit a4aae4
{
Packit a4aae4
    // Jose Garcia
Packit a4aae4
    // The same comment justifying throwing an Error in val2buf applies here.
Packit a4aae4
    if (!val)
Packit a4aae4
        throw InternalErr(__FILE__, __LINE__, "NULL pointer.");
Packit a4aae4
Packit a4aae4
    if (!*val)
Packit a4aae4
        *val = new dods_float64;
Packit a4aae4
Packit a4aae4
    *(dods_float64 *)*val = d_buf;
Packit a4aae4
Packit a4aae4
    return width();
Packit a4aae4
}
Packit a4aae4
Packit a4aae4
/** Return the value of the Float64 held by this instance. This is more
Packit a4aae4
    convenient than the general interface provided by buf2val, but its use
Packit a4aae4
    requires a downcase from BaseType to Float64.
Packit a4aae4
Packit a4aae4
    @return The dods_float32 value. */
Packit a4aae4
dods_float64
Packit a4aae4
Float64::value() const
Packit a4aae4
{
Packit a4aae4
    return d_buf;
Packit a4aae4
}
Packit a4aae4
Packit a4aae4
bool
Packit a4aae4
Float64::set_value(dods_float64 val)
Packit a4aae4
{
Packit a4aae4
    d_buf = val;
Packit a4aae4
    set_read_p(true);
Packit a4aae4
Packit a4aae4
    return true;
Packit a4aae4
}
Packit a4aae4
Packit a4aae4
void
Packit a4aae4
Float64::print_val(FILE *out, string space, bool print_decl_p)
Packit a4aae4
{
Packit a4aae4
    ostringstream oss;
Packit a4aae4
    print_val(oss, space, print_decl_p);
Packit a4aae4
    fwrite(oss.str().data(), sizeof(char), oss.str().length(), out);
Packit a4aae4
}
Packit a4aae4
Packit a4aae4
void
Packit a4aae4
Float64::print_val(ostream &out, string space, bool print_decl_p)
Packit a4aae4
{
Packit a4aae4
    // Set the precision to 15 digits
Packit a4aae4
    std::streamsize prec = out.precision(15);
Packit a4aae4
Packit a4aae4
    if (print_decl_p) {
Packit a4aae4
        print_decl(out, space, false);
Packit a4aae4
        out << " = " << d_buf << ";\n";
Packit a4aae4
    }
Packit a4aae4
    else
Packit a4aae4
        out << d_buf;
Packit a4aae4
Packit a4aae4
    // reset the precision
Packit a4aae4
    out.precision(prec);
Packit a4aae4
}
Packit a4aae4
Packit a4aae4
bool
Packit a4aae4
Float64::ops(BaseType *b, int op)
Packit a4aae4
{
Packit a4aae4
    // Extract the Byte arg's value.
Packit a4aae4
    if (!read_p() && !read()) {
Packit a4aae4
        // Jose Garcia
Packit a4aae4
        // Since the read method is virtual and implemented outside
Packit a4aae4
        // libdap++ if we cannot read the data that is the problem
Packit a4aae4
        // of the user or of whoever wrote the surrogate library
Packit a4aae4
        // implemeting read therefore it is an internal error.
Packit a4aae4
        throw InternalErr(__FILE__, __LINE__, "This value not read!");
Packit a4aae4
    }
Packit a4aae4
Packit a4aae4
    // Extract the second arg's value.
Packit a4aae4
    if (!b->read_p() && !b->read()) {
Packit a4aae4
        // Jose Garcia
Packit a4aae4
        // Since the read method is virtual and implemented outside
Packit a4aae4
        // libdap++ if we cannot read the data that is the problem
Packit a4aae4
        // of the user or of whoever wrote the surrogate library
Packit a4aae4
        // implemeting read therefore it is an internal error.
Packit a4aae4
        throw InternalErr(__FILE__, __LINE__, "This value not read!");
Packit a4aae4
    }
Packit a4aae4
Packit a4aae4
    return d4_ops(b, op);
Packit a4aae4
}
Packit a4aae4
Packit a4aae4
bool
Packit a4aae4
Float64::d4_ops(BaseType *b, int op)
Packit a4aae4
{
Packit a4aae4
    DBG(cerr << "b->typename(): " << b->type_name() << endl);
Packit a4aae4
Packit a4aae4
    switch (b->type()) {
Packit a4aae4
        case dods_int8_c:
Packit a4aae4
            return Cmp<dods_float64, dods_int8>(op, d_buf, static_cast<Int8*>(b)->value());
Packit a4aae4
        case dods_byte_c:
Packit a4aae4
            return SUCmp<dods_float64, dods_byte>(op, d_buf, static_cast<Byte*>(b)->value());
Packit a4aae4
        case dods_int16_c:
Packit a4aae4
            return Cmp<dods_float64, dods_int16>(op, d_buf, static_cast<Int16*>(b)->value());
Packit a4aae4
        case dods_uint16_c:
Packit a4aae4
            return SUCmp<dods_float64, dods_uint16>(op, d_buf, static_cast<UInt16*>(b)->value());
Packit a4aae4
        case dods_int32_c:
Packit a4aae4
            return Cmp<dods_float64, dods_int32>(op, d_buf, static_cast<Int32*>(b)->value());
Packit a4aae4
        case dods_uint32_c:
Packit a4aae4
            return SUCmp<dods_float64, dods_uint32>(op, d_buf, static_cast<UInt32*>(b)->value());
Packit a4aae4
        case dods_int64_c:
Packit a4aae4
            return Cmp<dods_float64, dods_int64>(op, d_buf, static_cast<Int64*>(b)->value());
Packit a4aae4
        case dods_uint64_c:
Packit a4aae4
            return SUCmp<dods_float64, dods_uint64>(op, d_buf, static_cast<UInt64*>(b)->value());
Packit a4aae4
        case dods_float32_c:
Packit a4aae4
            // Note that this code casts the double (dods_float64) to a float because when
Packit a4aae4
            // real numbers are approximated using float or double, errors are larger in the float
Packit a4aae4
            // case, making <= and == operators fail. By casting to the smaller type, the
Packit a4aae4
            // same values have the same error and we can avoid using a range compare and a
Packit a4aae4
            // delta value.
Packit a4aae4
            return Cmp<dods_float32, dods_float32>(op, (float)d_buf, static_cast<Float32*>(b)->value());
Packit a4aae4
        case dods_float64_c:
Packit a4aae4
            return Cmp<dods_float64, dods_float64>(op, d_buf, static_cast<Float64*>(b)->value());
Packit a4aae4
        case dods_str_c:
Packit a4aae4
        case dods_url_c:
Packit a4aae4
            throw Error(malformed_expr, "Relational operators can only compare compatible types (number, string).");
Packit a4aae4
        default:
Packit a4aae4
            throw Error(malformed_expr, "Relational operators only work with scalar types.");
Packit a4aae4
    }
Packit a4aae4
Packit a4aae4
}
Packit a4aae4
/** @brief dumps information about this object
Packit a4aae4
 *
Packit a4aae4
 * Displays the pointer value of this instance and information about this
Packit a4aae4
 * instance.
Packit a4aae4
 *
Packit a4aae4
 * @param strm C++ i/o stream to dump the information to
Packit a4aae4
 * @return void
Packit a4aae4
 */
Packit a4aae4
void
Packit a4aae4
Float64::dump(ostream &strm) const
Packit a4aae4
{
Packit a4aae4
    strm << DapIndent::LMarg << "Float64::dump - ("
Packit a4aae4
    << (void *)this << ")" << endl ;
Packit a4aae4
    DapIndent::Indent() ;
Packit a4aae4
    BaseType::dump(strm) ;
Packit a4aae4
    strm << DapIndent::LMarg << "value: " << d_buf << endl ;
Packit a4aae4
    DapIndent::UnIndent() ;
Packit a4aae4
}
Packit a4aae4
Packit a4aae4
} // namespace libdap
Packit a4aae4