void TPZSkylMatrix<TVar>::InitializeElem(const TPZVec<int> &skyline, TPZManVector<REAL> &storage, TPZVec<REAL *> &point) {
int dim = skyline.NElements();
int nel = NumElements(skyline);
storage.Resize(nel);
storage.Fill(0.);
int i;
point.Resize(dim+1);
if(dim) {
point[0] = &storage[0];
point[dim] = &storage[0]+nel;
} else {
point[0] = 0;
}
for(i=1; i<dim+1; i++) point[i] = point[i-1]+(i-1)-skyline[i-1]+1;
}
