Programming

typedef struct{
  int thread_no;
  int allthread;
  mat *G,*C,*B;
  mat *us, *Z;
  vec *samples;
  double *Control, *Info;
} parm;

 
 void sampleLU(int thread_no, int allthread, mat *G, mat *C, mat *B, mat &us, mat &Z, vec &samples, double Control[], double Info[])
{
  Real_Timer lu_symbolic_init, lu_symbolic_free, lu_numerical, lu_solve_time;
  Real_Timer sCpG_run_time;
  int np = samples.size();
  int nDim = B->m;
  int start, stop;
 
  start = thread_no * (int)(np/allthread);
  stop = start + (int)(np/allthread) - 1;
  if( thread_no == allthread-1 ) stop = np-1; 
  for (int i = start; i<=stop; i++ ){
    cs *A;
    if(thread_no == 0) {
      sCpG_run_time.start();
      A = G+sample(i)*C;
      sCpG_run_time.stop();
    }
    else{
      A = G+sample(i)*C;
    }
 
 /* LU decomposition */
...
 if(i == start){
   if(thread_no == 0) {
     lu_symbolic_init.start();
     A = LU;
     lu_symbolic_init.stop();
   }
   else {
    A = LU
   }
 }

 /* solve Az = b  */
 double* z = new double[nDim];
 vec b(nDim);
 b.zeros();
 if(thread_no == 0) lu_solve_time.start();
 LUz=b
 vec zz(z, nDim);
 delete [] z;
 Z.set_col(i, zz); 
 ...
  }
  if(thread_no == 0){
    std::cout << "Thread" << thread_no << " sC+G \t:" << sCpG_run_time.get_time() << std::endl;
    std::cout << "Thread" << thread_no << " symbolic initial time: \t"<<lu_symbolic_init.get_time()<<std::endl;
    std::cout << "Thread" << thread_no << " LU decomposition time \t:" << lu_numerical.get_time() << std::endl;
    std::cout << "Thread" << thread_no << " symbolic free time: \t"<<lu_symbolic_init.get_time()<<std::endl;
    std::cout << "Thread" << thread_no << " LU solve time   \t: " << lu_solve_time.get_time() <<std::endl;
    
  }
}
void * psampleLU(void *arg)
{
  parm *p = (parm *)arg;
  sampleLU(p->thread_no, p->allthread, p->G, p->C, p->B, *(p->us), *(p->Z), *(p->samples), p->Control, p->Info);
  return NULL;
}