/usr/lib/petscdir/3.7.7/x86_64-linux-gnu-complex-debug/include/petsc/private/dmdaimpl.h is in libpetsc-complex-3.7.7-dbg 3.7.7+dfsg1-2build5.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 | /*
Distributed arrays - communication tools for parallel, rectangular grids.
*/
#if !defined(_DAIMPL_H)
#define _DAIMPL_H
#include <petscdmda.h>
#include <petsc/private/dmimpl.h>
typedef struct {
PetscInt M,N,P; /* array dimensions */
PetscInt m,n,p; /* processor layout */
PetscInt w; /* degrees of freedom per node */
PetscInt s; /* stencil width */
PetscInt xs,xe,ys,ye,zs,ze; /* range of local values */
PetscInt Xs,Xe,Ys,Ye,Zs,Ze; /* range including ghost values
values above already scaled by w */
PetscInt base; /* global number of 1st local node, includes the * w term */
DMBoundaryType bx,by,bz; /* indicates type of ghost nodes at boundary */
VecScatter gtol,ltol; /* scatters, see below for details */
DMDAStencilType stencil_type; /* stencil, either box or star */
DMDAInterpolationType interptype;
PetscInt nlocal,Nlocal; /* local size of local vector and global vector, includes the * w term */
PetscInt xol,yol,zol; /* overlap of local subdomains */
PetscInt xo,yo,zo; /* offsets for the indices in x y and z */
PetscInt Mo,No,Po; /* the size of the problem the offset is in to */
PetscInt Nsub; /* number of local subdomains to decompose into */
PetscInt nonxs,nonys,nonzs; /* the nonoverlapping starts in the case of a subdomain da */
PetscInt nonxm,nonym,nonzm; /* the nonoverlapping sizes in the case of a subdomain da */
AO ao; /* application ordering context */
AOType aotype; /* type of application ordering */
char **fieldname; /* names of individual components in vectors */
char **coordinatename; /* names of coordinate directions, for example, x, y, z */
PetscInt *lx,*ly,*lz; /* number of nodes in each partition block along 3 axis */
Vec natural; /* global vector for storing items in natural order */
VecScatter gton; /* vector scatter from global to natural */
PetscMPIInt *neighbors; /* ranks of all neighbors and self */
ISColoring localcoloring; /* set by DMCreateColoring() */
ISColoring ghostedcoloring;
DMDAElementType elementtype;
PetscInt ne; /* number of elements */
PetscInt *e; /* the elements */
PetscInt refine_x,refine_y,refine_z; /* ratio used in refining */
PetscInt coarsen_x,coarsen_y,coarsen_z; /* ratio used for coarsening */
PetscBool negativeMNP; /* used in DMSetFromOptions_DA() to check if the initial values provided in code can be changed with options database */
#define DMDA_MAX_WORK_ARRAYS 2 /* work arrays for holding work via DMDAGetArray() */
void *arrayin[DMDA_MAX_WORK_ARRAYS],*arrayout[DMDA_MAX_WORK_ARRAYS];
void *arrayghostedin[DMDA_MAX_WORK_ARRAYS],*arrayghostedout[DMDA_MAX_WORK_ARRAYS];
void *startin[DMDA_MAX_WORK_ARRAYS],*startout[DMDA_MAX_WORK_ARRAYS];
void *startghostedin[DMDA_MAX_WORK_ARRAYS],*startghostedout[DMDA_MAX_WORK_ARRAYS];
PetscErrorCode (*lf)(DM, Vec, Vec, void *);
PetscErrorCode (*lj)(DM, Vec, Vec, void *);
/* used by DMDASetBlockFills() */
PetscInt *ofill,*dfill;
PetscInt *ofillcols;
/* used by DMDASetMatPreallocateOnly() */
PetscBool prealloc_only;
PetscInt preallocCenterDim; /* Dimension of the points which connect adjacent points for preallocation */
} DM_DA;
/*
Vectors:
Global has on each processor the interior degrees of freedom and
no ghost points. This vector is what the solvers usually see.
Local has on each processor the ghost points as well. This is
what code to calculate Jacobians, etc. usually sees.
Vector scatters:
gtol - Global representation to local
ltog - Local representation to global (involves no communication)
ltol - Local representation to local representation, updates the
ghostpoint values in the second vector from (correct) interior
values in the first vector. This is good for explicit
nearest neighbor timestepping.
*/
PETSC_INTERN PetscErrorCode VecView_MPI_DA(Vec,PetscViewer);
PETSC_INTERN PetscErrorCode VecLoad_Default_DA(Vec, PetscViewer);
PETSC_INTERN PetscErrorCode DMView_DA_Matlab(DM,PetscViewer);
PETSC_INTERN PetscErrorCode DMView_DA_Binary(DM,PetscViewer);
PETSC_INTERN PetscErrorCode DMView_DA_VTK(DM,PetscViewer);
PETSC_EXTERN PetscErrorCode DMDAVTKWriteAll(PetscObject,PetscViewer);
PETSC_EXTERN PetscErrorCode DMDASelectFields(DM,PetscInt*,PetscInt**);
PETSC_EXTERN PetscLogEvent DMDA_LocalADFunction;
#endif
|