opt-param.cc

Go to the documentation of this file.

/////////////////////////////////////////////////////////////////////////////
//! \file opt-param.cc  Optimization parameters
//
// Copyright (c) 2000-2004 Petr Peringer
//
// This library is licensed under GNU Library GPL. See the file COPYING.
//
 
// EXPERIMENTAL
// parameters for optimization purposes
 
#include "simlib.h"
#include "internal.h"
#include "optimize.h"
#include <cstring> // strcmp()
 
namespace simlib3 {
 
SIMLIB_IMPLEMENTATION;
 
void Param::Print() const
{
    ::Print("#Parameter %s: value=%g (min=%g, max=%g)\n", name, value, min,
            max);
}
 
#if 0
Param & Param::operator = (double x) {
    if (x > max)
        x = max;
    if (x < min)
        x = min;
    value = x;
}
#endif
 
////////////////////////////////////////////////////////////////////////////
//
//
int ParameterVector::search(const char *name)
{
    for (int i = 0; i < n; i++)
        if (std::strcmp(p[i].Name(), name) == 0)
            return i;
    return -1;
}
 
//  ParameterVector(): n(0), p(0) {}
//  void Add(const Param&p) { add to the end }
//    pp.Add(Param("x",1,5));
 
// initialize by array
ParameterVector::ParameterVector(int sz, Param * a):
n(sz), p(new Param[n])
{
    for (int i = 0; i < n; i++)
        p[i] = a[i];
}
 
// copy constructor
ParameterVector::ParameterVector(const ParameterVector & a):n(a.n),
p(new Param[n])
{
    for (int i = 0; i < n; i++)
        p[i] = a[i];
}
 
// assignment
ParameterVector & ParameterVector::operator = (const ParameterVector & a) {
    if (this == &a)
        return *this;           // a=a
    n = a.n;
    delete[] p;                 // TODO add exception-safety
    p = new Param[n];
    for (int i = 0; i < n; i++)
        p[i] = a[i];            // copy objects
    return *this;
}
 
ParameterVector::~ParameterVector()
{
    delete[]p;
}
 
// PROTOTYPE, values only !!!!!!!#########
bool operator == (const ParameterVector & p1, const ParameterVector & p2) {
    int n = p1.size();
    if (n != p2.size())
        return false;
    for (int i = 0; i < n; i++)
        if (p1[i].Value() != p2[i].Value())
            return false;
    return true;
}
 
// print
void ParameterVector::PrintValues() const
{
    for (int i = 0; i < n; i++)
        ::Print("%g ", p[i].Value());
}
 
void ParameterVector::Print() const
{
    for (int i = 0; i < n; i++)
        p[i].Print();
}
 
}
// end