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]; }
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; }