Example #1
0
int main(void)
{
  short vec[256]; 
  int i, n = 20;

  printf("\nold: ");
  for (i = 0; i < n; i++){
    vec[i] = (int) floor(rand());
    printf("%d ", vec[i]);
  }
  
  fillvec(vec, n);
  
  printf("\nnew: ");
  for (i = 0; i < n; i++){
    printf("%d ", vec[i]);
  }
  return vec[1];
}
Example #2
0
double Quad::get_1diag_massmat(int ID){
  double *wa, *wb;
  double *wvec = dvector(0, qtot-1), vmmat;
  Mode mw,*m;

#ifndef PCONTBASE
  double **ba, **bb;
  Mode m1;
  get_moda_GL (qa, &ba);
  get_moda_GL (qb, &bb);
  m1.a = ba[ID];
  m1.b = bb[ID];
  m = &m1;
#else
  Basis *b = getbasis();
  m = b->vert+ID;
#endif

  getzw(qa,&wa,&wa,'a');
  getzw(qb,&wb,&wb,'a');

  mw.a = wa;  mw.b = wb;

  fillvec(&mw, wvec);

  if(curvX)
    dvmul(qtot, wvec, 1, geom->jac.p, 1, wvec, 1);
  else
    dscal(qtot, geom->jac.d, wvec, 1);

  vmmat = Quad_mass_mprod(this, m, wvec);

  free(wvec);

  return vmmat;
}