void RefinElemComp(TPZCompMesh *cMesh, int indexEl)
{
TPZVec<int64_t > subindex;
int64_t nel = cMesh->ElementVec().NElements();
for(int64_t el=0; el < nel; el++){
TPZCompEl * compEl = cMesh->ElementVec()[el];
if(!compEl) continue;
int64_t ind = compEl->Index();
if(ind==indexEl){
compEl->Divide(indexEl, subindex, 1);
}
}
}
void RefinUniformElemComp(TPZCompMesh *cMesh, int ndiv)
{
TPZVec<int64_t > subindex;
for (int64_t iref = 0; iref < ndiv; iref++) {
TPZAdmChunkVector<TPZCompEl *> elvec = cMesh->ElementVec();
int64_t nel = elvec.NElements();
for(int64_t el=0; el < nel; el++){
TPZCompEl * compEl = elvec[el];
if(!compEl) continue;
int64_t ind = compEl->Index();
compEl->Divide(ind, subindex, 0);
}
}
}