1 /* 2 Code for manipulating distributed regular arrays in parallel. 3 */ 4 5 #include <petsc/private/dmdaimpl.h> /*I "petscdmda.h" I*/ 6 7 PetscErrorCode DMDestroy_DA(DM da) 8 { 9 DM_DA *dd = (DM_DA *)da->data; 10 11 PetscFunctionBegin; 12 /* destroy the external/common part */ 13 for (PetscInt i = 0; i < DMDA_MAX_WORK_ARRAYS; i++) { 14 PetscCall(PetscFree(dd->startghostedout[i])); 15 PetscCall(PetscFree(dd->startghostedin[i])); 16 PetscCall(PetscFree(dd->startout[i])); 17 PetscCall(PetscFree(dd->startin[i])); 18 } 19 20 PetscCall(VecScatterDestroy(&dd->gtol)); 21 PetscCall(VecScatterDestroy(&dd->ltol)); 22 PetscCall(VecDestroy(&dd->natural)); 23 PetscCall(VecScatterDestroy(&dd->gton)); 24 PetscCall(AODestroy(&dd->ao)); 25 PetscCall(PetscFree(dd->aotype)); 26 27 PetscCall(PetscFree(dd->lx)); 28 PetscCall(PetscFree(dd->ly)); 29 PetscCall(PetscFree(dd->lz)); 30 31 PetscCall(PetscFree(dd->refine_x_hier)); 32 PetscCall(PetscFree(dd->refine_y_hier)); 33 PetscCall(PetscFree(dd->refine_z_hier)); 34 35 if (dd->fieldname) { 36 for (PetscInt i = 0; i < dd->w; i++) PetscCall(PetscFree(dd->fieldname[i])); 37 PetscCall(PetscFree(dd->fieldname)); 38 } 39 if (dd->coordinatename) { 40 for (PetscInt i = 0; i < da->dim; i++) PetscCall(PetscFree(dd->coordinatename[i])); 41 PetscCall(PetscFree(dd->coordinatename)); 42 } 43 PetscCall(ISColoringDestroy(&dd->localcoloring)); 44 PetscCall(ISColoringDestroy(&dd->ghostedcoloring)); 45 46 PetscCall(PetscFree(dd->neighbors)); 47 PetscCall(PetscFree(dd->dfill)); 48 PetscCall(PetscFree(dd->ofill)); 49 PetscCall(PetscFree(dd->ofillcols)); 50 PetscCall(PetscFree(dd->e)); 51 PetscCall(ISDestroy(&dd->ecorners)); 52 53 PetscCall(PetscObjectComposeFunction((PetscObject)da, "DMSetUpGLVisViewer_C", NULL)); 54 55 PetscCall(PetscFree(dd)); 56 PetscFunctionReturn(PETSC_SUCCESS); 57 } 58