density.cpp

Go to the documentation of this file.

/*
 * This work is licensed under a Creative Commons 
 * Attribution-Noncommercial-Share Alike 3.0 United States License.
 * 
 *    http://creativecommons.org/licenses/by-nc-sa/3.0/us/
 * 
 * You are free:
 * 
 *    to Share - to copy, distribute, display, and perform the work
 *    to Remix - to make derivative works
 * 
 * Under the following conditions:
 * 
 *    Attribution. You must attribute the work in the manner specified by the
 *    author or licensor (but not in any way that suggests that they endorse you
 *    or your use of the work).
 * 
 *    Noncommercial. You may not use this work for commercial purposes.
 * 
 *    Share Alike. If you alter, transform, or build upon this work, you may
 *    distribute the resulting work only under the same or similar license to
 *    this one.
 * 
 * For any reuse or distribution, you must make clear to others the license
 * terms of this work. The best way to do this is by including this header.
 * 
 * Any of the above conditions can be waived if you get permission from the
 * copyright holder.
 * 
 * Apart from the remix rights granted under this license, nothing in this
 * license impairs or restricts the author's moral rights.
 */


#include <config.h>

#include <iostream>

#include <jwsc++/base/exception.h>

#include "density.h"


using std::ostream;
using std::istream;
using jwscxx::base::Arg_error;


Density::Uniform_params::Uniform_params(float min, float max) throw (Arg_error)
{
    set(min, max);
}


void Density::Uniform_params::set(float min, float max) throw (Arg_error)
{
    if (min > max)
    {
        throw Arg_error("Minimum larger than maximum");
    }

    this->min = min;
    this->max = max;
}


ostream& Density::operator<<
(
     ostream&                       out, 
     const Density::Uniform_params& params
)
{
    return out << params.min << ' ' << params.max;
}


istream& Density::operator>>
(
     istream&                 in, 
     Density::Uniform_params& params
)
{
    in >> params.min;
    in >> params.max;
    return in;
}


Density::Gaussian_params::Gaussian_params
(
    float mu, 
    float sigma, 
    float min, 
    float max
) 
throw (Arg_error)
{
    set(mu, sigma, min, max);
}


void Density::Gaussian_params::set
(
    float mu, 
    float sigma, 
    float min, 
    float max
) 
throw (Arg_error)
{
    if (min > max)
    {
        throw Arg_error("Minimum larger than maximum");
    }
    if (mu < min || mu > max)
    {
        throw Arg_error("Mean not in min/max range");
    }
    if (sigma < 1.0e-8)
    {
        throw Arg_error("Sigma too small");
    }

    this->mu    = mu;
    this->sigma = sigma;
    this->min   = min;
    this->max   = max;
}


ostream& Density::operator<<
(
     ostream&                        out, 
     const Density::Gaussian_params& params
)
{
    return out << params.mu << ' ' << params.sigma << ' '
               << params.min << ' ' << params.max;
}


istream& Density::operator>>
(
     istream&                  in, 
     Density::Gaussian_params& params
)
{
    in >> params.mu;
    in >> params.sigma;
    in >> params.min;
    in >> params.max;
    return in;
}