convexity< T > Class Template Reference

#include <mathlib.h>

Collaboration diagram for convexity< T >:

List of all members.

Public Member Functions
	convexity (T a_, T b_, T af_, T bf_)
	convexity ()
void	operator() (bool &result, vector< T > const &x, vector< T > const &xf)
template<typename F >
void	operator() (bool &result, F f, uint n, T const a_, T const b_)
Public Attributes
T	a
T	b
T	af
T	bf

---

Detailed Description

template<typename T>
class convexity< T >

Definition at line 208 of file mathlib.h.

---

Constructor & Destructor Documentation

template<typename T>

convexity< T >::convexity	(	T	a_,
		T	b_,
		T	af_,
		T	bf_
	)			[inline]

Definition at line 220 of file mathlib.h.

    : a(a_), b(b_), af(af_), bf(bf_) {}

template<class T >

convexity< T >::convexity

(

)

[inline]

Definition at line 819 of file mathlib.h.

{}

---

Member Function Documentation

template<class T >

template<typename F >

void convexity< T >::operator()	(	bool &	result,
		F	f,
		uint	n,
		T const	a_,
		T const	b_
	)			[inline]

Definition at line 758 of file mathlib.h.

{
  if (n==0)
  {
    result=false;
    return;
  }

  a = a_;
  b = b_;

  f(af,a);
  f(bf,b);

//  cout << "af=" << af << ", bf=" << bf << endl;

  T one(1.0); //Unbelievabley  (j+1)/(n+1)*(b-a) was somewhere converting to int.

/*
// Building containers and forwarding.
  vector<T> x(n);
  vector<T> xf(n);


  for (uint j=0; j<n; ++j)
  {
    x[j] = a+(j+one)/(n+one)*(b-a);
    f(xf[j],x[j]);
    //cout << "x[" << j << "]=" << x[j] << ", xf[" << j << "]=" << xf[j] << endl;
  }
  operator()(result,x,xf);
*/

  T x;
  T xf;
  T lambda;

  for (uint j=0; j<n; ++j)
  {
    x = a+(j+one)/(n+one)*(b-a);
    f(xf,x);
//    cout << "x[" << j << "]=" << x << ", xf[" << j << "]=" << xf << endl;
    lambda = (x-b)/(a-b);
    if ( lambda*af+bf*(one-lambda) < xf )
    {
      result = false;
      return;
    }
  }
 
}

template<class T >

void convexity< T >::operator()	(	bool &	result,
		vector< T > const &	x,
		vector< T > const &	xf
	)			[inline]

Definition at line 824 of file mathlib.h.

{
  result = true;

  T one(1.0);

  T lambda;

  for (uint i=0; i<x.size(); ++i)
  {
    lambda = (x[i]-b)/(a-b);
    if ( lambda*af+bf*(one-lambda) < xf[i] )
    {
      result = false;
      return;
    }
  }
}

---

Member Data Documentation

template<typename T>

T convexity< T >::a

Definition at line 213 of file mathlib.h.

template<typename T>

T convexity< T >::af

Definition at line 216 of file mathlib.h.

template<typename T>

T convexity< T >::b

Definition at line 214 of file mathlib.h.

template<typename T>

T convexity< T >::bf

Definition at line 217 of file mathlib.h.

---

The documentation for this class was generated from the following file:

• mathlib/mathlib.h