*** empty log message ***

infer.cc

-// $Id $
+// $Id$
 
-#include <cstdlib>
-#include <cmath>
 #include <iostream>
+#include <cstdlib>
 #include <cmath>
 
+#include "infer.h"
+
 extern "C" { 
 #include "libprob.h"
 };
 
 #define POSIX_SOURCE 1
 
-#define USAGE "usage: infer <W> <u> <D> <S> <k1> [<k2> <k3>...<kN>]"
-
 #if !defined(M_PI)
 #define M_PI (4 * atan(1.0))
 #endif
 
-int main(unsigned int argc, char **argv)
-{
-  if ( argc < 6 )
-    {
-      cerr << USAGE << endl;
-      return EXIT_FAILURE;
-    }
-
-  argc--;
-  float W = atof(argv[1]); argc--;
-  float u = atof(argv[2]); argc--;
-  int D = atoi(argv[3]); argc--;
-  int S = atoi(argv[4]); argc--;
-
-  float Em = S * pow(2*M_PI,-0.5*(D+1));
-  cout << "S * pow(2*M_PI,-0.5*(D+1)) = " << Em << endl;
-  Em = pow(W,-D);
-  cout << "pow(W,-D) = " << Em << endl;
-  Em = pow(u,D-1);
-  cout << "pow(u,D-1) = " << Em << endl;
-  Em = exp(-0.5*(u*u));
-  cout << "exp(-0.5*(u*u)) = " << Em << endl;
-  Em = S * pow(2*M_PI,-0.5*(D+1))*pow(W,-D)*pow(u,D-1)*exp(-0.5*(u*u));
-  cout << "Em = " << Em << endl;
+template <class T> T sqr(const T& value) { return value * value; }
+
+Infer::Infer(float dLh, float t, unsigned int V) {
+  Em_ = V * pow(2*M_PI,-2) * dLh * (sqr(t) - 1) * exp(-sqr(t)/2.0);
+  B_ = pow((gamma(2.5)*Em_)/(V*(0.5 + 0.5*erf(-t/sqrt(2)))),(2.0/3.0));
+
+  cout << "E{m} " << Em_ << endl;
+  cout << "Beta = " << B_ << endl;
+}
   
-  cout << gamma(0.5*D+1)*Em << endl;
-  cout << S*(0.5 + 0.5*erf(-u)) << endl;
-  cout << (gamma(0.5*D+1)*Em)/(S*erf(-u)) << endl;
-  float B = pow((gamma(0.5*D+1)*Em)/(S*(0.5 + 0.5*erf(-u/sqrt(2)))),(2.0/(float)D));
-  cout << "Beta = " << B << endl;
-
-  for ( unsigned int n = 0; n < argc; n++ ) {
-    unsigned int k = atoi(argv[5 + n]);
-    cout << 1 - exp(-Em * exp(-B * pow( k * k, 2.0 / (float)D))) << endl;
-  }
-  return EXIT_SUCCESS;
+float Infer::operator() (unsigned int k) {
+  return 1 - exp(-Em_ * exp(-B_ * pow( k , 2.0 / 3.0)));
 }
-

infer.h

 #if !defined(Infer_h)
 #define Infer_h
 
-// $Id $
+// $Id$
 
 class Infer {
 public: