#include <simplex.hxx>

Collaboration diagram for Simplex:

Public Member Functions
	Simplex (long, long, Maestro &)

	~Simplex (void)

void	setStop (long)

void	start (void)

int	next (void)

void	finish (void)

void	solve (void)

Solution *	solution (void)

Protected Member Functions
Solution *	add (Solution *)

std::vector< double > *	minimum (void)

std::vector< double > *	barycentre (void)

Protected Attributes
long	size

long	nbin

long	maxe

long	nbeval

std::vector< Solution * >	simplx

std::vector< Point * >	work

Maestro *	calc

Detailed Description

Definition at line 29 of file simplex.hxx.

Constructor & Destructor Documentation

Simplex::Simplex	(	long	sz,
		long	nbgen,
		Maestro &	maest
	)

Definition at line 26 of file simplex.cxx.

References calc, maxe, nbeval, nbin, and size.

 {
     size = sz;
     nbin = nbgen;
     calc = &maest;
     maxe = 40000;
     nbeval = 0;
 }

Simplex::~Simplex ( void )

Definition at line 35 of file simplex.cxx.

References CORBAEngineTest::i, nbin, simplx, size, and work.

 {
     int     i;
     
     if (work.size() == size) {
         for (i=0; i<size; i++)
             delete work[i];
         for (i=0; i<nbin; i++)
             delete simplx[i];
     }
 }

Member Function Documentation

Solution * Simplex::add ( Solution * sol )

protected

Definition at line 101 of file simplex.cxx.

References CORBAEngineTest::i, nbin, Solution::obj, PMMLBasicsTestLauncher::ret, and simplx.

Referenced by next().

 {
     int     i;
     Solution    *ret, *swp;
     
     if (simplx.size() < nbin) {
         ret = (Solution *) NULL;
         i = simplx.size();
         simplx.push_back(sol);
         while (i && simplx[i]->obj[0] < simplx[i-1]->obj[0]) {
             // TODO swap interne
             swp = simplx[i];
             simplx[i] = simplx[i-1];
             simplx[i-1] = swp;
             i--;
         }
     }
     else if (sol->obj[0] > simplx[nbin-1]->obj[0])
         ret = sol;
     else {
         i = nbin -1;
         ret = simplx[i];
         simplx[i] = sol;
         while (i && simplx[i]->obj[0] < simplx[i-1]->obj[0]) {
             swp = simplx[i];
             simplx[i] = simplx[i-1];
             simplx[i-1] = swp;
             i--;
         }
     }
     return ret;
 }

std::vector< double > * Simplex::barycentre ( void )

protected

Definition at line 139 of file simplex.cxx.

References CORBAEngineTest::i, nbin, PMMLBasicsTestLauncher::ret, and simplx.

Referenced by next().

 {
     int     i,j;
     std::vector<double>     *ret;
 
     ret = new std::vector<double>(nbin);
     for (i=0; i<nbin; i++) {
         (*ret)[i] = 0.0;
         for (j=0; j<nbin; j++)
             (*ret)[i] += (*simplx[j]->param)[i];
         (*ret)[i] /= nbin;
     }
     return ret;
 }

void Simplex::finish ( void )

Definition at line 90 of file simplex.cxx.

References maxe, and nbeval.

Referenced by YACS::ENGINE::PluginSimplex::finish().

 {
     std::cout << maxe - nbeval << std::endl;
     return;
 }

std::vector< double > * Simplex::minimum ( void )

protected

Definition at line 134 of file simplex.cxx.

References simplx.

Referenced by next().

 {
     return new std::vector<double>(*simplx[0]->param);
 }

int Simplex::next ( void )

Definition at line 65 of file simplex.cxx.

References add(), barycentre(), calc, Maestro::get(), minimum(), nbeval, Maestro::put(), size, and work.

Referenced by start(), and YACS::ENGINE::PluginSimplex::takeDecision().

 {
     long        id;
     Solution    *pt, *pnt;
     std::vector<double>     *next, *res;
     
     res = calc->get(&id);
     nbeval--;
     pt = work[id]->inform(*res);
     if (pt) {
         pnt = add(pt);
         if (pnt)
             work[id]->mute(*pnt, *barycentre(), *minimum());
         else
             work[id]->reinit();
     }
     next = work[id]->next();
     if (nbeval > size)
         if (next) 
             calc->put(id, *next);
         //else
         //    nbeval = size - 1;
     return (nbeval > 0);
 }

void Simplex::setStop ( long max )

Definition at line 47 of file simplex.cxx.

References maxe.

Referenced by YACS::ENGINE::PluginSimplex::parseFileToInit().

 {
     maxe = max;
 }

Solution * Simplex::solution ( void )

Definition at line 96 of file simplex.cxx.

References simplx.

Referenced by YACS::ENGINE::PluginSimplex::finish().

 {
     return simplx[0];
 }

void Simplex::solve ( void )

void Simplex::start ( void )

Definition at line 52 of file simplex.cxx.

References calc, CORBAEngineTest::i, maxe, nbeval, nbin, next(), Maestro::put(), size, and work.

Referenced by YACS::ENGINE::PluginSimplex::start().

 {
     long    i;
 
     work.resize(size);
     for (i=0; i<size; i++) {
         work[i] = new Point(nbin);
         calc->put(i, *(work[i]->next()));
     }
     nbeval = maxe;
 }