mod_mc1.cpp

Go to the documentation of this file.

#include <admodel.h>

void initial_params::add_random_vector(const dvector& x)
{
  int ii=1;
  for (int i=0;i<num_initial_params;i++)
  {
    if (withinbound(0,(varsptr[i])->phase_start,current_phase))
    {
      (varsptr[i])->add_value((const dvector&)(x),ii);
    }
  }
}

void param_init_number::add_value(const dvector& ndev, const int& _ii)
{
  int& ii=(int&) _ii;
  (*this)+=ndev(ii);
  ii++;
}

void param_init_bounded_number::add_value(const dvector& ndev, const int& _ii)
{
  int& ii=(int&) _ii;
  (*this)+=ndev(ii);
  ii++;
}

void param_init_vector::add_value(const dvector& ndev, const int& _ii)
{
  int& ii=(int&) _ii;
  int mmin=indexmin();
  int mmax=indexmax();
  double* pndevii = ndev.get_v() + ii;
  for (int i=mmin;i<=mmax;i++)
  {
    (*this)(i) += *pndevii;

    ++pndevii;
    ++ii;
  }
}

void param_init_bounded_vector::add_value(const dvector& ndev, const int& _ii)
{
  int& ii=(int&) _ii;
  int mmin=indexmin();
  int mmax=indexmax();
  for (int i=mmin;i<=mmax;i++)
  {
    (*this)(i)+=ndev(ii);
    ii++;
  }
}

void param_init_matrix::add_value(const dvector& ndev, const int& _ii)
{
  int& ii=(int&) _ii;
  int rmin=rowmin();
  int rmax=rowmax();
  for (int i=rmin;i<=rmax;i++)
  {
    int cmin=(*this)(i).indexmin();
    int cmax=(*this)(i).indexmax();
    for (int j=cmin;j<=cmax;j++)
    {
      (*this)(i,j)+=ndev(ii);
      ii++;
    }
  }
}

void param_init_bounded_matrix::add_value(const dvector& ndev, const int& _ii)
{
  int& ii=(int&) _ii;
  int rmin=rowmin();
  int rmax=rowmax();
  for (int i=rmin;i<=rmax;i++)
  {
    int cmin=(*this)(i).indexmin();
    int cmax=(*this)(i).indexmax();
    for (int j=cmin;j<=cmax;j++)
    {
      (*this)(i,j)+=ndev(ii);
      ii++;
    }
  }
}