gem5/base/statistics.cc
Nathan Binkert 9dcfdce732 Improve python output support
fix some bugs

base/statistics.cc:
    -  Move the python dump stream stuff into the statistics package
    because it needs more control over the file when there
    are multiple dumps
    -  add a subname to the output
    -  make all of the stats a list instead of variable arguments to
    get around the 255 argument limit
    -  the description needs to be triple quoted
    -  avoid errors for formulas that don't have any root element
    set up.
base/statistics.hh:
    -  get rid of mode_python and just separate python dumping from
    regular stats dumping
    -  fix a huge bug that made a Formula use a VectorData instead
    of a FormulaData

--HG--
extra : convert_revision : 7303cff3ccdcc3d306ab17375219fc7fecac7e5e
2003-12-28 12:48:13 -05:00

1308 lines
29 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 <fstream>
#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;
ofstream *stream;
Python *py;
public:
Data();
~Data();
void dump(ostream &stream, DisplayMode mode);
void display(ostream &stream, DisplayMode mode);
void python_start(const string &file);
void python_dump(const string &name, const string &subname);
void python(const string &name, const string &subname,
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"; }
};
Data::Data()
: stream(0), py(0)
{
}
Data::~Data()
{
if (stream) {
delete py;
ccprintf(*stream, "if __name__ == '__main__':\n");
ccprintf(*stream, " program_display()\n");
stream->close();
delete stream;
}
}
void
Data::dump(ostream &stream, DisplayMode mode)
{
MainBin *orig = MainBin::curBin();
switch (mode) {
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_start(const string &file)
{
if (stream)
panic("can't start python twice!");
stream = new ofstream(file.c_str(), ios::trunc);
py = new Python(*stream);
ccprintf(*stream, "import sys\n");
ccprintf(*stream, "sys.path.append('.')\n");
ccprintf(*stream, "from m5stats import *\n\n");
}
void
Data::python_dump(const string &name, const string &subname)
{
if (!py)
panic("Can't dump python without first opening the file");
if (bins.empty()) {
python(name, subname, "");
} else {
list<MainBin *>::iterator i = bins.begin();
list<MainBin *>::iterator end = bins.end();
while (i != end) {
(*i)->activate();
python(name, subname, (*i)->name());
++i;
}
}
py->next();
}
void
Data::python(const string &name, const string &subname, const string &bin)
{
py->start("collections.append");
py->start("Collection");
py->qarg(name);
py->qarg(subname);
py->qarg(bin);
py->qarg(hostname());
py->qarg(Time::start.date());
py->startList();
list_t::iterator i = allStats.begin();
list_t::iterator end = allStats.end();
while (i != end) {
StatData *stat = *i;
stat->python(*py);
++i;
}
py->endList();
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.qqqarg(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.qqqarg(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.qqqarg(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.qqqarg(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.qqqarg(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.qqqarg(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
{
if (root)
vec = root->val();
}
result_t
FormulaBase::total() const
{
return root ? root->total() : 0.0;
}
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, DisplayMode mode)
{
Database::StatDB().dump(stream, mode);
}
void
python_start(const string &file)
{
Database::StatDB().python_start(file);
}
void
python_dump(const string &name, const string &subname)
{
Database::StatDB().python_dump(name, subname);
}
CallbackQueue resetQueue;
void
registerResetCallback(Callback *cb)
{
resetQueue.add(cb);
}
void
reset()
{
Database::StatDB().reset();
resetQueue.process();
}
} // namespace Statistics