Biomechanical Joint Model  Author: Anderson Maciel

---

vecarray.cpp

Go to the documentation of this file.

#include "vecarray.h"

namespace LinAlg{



VEC& VECArray::specialVEC::operator=(const VEC &A)
{ 
   assert(n_==A.dim()); // dimension MAY NOT be changed
   return VEC::operator=(A);
}


VECArray::~VECArray() 
{ 
   destroy(); 
}

void VECArray::init_vv()
{
   if (vv == NULL)
   {
      vv = new specialVEC[n_/k_];
      for (int k=0; k<n_/k_; k++)
         vv[k].resize(k_);
   }
   
   for (int i=0; i<n_/k_; i++)
   {
      vv[i].redirect(v_ + i*k_);
  }
}

void VECArray::initialize(unsigned N)
{
   VEC::initialize(N);
   init_vv();
}

void VECArray::destroy()
{  
   VEC::destroy();
   delete [] vv;
   vv = NULL;
}

}

---