fill.cpp

Go to the documentation of this file.

#include "statsLib.h"

void fill(const dmatrix& _d,const dvector& _v)
{
 ADUNCONST(dmatrix,d) 
 ADUNCONST(dvector,v)
 int mmin =d.indexmin();
 int mmax =d.indexmax();
 int sz=0;
 int i;
 ivector sz1(mmin,mmax);
 for (i=mmin;i<=mmax;i++)
 {
   sz1(i)=d(i).indexmax()-d(i).indexmin()+1;
   sz+=sz1(i);
 }
 
 if (sz != v.indexmax()-v.indexmin()+1)
 {
   cerr << "Error in sizes in function"
     " void fill(const dmatrix& _d,const dvector& v)" << endl;
   ad_exit(1);
 }
 int offset=v.indexmin();
 for (i=mmin;i<=mmax;i++)
 {
   d(i)=v(offset,offset+sz1(i)-1).shift(d(i).indexmin());
 }
}

void fill(const dvar_matrix& _d,const dvar_vector& _v)
{
 ADUNCONST(dvar_matrix,d)
 ADUNCONST(dvar_vector,v)
 int mmin =d.indexmin();
 int mmax =d.indexmax();
 int sz=0;
 int i;
 ivector sz1(mmin,mmax);
 for (i=mmin;i<=mmax;i++)
 {
   sz1(i)=d(i).indexmax()-d(i).indexmin()+1;
   sz+=sz1(i);
 }
 
 if (sz != v.indexmax()-v.indexmin()+1)
 {
   cerr << "Error in sizes in function"
     " void fill(const dmatrix& _d,const dvector& v)" << endl;
   ad_exit(1);
 }
 int offset=v.indexmin();
 for (i=mmin;i<=mmax;i++)
 {
   d(i)=v(offset,offset+sz1(i)-1).shift(d(i).indexmin());
 }
}