00001 /* -*- Mode:C++; c-basic-offset:8; tab-width:8; indent-tabs-mode:t -*- */ 00002 /* 00003 * Copyright (c) 1997 Regents of the University of California. 00004 * All rights reserved. 00005 * 00006 * Redistribution and use in source and binary forms, with or without 00007 * modification, are permitted provided that the following conditions 00008 * are met: 00009 * 1. Redistributions of source code must retain the above copyright 00010 * notice, this list of conditions and the following disclaimer. 00011 * 2. Redistributions in binary form must reproduce the above copyright 00012 * notice, this list of conditions and the following disclaimer in the 00013 * documentation and/or other materials provided with the distribution. 00014 * 3. All advertising materials mentioning features or use of this software 00015 * must display the following acknowledgement: 00016 * This product includes software developed by the MASH Research 00017 * Group at the University of California Berkeley. 00018 * 4. Neither the name of the University nor of the Research Group may be 00019 * used to endorse or promote products derived from this software without 00020 * specific prior written permission. 00021 * 00022 * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND 00023 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 00024 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE 00025 * ARE DISCLAIMED. IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE 00026 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL 00027 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS 00028 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) 00029 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT 00030 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY 00031 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF 00032 * SUCH DAMAGE. 00033 * 00034 * @(#) $Header: /cvsroot/nsnam/ns-2/tools/integrator.h,v 1.9 1999/09/28 03:46:29 heideman Exp $ (LBL) 00035 */ 00036 00037 #ifndef ns_integrator_h 00038 #define ns_integrator_h 00039 00040 #include "config.h" 00041 00042 class Integrator : public TclObject { 00043 public: 00044 Integrator(); 00045 void set(double x, double y); 00046 void newPoint(double x, double y); 00047 int command(int argc, const char*const* argv); 00048 protected: 00049 double lastx_; 00050 double lasty_; 00051 double sum_; 00052 }; 00053 00054 // a set of statistical samples 00055 class Samples : public TclObject { 00056 public: 00057 Samples() : cnt_(0), sum_(0.0), sqsum_(0.0) { } 00058 void newPoint(double val) { 00059 cnt_++; 00060 sum_ += val; 00061 val *= val; 00062 sqsum_ += val; 00063 } 00064 int cnt() const { return (cnt_); } 00065 double sum() const { return (sum_); } 00066 double mean() const { 00067 if (cnt_) 00068 return (sum_ / cnt_); 00069 return sum_; // yes, that is 0.0... 00070 } 00071 double variance() const { 00072 // use cnt_-1 degrees of freedom 00073 if (cnt_ > 1) 00074 return ((sqsum_ - mean() * sum_) / (cnt_ - 1)); 00075 return 0.0; 00076 } 00077 void reset() { cnt_ = 0; sum_ = sqsum_ = 0.0; } 00078 int command(int argc, const char*const* argv); 00079 protected: 00080 int cnt_; // count of samples 00081 double sum_; // sum of x_i 00082 double sqsum_; // sum of (x_i)^2 00083 }; 00084 #endif