gem5/base/statistics.cc
Nathan Binkert 55d94ba2e8 Add python output support to the statistics package!
base/statistics.cc:
base/statistics.hh:
    -  add python output support to the statistics package
    -  each statistic type has a python() member function that takes
    a Python object to which the stat will output it's python
    representation
    -  add getStatData hack so that the StatData pointer can be looked
    up by the proxies with their opaque pointer to the stat they're
    proxying for.  This is necessary because the proxy really proxies
    for the bin and not the stat.  Be nice to figure out how to get
    rid of it.  The hack is used so that the str() function of a
    proxy can properly name itself.
    -  To print formula stats, every stat has a str() function that
    converts that stat to a string that python can execute to get
    a value.
test/Makefile:
    add python stuff
test/stattest.cc:
    add more tests and test python support

--HG--
extra : convert_revision : 513814ab0a125606897f2c57dccdf22879032ef9
2003-12-24 03:25:36 -05:00

1263 lines
28 KiB
C++

/*
* Copyright (c) 2003 The Regents of The University of Michigan
* All rights reserved.
*
* Redistribution and use in source and binary forms, with or without
* modification, are permitted provided that the following conditions are
* met: redistributions of source code must retain the above copyright
* notice, this list of conditions and the following disclaimer;
* redistributions in binary form must reproduce the above copyright
* notice, this list of conditions and the following disclaimer in the
* documentation and/or other materials provided with the distribution;
* neither the name of the copyright holders nor the names of its
* contributors may be used to endorse or promote products derived from
* this software without specific prior written permission.
*
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
* A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
* OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
* SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
* LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
* DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
* THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
* (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
* OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*/
#include <iomanip>
#include <iostream>
#include <list>
#include <map>
#include <string>
#include <sstream>
#include "base/callback.hh"
#include "base/cprintf.hh"
#include "base/hostinfo.hh"
#include "base/misc.hh"
#include "base/python.hh"
#include "base/statistics.hh"
#include "base/str.hh"
#include "base/time.hh"
#include "base/trace.hh"
#ifdef __M5_NAN
float
__nan()
{
union {
uint32_t ui;
float f;
} nan;
nan.ui = 0x7fc00000;
return nan.f;
}
#endif
#ifdef STAT_DEBUG
static int total_stats = 0;
#endif
using namespace std;
// This is a hack to get this parameter from the old stats package.
namespace Statistics {
bool PrintDescriptions = true;
DisplayMode DefaultMode = mode_simplescalar;
namespace Database
{
class Data
{
private:
typedef list<StatData *> list_t;
typedef map<void *, StatData *> map_t;
list<MainBin *> bins;
list_t allStats;
list_t printStats;
map_t statMap;
public:
void dump(ostream &stream, const string &name, DisplayMode mode);
void display(ostream &stream, DisplayMode mode);
void python(ostream &stream, const string &name);
void python(Python &py, const string &name, const string &bin);
StatData *find(void *stat);
void mapStat(void *stat, StatData *data);
void check();
void reset();
void regBin(MainBin *bin, string name);
void regPrint(void *stat);
static std::string name() { return "Statistics Database"; }
};
void
Data::dump(ostream &stream, const string &name, DisplayMode mode)
{
MainBin *orig = MainBin::curBin();
switch (mode) {
case mode_python:
python(stream, name);
break;
case mode_m5:
case mode_simplescalar:
display(stream, mode);
break;
default:
warn("invalid display mode!\n");
display(stream, mode_m5);
break;
}
if (orig)
orig->activate();
}
void
Data::display(ostream &stream, DisplayMode mode)
{
if (!bins.empty()) {
list<MainBin *>::iterator i = bins.begin();
list<MainBin *>::iterator bins_end = bins.end();
ccprintf(stream, "PRINTING BINNED STATS\n");
while (i != bins_end) {
(*i)->activate();
ccprintf(stream,"---%s Bin------------\n", (*i)->name());
list_t::iterator j = printStats.begin();
list_t::iterator end = printStats.end();
while (j != end) {
StatData *stat = *j;
if (stat->dodisplay())
stat->display(stream, mode);
++j;
}
++i;
ccprintf(stream, "---------------------------------\n");
}
} else {
list_t::iterator i = printStats.begin();
list_t::iterator end = printStats.end();
while (i != end) {
StatData *stat = *i;
if (stat->dodisplay() && !stat->binned())
stat->display(stream, mode);
++i;
}
}
}
void
Data::python(ostream &stream, const string &name)
{
Python py(stream);
ccprintf(stream, "import sys\n");
ccprintf(stream, "sys.path.append('.')\n");
ccprintf(stream, "from m5stats import *\n\n");
if (bins.empty()) {
python(py, name, "");
} else {
list<MainBin *>::iterator i = bins.begin();
list<MainBin *>::iterator end = bins.end();
while (i != end) {
(*i)->activate();
python(py, name, (*i)->name());
++i;
}
}
py.next();
ccprintf(stream, "if __name__ == '__main__':\n");
ccprintf(stream, " program_display()\n");
}
void
Data::python(Python &py, const string &name, const string &bin)
{
py.start("collections.append");
py.start("Collection");
py.qarg(name);
py.qarg(bin);
py.qarg(hostname());
py.qarg(Time::start.date());
list_t::iterator i = allStats.begin();
list_t::iterator end = allStats.end();
while (i != end) {
StatData *stat = *i;
stat->python(py);
++i;
}
py.end();
py.end();
}
StatData *
Data::find(void *stat)
{
map_t::const_iterator i = statMap.find(stat);
if (i == statMap.end())
return NULL;
return (*i).second;
}
void
Data::check()
{
list_t::iterator i = allStats.begin();
list_t::iterator end = allStats.end();
while (i != end) {
StatData *data = *i;
assert(data);
data->check();
++i;
}
i = allStats.begin();
int j = 0;
while (i != end) {
StatData *data = *i;
if (!(data->flags & print))
data->name = "__Stat" + to_string(j++);
++i;
}
}
void
Data::reset()
{
// reset non-binned stats
list_t::iterator i = allStats.begin();
list_t::iterator end = allStats.end();
while (i != end) {
StatData *data = *i;
if (!data->binned())
data->reset();
++i;
}
// save the bin so we can go back to where we were
MainBin *orig = MainBin::curBin();
// reset binned stats
list<MainBin *>::iterator bi = bins.begin();
list<MainBin *>::iterator be = bins.end();
while (bi != be) {
MainBin *bin = *bi;
bin->activate();
i = allStats.begin();
while (i != end) {
StatData *data = *i;
if (data->binned())
data->reset();
++i;
}
++bi;
}
// restore bin
MainBin::curBin() = orig;
}
void
Data::mapStat(void *stat, StatData *data)
{
if (statMap.find(stat) != statMap.end())
panic("shouldn't register stat twice!");
allStats.push_back(data);
#ifndef NDEBUG
bool success =
#endif
(statMap.insert(make_pair(stat, data))).second;
assert(statMap.find(stat) != statMap.end());
assert(success && "this should never fail");
}
void
Data::regBin(MainBin *bin, string _name)
{
bins.push_back(bin);
DPRINTF(Stats, "registering %s\n", _name);
}
void
Data::regPrint(void *stat)
{
StatData *data = find(stat);
if (data->flags & print)
return;
data->flags |= print;
list_t::iterator j = printStats.insert(printStats.end(), data);
inplace_merge(printStats.begin(), j, printStats.end(), StatData::less);
}
Data &
StatDB()
{
static Data db;
return db;
}
}
StatData *
DataAccess::find() const
{
return Database::StatDB().find(const_cast<void *>((const void *)this));
}
const StatData *
getStatData(const void *stat)
{
return Database::StatDB().find(const_cast<void *>(stat));
}
void
DataAccess::map(StatData *data)
{
Database::StatDB().mapStat(this, data);
}
StatData *
DataAccess::statData()
{
StatData *ptr = find();
assert(ptr);
return ptr;
}
const StatData *
DataAccess::statData() const
{
const StatData *ptr = find();
assert(ptr);
return ptr;
}
void
DataAccess::setInit()
{
statData()->flags |= init;
}
void
DataAccess::setPrint()
{
Database::StatDB().regPrint(this);
}
StatData::~StatData()
{
}
bool
StatData::less(StatData *stat1, StatData *stat2)
{
const string &name1 = stat1->name;
const string &name2 = stat2->name;
vector<string> v1;
vector<string> v2;
tokenize(v1, name1, '.');
tokenize(v2, name2, '.');
int last = min(v1.size(), v2.size()) - 1;
for (int i = 0; i < last; ++i)
if (v1[i] != v2[i])
return v1[i] < v2[i];
// Special compare for last element.
if (v1[last] == v2[last])
return v1.size() < v2.size();
else
return v1[last] < v2[last];
return false;
}
bool
StatData::baseCheck() const
{
if (!(flags & init)) {
#ifdef STAT_DEBUG
cprintf("this is stat number %d\n",(*i)->number);
#endif
panic("Not all stats have been initialized");
return false;
}
if ((flags & print) && name.empty()) {
panic("all printable stats must be named");
return false;
}
return true;
}
string
ValueToString(result_t value, DisplayMode mode, int precision)
{
stringstream val;
if (!isnan(value)) {
if (precision != -1)
val.precision(precision);
else if (value == rint(value))
val.precision(0);
val.unsetf(ios::showpoint);
val.setf(ios::fixed);
val << value;
} else {
val << (mode == mode_m5 ? "no value" : "<err: div-0>");
}
return val.str();
}
struct ScalarPrint
{
result_t value;
string name;
string desc;
StatFlags flags;
DisplayMode mode;
int precision;
result_t pdf;
result_t cdf;
ScalarPrint()
: value(0.0), flags(0), mode(DefaultMode), precision(0),
pdf(NAN), cdf(NAN)
{}
void operator()(ostream &stream) const;
};
void
ScalarPrint::operator()(ostream &stream) const
{
if (flags & nozero && value == 0.0 ||
flags & nonan && isnan(value))
return;
stringstream pdfstr, cdfstr;
if (!isnan(pdf))
ccprintf(pdfstr, "%.2f%%", pdf * 100.0);
if (!isnan(cdf))
ccprintf(cdfstr, "%.2f%%", cdf * 100.0);
if (mode == mode_simplescalar && flags & __substat) {
ccprintf(stream, "%32s %12s %10s %10s", name,
ValueToString(value, mode, precision), pdfstr, cdfstr);
} else {
ccprintf(stream, "%-40s %12s %10s %10s", name,
ValueToString(value, mode, precision), pdfstr, cdfstr);
}
if (PrintDescriptions) {
if (!desc.empty())
ccprintf(stream, " # %s", desc);
}
stream << endl;
}
struct VectorPrint
{
string name;
string desc;
vector<string> subnames;
vector<string> subdescs;
StatFlags flags;
DisplayMode mode;
int precision;
rvec_t vec;
result_t total;
VectorPrint()
: subnames(0), subdescs(0), flags(0), mode(DefaultMode),
precision(-1), total(NAN)
{}
void operator()(ostream &stream) const;
};
void
VectorPrint::operator()(std::ostream &stream) const
{
int _size = vec.size();
result_t _total = 0.0;
if (flags & (pdf | cdf)) {
for (int i = 0; i < _size; ++i) {
_total += vec[i];
}
}
string base = name + ((mode == mode_simplescalar) ? "_" : "::");
ScalarPrint print;
print.name = name;
print.desc = desc;
print.precision = precision;
print.flags = flags;
bool havesub = !subnames.empty();
if (_size == 1) {
print.value = vec[0];
print(stream);
} else if (mode == mode_m5) {
for (int i = 0; i < _size; ++i) {
if (havesub && (i >= subnames.size() || subnames[i].empty()))
continue;
print.name = base + (havesub ? subnames[i] : to_string(i));
print.desc = subdescs.empty() ? desc : subdescs[i];
print.value = vec[i];
if (_total && (flags & pdf)) {
print.pdf = vec[i] / _total;
print.cdf += print.pdf;
}
print(stream);
}
if (flags & ::Statistics::total) {
print.name = base + "total";
print.desc = desc;
print.value = total;
print(stream);
}
} else {
if (flags & ::Statistics::total) {
print.value = total;
print(stream);
}
result_t _pdf = 0.0;
result_t _cdf = 0.0;
if (flags & dist) {
ccprintf(stream, "%s.start_dist\n", name);
for (int i = 0; i < _size; ++i) {
print.name = havesub ? subnames[i] : to_string(i);
print.desc = subdescs.empty() ? desc : subdescs[i];
print.flags |= __substat;
print.value = vec[i];
if (_total) {
_pdf = vec[i] / _total;
_cdf += _pdf;
}
if (flags & pdf)
print.pdf = _pdf;
if (flags & cdf)
print.cdf = _cdf;
print(stream);
}
ccprintf(stream, "%s.end_dist\n", name);
} else {
for (int i = 0; i < _size; ++i) {
if (havesub && subnames[i].empty())
continue;
print.name = base;
print.name += havesub ? subnames[i] : to_string(i);
print.desc = subdescs.empty() ? desc : subdescs[i];
print.value = vec[i];
if (_total) {
_pdf = vec[i] / _total;
_cdf += _pdf;
} else {
_pdf = _cdf = NAN;
}
if (flags & pdf) {
print.pdf = _pdf;
print.cdf = _cdf;
}
print(stream);
}
}
}
}
struct DistPrint
{
string name;
string desc;
StatFlags flags;
DisplayMode mode;
int precision;
result_t min_val;
result_t max_val;
result_t underflow;
result_t overflow;
rvec_t vec;
result_t sum;
result_t squares;
result_t samples;
int min;
int max;
int bucket_size;
int size;
bool fancy;
void operator()(ostream &stream) const;
};
void
DistPrint::operator()(ostream &stream) const
{
if (fancy) {
ScalarPrint print;
string base = name + ((mode == mode_m5) ? "::" : "_");
print.precision = precision;
print.flags = flags;
print.mode = mode;
print.desc = desc;
print.name = base + "mean";
print.value = samples ? sum / samples : NAN;
print(stream);
print.name = base + "stdev";
print.value = samples ? sqrt((samples * squares - sum * sum) /
(samples * (samples - 1.0))) : NAN;
print(stream);
print.name = "**Ignore: " + base + "TOT";
print.value = samples;
print(stream);
return;
}
assert(size == vec.size());
result_t total = 0.0;
total += underflow;
for (int i = 0; i < size; ++i)
total += vec[i];
total += overflow;
string base = name + (mode == mode_m5 ? "::" : ".");
ScalarPrint print;
print.desc = (mode == mode_m5) ? desc : "";
print.flags = flags;
print.mode = mode;
print.precision = precision;
if (mode == mode_simplescalar) {
ccprintf(stream, "%-42s", base + "start_dist");
if (PrintDescriptions && !desc.empty())
ccprintf(stream, " # %s", desc);
stream << endl;
}
print.name = base + "samples";
print.value = samples;
print(stream);
print.name = base + "min_value";
print.value = min_val;
print(stream);
if (mode == mode_m5 || underflow > 0.0) {
print.name = base + "underflows";
print.value = underflow;
if (mode == mode_m5 && total) {
print.pdf = underflow / total;
print.cdf += print.pdf;
}
print(stream);
}
if (mode == mode_m5) {
for (int i = 0; i < size; ++i) {
stringstream namestr;
namestr << name;
int low = i * bucket_size + min;
int high = ::min((i + 1) * bucket_size + min - 1, max);
namestr << low;
if (low < high)
namestr << "-" << high;
print.name = namestr.str();
print.value = vec[i];
if (total) {
print.pdf = vec[i] / total;
print.cdf += print.pdf;
}
print(stream);
}
} else {
int _min;
result_t _pdf;
result_t _cdf = 0.0;
print.flags = flags | __substat;
for (int i = 0; i < size; ++i) {
if (flags & nozero && vec[i] == 0.0 ||
flags & nonan && isnan(vec[i]))
continue;
_min = i * bucket_size + min;
_pdf = vec[i] / total * 100.0;
_cdf += _pdf;
print.name = ValueToString(_min, mode, 0);
print.value = vec[i];
print.pdf = (flags & pdf) ? _pdf : NAN;
print.cdf = (flags & cdf) ? _cdf : NAN;
print(stream);
}
print.flags = flags;
}
if (mode == mode_m5 || overflow > 0.0) {
print.name = base + "overflows";
print.value = overflow;
if (mode == mode_m5 && total) {
print.pdf = overflow / total;
print.cdf += print.pdf;
} else {
print.pdf = NAN;
print.cdf = NAN;
}
print(stream);
}
print.pdf = NAN;
print.cdf = NAN;
if (mode != mode_simplescalar) {
print.name = base + "total";
print.value = total;
print(stream);
}
print.name = base + "max_value";
print.value = max_val;
print(stream);
if (mode != mode_simplescalar && samples != 0) {
print.name = base + "mean";
print.value = sum / samples;
print(stream);
print.name = base + "stdev";
print.value = sqrt((samples * squares - sum * sum) /
(samples * (samples - 1.0)));
print(stream);
}
if (mode == mode_simplescalar)
ccprintf(stream, "%send_dist\n\n", base);
}
void
ScalarDataBase::display(ostream &stream, DisplayMode mode) const
{
ScalarPrint print;
print.value = val();
print.name = name;
print.desc = desc;
print.flags = flags;
print.mode = mode;
print.precision = precision;
print(stream);
}
void
VectorDataBase::display(ostream &stream, DisplayMode mode) const
{
int size = this->size();
const_cast<VectorDataBase *>(this)->update();
VectorPrint print;
print.name = name;
print.desc = desc;
print.flags = flags;
print.mode = mode;
print.precision = precision;
print.vec = val();
print.total = total();
if (!subnames.empty()) {
for (int i = 0; i < size; ++i) {
if (!subnames[i].empty()) {
print.subnames = subnames;
print.subnames.resize(size);
for (int i = 0; i < size; ++i) {
if (!subnames[i].empty() && !subdescs[i].empty()) {
print.subdescs = subdescs;
print.subdescs.resize(size);
break;
}
}
break;
}
}
}
print(stream);
}
void
Vector2dDataBase::display(ostream &stream, DisplayMode mode) const
{
const_cast<Vector2dDataBase *>(this)->update();
bool havesub = false;
VectorPrint print;
print.subnames = y_subnames;
print.flags = flags;
print.mode = mode;
print.precision = precision;
if (!subnames.empty()) {
for (int i = 0; i < x; ++i)
if (!subnames[i].empty())
havesub = true;
}
rvec_t tot_vec(y);
result_t super_total = 0.0;
for (int i = 0; i < x; ++i) {
if (havesub && (i >= subnames.size() || subnames[i].empty()))
continue;
int iy = i * y;
rvec_t yvec(y);
result_t total = 0.0;
for (int j = 0; j < y; ++j) {
yvec[j] = vec[iy + j];
tot_vec[j] += yvec[j];
total += yvec[j];
super_total += yvec[j];
}
print.name = name + "_" + (havesub ? subnames[i] : to_string(i));
print.desc = desc;
print.vec = yvec;
print.total = total;
print(stream);
}
if ((flags & ::Statistics::total) && (x > 1)) {
print.name = name;
print.desc = desc;
print.vec = tot_vec;
print.total = super_total;
print(stream);
}
}
void
DistDataBase::display(ostream &stream, DisplayMode mode) const
{
const_cast<DistDataBase *>(this)->update();
DistPrint print;
print.name = name;
print.desc = desc;
print.flags = flags;
print.mode = mode;
print.precision = precision;
print.min_val = data.min_val;
print.max_val = data.max_val;
print.underflow = data.underflow;
print.overflow = data.overflow;
print.vec = data.vec;
print.sum = data.sum;
print.squares = data.squares;
print.samples = data.samples;
print.min = data.min;
print.max = data.max;
print.bucket_size = data.bucket_size;
print.size = data.size;
print.fancy = data.fancy;
print(stream);
}
void
VectorDistDataBase::display(ostream &stream, DisplayMode mode) const
{
const_cast<VectorDistDataBase *>(this)->update();
for (int i = 0; i < size(); ++i) {
DistPrint print;
print.name = name +
(subnames[i].empty() ? ("_" + to_string(i)) : subnames[i]);
print.desc = subdescs[i].empty() ? desc : subdescs[i];
print.flags = flags;
print.mode = mode;
print.precision = precision;
print.min_val = data[i].min_val;
print.max_val = data[i].max_val;
print.underflow = data[i].underflow;
print.overflow = data[i].overflow;
print.vec = data[i].vec;
print.sum = data[i].sum;
print.squares = data[i].squares;
print.samples = data[i].samples;
print.min = data[i].min;
print.max = data[i].max;
print.bucket_size = data[i].bucket_size;
print.size = data[i].size;
print.fancy = data[i].fancy;
print(stream);
}
}
void
ScalarDataBase::python(Python &py) const
{
py.start("Scalar");
py.qarg(name);
py.qarg(desc);
py.kwarg("binned", binned());
py.kwarg("precision", precision);
py.kwarg("flags", flags);
if (prereq)
py.qkwarg("prereq", prereq->name);
py.kwarg("value", val());
py.end();
}
void
VectorDataBase::python(Python &py) const
{
const_cast<VectorDataBase *>(this)->update();
py.start("Vector");
py.qarg(name);
py.qarg(desc);
py.kwarg("binned", binned());
py.kwarg("precision", precision);
py.kwarg("flags", flags);
if (prereq)
py.qkwarg("prereq", prereq->name);
py.kwarg("value", val());
if (!subnames.empty())
py.qkwarg("subnames", subnames);
if (!subdescs.empty())
py.qkwarg("subdescs", subdescs);
py.end();
}
void
DistDataData::python(Python &py, const string &name) const
{
string s = name.empty() ? "" : name + "=";
if (samples == 0 || fancy)
s += "SimpleDist";
else
s += "FullDist";
py.start(s);
py.arg(sum);
py.arg(squares);
py.arg(samples);
if (samples && !fancy) {
py.arg(min_val);
py.arg(min_val);
py.arg(underflow);
py.arg(vec);
py.arg(overflow);
py.arg(min);
py.arg(max);
py.arg(bucket_size);
py.arg(size);
}
py.end();
}
void
FormulaDataBase::python(Python &py) const
{
const_cast<FormulaDataBase *>(this)->update();
py.start("Formula");
py.qarg(name);
py.qarg(desc);
py.kwarg("binned", binned());
py.kwarg("precision", precision);
py.kwarg("flags", flags);
if (prereq)
py.qkwarg("prereq", prereq->name);
py.qkwarg("formula", str());
if (!subnames.empty())
py.qkwarg("subnames", subnames);
if (!subdescs.empty())
py.qkwarg("subdescs", subdescs);
py.end();
}
void
DistDataBase::python(Python &py) const
{
const_cast<DistDataBase *>(this)->update();
py.start("Dist");
py.qarg(name);
py.qarg(desc);
py.kwarg("binned", binned());
py.kwarg("precision", precision);
py.kwarg("flags", flags);
if (prereq)
py.qkwarg("prereq", prereq->name);
data.python(py, "dist");
py.end();
}
void
VectorDistDataBase::python(Python &py) const
{
const_cast<VectorDistDataBase *>(this)->update();
py.start("VectorDist");
py.qarg(name);
py.qarg(desc);
py.kwarg("binned", binned());
py.kwarg("precision", precision);
py.kwarg("flags", flags);
if (prereq)
py.qkwarg("prereq", prereq->name);
if (!subnames.empty())
py.qkwarg("subnames", subnames);
if (!subdescs.empty())
py.qkwarg("subdescs", subdescs);
py.tuple("dist");
typedef std::vector<DistDataData>::const_iterator iter;
iter i = data.begin();
iter end = data.end();
while (i != end) {
i->python(py, "");
++i;
}
py.endTuple();
py.end();
}
void
Vector2dDataBase::python(Python &py) const
{
const_cast<Vector2dDataBase *>(this)->update();
py.start("Vector2d");
py.qarg(name);
py.qarg(desc);
py.kwarg("binned", binned());
py.kwarg("precision", precision);
py.kwarg("flags", flags);
if (prereq)
py.qkwarg("prereq", prereq->name);
py.kwarg("value", vec);
if (!subnames.empty())
py.qkwarg("subnames", subnames);
if (!subdescs.empty())
py.qkwarg("subdescs", subdescs);
if (!y_subnames.empty())
py.qkwarg("ysubnames", y_subnames);
py.kwarg("x", x);
py.kwarg("y", y);
py.end();
}
void
FormulaBase::val(rvec_t &vec) const
{
vec = root->val();
}
result_t
FormulaBase::total() const
{
return root->total();
}
size_t
FormulaBase::size() const
{
if (!root)
return 0;
else
return root->size();
}
bool
FormulaBase::binned() const
{
return root && root->binned();
}
void
FormulaBase::reset()
{
}
bool
FormulaBase::zero() const
{
rvec_t vec;
val(vec);
for (int i = 0; i < vec.size(); ++i)
if (vec[i] != 0.0)
return false;
return true;
}
void
FormulaBase::update(StatData *)
{
}
string
FormulaBase::str() const
{
return root ? root->str() : "";
}
Formula::Formula()
{
setInit();
}
Formula::Formula(Temp r)
{
root = r;
assert(size());
}
const Formula &
Formula::operator=(Temp r)
{
assert(!root && "Can't change formulas");
root = r;
assert(size());
return *this;
}
const Formula &
Formula::operator+=(Temp r)
{
if (root)
root = NodePtr(new BinaryNode<std::plus<result_t> >(root, r));
else
root = r;
assert(size());
return *this;
}
MainBin::MainBin(const string &name)
: _name(name), mem(NULL), memsize(-1)
{
Database::StatDB().regBin(this, name);
}
MainBin::~MainBin()
{
if (mem)
delete [] mem;
}
char *
MainBin::memory(off_t off)
{
if (memsize == -1)
memsize = CeilPow2((size_t) offset());
if (!mem) {
mem = new char[memsize];
memset(mem, 0, memsize);
}
assert(offset() <= size());
return mem + off;
}
void
check()
{
Database::StatDB().check();
}
void
dump(ostream &stream, const string &name, DisplayMode mode)
{
Database::StatDB().dump(stream, name, mode);
}
CallbackQueue resetQueue;
void
registerResetCallback(Callback *cb)
{
resetQueue.add(cb);
}
void
reset()
{
Database::StatDB().reset();
resetQueue.process();
}
} // namespace Statistics