Bug Summary

File:dm/impls/da/dalocal.c
Warning:line 346, column 9
Dereference of null pointer (loaded from variable 'cone')

Annotated Source Code

[?] Use j/k keys for keyboard navigation

1
2/*
3 Code for manipulating distributed regular arrays in parallel.
4*/
5
6#include <petsc/private/dmdaimpl.h> /*I "petscdmda.h" I*/
7#include <petscbt.h>
8#include <petscsf.h>
9#include <petscds.h>
10#include <petscfe.h>
11
12/*
13 This allows the DMDA vectors to properly tell MATLAB their dimensions
14*/
15#if defined(PETSC_HAVE_MATLAB_ENGINE)
16#include <engine.h> /* MATLAB include file */
17#include <mex.h> /* MATLAB include file */
18static PetscErrorCode VecMatlabEnginePut_DA2d(PetscObject obj,void *mengine)
19{
20 PetscErrorCode ierr;
21 PetscInt n,m;
22 Vec vec = (Vec)obj;
23 PetscScalar *array;
24 mxArray *mat;
25 DM da;
26
27 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 27; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
28 ierr = VecGetDM(vec, &da);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),28,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
29 if (!da) SETERRQ(PetscObjectComm((PetscObject)vec),PETSC_ERR_ARG_WRONGSTATE,"Vector not associated with a DMDA")return PetscError(PetscObjectComm((PetscObject)vec),29,__func__
,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c",73
,PETSC_ERROR_INITIAL,"Vector not associated with a DMDA");
30 ierr = DMDAGetGhostCorners(da,0,0,0,&m,&n,0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),30,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
31
32 ierr = VecGetArray(vec,&array);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),32,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
33#if !defined(PETSC_USE_COMPLEX)
34 mat = mxCreateDoubleMatrix(m,n,mxREAL);
35#else
36 mat = mxCreateDoubleMatrix(m,n,mxCOMPLEX);
37#endif
38 ierr = PetscMemcpy(mxGetPr(mat),array,n*m*sizeof(PetscScalar));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),38,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
39 ierr = PetscObjectName(obj);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),39,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
40 engPutVariable((Engine*)mengine,obj->name,mat);
41
42 ierr = VecRestoreArray(vec,&array);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),42,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
43 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
44}
45#endif
46
47
48PetscErrorCode DMCreateLocalVector_DA(DM da,Vec *g)
49{
50 PetscErrorCode ierr;
51 DM_DA *dd = (DM_DA*)da->data;
52
53 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 53; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
54 PetscValidHeaderSpecific(da,DM_CLASSID,1)do { if (!da) return PetscError(((MPI_Comm)0x44000001),54,__func__
,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c",85
,PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",1); if (!PetscCheckPointer
(da,PETSC_OBJECT)) return PetscError(((MPI_Comm)0x44000001),54
,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,1); if (((PetscObject)(da))->classid != DM_CLASSID) { if (
((PetscObject)(da))->classid == -1) return PetscError(((MPI_Comm
)0x44000001),54,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,1); else return PetscError(((MPI_Comm)0x44000001),54,__func__
,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c",62
,PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d",1
); } } while (0);
55 PetscValidPointer(g,2)do { if (!g) return PetscError(((MPI_Comm)0x44000001),55,__func__
,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c",85
,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",2); if (!
PetscCheckPointer(g,PETSC_CHAR)) return PetscError(((MPI_Comm
)0x44000001),55,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer: Parameter # %d",2);
} while (0);
56 if (da->defaultSection) {
57 ierr = DMCreateLocalVector_Section_Private(da,g);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),57,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
58 } else {
59 ierr = VecCreate(PETSC_COMM_SELF((MPI_Comm)0x44000001),g);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),59,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
60 ierr = VecSetSizes(*g,dd->nlocal,PETSC_DETERMINE-1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),60,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
61 ierr = VecSetBlockSize(*g,dd->w);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),61,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
62 ierr = VecSetType(*g,da->vectype);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),62,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
63 ierr = VecSetDM(*g, da);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),63,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
64#if defined(PETSC_HAVE_MATLAB_ENGINE)
65 if (dd->w == 1 && da->dim == 2) {
66 ierr = PetscObjectComposeFunction((PetscObject)*g,"PetscMatlabEnginePut_C",VecMatlabEnginePut_DA2d)PetscObjectComposeFunction_Private((PetscObject)*g,"PetscMatlabEnginePut_C"
,(PetscVoidFunction)(VecMatlabEnginePut_DA2d));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),66,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
67 }
68#endif
69 }
70 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
71}
72
73/*@
74 DMDAGetNumCells - Get the number of cells in the local piece of the DMDA. This includes ghost cells.
75
76 Input Parameter:
77. dm - The DM object
78
79 Output Parameters:
80+ numCellsX - The number of local cells in the x-direction
81. numCellsY - The number of local cells in the y-direction
82. numCellsZ - The number of local cells in the z-direction
83- numCells - The number of local cells
84
85 Level: developer
86
87.seealso: DMDAGetCellPoint()
88@*/
89PetscErrorCode DMDAGetNumCells(DM dm, PetscInt *numCellsX, PetscInt *numCellsY, PetscInt *numCellsZ, PetscInt *numCells)
90{
91 DM_DA *da = (DM_DA*) dm->data;
92 const PetscInt dim = dm->dim;
93 const PetscInt mx = (da->Xe - da->Xs)/da->w, my = da->Ye - da->Ys, mz = da->Ze - da->Zs;
94 const PetscInt nC = (mx)*(dim > 1 ? (my)*(dim > 2 ? (mz) : 1) : 1);
95
96 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 96; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
97 PetscValidHeaderSpecificType(dm, DM_CLASSID, 1,DMDA)do { PetscErrorCode __ierr; PetscBool same; do { if (!dm) return
PetscError(((MPI_Comm)0x44000001),97,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,85,PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",1); if (
!PetscCheckPointer(dm,PETSC_OBJECT)) return PetscError(((MPI_Comm
)0x44000001),97,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,1); if (((PetscObject)(dm))->classid != DM_CLASSID) { if (
((PetscObject)(dm))->classid == -1) return PetscError(((MPI_Comm
)0x44000001),97,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,1); else return PetscError(((MPI_Comm)0x44000001),97,__func__
,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c",62
,PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d",1
); } } while (0); __ierr = PetscObjectTypeCompare((PetscObject
)dm,"da",&same);do {if (__builtin_expect(!!(__ierr),0)) return
PetscError(((MPI_Comm)0x44000001),97,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0); if (!same) return
PetscError(((MPI_Comm)0x44000001),97,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,62,PETSC_ERROR_INITIAL,"Wrong subtype object:Parameter # %d must have implementation %s it is %s"
,1,"da",((PetscObject)dm)->type_name); } while (0);
98 if (numCellsX) {
99 PetscValidIntPointer(numCellsX,2)do { if (!numCellsX) return PetscError(((MPI_Comm)0x44000001)
,99,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",2); if
(!PetscCheckPointer(numCellsX,PETSC_INT)) return PetscError(
((MPI_Comm)0x44000001),99,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,2); } while (0);
100 *numCellsX = mx;
101 }
102 if (numCellsY) {
103 PetscValidIntPointer(numCellsY,3)do { if (!numCellsY) return PetscError(((MPI_Comm)0x44000001)
,103,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",3); if
(!PetscCheckPointer(numCellsY,PETSC_INT)) return PetscError(
((MPI_Comm)0x44000001),103,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,3); } while (0);
104 *numCellsY = my;
105 }
106 if (numCellsZ) {
107 PetscValidIntPointer(numCellsZ,4)do { if (!numCellsZ) return PetscError(((MPI_Comm)0x44000001)
,107,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",4); if
(!PetscCheckPointer(numCellsZ,PETSC_INT)) return PetscError(
((MPI_Comm)0x44000001),107,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,4); } while (0);
108 *numCellsZ = mz;
109 }
110 if (numCells) {
111 PetscValidIntPointer(numCells,5)do { if (!numCells) return PetscError(((MPI_Comm)0x44000001),
111,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",5); if
(!PetscCheckPointer(numCells,PETSC_INT)) return PetscError((
(MPI_Comm)0x44000001),111,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,5); } while (0);
112 *numCells = nC;
113 }
114 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
115}
116
117/*@
118 DMDAGetCellPoint - Get the DM point corresponding to the tuple (i, j, k) in the DMDA
119
120 Input Parameters:
121+ dm - The DM object
122- i,j,k - The global indices for the cell
123
124 Output Parameters:
125. point - The local DM point
126
127 Level: developer
128
129.seealso: DMDAGetNumCells()
130@*/
131PetscErrorCode DMDAGetCellPoint(DM dm, PetscInt i, PetscInt j, PetscInt k, PetscInt *point)
132{
133 const PetscInt dim = dm->dim;
134 DMDALocalInfo info;
135 PetscErrorCode ierr;
136
137 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 137; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
138 PetscValidHeaderSpecificType(dm, DM_CLASSID, 1,DMDA)do { PetscErrorCode __ierr; PetscBool same; do { if (!dm) return
PetscError(((MPI_Comm)0x44000001),138,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,85,PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",1); if (
!PetscCheckPointer(dm,PETSC_OBJECT)) return PetscError(((MPI_Comm
)0x44000001),138,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,1); if (((PetscObject)(dm))->classid != DM_CLASSID) { if (
((PetscObject)(dm))->classid == -1) return PetscError(((MPI_Comm
)0x44000001),138,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,1); else return PetscError(((MPI_Comm)0x44000001),138,__func__
,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c",62
,PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d",1
); } } while (0); __ierr = PetscObjectTypeCompare((PetscObject
)dm,"da",&same);do {if (__builtin_expect(!!(__ierr),0)) return
PetscError(((MPI_Comm)0x44000001),138,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0); if (!same) return
PetscError(((MPI_Comm)0x44000001),138,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,62,PETSC_ERROR_INITIAL,"Wrong subtype object:Parameter # %d must have implementation %s it is %s"
,1,"da",((PetscObject)dm)->type_name); } while (0);
139 PetscValidIntPointer(point,5)do { if (!point) return PetscError(((MPI_Comm)0x44000001),139
,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",5); if
(!PetscCheckPointer(point,PETSC_INT)) return PetscError(((MPI_Comm
)0x44000001),139,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,5); } while (0);
140 ierr = DMDAGetLocalInfo(dm, &info);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),140,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
141 if (dim > 0) {if ((i < info.gxs) || (i >= info.gxs+info.gxm)) SETERRQ3(PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "X index %d not in [%d, %d)", i, info.gxs, info.gxs+info.gxm)return PetscError(((MPI_Comm)0x44000001),141,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,63,PETSC_ERROR_INITIAL,"X index %d not in [%d, %d)",i,info.gxs
,info.gxs+info.gxm);}
142 if (dim > 1) {if ((j < info.gys) || (j >= info.gys+info.gym)) SETERRQ3(PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Y index %d not in [%d, %d)", j, info.gys, info.gys+info.gym)return PetscError(((MPI_Comm)0x44000001),142,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,63,PETSC_ERROR_INITIAL,"Y index %d not in [%d, %d)",j,info.gys
,info.gys+info.gym);}
143 if (dim > 2) {if ((k < info.gzs) || (k >= info.gzs+info.gzm)) SETERRQ3(PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Z index %d not in [%d, %d)", k, info.gzs, info.gzs+info.gzm)return PetscError(((MPI_Comm)0x44000001),143,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,63,PETSC_ERROR_INITIAL,"Z index %d not in [%d, %d)",k,info.gzs
,info.gzs+info.gzm);}
144 *point = i + (dim > 1 ? (j + (dim > 2 ? k*info.gym : 0))*info.gxm : 0);
145 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
146}
147
148PetscErrorCode DMDAGetNumVertices(DM dm, PetscInt *numVerticesX, PetscInt *numVerticesY, PetscInt *numVerticesZ, PetscInt *numVertices)
149{
150 DM_DA *da = (DM_DA*) dm->data;
151 const PetscInt dim = dm->dim;
152 const PetscInt mx = (da->Xe - da->Xs)/da->w, my = da->Ye - da->Ys, mz = da->Ze - da->Zs;
153 const PetscInt nVx = mx+1;
154 const PetscInt nVy = dim > 1 ? (my+1) : 1;
155 const PetscInt nVz = dim > 2 ? (mz+1) : 1;
156 const PetscInt nV = nVx*nVy*nVz;
157
158 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 158; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
159 if (numVerticesX) {
160 PetscValidIntPointer(numVerticesX,2)do { if (!numVerticesX) return PetscError(((MPI_Comm)0x44000001
),160,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",2); if
(!PetscCheckPointer(numVerticesX,PETSC_INT)) return PetscError
(((MPI_Comm)0x44000001),160,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,2); } while (0);
161 *numVerticesX = nVx;
162 }
163 if (numVerticesY) {
164 PetscValidIntPointer(numVerticesY,3)do { if (!numVerticesY) return PetscError(((MPI_Comm)0x44000001
),164,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",3); if
(!PetscCheckPointer(numVerticesY,PETSC_INT)) return PetscError
(((MPI_Comm)0x44000001),164,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,3); } while (0);
165 *numVerticesY = nVy;
166 }
167 if (numVerticesZ) {
168 PetscValidIntPointer(numVerticesZ,4)do { if (!numVerticesZ) return PetscError(((MPI_Comm)0x44000001
),168,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",4); if
(!PetscCheckPointer(numVerticesZ,PETSC_INT)) return PetscError
(((MPI_Comm)0x44000001),168,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,4); } while (0);
169 *numVerticesZ = nVz;
170 }
171 if (numVertices) {
172 PetscValidIntPointer(numVertices,5)do { if (!numVertices) return PetscError(((MPI_Comm)0x44000001
),172,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",5); if
(!PetscCheckPointer(numVertices,PETSC_INT)) return PetscError
(((MPI_Comm)0x44000001),172,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,5); } while (0);
173 *numVertices = nV;
174 }
175 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
176}
177
178PetscErrorCode DMDAGetNumFaces(DM dm, PetscInt *numXFacesX, PetscInt *numXFaces, PetscInt *numYFacesY, PetscInt *numYFaces, PetscInt *numZFacesZ, PetscInt *numZFaces)
179{
180 DM_DA *da = (DM_DA*) dm->data;
181 const PetscInt dim = dm->dim;
182 const PetscInt mx = (da->Xe - da->Xs)/da->w, my = da->Ye - da->Ys, mz = da->Ze - da->Zs;
183 const PetscInt nxF = (dim > 1 ? (my)*(dim > 2 ? (mz) : 1) : 1);
184 const PetscInt nXF = (mx+1)*nxF;
185 const PetscInt nyF = mx*(dim > 2 ? mz : 1);
186 const PetscInt nYF = dim > 1 ? (my+1)*nyF : 0;
187 const PetscInt nzF = mx*(dim > 1 ? my : 0);
188 const PetscInt nZF = dim > 2 ? (mz+1)*nzF : 0;
189
190 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 190; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
191 if (numXFacesX) {
192 PetscValidIntPointer(numXFacesX,2)do { if (!numXFacesX) return PetscError(((MPI_Comm)0x44000001
),192,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",2); if
(!PetscCheckPointer(numXFacesX,PETSC_INT)) return PetscError
(((MPI_Comm)0x44000001),192,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,2); } while (0);
193 *numXFacesX = nxF;
194 }
195 if (numXFaces) {
196 PetscValidIntPointer(numXFaces,3)do { if (!numXFaces) return PetscError(((MPI_Comm)0x44000001)
,196,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",3); if
(!PetscCheckPointer(numXFaces,PETSC_INT)) return PetscError(
((MPI_Comm)0x44000001),196,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,3); } while (0);
197 *numXFaces = nXF;
198 }
199 if (numYFacesY) {
200 PetscValidIntPointer(numYFacesY,4)do { if (!numYFacesY) return PetscError(((MPI_Comm)0x44000001
),200,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",4); if
(!PetscCheckPointer(numYFacesY,PETSC_INT)) return PetscError
(((MPI_Comm)0x44000001),200,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,4); } while (0);
201 *numYFacesY = nyF;
202 }
203 if (numYFaces) {
204 PetscValidIntPointer(numYFaces,5)do { if (!numYFaces) return PetscError(((MPI_Comm)0x44000001)
,204,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",5); if
(!PetscCheckPointer(numYFaces,PETSC_INT)) return PetscError(
((MPI_Comm)0x44000001),204,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,5); } while (0);
205 *numYFaces = nYF;
206 }
207 if (numZFacesZ) {
208 PetscValidIntPointer(numZFacesZ,6)do { if (!numZFacesZ) return PetscError(((MPI_Comm)0x44000001
),208,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",6); if
(!PetscCheckPointer(numZFacesZ,PETSC_INT)) return PetscError
(((MPI_Comm)0x44000001),208,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,6); } while (0);
209 *numZFacesZ = nzF;
210 }
211 if (numZFaces) {
212 PetscValidIntPointer(numZFaces,7)do { if (!numZFaces) return PetscError(((MPI_Comm)0x44000001)
,212,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",7); if
(!PetscCheckPointer(numZFaces,PETSC_INT)) return PetscError(
((MPI_Comm)0x44000001),212,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,7); } while (0);
213 *numZFaces = nZF;
214 }
215 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
216}
217
218PetscErrorCode DMDAGetHeightStratum(DM dm, PetscInt height, PetscInt *pStart, PetscInt *pEnd)
219{
220 const PetscInt dim = dm->dim;
221 PetscInt nC, nV, nXF, nYF, nZF;
222 PetscErrorCode ierr;
223
224 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 224; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
225 if (pStart) PetscValidIntPointer(pStart,3)do { if (!pStart) return PetscError(((MPI_Comm)0x44000001),225
,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",3); if
(!PetscCheckPointer(pStart,PETSC_INT)) return PetscError(((MPI_Comm
)0x44000001),225,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,3); } while (0);
226 if (pEnd) PetscValidIntPointer(pEnd,4)do { if (!pEnd) return PetscError(((MPI_Comm)0x44000001),226,
__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",4); if
(!PetscCheckPointer(pEnd,PETSC_INT)) return PetscError(((MPI_Comm
)0x44000001),226,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,4); } while (0);
227 ierr = DMDAGetNumCells(dm, NULL((void*)0), NULL((void*)0), NULL((void*)0), &nC);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),227,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
228 ierr = DMDAGetNumVertices(dm, NULL((void*)0), NULL((void*)0), NULL((void*)0), &nV);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),228,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
229 ierr = DMDAGetNumFaces(dm, NULL((void*)0), &nXF, NULL((void*)0), &nYF, NULL((void*)0), &nZF);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),229,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
230 if (height == 0) {
231 /* Cells */
232 if (pStart) *pStart = 0;
233 if (pEnd) *pEnd = nC;
234 } else if (height == 1) {
235 /* Faces */
236 if (pStart) *pStart = nC+nV;
237 if (pEnd) *pEnd = nC+nV+nXF+nYF+nZF;
238 } else if (height == dim) {
239 /* Vertices */
240 if (pStart) *pStart = nC;
241 if (pEnd) *pEnd = nC+nV;
242 } else if (height < 0) {
243 /* All points */
244 if (pStart) *pStart = 0;
245 if (pEnd) *pEnd = nC+nV+nXF+nYF+nZF;
246 } else SETERRQ1(PetscObjectComm((PetscObject)dm), PETSC_ERR_ARG_OUTOFRANGE, "No points of height %d in the DA", height)return PetscError(PetscObjectComm((PetscObject)dm),246,__func__
,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c",63
,PETSC_ERROR_INITIAL,"No points of height %d in the DA",height
);
247 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
248}
249
250PetscErrorCode DMDAGetDepthStratum(DM dm, PetscInt depth, PetscInt *pStart, PetscInt *pEnd)
251{
252 const PetscInt dim = dm->dim;
253 PetscInt nC, nV, nXF, nYF, nZF;
254 PetscErrorCode ierr;
255
256 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 256; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
257 if (pStart) PetscValidIntPointer(pStart,3)do { if (!pStart) return PetscError(((MPI_Comm)0x44000001),257
,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",3); if
(!PetscCheckPointer(pStart,PETSC_INT)) return PetscError(((MPI_Comm
)0x44000001),257,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,3); } while (0);
258 if (pEnd) PetscValidIntPointer(pEnd,4)do { if (!pEnd) return PetscError(((MPI_Comm)0x44000001),258,
__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",4); if
(!PetscCheckPointer(pEnd,PETSC_INT)) return PetscError(((MPI_Comm
)0x44000001),258,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to Int: Parameter # %d"
,4); } while (0);
259 ierr = DMDAGetNumCells(dm, NULL((void*)0), NULL((void*)0), NULL((void*)0), &nC);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),259,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
260 ierr = DMDAGetNumVertices(dm, NULL((void*)0), NULL((void*)0), NULL((void*)0), &nV);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),260,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
261 ierr = DMDAGetNumFaces(dm, NULL((void*)0), &nXF, NULL((void*)0), &nYF, NULL((void*)0), &nZF);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),261,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
262 if (depth == dim) {
263 /* Cells */
264 if (pStart) *pStart = 0;
265 if (pEnd) *pEnd = nC;
266 } else if (depth == dim-1) {
267 /* Faces */
268 if (pStart) *pStart = nC+nV;
269 if (pEnd) *pEnd = nC+nV+nXF+nYF+nZF;
270 } else if (depth == 0) {
271 /* Vertices */
272 if (pStart) *pStart = nC;
273 if (pEnd) *pEnd = nC+nV;
274 } else if (depth < 0) {
275 /* All points */
276 if (pStart) *pStart = 0;
277 if (pEnd) *pEnd = nC+nV+nXF+nYF+nZF;
278 } else SETERRQ1(PetscObjectComm((PetscObject)dm), PETSC_ERR_ARG_OUTOFRANGE, "No points of depth %d in the DA", depth)return PetscError(PetscObjectComm((PetscObject)dm),278,__func__
,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c",63
,PETSC_ERROR_INITIAL,"No points of depth %d in the DA",depth);
279 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
280}
281
282PetscErrorCode DMDAGetConeSize(DM dm, PetscInt p, PetscInt *coneSize)
283{
284 const PetscInt dim = dm->dim;
285 PetscInt nC, nV, nXF, nYF, nZF;
286 PetscErrorCode ierr;
287
288 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 288; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
289 *coneSize = 0;
290 ierr = DMDAGetNumCells(dm, NULL((void*)0), NULL((void*)0), NULL((void*)0), &nC);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),290,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
291 ierr = DMDAGetNumVertices(dm, NULL((void*)0), NULL((void*)0), NULL((void*)0), &nV);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),291,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
292 ierr = DMDAGetNumFaces(dm, NULL((void*)0), &nXF, NULL((void*)0), &nYF, NULL((void*)0), &nZF);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),292,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
293 switch (dim) {
294 case 2:
295 if (p >= 0) {
296 if (p < nC) {
297 *coneSize = 4;
298 } else if (p < nC+nV) {
299 *coneSize = 0;
300 } else if (p < nC+nV+nXF+nYF+nZF) {
301 *coneSize = 2;
302 } else SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Point %d should be in [0, %d)", p, nC+nV+nXF+nYF+nZF)return PetscError(((MPI_Comm)0x44000001),302,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,63,PETSC_ERROR_INITIAL,"Point %d should be in [0, %d)",p,nC+
nV+nXF+nYF+nZF);
303 } else SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Negative point %d is invalid", p)return PetscError(((MPI_Comm)0x44000001),303,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,63,PETSC_ERROR_INITIAL,"Negative point %d is invalid",p);
304 break;
305 case 3:
306 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_SUP, "Too lazy to do 3D")return PetscError(((MPI_Comm)0x44000001),306,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,56,PETSC_ERROR_INITIAL,"Too lazy to do 3D");
307 break;
308 }
309 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
310}
311
312PetscErrorCode DMDAGetCone(DM dm, PetscInt p, PetscInt *cone[])
313{
314 const PetscInt dim = dm->dim;
315 PetscInt nCx, nCy, nCz, nC, nVx, nVy, nVz, nV, nxF, nyF, nzF, nXF, nYF, nZF;
316 PetscErrorCode ierr;
317
318 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 318; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
319 if (!cone) {ierr = DMGetWorkArray(dm, 6, MPIU_INT((MPI_Datatype)0x4c000405), cone);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),319,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}
1
Assuming 'cone' is null
→
2
←
Taking true branch
→
320 ierr = DMDAGetNumCells(dm, &nCx, &nCy, &nCz, &nC);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),320,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
321 ierr = DMDAGetNumVertices(dm, &nVx, &nVy, &nVz, &nV);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),321,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
322 ierr = DMDAGetNumFaces(dm, &nxF, &nXF, &nyF, &nYF, &nzF, &nZF);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),322,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
323 switch (dim) {
3
←
Control jumps to 'case 2:' at line 324
→
324 case 2:
325 if (p >= 0) {
4
←
Assuming 'p' is >= 0
→
5
←
Taking true branch
→
326 if (p < nC) {
6
←
Taking false branch
→
327 const PetscInt cy = p / nCx;
328 const PetscInt cx = p % nCx;
329
330 (*cone)[0] = cy*nxF + cx + nC+nV;
331 (*cone)[1] = cx*nyF + cy + nyF + nC+nV+nXF;
332 (*cone)[2] = cy*nxF + cx + nxF + nC+nV;
333 (*cone)[3] = cx*nyF + cy + nC+nV+nXF;
334 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_SUP, "Too lazy to do cell cones")return PetscError(((MPI_Comm)0x44000001),334,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,56,PETSC_ERROR_INITIAL,"Too lazy to do cell cones");
335 } else if (p < nC+nV) {
7
←
Taking false branch
→
336 } else if (p < nC+nV+nXF) {
8
←
Taking false branch
→
337 const PetscInt fy = (p - nC+nV) / nxF;
338 const PetscInt fx = (p - nC+nV) % nxF;
339
340 (*cone)[0] = fy*nVx + fx + nC;
341 (*cone)[1] = fy*nVx + fx + 1 + nC;
342 } else if (p < nC+nV+nXF+nYF) {
9
←
Taking true branch
→
343 const PetscInt fx = (p - nC+nV+nXF) / nyF;
344 const PetscInt fy = (p - nC+nV+nXF) % nyF;
345
346 (*cone)[0] = fy*nVx + fx + nC;
10
←
Dereference of null pointer (loaded from variable 'cone')
347 (*cone)[1] = fy*nVx + fx + nVx + nC;
348 } else SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Point %d should be in [0, %d)", p, nC+nV+nXF+nYF+nZF)return PetscError(((MPI_Comm)0x44000001),348,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,63,PETSC_ERROR_INITIAL,"Point %d should be in [0, %d)",p,nC+
nV+nXF+nYF+nZF);
349 } else SETERRQ1(PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Negative point %d is invalid", p)return PetscError(((MPI_Comm)0x44000001),349,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,63,PETSC_ERROR_INITIAL,"Negative point %d is invalid",p);
350 break;
351 case 3:
352 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_SUP, "Too lazy to do 3D")return PetscError(((MPI_Comm)0x44000001),352,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,56,PETSC_ERROR_INITIAL,"Too lazy to do 3D");
353 break;
354 }
355 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
356}
357
358PetscErrorCode DMDARestoreCone(DM dm, PetscInt p, PetscInt *cone[])
359{
360 PetscErrorCode ierr;
361
362 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 362; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
363 ierr = DMGetWorkArray(dm, 6, MPIU_INT((MPI_Datatype)0x4c000405), cone);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),363,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
364 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
365}
366
367/*@C
368 DMDACreateSection - Create a PetscSection inside the DMDA that describes data layout. This allows multiple fields with
369 different numbers of dofs on vertices, cells, and faces in each direction.
370
371 Input Parameters:
372+ dm- The DMDA
373. numFields - The number of fields
374. numComp - The number of components in each field
375. numDof - The number of dofs per dimension for each field
376. numFaceDof - The number of dofs per face for each field and direction, or NULL
377
378 Level: developer
379
380 Note:
381 The default DMDA numbering is as follows:
382
383 - Cells: [0, nC)
384 - Vertices: [nC, nC+nV)
385 - X-Faces: [nC+nV, nC+nV+nXF) normal is +- x-dir
386 - Y-Faces: [nC+nV+nXF, nC+nV+nXF+nYF) normal is +- y-dir
387 - Z-Faces: [nC+nV+nXF+nYF, nC+nV+nXF+nYF+nZF) normal is +- z-dir
388
389 We interpret the default DMDA partition as a cell partition, and the data assignment as a cell assignment.
390@*/
391PetscErrorCode DMDACreateSection(DM dm, const PetscInt numComp[], const PetscInt numDof[], const PetscInt numFaceDof[], PetscSection *s)
392{
393 PetscSection section;
394 const PetscInt dim = dm->dim;
395 PetscInt numFields, numVertexTotDof = 0, numCellTotDof = 0, numFaceTotDof[3] = {0, 0, 0};
396 PetscBT isLeaf;
397 PetscSF sf;
398 PetscMPIInt rank;
399 const PetscMPIInt *neighbors;
400 PetscInt *localPoints;
401 PetscSFNode *remotePoints;
402 PetscInt nleaves = 0, nleavesCheck = 0, nL = 0;
403 PetscInt nC, nVx, nVy, nVz, nV, nxF, nXF, nyF, nYF, nzF, nZF;
404 PetscInt pStart, pEnd, cStart, cEnd, vStart, vEnd, fStart, fEnd, xfStart, xfEnd, yfStart, yfEnd, zfStart, zfEnd;
405 PetscInt f, v, c, xf, yf, zf, xn, yn, zn;
406 PetscErrorCode ierr;
407
408 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 408; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
409 PetscValidHeaderSpecificType(dm, DM_CLASSID, 1,DMDA)do { PetscErrorCode __ierr; PetscBool same; do { if (!dm) return
PetscError(((MPI_Comm)0x44000001),409,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,85,PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",1); if (
!PetscCheckPointer(dm,PETSC_OBJECT)) return PetscError(((MPI_Comm
)0x44000001),409,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,1); if (((PetscObject)(dm))->classid != DM_CLASSID) { if (
((PetscObject)(dm))->classid == -1) return PetscError(((MPI_Comm
)0x44000001),409,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,1); else return PetscError(((MPI_Comm)0x44000001),409,__func__
,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c",62
,PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d",1
); } } while (0); __ierr = PetscObjectTypeCompare((PetscObject
)dm,"da",&same);do {if (__builtin_expect(!!(__ierr),0)) return
PetscError(((MPI_Comm)0x44000001),409,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0); if (!same) return
PetscError(((MPI_Comm)0x44000001),409,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,62,PETSC_ERROR_INITIAL,"Wrong subtype object:Parameter # %d must have implementation %s it is %s"
,1,"da",((PetscObject)dm)->type_name); } while (0);
410 PetscValidPointer(s, 4)do { if (!s) return PetscError(((MPI_Comm)0x44000001),410,__func__
,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c",85
,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",4); if (!
PetscCheckPointer(s,PETSC_CHAR)) return PetscError(((MPI_Comm
)0x44000001),410,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer: Parameter # %d",4);
} while (0);
411 ierr = MPI_Comm_rank(PetscObjectComm((PetscObject)dm), &rank);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),411,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
412 ierr = DMDAGetNumCells(dm, NULL((void*)0), NULL((void*)0), NULL((void*)0), &nC);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),412,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
413 ierr = DMDAGetNumVertices(dm, &nVx, &nVy, &nVz, &nV);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),413,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
414 ierr = DMDAGetNumFaces(dm, &nxF, &nXF, &nyF, &nYF, &nzF, &nZF);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),414,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
415 ierr = DMDAGetHeightStratum(dm, -1, &pStart, &pEnd);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),415,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
416 ierr = DMDAGetHeightStratum(dm, 0, &cStart, &cEnd);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),416,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
417 ierr = DMDAGetHeightStratum(dm, 1, &fStart, &fEnd);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),417,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
418 ierr = DMDAGetHeightStratum(dm, dim, &vStart, &vEnd);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),418,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
419 xfStart = vEnd; xfEnd = xfStart+nXF;
420 yfStart = xfEnd; yfEnd = yfStart+nYF;
421 zfStart = yfEnd; zfEnd = zfStart+nZF;
422 if (zfEnd != fEnd) SETERRQ2(PetscObjectComm((PetscObject)dm), PETSC_ERR_PLIB, "Invalid face end %d, should be %d", zfEnd, fEnd)return PetscError(PetscObjectComm((PetscObject)dm),422,__func__
,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c",77
,PETSC_ERROR_INITIAL,"Invalid face end %d, should be %d",zfEnd
,fEnd);
423 /* Create local section */
424 ierr = DMDAGetInfo(dm, 0,0,0,0,0,0,0, &numFields, 0,0,0,0,0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),424,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
425 for (f = 0; f < numFields; ++f) {
426 numVertexTotDof += numDof[f*(dim+1)+0];
427 numCellTotDof += numDof[f*(dim+1)+dim];
428 numFaceTotDof[0] += dim > 0 ? (numFaceDof ? numFaceDof[f*dim+0] : numDof[f*(dim+1)+dim-1]) : 0;
429 numFaceTotDof[1] += dim > 1 ? (numFaceDof ? numFaceDof[f*dim+1] : numDof[f*(dim+1)+dim-1]) : 0;
430 numFaceTotDof[2] += dim > 2 ? (numFaceDof ? numFaceDof[f*dim+2] : numDof[f*(dim+1)+dim-1]) : 0;
431 }
432 ierr = PetscSectionCreate(PetscObjectComm((PetscObject)dm), &section);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),432,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
433 if (numFields > 0) {
434 ierr = PetscSectionSetNumFields(section, numFields);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),434,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
435 for (f = 0; f < numFields; ++f) {
436 const char *name;
437
438 ierr = DMDAGetFieldName(dm, f, &name);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),438,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
439 ierr = PetscSectionSetFieldName(section, f, name ? name : "Field");CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),439,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
440 if (numComp) {
441 ierr = PetscSectionSetFieldComponents(section, f, numComp[f]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),441,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
442 }
443 }
444 }
445 ierr = PetscSectionSetChart(section, pStart, pEnd);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),445,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
446 for (v = vStart; v < vEnd; ++v) {
447 for (f = 0; f < numFields; ++f) {
448 ierr = PetscSectionSetFieldDof(section, v, f, numDof[f*(dim+1)+0]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),448,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
449 }
450 ierr = PetscSectionSetDof(section, v, numVertexTotDof);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),450,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
451 }
452 for (xf = xfStart; xf < xfEnd; ++xf) {
453 for (f = 0; f < numFields; ++f) {
454 ierr = PetscSectionSetFieldDof(section, xf, f, numFaceDof ? numFaceDof[f*dim+0] : numDof[f*(dim+1)+dim-1]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),454,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
455 }
456 ierr = PetscSectionSetDof(section, xf, numFaceTotDof[0]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),456,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
457 }
458 for (yf = yfStart; yf < yfEnd; ++yf) {
459 for (f = 0; f < numFields; ++f) {
460 ierr = PetscSectionSetFieldDof(section, yf, f, numFaceDof ? numFaceDof[f*dim+1] : numDof[f*(dim+1)+dim-1]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),460,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
461 }
462 ierr = PetscSectionSetDof(section, yf, numFaceTotDof[1]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),462,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
463 }
464 for (zf = zfStart; zf < zfEnd; ++zf) {
465 for (f = 0; f < numFields; ++f) {
466 ierr = PetscSectionSetFieldDof(section, zf, f, numFaceDof ? numFaceDof[f*dim+2] : numDof[f*(dim+1)+dim-1]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),466,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
467 }
468 ierr = PetscSectionSetDof(section, zf, numFaceTotDof[2]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),468,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
469 }
470 for (c = cStart; c < cEnd; ++c) {
471 for (f = 0; f < numFields; ++f) {
472 ierr = PetscSectionSetFieldDof(section, c, f, numDof[f*(dim+1)+dim]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),472,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
473 }
474 ierr = PetscSectionSetDof(section, c, numCellTotDof);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),474,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
475 }
476 ierr = PetscSectionSetUp(section);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),476,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
477 /* Create mesh point SF */
478 ierr = PetscBTCreate(pEnd-pStart, &isLeaf);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),478,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
479 ierr = DMDAGetNeighbors(dm, &neighbors);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),479,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
480 for (zn = 0; zn < (dim > 2 ? 3 : 1); ++zn) {
481 for (yn = 0; yn < (dim > 1 ? 3 : 1); ++yn) {
482 for (xn = 0; xn < 3; ++xn) {
483 const PetscInt xp = xn-1, yp = dim > 1 ? yn-1 : 0, zp = dim > 2 ? zn-1 : 0;
484 const PetscInt neighbor = neighbors[(zn*3+yn)*3+xn];
485 PetscInt xv, yv, zv;
486
487 if (neighbor >= 0 && neighbor < rank) {
488 if (xp < 0) { /* left */
489 if (yp < 0) { /* bottom */
490 if (zp < 0) { /* back */
491 const PetscInt localVertex = ( 0*nVy + 0)*nVx + 0 + nC; /* left bottom back vertex */
492 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
493 } else if (zp > 0) { /* front */
494 const PetscInt localVertex = ((nVz-1)*nVy + 0)*nVx + 0 + nC; /* left bottom front vertex */
495 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
496 } else {
497 for (zv = 0; zv < nVz; ++zv) {
498 const PetscInt localVertex = (zv*nVy + 0)*nVx + 0 + nC; /* left bottom vertices */
499 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
500 }
501 }
502 } else if (yp > 0) { /* top */
503 if (zp < 0) { /* back */
504 const PetscInt localVertex = ( 0*nVy + nVy-1)*nVx + 0 + nC; /* left top back vertex */
505 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
506 } else if (zp > 0) { /* front */
507 const PetscInt localVertex = ((nVz-1)*nVy + nVy-1)*nVx + 0 + nC; /* left top front vertex */
508 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
509 } else {
510 for (zv = 0; zv < nVz; ++zv) {
511 const PetscInt localVertex = (zv*nVy + nVy-1)*nVx + 0 + nC; /* left top vertices */
512 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
513 }
514 }
515 } else {
516 if (zp < 0) { /* back */
517 for (yv = 0; yv < nVy; ++yv) {
518 const PetscInt localVertex = ( 0*nVy + yv)*nVx + 0 + nC; /* left back vertices */
519 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
520 }
521 } else if (zp > 0) { /* front */
522 for (yv = 0; yv < nVy; ++yv) {
523 const PetscInt localVertex = ((nVz-1)*nVy + yv)*nVx + 0 + nC; /* left front vertices */
524 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
525 }
526 } else {
527 for (zv = 0; zv < nVz; ++zv) {
528 for (yv = 0; yv < nVy; ++yv) {
529 const PetscInt localVertex = (zv*nVy + yv)*nVx + 0 + nC; /* left vertices */
530 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
531 }
532 }
533#if 0
534 for (xf = 0; xf < nxF; ++xf) {
535 /* THIS IS WRONG */
536 const PetscInt localFace = 0 + nC+nV; /* left faces */
537 if (!PetscBTLookupSet(isLeaf, localFace)) ++nleaves;
538 }
539#endif
540 }
541 }
542 } else if (xp > 0) { /* right */
543 if (yp < 0) { /* bottom */
544 if (zp < 0) { /* back */
545 const PetscInt localVertex = ( 0*nVy + 0)*nVx + nVx-1 + nC; /* right bottom back vertex */
546 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
547 } else if (zp > 0) { /* front */
548 const PetscInt localVertex = ((nVz-1)*nVy + 0)*nVx + nVx-1 + nC; /* right bottom front vertex */
549 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
550 } else {
551 for (zv = 0; zv < nVz; ++zv) {
552 const PetscInt localVertex = (zv*nVy + 0)*nVx + nVx-1 + nC; /* right bottom vertices */
553 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
554 }
555 }
556 } else if (yp > 0) { /* top */
557 if (zp < 0) { /* back */
558 const PetscInt localVertex = ( 0*nVy + nVy-1)*nVx + nVx-1 + nC; /* right top back vertex */
559 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
560 } else if (zp > 0) { /* front */
561 const PetscInt localVertex = ((nVz-1)*nVy + nVy-1)*nVx + nVx-1 + nC; /* right top front vertex */
562 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
563 } else {
564 for (zv = 0; zv < nVz; ++zv) {
565 const PetscInt localVertex = (zv*nVy + nVy-1)*nVx + nVx-1 + nC; /* right top vertices */
566 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
567 }
568 }
569 } else {
570 if (zp < 0) { /* back */
571 for (yv = 0; yv < nVy; ++yv) {
572 const PetscInt localVertex = ( 0*nVy + yv)*nVx + nVx-1 + nC; /* right back vertices */
573 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
574 }
575 } else if (zp > 0) { /* front */
576 for (yv = 0; yv < nVy; ++yv) {
577 const PetscInt localVertex = ((nVz-1)*nVy + yv)*nVx + nVx-1 + nC; /* right front vertices */
578 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
579 }
580 } else {
581 for (zv = 0; zv < nVz; ++zv) {
582 for (yv = 0; yv < nVy; ++yv) {
583 const PetscInt localVertex = (zv*nVy + yv)*nVx + nVx-1 + nC; /* right vertices */
584 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
585 }
586 }
587#if 0
588 for (xf = 0; xf < nxF; ++xf) {
589 /* THIS IS WRONG */
590 const PetscInt localFace = 0 + nC+nV; /* right faces */
591 if (!PetscBTLookupSet(isLeaf, localFace)) ++nleaves;
592 }
593#endif
594 }
595 }
596 } else {
597 if (yp < 0) { /* bottom */
598 if (zp < 0) { /* back */
599 for (xv = 0; xv < nVx; ++xv) {
600 const PetscInt localVertex = ( 0*nVy + 0)*nVx + xv + nC; /* bottom back vertices */
601 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
602 }
603 } else if (zp > 0) { /* front */
604 for (xv = 0; xv < nVx; ++xv) {
605 const PetscInt localVertex = ((nVz-1)*nVy + 0)*nVx + xv + nC; /* bottom front vertices */
606 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
607 }
608 } else {
609 for (zv = 0; zv < nVz; ++zv) {
610 for (xv = 0; xv < nVx; ++xv) {
611 const PetscInt localVertex = (zv*nVy + 0)*nVx + xv + nC; /* bottom vertices */
612 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
613 }
614 }
615#if 0
616 for (yf = 0; yf < nyF; ++yf) {
617 /* THIS IS WRONG */
618 const PetscInt localFace = 0 + nC+nV; /* bottom faces */
619 if (!PetscBTLookupSet(isLeaf, localVFace)) ++nleaves;
620 }
621#endif
622 }
623 } else if (yp > 0) { /* top */
624 if (zp < 0) { /* back */
625 for (xv = 0; xv < nVx; ++xv) {
626 const PetscInt localVertex = ( 0*nVy + nVy-1)*nVx + xv + nC; /* top back vertices */
627 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
628 }
629 } else if (zp > 0) { /* front */
630 for (xv = 0; xv < nVx; ++xv) {
631 const PetscInt localVertex = ((nVz-1)*nVy + nVy-1)*nVx + xv + nC; /* top front vertices */
632 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
633 }
634 } else {
635 for (zv = 0; zv < nVz; ++zv) {
636 for (xv = 0; xv < nVx; ++xv) {
637 const PetscInt localVertex = (zv*nVy + nVy-1)*nVx + xv + nC; /* top vertices */
638 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
639 }
640 }
641#if 0
642 for (yf = 0; yf < nyF; ++yf) {
643 /* THIS IS WRONG */
644 const PetscInt localFace = 0 + nC+nV; /* top faces */
645 if (!PetscBTLookupSet(isLeaf, localVFace)) ++nleaves;
646 }
647#endif
648 }
649 } else {
650 if (zp < 0) { /* back */
651 for (yv = 0; yv < nVy; ++yv) {
652 for (xv = 0; xv < nVx; ++xv) {
653 const PetscInt localVertex = ( 0*nVy + yv)*nVx + xv + nC; /* back vertices */
654 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
655 }
656 }
657#if 0
658 for (zf = 0; zf < nzF; ++zf) {
659 /* THIS IS WRONG */
660 const PetscInt localFace = 0 + nC+nV; /* back faces */
661 if (!PetscBTLookupSet(isLeaf, localFace)) ++nleaves;
662 }
663#endif
664 } else if (zp > 0) { /* front */
665 for (yv = 0; yv < nVy; ++yv) {
666 for (xv = 0; xv < nVx; ++xv) {
667 const PetscInt localVertex = ((nVz-1)*nVy + yv)*nVx + xv + nC; /* front vertices */
668 if (!PetscBTLookupSet(isLeaf, localVertex)) ++nleaves;
669 }
670 }
671#if 0
672 for (zf = 0; zf < nzF; ++zf) {
673 /* THIS IS WRONG */
674 const PetscInt localFace = 0 + nC+nV; /* front faces */
675 if (!PetscBTLookupSet(isLeaf, localFace)) ++nleaves;
676 }
677#endif
678 } else {
679 /* Nothing is shared from the interior */
680 }
681 }
682 }
683 }
684 }
685 }
686 }
687 ierr = PetscBTMemzero(pEnd-pStart, isLeaf);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),687,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
688 ierr = PetscMalloc1(nleaves,&localPoints)PetscMallocA(1,PETSC_FALSE,688,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,(size_t)(nleaves)*sizeof(**(&localPoints)),(&localPoints
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),688,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
689 ierr = PetscMalloc1(nleaves,&remotePoints)PetscMallocA(1,PETSC_FALSE,689,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,(size_t)(nleaves)*sizeof(**(&remotePoints)),(&remotePoints
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),689,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
690 for (zn = 0; zn < (dim > 2 ? 3 : 1); ++zn) {
691 for (yn = 0; yn < (dim > 1 ? 3 : 1); ++yn) {
692 for (xn = 0; xn < 3; ++xn) {
693 const PetscInt xp = xn-1, yp = dim > 1 ? yn-1 : 0, zp = dim > 2 ? zn-1 : 0;
694 const PetscInt neighbor = neighbors[(zn*3+yn)*3+xn];
695 PetscInt xv, yv, zv;
696
697 if (neighbor >= 0 && neighbor < rank) {
698 if (xp < 0) { /* left */
699 if (yp < 0) { /* bottom */
700 if (zp < 0) { /* back */
701 const PetscInt localVertex = ( 0*nVy + 0)*nVx + 0 + nC; /* left bottom back vertex */
702 const PetscInt remoteVertex = ((nVz-1)*nVy + nVy-1)*nVx + nVx-1 + nC; /* TODO: Correct this for neighbor sizes */
703
704 if (!PetscBTLookupSet(isLeaf, localVertex)) {
705 localPoints[nL] = localVertex;
706 remotePoints[nL].rank = neighbor;
707 remotePoints[nL].index = remoteVertex;
708 ++nL;
709 }
710 } else if (zp > 0) { /* front */
711 const PetscInt localVertex = ((nVz-1)*nVy + 0)*nVx + 0 + nC; /* left bottom front vertex */
712 const PetscInt remoteVertex = ( 0*nVy + nVy-1)*nVx + nVx-1 + nC; /* TODO: Correct this for neighbor sizes */
713
714 if (!PetscBTLookupSet(isLeaf, localVertex)) {
715 localPoints[nL] = localVertex;
716 remotePoints[nL].rank = neighbor;
717 remotePoints[nL].index = remoteVertex;
718 ++nL;
719 }
720 } else {
721 for (zv = 0; zv < nVz; ++zv) {
722 const PetscInt localVertex = (zv*nVy + 0)*nVx + 0 + nC; /* left bottom vertices */
723 const PetscInt remoteVertex = (zv*nVy + nVy-1)*nVx + nVx-1 + nC; /* TODO: Correct this for neighbor sizes */
724
725 if (!PetscBTLookupSet(isLeaf, localVertex)) {
726 localPoints[nL] = localVertex;
727 remotePoints[nL].rank = neighbor;
728 remotePoints[nL].index = remoteVertex;
729 ++nL;
730 }
731 }
732 }
733 } else if (yp > 0) { /* top */
734 if (zp < 0) { /* back */
735 const PetscInt localVertex = ( 0*nVy + nVy-1)*nVx + 0 + nC; /* left top back vertex */
736 const PetscInt remoteVertex = ((nVz-1)*nVy + 0)*nVx + nVx-1 + nC; /* TODO: Correct this for neighbor sizes */
737
738 if (!PetscBTLookupSet(isLeaf, localVertex)) {
739 localPoints[nL] = localVertex;
740 remotePoints[nL].rank = neighbor;
741 remotePoints[nL].index = remoteVertex;
742 ++nL;
743 }
744 } else if (zp > 0) { /* front */
745 const PetscInt localVertex = ((nVz-1)*nVy + nVy-1)*nVx + 0 + nC; /* left top front vertex */
746 const PetscInt remoteVertex = ( 0*nVy + 0)*nVx + nVx-1 + nC; /* TODO: Correct this for neighbor sizes */
747
748 if (!PetscBTLookupSet(isLeaf, localVertex)) {
749 localPoints[nL] = localVertex;
750 remotePoints[nL].rank = neighbor;
751 remotePoints[nL].index = remoteVertex;
752 ++nL;
753 }
754 } else {
755 for (zv = 0; zv < nVz; ++zv) {
756 const PetscInt localVertex = (zv*nVy + nVy-1)*nVx + 0 + nC; /* left top vertices */
757 const PetscInt remoteVertex = (zv*nVy + 0)*nVx + nVx-1 + nC; /* TODO: Correct this for neighbor sizes */
758
759 if (!PetscBTLookupSet(isLeaf, localVertex)) {
760 localPoints[nL] = localVertex;
761 remotePoints[nL].rank = neighbor;
762 remotePoints[nL].index = remoteVertex;
763 ++nL;
764 }
765 }
766 }
767 } else {
768 if (zp < 0) { /* back */
769 for (yv = 0; yv < nVy; ++yv) {
770 const PetscInt localVertex = ( 0*nVy + yv)*nVx + 0 + nC; /* left back vertices */
771 const PetscInt remoteVertex = ((nVz-1)*nVy + yv)*nVx + nVx-1 + nC; /* TODO: Correct this for neighbor sizes */
772
773 if (!PetscBTLookupSet(isLeaf, localVertex)) {
774 localPoints[nL] = localVertex;
775 remotePoints[nL].rank = neighbor;
776 remotePoints[nL].index = remoteVertex;
777 ++nL;
778 }
779 }
780 } else if (zp > 0) { /* front */
781 for (yv = 0; yv < nVy; ++yv) {
782 const PetscInt localVertex = ((nVz-1)*nVy + yv)*nVx + 0 + nC; /* left front vertices */
783 const PetscInt remoteVertex = ( 0*nVy + yv)*nVx + nVx-1 + nC; /* TODO: Correct this for neighbor sizes */
784
785 if (!PetscBTLookupSet(isLeaf, localVertex)) {
786 localPoints[nL] = localVertex;
787 remotePoints[nL].rank = neighbor;
788 remotePoints[nL].index = remoteVertex;
789 ++nL;
790 }
791 }
792 } else {
793 for (zv = 0; zv < nVz; ++zv) {
794 for (yv = 0; yv < nVy; ++yv) {
795 const PetscInt localVertex = (zv*nVy + yv)*nVx + 0 + nC; /* left vertices */
796 const PetscInt remoteVertex = (zv*nVy + yv)*nVx + nVx-1 + nC; /* TODO: Correct this for neighbor sizes */
797
798 if (!PetscBTLookupSet(isLeaf, localVertex)) {
799 localPoints[nL] = localVertex;
800 remotePoints[nL].rank = neighbor;
801 remotePoints[nL].index = remoteVertex;
802 ++nL;
803 }
804 }
805 }
806#if 0
807 for (xf = 0; xf < nxF; ++xf) {
808 /* THIS IS WRONG */
809 const PetscInt localFace = 0 + nC+nV; /* left faces */
810 const PetscInt remoteFace = 0 + nC+nV;
811
812 if (!PetscBTLookupSet(isLeaf, localFace)) {
813 localPoints[nL] = localFace;
814 remotePoints[nL].rank = neighbor;
815 remotePoints[nL].index = remoteFace;
816 }
817 }
818#endif
819 }
820 }
821 } else if (xp > 0) { /* right */
822 if (yp < 0) { /* bottom */
823 if (zp < 0) { /* back */
824 const PetscInt localVertex = ( 0*nVy + 0)*nVx + nVx-1 + nC; /* right bottom back vertex */
825 const PetscInt remoteVertex = ((nVz-1)*nVy + nVy-1)*nVx + 0 + nC; /* TODO: Correct this for neighbor sizes */
826
827 if (!PetscBTLookupSet(isLeaf, localVertex)) {
828 localPoints[nL] = localVertex;
829 remotePoints[nL].rank = neighbor;
830 remotePoints[nL].index = remoteVertex;
831 ++nL;
832 }
833 } else if (zp > 0) { /* front */
834 const PetscInt localVertex = ((nVz-1)*nVy + 0)*nVx + nVx-1 + nC; /* right bottom front vertex */
835 const PetscInt remoteVertex = ( 0*nVy + nVy-1)*nVx + 0 + nC; /* TODO: Correct this for neighbor sizes */
836
837 if (!PetscBTLookupSet(isLeaf, localVertex)) {
838 localPoints[nL] = localVertex;
839 remotePoints[nL].rank = neighbor;
840 remotePoints[nL].index = remoteVertex;
841 ++nL;
842 }
843 } else {
844 nleavesCheck += nVz;
845 for (zv = 0; zv < nVz; ++zv) {
846 const PetscInt localVertex = (zv*nVy + 0)*nVx + nVx-1 + nC; /* right bottom vertices */
847 const PetscInt remoteVertex = (zv*nVy + nVy-1)*nVx + 0 + nC; /* TODO: Correct this for neighbor sizes */
848
849 if (!PetscBTLookupSet(isLeaf, localVertex)) {
850 localPoints[nL] = localVertex;
851 remotePoints[nL].rank = neighbor;
852 remotePoints[nL].index = remoteVertex;
853 ++nL;
854 }
855 }
856 }
857 } else if (yp > 0) { /* top */
858 if (zp < 0) { /* back */
859 const PetscInt localVertex = ( 0*nVy + nVy-1)*nVx + nVx-1 + nC; /* right top back vertex */
860 const PetscInt remoteVertex = ((nVz-1)*nVy + 0)*nVx + 0 + nC; /* TODO: Correct this for neighbor sizes */
861
862 if (!PetscBTLookupSet(isLeaf, localVertex)) {
863 localPoints[nL] = localVertex;
864 remotePoints[nL].rank = neighbor;
865 remotePoints[nL].index = remoteVertex;
866 ++nL;
867 }
868 } else if (zp > 0) { /* front */
869 const PetscInt localVertex = ((nVz-1)*nVy + nVy-1)*nVx + nVx-1 + nC; /* right top front vertex */
870 const PetscInt remoteVertex = ( 0*nVy + 0)*nVx + 0 + nC; /* TODO: Correct this for neighbor sizes */
871
872 if (!PetscBTLookupSet(isLeaf, localVertex)) {
873 localPoints[nL] = localVertex;
874 remotePoints[nL].rank = neighbor;
875 remotePoints[nL].index = remoteVertex;
876 ++nL;
877 }
878 } else {
879 for (zv = 0; zv < nVz; ++zv) {
880 const PetscInt localVertex = (zv*nVy + nVy-1)*nVx + nVx-1 + nC; /* right top vertices */
881 const PetscInt remoteVertex = (zv*nVy + 0)*nVx + 0 + nC; /* TODO: Correct this for neighbor sizes */
882
883 if (!PetscBTLookupSet(isLeaf, localVertex)) {
884 localPoints[nL] = localVertex;
885 remotePoints[nL].rank = neighbor;
886 remotePoints[nL].index = remoteVertex;
887 ++nL;
888 }
889 }
890 }
891 } else {
892 if (zp < 0) { /* back */
893 for (yv = 0; yv < nVy; ++yv) {
894 const PetscInt localVertex = ( 0*nVy + yv)*nVx + nVx-1 + nC; /* right back vertices */
895 const PetscInt remoteVertex = ((nVz-1)*nVy + yv)*nVx + 0 + nC; /* TODO: Correct this for neighbor sizes */
896
897 if (!PetscBTLookupSet(isLeaf, localVertex)) {
898 localPoints[nL] = localVertex;
899 remotePoints[nL].rank = neighbor;
900 remotePoints[nL].index = remoteVertex;
901 ++nL;
902 }
903 }
904 } else if (zp > 0) { /* front */
905 for (yv = 0; yv < nVy; ++yv) {
906 const PetscInt localVertex = ((nVz-1)*nVy + yv)*nVx + nVx-1 + nC; /* right front vertices */
907 const PetscInt remoteVertex = ( 0*nVy + yv)*nVx + 0 + nC; /* TODO: Correct this for neighbor sizes */
908
909 if (!PetscBTLookupSet(isLeaf, localVertex)) {
910 localPoints[nL] = localVertex;
911 remotePoints[nL].rank = neighbor;
912 remotePoints[nL].index = remoteVertex;
913 ++nL;
914 }
915 }
916 } else {
917 for (zv = 0; zv < nVz; ++zv) {
918 for (yv = 0; yv < nVy; ++yv) {
919 const PetscInt localVertex = (zv*nVy + yv)*nVx + nVx-1 + nC; /* right vertices */
920 const PetscInt remoteVertex = (zv*nVy + yv)*nVx + 0 + nC; /* TODO: Correct this for neighbor sizes */
921
922 if (!PetscBTLookupSet(isLeaf, localVertex)) {
923 localPoints[nL] = localVertex;
924 remotePoints[nL].rank = neighbor;
925 remotePoints[nL].index = remoteVertex;
926 ++nL;
927 }
928 }
929 }
930#if 0
931 for (xf = 0; xf < nxF; ++xf) {
932 /* THIS IS WRONG */
933 const PetscInt localFace = 0 + nC+nV; /* right faces */
934 const PetscInt remoteFace = 0 + nC+nV;
935
936 if (!PetscBTLookupSet(isLeaf, localFace)) {
937 localPoints[nL] = localFace;
938 remotePoints[nL].rank = neighbor;
939 remotePoints[nL].index = remoteFace;
940 ++nL;
941 }
942 }
943#endif
944 }
945 }
946 } else {
947 if (yp < 0) { /* bottom */
948 if (zp < 0) { /* back */
949 for (xv = 0; xv < nVx; ++xv) {
950 const PetscInt localVertex = ( 0*nVy + 0)*nVx + xv + nC; /* bottom back vertices */
951 const PetscInt remoteVertex = ((nVz-1)*nVy + nVy-1)*nVx + xv + nC; /* TODO: Correct this for neighbor sizes */
952
953 if (!PetscBTLookupSet(isLeaf, localVertex)) {
954 localPoints[nL] = localVertex;
955 remotePoints[nL].rank = neighbor;
956 remotePoints[nL].index = remoteVertex;
957 ++nL;
958 }
959 }
960 } else if (zp > 0) { /* front */
961 for (xv = 0; xv < nVx; ++xv) {
962 const PetscInt localVertex = ((nVz-1)*nVy + 0)*nVx + xv + nC; /* bottom front vertices */
963 const PetscInt remoteVertex = ( 0*nVy + nVy-1)*nVx + xv + nC; /* TODO: Correct this for neighbor sizes */
964
965 if (!PetscBTLookupSet(isLeaf, localVertex)) {
966 localPoints[nL] = localVertex;
967 remotePoints[nL].rank = neighbor;
968 remotePoints[nL].index = remoteVertex;
969 ++nL;
970 }
971 }
972 } else {
973 for (zv = 0; zv < nVz; ++zv) {
974 for (xv = 0; xv < nVx; ++xv) {
975 const PetscInt localVertex = (zv*nVy + 0)*nVx + xv + nC; /* bottom vertices */
976 const PetscInt remoteVertex = (zv*nVy + nVy-1)*nVx + xv + nC; /* TODO: Correct this for neighbor sizes */
977
978 if (!PetscBTLookupSet(isLeaf, localVertex)) {
979 localPoints[nL] = localVertex;
980 remotePoints[nL].rank = neighbor;
981 remotePoints[nL].index = remoteVertex;
982 ++nL;
983 }
984 }
985 }
986#if 0
987 for (yf = 0; yf < nyF; ++yf) {
988 /* THIS IS WRONG */
989 const PetscInt localFace = 0 + nC+nV; /* bottom faces */
990 const PetscInt remoteFace = 0 + nC+nV;
991
992 if (!PetscBTLookupSet(isLeaf, localFace)) {
993 localPoints[nL] = localFace;
994 remotePoints[nL].rank = neighbor;
995 remotePoints[nL].index = remoteFace;
996 ++nL;
997 }
998 }
999#endif
1000 }
1001 } else if (yp > 0) { /* top */
1002 if (zp < 0) { /* back */
1003 for (xv = 0; xv < nVx; ++xv) {
1004 const PetscInt localVertex = ( 0*nVy + nVy-1)*nVx + xv + nC; /* top back vertices */
1005 const PetscInt remoteVertex = ((nVz-1)*nVy + 0)*nVx + xv + nC; /* TODO: Correct this for neighbor sizes */
1006
1007 if (!PetscBTLookupSet(isLeaf, localVertex)) {
1008 localPoints[nL] = localVertex;
1009 remotePoints[nL].rank = neighbor;
1010 remotePoints[nL].index = remoteVertex;
1011 ++nL;
1012 }
1013 }
1014 } else if (zp > 0) { /* front */
1015 for (xv = 0; xv < nVx; ++xv) {
1016 const PetscInt localVertex = ((nVz-1)*nVy + nVy-1)*nVx + xv + nC; /* top front vertices */
1017 const PetscInt remoteVertex = ( 0*nVy + 0)*nVx + xv + nC; /* TODO: Correct this for neighbor sizes */
1018
1019 if (!PetscBTLookupSet(isLeaf, localVertex)) {
1020 localPoints[nL] = localVertex;
1021 remotePoints[nL].rank = neighbor;
1022 remotePoints[nL].index = remoteVertex;
1023 ++nL;
1024 }
1025 }
1026 } else {
1027 for (zv = 0; zv < nVz; ++zv) {
1028 for (xv = 0; xv < nVx; ++xv) {
1029 const PetscInt localVertex = (zv*nVy + nVy-1)*nVx + xv + nC; /* top vertices */
1030 const PetscInt remoteVertex = (zv*nVy + 0)*nVx + xv + nC; /* TODO: Correct this for neighbor sizes */
1031
1032 if (!PetscBTLookupSet(isLeaf, localVertex)) {
1033 localPoints[nL] = localVertex;
1034 remotePoints[nL].rank = neighbor;
1035 remotePoints[nL].index = remoteVertex;
1036 ++nL;
1037 }
1038 }
1039 }
1040#if 0
1041 for (yf = 0; yf < nyF; ++yf) {
1042 /* THIS IS WRONG */
1043 const PetscInt localFace = 0 + nC+nV; /* top faces */
1044 const PetscInt remoteFace = 0 + nC+nV;
1045
1046 if (!PetscBTLookupSet(isLeaf, localFace)) {
1047 localPoints[nL] = localFace;
1048 remotePoints[nL].rank = neighbor;
1049 remotePoints[nL].index = remoteFace;
1050 ++nL;
1051 }
1052 }
1053#endif
1054 }
1055 } else {
1056 if (zp < 0) { /* back */
1057 for (yv = 0; yv < nVy; ++yv) {
1058 for (xv = 0; xv < nVx; ++xv) {
1059 const PetscInt localVertex = ( 0*nVy + yv)*nVx + xv + nC; /* back vertices */
1060 const PetscInt remoteVertex = ((nVz-1)*nVy + yv)*nVx + xv + nC; /* TODO: Correct this for neighbor sizes */
1061
1062 if (!PetscBTLookupSet(isLeaf, localVertex)) {
1063 localPoints[nL] = localVertex;
1064 remotePoints[nL].rank = neighbor;
1065 remotePoints[nL].index = remoteVertex;
1066 ++nL;
1067 }
1068 }
1069 }
1070#if 0
1071 for (zf = 0; zf < nzF; ++zf) {
1072 /* THIS IS WRONG */
1073 const PetscInt localFace = 0 + nC+nV; /* back faces */
1074 const PetscInt remoteFace = 0 + nC+nV;
1075
1076 if (!PetscBTLookupSet(isLeaf, localFace)) {
1077 localPoints[nL] = localFace;
1078 remotePoints[nL].rank = neighbor;
1079 remotePoints[nL].index = remoteFace;
1080 ++nL;
1081 }
1082 }
1083#endif
1084 } else if (zp > 0) { /* front */
1085 for (yv = 0; yv < nVy; ++yv) {
1086 for (xv = 0; xv < nVx; ++xv) {
1087 const PetscInt localVertex = ((nVz-1)*nVy + yv)*nVx + xv + nC; /* front vertices */
1088 const PetscInt remoteVertex = ( 0*nVy + yv)*nVx + xv + nC; /* TODO: Correct this for neighbor sizes */
1089
1090 if (!PetscBTLookupSet(isLeaf, localVertex)) {
1091 localPoints[nL] = localVertex;
1092 remotePoints[nL].rank = neighbor;
1093 remotePoints[nL].index = remoteVertex;
1094 ++nL;
1095 }
1096 }
1097 }
1098#if 0
1099 for (zf = 0; zf < nzF; ++zf) {
1100 /* THIS IS WRONG */
1101 const PetscInt localFace = 0 + nC+nV; /* front faces */
1102 const PetscInt remoteFace = 0 + nC+nV;
1103
1104 if (!PetscBTLookupSet(isLeaf, localFace)) {
1105 localPoints[nL] = localFace;
1106 remotePoints[nL].rank = neighbor;
1107 remotePoints[nL].index = remoteFace;
1108 ++nL;
1109 }
1110 }
1111#endif
1112 } else {
1113 /* Nothing is shared from the interior */
1114 }
1115 }
1116 }
1117 }
1118 }
1119 }
1120 }
1121 ierr = PetscBTDestroy(&isLeaf);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1121,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1122 /* Remove duplication in leaf determination */
1123 if (nleaves != nL) SETERRQ2(PetscObjectComm((PetscObject)dm), PETSC_ERR_PLIB, "The number of leaves %d did not match the number of remote leaves %d", nleaves, nleavesCheck)return PetscError(PetscObjectComm((PetscObject)dm),1123,__func__
,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c",77
,PETSC_ERROR_INITIAL,"The number of leaves %d did not match the number of remote leaves %d"
,nleaves,nleavesCheck);
1124 ierr = PetscSFCreate(PetscObjectComm((PetscObject)dm), &sf);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1124,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1125 ierr = PetscSFSetGraph(sf, pEnd, nleaves, localPoints, PETSC_OWN_POINTER, remotePoints, PETSC_OWN_POINTER);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1125,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1126 ierr = DMSetPointSF(dm, sf);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1126,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1127 ierr = PetscSFDestroy(&sf);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1127,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1128 *s = section;
1129 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
1130}
1131
1132PetscErrorCode DMDASetVertexCoordinates(DM dm, PetscReal xl, PetscReal xu, PetscReal yl, PetscReal yu, PetscReal zl, PetscReal zu)
1133{
1134 DM_DA *da = (DM_DA *) dm->data;
1135 Vec coordinates;
1136 PetscSection section;
1137 PetscScalar *coords;
1138 PetscReal h[3];
1139 PetscInt dim, size, M, N, P, nVx, nVy, nVz, nV, vStart, vEnd, v, i, j, k;
1140 PetscErrorCode ierr;
1141
1142 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 1142; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
1143 PetscValidHeaderSpecificType(dm, DM_CLASSID, 1,DMDA)do { PetscErrorCode __ierr; PetscBool same; do { if (!dm) return
PetscError(((MPI_Comm)0x44000001),1143,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,85,PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",1); if (
!PetscCheckPointer(dm,PETSC_OBJECT)) return PetscError(((MPI_Comm
)0x44000001),1143,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,1); if (((PetscObject)(dm))->classid != DM_CLASSID) { if (
((PetscObject)(dm))->classid == -1) return PetscError(((MPI_Comm
)0x44000001),1143,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,1); else return PetscError(((MPI_Comm)0x44000001),1143,__func__
,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c",62
,PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d",1
); } } while (0); __ierr = PetscObjectTypeCompare((PetscObject
)dm,"da",&same);do {if (__builtin_expect(!!(__ierr),0)) return
PetscError(((MPI_Comm)0x44000001),1143,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0); if (!same) return
PetscError(((MPI_Comm)0x44000001),1143,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,62,PETSC_ERROR_INITIAL,"Wrong subtype object:Parameter # %d must have implementation %s it is %s"
,1,"da",((PetscObject)dm)->type_name); } while (0);
1144 ierr = DMDAGetInfo(dm, &dim, &M, &N, &P, 0,0,0,0,0,0,0,0,0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1144,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1145 if (dim > 3) SETERRQ(PetscObjectComm((PetscObject)dm),PETSC_ERR_PLIB,"The following code only works for dim <= 3")return PetscError(PetscObjectComm((PetscObject)dm),1145,__func__
,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c",77
,PETSC_ERROR_INITIAL,"The following code only works for dim <= 3"
);
1146 h[0] = (xu - xl)/M;
1147 h[1] = (yu - yl)/N;
1148 h[2] = (zu - zl)/P;
1149 ierr = DMDAGetDepthStratum(dm, 0, &vStart, &vEnd);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1149,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1150 ierr = DMDAGetNumVertices(dm, &nVx, &nVy, &nVz, &nV);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1150,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1151 ierr = PetscSectionCreate(PetscObjectComm((PetscObject) dm), &section);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1151,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1152 ierr = PetscSectionSetNumFields(section, 1);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1152,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1153 ierr = PetscSectionSetFieldComponents(section, 0, dim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1153,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1154 ierr = PetscSectionSetChart(section, vStart, vEnd);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1154,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1155 for (v = vStart; v < vEnd; ++v) {
1156 ierr = PetscSectionSetDof(section, v, dim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1156,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1157 }
1158 ierr = PetscSectionSetUp(section);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1158,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1159 ierr = PetscSectionGetStorageSize(section, &size);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1159,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1160 ierr = VecCreateSeq(PETSC_COMM_SELF((MPI_Comm)0x44000001), size, &coordinates);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1160,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1161 ierr = PetscObjectSetName((PetscObject)coordinates,"coordinates");CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1161,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1162 ierr = VecGetArray(coordinates, &coords);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1162,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1163 for (k = 0; k < nVz; ++k) {
1164 PetscInt ind[3], d, off;
1165
1166 ind[0] = 0;
1167 ind[1] = 0;
1168 ind[2] = k + da->zs;
1169 for (j = 0; j < nVy; ++j) {
1170 ind[1] = j + da->ys;
1171 for (i = 0; i < nVx; ++i) {
1172 const PetscInt vertex = (k*nVy + j)*nVx + i + vStart;
1173
1174 ierr = PetscSectionGetOffset(section, vertex, &off);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1174,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1175 ind[0] = i + da->xs;
1176 for (d = 0; d < dim; ++d) {
1177 coords[off+d] = h[d]*ind[d];
1178 }
1179 }
1180 }
1181 }
1182 ierr = VecRestoreArray(coordinates, &coords);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1182,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1183 ierr = DMSetCoordinateSection(dm, PETSC_DETERMINE-1, section);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1183,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1184 ierr = DMSetCoordinatesLocal(dm, coordinates);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1184,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1185 ierr = PetscSectionDestroy(&section);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1185,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1186 ierr = VecDestroy(&coordinates);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1186,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1187 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
1188}
1189
1190PetscErrorCode DMProjectFunctionLocal_DA(DM dm, PetscReal time, PetscErrorCode (**funcs)(PetscInt, PetscReal, const PetscReal [], PetscInt, PetscScalar *, void *), void **ctxs, InsertMode mode, Vec localX)
1191{
1192 PetscDS prob;
1193 PetscFE fe;
1194 PetscDualSpace sp;
1195 PetscQuadrature q;
1196 PetscSection section;
1197 PetscScalar *values;
1198 PetscInt numFields, Nc, dim, dimEmbed, spDim, totDim, numValues, cStart, cEnd, f, c, v, d;
1199 PetscFEGeom *geom;
1200 PetscErrorCode ierr;
1201
1202 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 1202; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
1203 ierr = DMGetDS(dm, &prob);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1203,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1204 ierr = PetscDSGetTotalDimension(prob, &totDim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1204,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1205 ierr = PetscDSGetDiscretization(prob, 0, (PetscObject *) &fe);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1205,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1206 ierr = PetscFEGetQuadrature(fe, &q);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1206,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1207 ierr = DMGetSection(dm, &section);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1207,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1208 ierr = PetscSectionGetNumFields(section, &numFields);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1208,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1209 ierr = DMDAGetInfo(dm, &dim,0,0,0,0,0,0,0,0,0,0,0,0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1209,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1210 ierr = DMGetCoordinateDim(dm, &dimEmbed);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1210,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1211 ierr = DMDAGetHeightStratum(dm, 0, &cStart, &cEnd);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1211,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1212 ierr = DMDAVecGetClosure(dm, section, localX, cStart, &numValues, NULL((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1212,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1213 if (numValues != totDim) SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "The section cell closure size %d != dual space dimension %d", numValues, totDim)return PetscError(((MPI_Comm)0x44000001),1213,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,62,PETSC_ERROR_INITIAL,"The section cell closure size %d != dual space dimension %d"
,numValues,totDim);
1214 ierr = DMGetWorkArray(dm, numValues, MPIU_SCALAR((MPI_Datatype)0x4c00080b), &values);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1214,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1215 ierr = PetscFEGeomCreate(q,1,dimEmbed,PETSC_FALSE,&geom);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1215,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1216 for (c = cStart; c < cEnd; ++c) {
1217 ierr = DMDAComputeCellGeometryFEM(dm, c, q, geom->v, geom->J, NULL((void*)0), geom->detJ);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1217,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1218 for (f = 0, v = 0; f < numFields; ++f) {
1219 void * const ctx = ctxs ? ctxs[f] : NULL((void*)0);
1220
1221 ierr = PetscDSGetDiscretization(prob, f, (PetscObject *) &fe);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1221,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1222 ierr = PetscFEGetNumComponents(fe, &Nc);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1222,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1223 ierr = PetscFEGetDualSpace(fe, &sp);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1223,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1224 ierr = PetscDualSpaceGetDimension(sp, &spDim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1224,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1225 for (d = 0; d < spDim; ++d) {
1226 ierr = PetscDualSpaceApply(sp, d, time, geom, Nc, funcs[f], ctx, &values[v]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1226,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1227 }
1228 }
1229 ierr = DMDAVecSetClosure(dm, section, localX, c, values, mode);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1229,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1230 }
1231 ierr = PetscFEGeomDestroy(&geom);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1231,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1232 ierr = DMRestoreWorkArray(dm, numValues, MPIU_SCALAR((MPI_Datatype)0x4c00080b), &values);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1232,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1233 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
1234}
1235
1236PetscErrorCode DMComputeL2Diff_DA(DM dm, PetscReal time, PetscErrorCode (**funcs)(PetscInt, PetscReal, const PetscReal [], PetscInt, PetscScalar *, void *), void **ctxs, Vec X, PetscReal *diff)
1237{
1238 const PetscInt debug = 0;
1239 PetscDS prob;
1240 PetscFE fe;
1241 PetscQuadrature quad;
1242 PetscSection section;
1243 Vec localX;
1244 PetscScalar *funcVal;
1245 PetscReal *coords, *v0, *J, *invJ, detJ;
1246 PetscReal localDiff = 0.0;
1247 PetscInt dim, numFields, Nc, cStart, cEnd, c, field, fieldOffset, comp;
1248 PetscErrorCode ierr;
1249
1250 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 1250; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
1251 ierr = DMDAGetInfo(dm, &dim,0,0,0,0,0,0,0,0,0,0,0,0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1251,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1252 ierr = DMGetDS(dm, &prob);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1252,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1253 ierr = PetscDSGetTotalComponents(prob, &Nc);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1253,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1254 ierr = PetscDSGetDiscretization(prob, 0, (PetscObject *) &fe);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1254,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1255 ierr = PetscFEGetQuadrature(fe, &quad);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1255,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1256 ierr = DMGetSection(dm, &section);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1256,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1257 ierr = PetscSectionGetNumFields(section, &numFields);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1257,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1258 ierr = DMGetLocalVector(dm, &localX);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1258,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1259 ierr = DMGlobalToLocalBegin(dm, X, INSERT_VALUES, localX);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1259,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1260 ierr = DMGlobalToLocalEnd(dm, X, INSERT_VALUES, localX);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1260,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1261 /* There are no BC values in DAs right now: ierr = DMDAProjectFunctionLocal(dm, fe, funcs, INSERT_BC_VALUES, localX);CHKERRQ(ierr); */
1262 ierr = PetscMalloc5(Nc,&funcVal,dim,&coords,dim,&v0,dim*dim,&J,dim*dim,&invJ)PetscMallocA(5,PETSC_FALSE,1262,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,(size_t)(Nc)*sizeof(**(&funcVal)),(&funcVal),(size_t
)(dim)*sizeof(**(&coords)),(&coords),(size_t)(dim)*sizeof
(**(&v0)),(&v0),(size_t)(dim*dim)*sizeof(**(&J)),
(&J),(size_t)(dim*dim)*sizeof(**(&invJ)),(&invJ));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1262,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1263 ierr = DMDAGetHeightStratum(dm, 0, &cStart, &cEnd);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1263,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1264 for (c = cStart; c < cEnd; ++c) {
1265 PetscScalar *x = NULL((void*)0);
1266 PetscReal elemDiff = 0.0;
1267
1268 ierr = DMDAComputeCellGeometryFEM(dm, c, quad, v0, J, invJ, &detJ);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1268,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1269 if (detJ <= 0.0) SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid determinant %g for element %d", detJ, c)return PetscError(((MPI_Comm)0x44000001),1269,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,63,PETSC_ERROR_INITIAL,"Invalid determinant %g for element %d"
,detJ,c);
1270 ierr = DMDAVecGetClosure(dm, NULL((void*)0), localX, c, NULL((void*)0), &x);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1270,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1271
1272 for (field = 0, comp = 0, fieldOffset = 0; field < numFields; ++field) {
1273 void * const ctx = ctxs ? ctxs[field] : NULL((void*)0);
1274 const PetscReal *quadPoints, *quadWeights;
1275 PetscReal *basis;
1276 PetscInt qNc, Nq, Nb, Nc, q, d, e, fc, f;
1277
1278 ierr = PetscDSGetDiscretization(prob, field, (PetscObject *) &fe);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1278,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1279 ierr = PetscQuadratureGetData(quad, NULL((void*)0), &qNc, &Nq, &quadPoints, &quadWeights);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1279,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1280 ierr = PetscFEGetDimension(fe, &Nb);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1280,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1281 ierr = PetscFEGetNumComponents(fe, &Nc);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1281,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1282 if (qNc != Nc) SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_SUP, "Quadrature components %D != %D field components\n", qNc, Nc)return PetscError(((MPI_Comm)0x44000001),1282,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,56,PETSC_ERROR_INITIAL,"Quadrature components %D != %D field components\n"
,qNc,Nc);
1283 ierr = PetscFEGetDefaultTabulation(fe, &basis, NULL((void*)0), NULL((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1283,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1284 if (debug) {
1285 char title[1024];
1286 ierr = PetscSNPrintf(title, 1023, "Solution for Field %d", field);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1286,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1287 ierr = DMPrintCellVector(c, title, Nb, &x[fieldOffset]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1287,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1288 }
1289 for (q = 0; q < Nq; ++q) {
1290 for (d = 0; d < dim; d++) {
1291 coords[d] = v0[d];
1292 for (e = 0; e < dim; e++) {
1293 coords[d] += J[d*dim+e]*(quadPoints[q*dim+e] + 1.0);
1294 }
1295 }
1296 ierr = (*funcs[field])(dim, time, coords, numFields, funcVal, ctx);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1296,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1297 for (fc = 0; fc < Nc; ++fc) {
1298 PetscScalar interpolant = 0.0;
1299
1300 for (f = 0; f < Nb; ++f) {
1301 interpolant += x[fieldOffset+f]*basis[(q*Nb+f)*Nc+fc];
1302 }
1303 if (debug) {ierr = PetscPrintf(PETSC_COMM_SELF((MPI_Comm)0x44000001), " elem %d field %d diff %g\n", c, field, PetscSqr(PetscRealPart(interpolant - funcVal[fc]))(((interpolant - funcVal[fc]))*((interpolant - funcVal[fc])))*quadWeights[q*Nc+fc]*detJ);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1303,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}
1304 elemDiff += PetscSqr(PetscRealPart(interpolant - funcVal[fc]))(((interpolant - funcVal[fc]))*((interpolant - funcVal[fc])))*quadWeights[q*Nc+c]*detJ;
1305 }
1306 }
1307 comp += Nc;
1308 fieldOffset += Nb;
1309 }
1310 ierr = DMDAVecRestoreClosure(dm, NULL((void*)0), localX, c, NULL((void*)0), &x);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1310,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1311 if (debug) {ierr = PetscPrintf(PETSC_COMM_SELF((MPI_Comm)0x44000001), " elem %d diff %g\n", c, elemDiff);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1311,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}
1312 localDiff += elemDiff;
1313 }
1314 ierr = PetscFree5(funcVal,coords,v0,J,invJ)PetscFreeA(5,1314,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,&(funcVal),&(coords),&(v0),&(J),&(invJ));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1314,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1315 ierr = DMRestoreLocalVector(dm, &localX);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1315,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1316 ierr = MPIU_Allreduce(&localDiff, diff, 1, MPIU_REAL, MPIU_SUM, PetscObjectComm((PetscObject)dm))(PetscAllreduceBarrierCheck(PetscObjectComm((PetscObject)dm),
1,1316,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
) || ((petsc_allreduce_ct += PetscMPIParallelComm((PetscObjectComm
((PetscObject)dm))),0) || MPI_Allreduce((&localDiff),(diff
),(1),(((MPI_Datatype)0x4c00080b)),((MPI_Op)(0x58000003)),(PetscObjectComm
((PetscObject)dm)))));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1316,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1317 *diff = PetscSqrtReal(*diff)sqrt(*diff);
1318 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
1319}
1320
1321PetscErrorCode DMComputeL2GradientDiff_DA(DM dm, PetscReal time, PetscErrorCode (**funcs)(PetscInt, PetscReal, const PetscReal [], const PetscReal [], PetscInt, PetscScalar *, void *), void **ctxs, Vec X, const PetscReal n[], PetscReal *diff)
1322{
1323 const PetscInt debug = 0;
1324 PetscDS prob;
1325 PetscFE fe;
1326 PetscQuadrature quad;
1327 PetscSection section;
1328 Vec localX;
1329 PetscScalar *funcVal, *interpolantVec;
1330 PetscReal *coords, *realSpaceDer, *v0, *J, *invJ, detJ;
1331 PetscReal localDiff = 0.0;
1332 PetscInt dim, numFields, Nc, cStart, cEnd, c, field, fieldOffset, comp;
1333 PetscErrorCode ierr;
1334
1335 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 1335; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
1336 ierr = DMDAGetInfo(dm, &dim,0,0,0,0,0,0,0,0,0,0,0,0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1336,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1337 ierr = DMGetDS(dm, &prob);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1337,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1338 ierr = PetscDSGetTotalComponents(prob, &Nc);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1338,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1339 ierr = PetscDSGetDiscretization(prob, 0, (PetscObject *) &fe);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1339,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1340 ierr = PetscFEGetQuadrature(fe, &quad);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1340,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1341 ierr = DMGetSection(dm, &section);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1341,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1342 ierr = PetscSectionGetNumFields(section, &numFields);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1342,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1343 ierr = DMGetLocalVector(dm, &localX);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1343,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1344 ierr = DMGlobalToLocalBegin(dm, X, INSERT_VALUES, localX);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1344,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1345 ierr = DMGlobalToLocalEnd(dm, X, INSERT_VALUES, localX);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1345,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1346 /* There are no BC values in DAs right now: ierr = DMDAProjectFunctionLocal(dm, fe, funcs, INSERT_BC_VALUES, localX);CHKERRQ(ierr); */
1347 ierr = PetscMalloc7(Nc,&funcVal,dim,&coords,dim,&realSpaceDer,dim,&v0,dim*dim,&J,dim*dim,&invJ,dim,&interpolantVec)PetscMallocA(7,PETSC_FALSE,1347,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,(size_t)(Nc)*sizeof(**(&funcVal)),(&funcVal),(size_t
)(dim)*sizeof(**(&coords)),(&coords),(size_t)(dim)*sizeof
(**(&realSpaceDer)),(&realSpaceDer),(size_t)(dim)*sizeof
(**(&v0)),(&v0),(size_t)(dim*dim)*sizeof(**(&J)),
(&J),(size_t)(dim*dim)*sizeof(**(&invJ)),(&invJ),
(size_t)(dim)*sizeof(**(&interpolantVec)),(&interpolantVec
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1347,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1348 ierr = DMDAGetHeightStratum(dm, 0, &cStart, &cEnd);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1348,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1349 for (c = cStart; c < cEnd; ++c) {
1350 PetscScalar *x = NULL((void*)0);
1351 PetscReal elemDiff = 0.0;
1352
1353 ierr = DMDAComputeCellGeometryFEM(dm, c, quad, v0, J, invJ, &detJ);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1353,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1354 if (detJ <= 0.0) SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid determinant %g for element %d", detJ, c)return PetscError(((MPI_Comm)0x44000001),1354,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,63,PETSC_ERROR_INITIAL,"Invalid determinant %g for element %d"
,detJ,c);
1355 ierr = DMDAVecGetClosure(dm, NULL((void*)0), localX, c, NULL((void*)0), &x);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1355,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1356
1357 for (field = 0, comp = 0, fieldOffset = 0; field < numFields; ++field) {
1358 void * const ctx = ctxs ? ctxs[field] : NULL((void*)0);
1359 const PetscReal *quadPoints, *quadWeights;
1360 PetscReal *basisDer;
1361 PetscInt qNc, Nq, Nb, Nc, q, d, e, fc, f, g;
1362
1363 ierr = PetscDSGetDiscretization(prob, field, (PetscObject *) &fe);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1363,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1364 ierr = PetscQuadratureGetData(quad, NULL((void*)0), &qNc, &Nq, &quadPoints, &quadWeights);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1364,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1365 ierr = PetscFEGetDimension(fe, &Nb);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1365,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1366 ierr = PetscFEGetNumComponents(fe, &Nc);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1366,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1367 if (qNc != Nc) SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_SUP, "Quadrature components %D != %D field components\n", qNc, Nc)return PetscError(((MPI_Comm)0x44000001),1367,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,56,PETSC_ERROR_INITIAL,"Quadrature components %D != %D field components\n"
,qNc,Nc);
1368 ierr = PetscFEGetDefaultTabulation(fe, NULL((void*)0), &basisDer, NULL((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1368,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1369 if (debug) {
1370 char title[1024];
1371 ierr = PetscSNPrintf(title, 1023, "Solution for Field %d", field);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1371,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1372 ierr = DMPrintCellVector(c, title, Nb*Nc, &x[fieldOffset]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1372,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1373 }
1374 for (q = 0; q < Nq; ++q) {
1375 for (d = 0; d < dim; d++) {
1376 coords[d] = v0[d];
1377 for (e = 0; e < dim; e++) {
1378 coords[d] += J[d*dim+e]*(quadPoints[q*dim+e] + 1.0);
1379 }
1380 }
1381 ierr = (*funcs[field])(dim, time, coords, n, numFields, funcVal, ctx);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1381,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1382 for (fc = 0; fc < Nc; ++fc) {
1383 PetscScalar interpolant = 0.0;
1384
1385 for (d = 0; d < dim; ++d) interpolantVec[d] = 0.0;
1386 for (f = 0; f < Nb; ++f) {
1387 for (d = 0; d < dim; ++d) {
1388 realSpaceDer[d] = 0.0;
1389 for (g = 0; g < dim; ++g) {
1390 realSpaceDer[d] += invJ[g*dim+d]*basisDer[(q*Nb+f)*dim+g];
1391 }
1392 interpolantVec[d] += x[fieldOffset+f]*realSpaceDer[d];
1393 }
1394 }
1395 for (d = 0; d < dim; ++d) interpolant += interpolantVec[d]*n[d];
1396 if (debug) {ierr = PetscPrintf(PETSC_COMM_SELF((MPI_Comm)0x44000001), " elem %d fieldDer %d diff %g\n", c, field, PetscSqr(PetscRealPart(interpolant - funcVal[fc]))(((interpolant - funcVal[fc]))*((interpolant - funcVal[fc])))*quadWeights[q*Nc+c]*detJ);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1396,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}
1397 elemDiff += PetscSqr(PetscRealPart(interpolant - funcVal[fc]))(((interpolant - funcVal[fc]))*((interpolant - funcVal[fc])))*quadWeights[q*Nc+c]*detJ;
1398 }
1399 }
1400 comp += Nc;
1401 fieldOffset += Nb*Nc;
1402 }
1403 ierr = DMDAVecRestoreClosure(dm, NULL((void*)0), localX, c, NULL((void*)0), &x);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1403,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1404 if (debug) {ierr = PetscPrintf(PETSC_COMM_SELF((MPI_Comm)0x44000001), " elem %d diff %g\n", c, elemDiff);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1404,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}
1405 localDiff += elemDiff;
1406 }
1407 ierr = PetscFree7(funcVal,coords,realSpaceDer,v0,J,invJ,interpolantVec)PetscFreeA(7,1407,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,&(funcVal),&(coords),&(realSpaceDer),&(v0),&
(J),&(invJ),&(interpolantVec));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1407,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1408 ierr = DMRestoreLocalVector(dm, &localX);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1408,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1409 ierr = MPIU_Allreduce(&localDiff, diff, 1, MPIU_REAL, MPIU_SUM, PetscObjectComm((PetscObject)dm))(PetscAllreduceBarrierCheck(PetscObjectComm((PetscObject)dm),
1,1409,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
) || ((petsc_allreduce_ct += PetscMPIParallelComm((PetscObjectComm
((PetscObject)dm))),0) || MPI_Allreduce((&localDiff),(diff
),(1),(((MPI_Datatype)0x4c00080b)),((MPI_Op)(0x58000003)),(PetscObjectComm
((PetscObject)dm)))));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1409,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1410 *diff = PetscSqrtReal(*diff)sqrt(*diff);
1411 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
1412}
1413
1414/* ------------------------------------------------------------------- */
1415
1416/*@C
1417 DMDAGetArray - Gets a work array for a DMDA
1418
1419 Input Parameter:
1420+ da - information about my local patch
1421- ghosted - do you want arrays for the ghosted or nonghosted patch
1422
1423 Output Parameters:
1424. vptr - array data structured
1425
1426 Note: The vector values are NOT initialized and may have garbage in them, so you may need
1427 to zero them.
1428
1429 Level: advanced
1430
1431.seealso: DMDARestoreArray()
1432
1433@*/
1434PetscErrorCode DMDAGetArray(DM da,PetscBool ghosted,void *vptr)
1435{
1436 PetscErrorCode ierr;
1437 PetscInt j,i,xs,ys,xm,ym,zs,zm;
1438 char *iarray_start;
1439 void **iptr = (void**)vptr;
1440 DM_DA *dd = (DM_DA*)da->data;
1441
1442 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 1442; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
1443 PetscValidHeaderSpecificType(da,DM_CLASSID,1,DMDA)do { PetscErrorCode __ierr; PetscBool same; do { if (!da) return
PetscError(((MPI_Comm)0x44000001),1443,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,85,PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",1); if (
!PetscCheckPointer(da,PETSC_OBJECT)) return PetscError(((MPI_Comm
)0x44000001),1443,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,1); if (((PetscObject)(da))->classid != DM_CLASSID) { if (
((PetscObject)(da))->classid == -1) return PetscError(((MPI_Comm
)0x44000001),1443,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,1); else return PetscError(((MPI_Comm)0x44000001),1443,__func__
,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c",62
,PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d",1
); } } while (0); __ierr = PetscObjectTypeCompare((PetscObject
)da,"da",&same);do {if (__builtin_expect(!!(__ierr),0)) return
PetscError(((MPI_Comm)0x44000001),1443,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0); if (!same) return
PetscError(((MPI_Comm)0x44000001),1443,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,62,PETSC_ERROR_INITIAL,"Wrong subtype object:Parameter # %d must have implementation %s it is %s"
,1,"da",((PetscObject)da)->type_name); } while (0);
1444 if (ghosted) {
1445 for (i=0; i<DMDA_MAX_WORK_ARRAYS2; i++) {
1446 if (dd->arrayghostedin[i]) {
1447 *iptr = dd->arrayghostedin[i];
1448 iarray_start = (char*)dd->startghostedin[i];
1449 dd->arrayghostedin[i] = NULL((void*)0);
1450 dd->startghostedin[i] = NULL((void*)0);
1451
1452 goto done;
1453 }
1454 }
1455 xs = dd->Xs;
1456 ys = dd->Ys;
1457 zs = dd->Zs;
1458 xm = dd->Xe-dd->Xs;
1459 ym = dd->Ye-dd->Ys;
1460 zm = dd->Ze-dd->Zs;
1461 } else {
1462 for (i=0; i<DMDA_MAX_WORK_ARRAYS2; i++) {
1463 if (dd->arrayin[i]) {
1464 *iptr = dd->arrayin[i];
1465 iarray_start = (char*)dd->startin[i];
1466 dd->arrayin[i] = NULL((void*)0);
1467 dd->startin[i] = NULL((void*)0);
1468
1469 goto done;
1470 }
1471 }
1472 xs = dd->xs;
1473 ys = dd->ys;
1474 zs = dd->zs;
1475 xm = dd->xe-dd->xs;
1476 ym = dd->ye-dd->ys;
1477 zm = dd->ze-dd->zs;
1478 }
1479
1480 switch (da->dim) {
1481 case 1: {
1482 void *ptr;
1483
1484 ierr = PetscMalloc(xm*sizeof(PetscScalar),&iarray_start)((*PetscTrMalloc)((xm*sizeof(PetscScalar)),1484,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,(void**)(&iarray_start)));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1484,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1485
1486 ptr = (void*)(iarray_start - xs*sizeof(PetscScalar));
1487 *iptr = (void*)ptr;
1488 break;
1489 }
1490 case 2: {
1491 void **ptr;
1492
1493 ierr = PetscMalloc((ym+1)*sizeof(void*)+xm*ym*sizeof(PetscScalar),&iarray_start)((*PetscTrMalloc)(((ym+1)*sizeof(void*)+xm*ym*sizeof(PetscScalar
)),1493,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,(void**)(&iarray_start)));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1493,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1494
1495 ptr = (void**)(iarray_start + xm*ym*sizeof(PetscScalar) - ys*sizeof(void*));
1496 for (j=ys; j<ys+ym; j++) ptr[j] = iarray_start + sizeof(PetscScalar)*(xm*(j-ys) - xs);
1497 *iptr = (void*)ptr;
1498 break;
1499 }
1500 case 3: {
1501 void ***ptr,**bptr;
1502
1503 ierr = PetscMalloc((zm+1)*sizeof(void**)+(ym*zm+1)*sizeof(void*)+xm*ym*zm*sizeof(PetscScalar),&iarray_start)((*PetscTrMalloc)(((zm+1)*sizeof(void**)+(ym*zm+1)*sizeof(void
*)+xm*ym*zm*sizeof(PetscScalar)),1503,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,(void**)(&iarray_start)));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1503,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
1504
1505 ptr = (void***)(iarray_start + xm*ym*zm*sizeof(PetscScalar) - zs*sizeof(void*));
1506 bptr = (void**)(iarray_start + xm*ym*zm*sizeof(PetscScalar) + zm*sizeof(void**));
1507 for (i=zs; i<zs+zm; i++) ptr[i] = bptr + ((i-zs)*ym - ys);
1508 for (i=zs; i<zs+zm; i++) {
1509 for (j=ys; j<ys+ym; j++) {
1510 ptr[i][j] = iarray_start + sizeof(PetscScalar)*(xm*ym*(i-zs) + xm*(j-ys) - xs);
1511 }
1512 }
1513 *iptr = (void*)ptr;
1514 break;
1515 }
1516 default:
1517 SETERRQ1(PetscObjectComm((PetscObject)da),PETSC_ERR_SUP,"Dimension %D not supported",da->dim)return PetscError(PetscObjectComm((PetscObject)da),1517,__func__
,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c",56
,PETSC_ERROR_INITIAL,"Dimension %D not supported",da->dim);
1518 }
1519
1520done:
1521 /* add arrays to the checked out list */
1522 if (ghosted) {
1523 for (i=0; i<DMDA_MAX_WORK_ARRAYS2; i++) {
1524 if (!dd->arrayghostedout[i]) {
1525 dd->arrayghostedout[i] = *iptr;
1526 dd->startghostedout[i] = iarray_start;
1527 break;
1528 }
1529 }
1530 } else {
1531 for (i=0; i<DMDA_MAX_WORK_ARRAYS2; i++) {
1532 if (!dd->arrayout[i]) {
1533 dd->arrayout[i] = *iptr;
1534 dd->startout[i] = iarray_start;
1535 break;
1536 }
1537 }
1538 }
1539 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
1540}
1541
1542/*@C
1543 DMDARestoreArray - Restores an array of derivative types for a DMDA
1544
1545 Input Parameter:
1546+ da - information about my local patch
1547. ghosted - do you want arrays for the ghosted or nonghosted patch
1548- vptr - array data structured to be passed to ad_FormFunctionLocal()
1549
1550 Level: advanced
1551
1552.seealso: DMDAGetArray()
1553
1554@*/
1555PetscErrorCode DMDARestoreArray(DM da,PetscBool ghosted,void *vptr)
1556{
1557 PetscInt i;
1558 void **iptr = (void**)vptr,*iarray_start = 0;
1559 DM_DA *dd = (DM_DA*)da->data;
1560
1561 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c";
petscstack->line[petscstack->currentsize] = 1561; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
1562 PetscValidHeaderSpecificType(da,DM_CLASSID,1,DMDA)do { PetscErrorCode __ierr; PetscBool same; do { if (!da) return
PetscError(((MPI_Comm)0x44000001),1562,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,85,PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",1); if (
!PetscCheckPointer(da,PETSC_OBJECT)) return PetscError(((MPI_Comm
)0x44000001),1562,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,1); if (((PetscObject)(da))->classid != DM_CLASSID) { if (
((PetscObject)(da))->classid == -1) return PetscError(((MPI_Comm
)0x44000001),1562,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,1); else return PetscError(((MPI_Comm)0x44000001),1562,__func__
,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c",62
,PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d",1
); } } while (0); __ierr = PetscObjectTypeCompare((PetscObject
)da,"da",&same);do {if (__builtin_expect(!!(__ierr),0)) return
PetscError(((MPI_Comm)0x44000001),1562,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0); if (!same) return
PetscError(((MPI_Comm)0x44000001),1562,__func__,"/sandbox/petsc/petsc.next-tmp/src/dm/impls/da/dalocal.c"
,62,PETSC_ERROR_INITIAL,"Wrong subtype object:Parameter # %d must have implementation %s it is %s"
,1,"da",((PetscObject)da)->type_name); } while (0);
1563 if (ghosted) {
1564 for (i=0; i<DMDA_MAX_WORK_ARRAYS2; i++) {
1565 if (dd->arrayghostedout[i] == *iptr) {
1566 iarray_start = dd->startghostedout[i];
1567 dd->arrayghostedout[i] = NULL((void*)0);
1568 dd->startghostedout[i] = NULL((void*)0);
1569 break;
1570 }
1571 }
1572 for (i=0; i<DMDA_MAX_WORK_ARRAYS2; i++) {
1573 if (!dd->arrayghostedin[i]) {
1574 dd->arrayghostedin[i] = *iptr;
1575 dd->startghostedin[i] = iarray_start;
1576 break;
1577 }
1578 }
1579 } else {
1580 for (i=0; i<DMDA_MAX_WORK_ARRAYS2; i++) {
1581 if (dd->arrayout[i] == *iptr) {
1582 iarray_start = dd->startout[i];
1583 dd->arrayout[i] = NULL((void*)0);
1584 dd->startout[i] = NULL((void*)0);
1585 break;
1586 }
1587 }
1588 for (i=0; i<DMDA_MAX_WORK_ARRAYS2; i++) {
1589 if (!dd->arrayin[i]) {
1590 dd->arrayin[i] = *iptr;
1591 dd->startin[i] = iarray_start;
1592 break;
1593 }
1594 }
1595 }
1596 PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);
1597}
1598