Bug Summary

File:dm/dt/dualspace/impls/lagrange/dspacelagrange.c
Warning:line 71, column 10
5th function call argument is an uninitialized value

Annotated Source Code

[?] Use j/k keys for keyboard navigation

1#include <petsc/private/petscfeimpl.h> /*I "petscfe.h" I*/
2#include <petscdmplex.h>
3
4static PetscErrorCode PetscDualSpaceDestroy_Lagrange(PetscDualSpace sp)
5{
6 PetscDualSpace_Lag *lag = (PetscDualSpace_Lag *) sp->data;
7 PetscInt i;
8 PetscErrorCode ierr;
9
10 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 10; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
11 if (lag->symmetries) {
12 PetscInt **selfSyms = lag->symmetries[0];
13
14 if (selfSyms) {
15 PetscInt i, **allocated = &selfSyms[-lag->selfSymOff];
16
17 for (i = 0; i < lag->numSelfSym; i++) {
18 ierr = PetscFree(allocated[i])((*PetscTrFree)((void*)(allocated[i]),18,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
) || ((allocated[i]) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),18,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
19 }
20 ierr = PetscFree(allocated)((*PetscTrFree)((void*)(allocated),20,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
) || ((allocated) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),20,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
21 }
22 ierr = PetscFree(lag->symmetries)((*PetscTrFree)((void*)(lag->symmetries),22,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
) || ((lag->symmetries) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),22,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
23 }
24 for (i = 0; i < lag->height; i++) {
25 ierr = PetscDualSpaceDestroy(&lag->subspaces[i]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),25,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
26 }
27 ierr = PetscFree(lag->subspaces)((*PetscTrFree)((void*)(lag->subspaces),27,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
) || ((lag->subspaces) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),27,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
28 ierr = PetscFree(lag->numDof)((*PetscTrFree)((void*)(lag->numDof),28,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
) || ((lag->numDof) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),28,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
29 ierr = PetscFree(lag)((*PetscTrFree)((void*)(lag),29,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
) || ((lag) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),29,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
30 ierr = PetscObjectComposeFunction((PetscObject) sp, "PetscDualSpaceLagrangeGetContinuity_C", NULL)PetscObjectComposeFunction_Private((PetscObject) sp,"PetscDualSpaceLagrangeGetContinuity_C"
,(PetscVoidFunction)(((void*)0)));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),30,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
31 ierr = PetscObjectComposeFunction((PetscObject) sp, "PetscDualSpaceLagrangeSetContinuity_C", NULL)PetscObjectComposeFunction_Private((PetscObject) sp,"PetscDualSpaceLagrangeSetContinuity_C"
,(PetscVoidFunction)(((void*)0)));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),31,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
32 ierr = PetscObjectComposeFunction((PetscObject) sp, "PetscDualSpaceLagrangeGetTensor_C", NULL)PetscObjectComposeFunction_Private((PetscObject) sp,"PetscDualSpaceLagrangeGetTensor_C"
,(PetscVoidFunction)(((void*)0)));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),32,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
33 ierr = PetscObjectComposeFunction((PetscObject) sp, "PetscDualSpaceLagrangeSetTensor_C", NULL)PetscObjectComposeFunction_Private((PetscObject) sp,"PetscDualSpaceLagrangeSetTensor_C"
,(PetscVoidFunction)(((void*)0)));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),33,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
34 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);
35}
36
37static PetscErrorCode PetscDualSpaceLagrangeView_Ascii(PetscDualSpace sp, PetscViewer viewer)
38{
39 PetscDualSpace_Lag *lag = (PetscDualSpace_Lag *) sp->data;
40 PetscErrorCode ierr;
41
42 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 42; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
43 ierr = PetscViewerASCIIPrintf(viewer, "%s %sLagrange dual space\n", lag->continuous ? "Continuous" : "Discontinuous", lag->tensorSpace ? "tensor " : "");CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),43,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
44 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);
45}
46
47static PetscErrorCode PetscDualSpaceView_Lagrange(PetscDualSpace sp, PetscViewer viewer)
48{
49 PetscBool iascii;
50 PetscErrorCode ierr;
51
52 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 52; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
53 PetscValidHeaderSpecific(sp, PETSCDUALSPACE_CLASSID, 1)do { if (!sp) return PetscError(((MPI_Comm)0x44000001),53,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,85,PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",1); if (
!PetscCheckPointer(sp,PETSC_OBJECT)) return PetscError(((MPI_Comm
)0x44000001),53,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,1); if (((PetscObject)(sp))->classid != PETSCDUALSPACE_CLASSID
) { if (((PetscObject)(sp))->classid == -1) return PetscError
(((MPI_Comm)0x44000001),53,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,1); else return PetscError(((MPI_Comm)0x44000001),53,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,62,PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d"
,1); } } while (0);
54 PetscValidHeaderSpecific(viewer, PETSC_VIEWER_CLASSID, 2)do { if (!viewer) return PetscError(((MPI_Comm)0x44000001),54
,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,85,PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",2); if (
!PetscCheckPointer(viewer,PETSC_OBJECT)) return PetscError(((
MPI_Comm)0x44000001),54,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,2); if (((PetscObject)(viewer))->classid != PETSC_VIEWER_CLASSID
) { if (((PetscObject)(viewer))->classid == -1) return PetscError
(((MPI_Comm)0x44000001),54,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,2); else return PetscError(((MPI_Comm)0x44000001),54,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,62,PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d"
,2); } } while (0);
55 ierr = PetscObjectTypeCompare((PetscObject) viewer, PETSCVIEWERASCII"ascii", &iascii);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),55,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
56 if (iascii) {ierr = PetscDualSpaceLagrangeView_Ascii(sp, viewer);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),56,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}
57 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);
58}
59
60static PetscErrorCode PetscDualSpaceSetFromOptions_Lagrange(PetscOptionItems *PetscOptionsObject,PetscDualSpace sp)
61{
62 PetscBool continuous, tensor, flg;
1
'tensor' declared without an initial value
→
63 PetscErrorCode ierr;
64
65 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 65; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
66 ierr = PetscDualSpaceLagrangeGetContinuity(sp, &continuous);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),66,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
67 ierr = PetscDualSpaceLagrangeGetTensor(sp, &tensor);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),67,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
2
←
Calling 'PetscDualSpaceLagrangeGetTensor'
→
9
←
Returning from 'PetscDualSpaceLagrangeGetTensor'
→
68 ierr = PetscOptionsHead(PetscOptionsObject,"PetscDualSpace Lagrange Options");CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),68,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
69 ierr = PetscOptionsBool("-petscdualspace_lagrange_continuity", "Flag for continuous element", "PetscDualSpaceLagrangeSetContinuity", continuous, &continuous, &flg)PetscOptionsBool_Private(PetscOptionsObject,"-petscdualspace_lagrange_continuity"
,"Flag for continuous element","PetscDualSpaceLagrangeSetContinuity"
,continuous,&continuous,&flg);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),69,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
70 if (flg) {ierr = PetscDualSpaceLagrangeSetContinuity(sp, continuous);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),70,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}
10
←
Assuming 'flg' is 0
→
11
←
Taking false branch
→
71 ierr = PetscOptionsBool("-petscdualspace_lagrange_tensor", "Flag for tensor dual space", "PetscDualSpaceLagrangeSetContinuity", tensor, &tensor, &flg)PetscOptionsBool_Private(PetscOptionsObject,"-petscdualspace_lagrange_tensor"
,"Flag for tensor dual space","PetscDualSpaceLagrangeSetContinuity"
,tensor,&tensor,&flg);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),71,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
12
←
Within the expansion of the macro 'PetscOptionsBool':
a
5th function call argument is an uninitialized value
72 if (flg) {ierr = PetscDualSpaceLagrangeSetTensor(sp, tensor);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),72,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}
73 ierr = PetscOptionsTail()0; do {if (PetscOptionsObject->count != 1) 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);} while(0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),73,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
74 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);
75}
76
77static PetscErrorCode PetscDualSpaceDuplicate_Lagrange(PetscDualSpace sp, PetscDualSpace *spNew)
78{
79 PetscInt order, Nc;
80 PetscBool cont, tensor;
81 const char *name;
82 PetscErrorCode ierr;
83
84 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 84; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
85 ierr = PetscDualSpaceCreate(PetscObjectComm((PetscObject) sp), spNew);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),85,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
86 ierr = PetscObjectGetName((PetscObject) sp, &name);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),86,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
87 ierr = PetscObjectSetName((PetscObject) *spNew, name);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),87,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
88 ierr = PetscDualSpaceSetType(*spNew, PETSCDUALSPACELAGRANGE"lagrange");CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),88,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
89 ierr = PetscDualSpaceGetOrder(sp, &order);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),89,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
90 ierr = PetscDualSpaceSetOrder(*spNew, order);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),90,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
91 ierr = PetscDualSpaceGetNumComponents(sp, &Nc);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),91,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
92 ierr = PetscDualSpaceSetNumComponents(*spNew, Nc);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),92,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
93 ierr = PetscDualSpaceLagrangeGetContinuity(sp, &cont);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),93,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
94 ierr = PetscDualSpaceLagrangeSetContinuity(*spNew, cont);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),94,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
95 ierr = PetscDualSpaceLagrangeGetTensor(sp, &tensor);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),95,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
96 ierr = PetscDualSpaceLagrangeSetTensor(*spNew, tensor);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),96,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
97 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);
98}
99
100static PetscErrorCode PetscDualSpaceGetDimension_SingleCell_Lagrange(PetscDualSpace sp, PetscInt order, PetscInt *dim)
101{
102 PetscDualSpace_Lag *lag = (PetscDualSpace_Lag *) sp->data;
103 PetscReal D = 1.0;
104 PetscInt n, d;
105 PetscErrorCode ierr;
106
107 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 107; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
108 *dim = -1;
109 ierr = DMGetDimension(sp->dm, &n);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),109,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
110 if (!lag->tensorSpace) {
111 for (d = 1; d <= n; ++d) {
112 D *= ((PetscReal) (order+d))/d;
113 }
114 *dim = (PetscInt) (D + 0.5);
115 } else {
116 *dim = 1;
117 for (d = 0; d < n; ++d) *dim *= (order+1);
118 }
119 *dim *= sp->Nc;
120 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);
121}
122
123static PetscErrorCode PetscDualSpaceCreateHeightSubspace_Lagrange(PetscDualSpace sp, PetscInt height, PetscDualSpace *bdsp)
124{
125 PetscDualSpace_Lag *lag = (PetscDualSpace_Lag *) sp->data;
126 PetscBool continuous, tensor;
127 PetscInt order;
128 PetscErrorCode ierr;
129
130 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 130; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
131 ierr = PetscDualSpaceLagrangeGetContinuity(sp,&continuous);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),131,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
132 ierr = PetscDualSpaceGetOrder(sp,&order);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),132,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
133 if (height == 0) {
134 ierr = PetscObjectReference((PetscObject)sp);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),134,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
135 *bdsp = sp;
136 } else if (continuous == PETSC_FALSE || !order) {
137 *bdsp = NULL((void*)0);
138 } else {
139 DM dm, K;
140 PetscInt dim;
141
142 ierr = PetscDualSpaceGetDM(sp,&dm);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),142,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
143 ierr = DMGetDimension(dm,&dim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),143,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
144 if (height > dim || height < 0) SETERRQ2(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE,"Asked for dual space at height %d for dimension %d reference element\n",height,dim)return PetscError(((MPI_Comm)0x44000001),144,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,63,PETSC_ERROR_INITIAL,"Asked for dual space at height %d for dimension %d reference element\n"
,height,dim);
145 ierr = PetscDualSpaceDuplicate(sp,bdsp);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),145,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
146 ierr = PetscDualSpaceCreateReferenceCell(*bdsp, dim-height, lag->simplexCell, &K);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),146,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
147 ierr = PetscDualSpaceSetDM(*bdsp, K);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),147,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
148 ierr = DMDestroy(&K);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),148,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
149 ierr = PetscDualSpaceLagrangeGetTensor(sp,&tensor);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),149,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
150 ierr = PetscDualSpaceLagrangeSetTensor(*bdsp,tensor);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),150,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
151 ierr = PetscDualSpaceSetUp(*bdsp);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),151,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
152 }
153 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);
154}
155
156static PetscErrorCode PetscDualSpaceSetUp_Lagrange(PetscDualSpace sp)
157{
158 PetscDualSpace_Lag *lag = (PetscDualSpace_Lag *) sp->data;
159 DM dm = sp->dm;
160 PetscInt order = sp->order;
161 PetscInt Nc = sp->Nc;
162 MPI_Comm comm;
163 PetscBool continuous;
164 PetscSection csection;
165 Vec coordinates;
166 PetscReal *qpoints, *qweights;
167 PetscInt depth, dim, pdimMax, pStart, pEnd, p, *pStratStart, *pStratEnd, coneSize, d, f = 0, c;
168 PetscBool simplex, tensorSpace;
169 PetscErrorCode ierr;
170
171 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 171; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
172 ierr = PetscObjectGetComm((PetscObject) sp, &comm);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),172,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
173 if (!order) lag->continuous = PETSC_FALSE;
174 continuous = lag->continuous;
175 ierr = DMGetDimension(dm, &dim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),175,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
176 ierr = DMPlexGetDepth(dm, &depth);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),176,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
177 ierr = DMPlexGetChart(dm, &pStart, &pEnd);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),177,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
178 ierr = PetscCalloc1(dim+1, &lag->numDof)PetscMallocA(1,PETSC_TRUE,178,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(dim+1)*sizeof(**(&lag->numDof)),(&lag->
numDof));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),178,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
179 ierr = PetscMalloc2(depth+1,&pStratStart,depth+1,&pStratEnd)PetscMallocA(2,PETSC_FALSE,179,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(depth+1)*sizeof(**(&pStratStart)),(&pStratStart
),(size_t)(depth+1)*sizeof(**(&pStratEnd)),(&pStratEnd
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),179,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
180 for (d = 0; d <= depth; ++d) {ierr = DMPlexGetDepthStratum(dm, d, &pStratStart[d], &pStratEnd[d]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),180,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}
181 ierr = DMPlexGetConeSize(dm, pStratStart[depth], &coneSize);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),181,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
182 ierr = DMGetCoordinateSection(dm, &csection);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),182,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
183 ierr = DMGetCoordinatesLocal(dm, &coordinates);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),183,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
184 if (depth == 1) {
185 if (coneSize == dim+1) simplex = PETSC_TRUE;
186 else if (coneSize == 1 << dim) simplex = PETSC_FALSE;
187 else SETERRQ(comm, PETSC_ERR_SUP, "Only support simplices and tensor product cells")return PetscError(comm,187,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,56,PETSC_ERROR_INITIAL,"Only support simplices and tensor product cells"
);
188 } else if (depth == dim) {
189 if (coneSize == dim+1) simplex = PETSC_TRUE;
190 else if (coneSize == 2 * dim) simplex = PETSC_FALSE;
191 else SETERRQ(comm, PETSC_ERR_SUP, "Only support simplices and tensor product cells")return PetscError(comm,191,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,56,PETSC_ERROR_INITIAL,"Only support simplices and tensor product cells"
);
192 } else SETERRQ(comm, PETSC_ERR_SUP, "Only support cell-vertex meshes or fully interpolated meshes")return PetscError(comm,192,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,56,PETSC_ERROR_INITIAL,"Only support cell-vertex meshes or fully interpolated meshes"
);
193 lag->simplexCell = simplex;
194 if (dim > 1 && continuous && lag->simplexCell == lag->tensorSpace) SETERRQ(PETSC_COMM_SELF,PETSC_ERR_SUP, "Mismatching simplex/tensor cells and spaces only allowed for discontinuous elements")return PetscError(((MPI_Comm)0x44000001),194,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,56,PETSC_ERROR_INITIAL,"Mismatching simplex/tensor cells and spaces only allowed for discontinuous elements"
);
195 tensorSpace = lag->tensorSpace;
196 lag->height = 0;
197 lag->subspaces = NULL((void*)0);
198 if (continuous && order > 0 && dim > 0) {
199 PetscInt i;
200
201 lag->height = dim;
202 ierr = PetscMalloc1(dim,&lag->subspaces)PetscMallocA(1,PETSC_FALSE,202,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(dim)*sizeof(**(&lag->subspaces)),(&lag->
subspaces));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),202,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
203 ierr = PetscDualSpaceCreateHeightSubspace_Lagrange(sp,1,&lag->subspaces[0]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),203,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
204 ierr = PetscDualSpaceSetUp(lag->subspaces[0]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),204,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
205 for (i = 1; i < dim; i++) {
206 ierr = PetscDualSpaceGetHeightSubspace(lag->subspaces[i-1],1,&lag->subspaces[i]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),206,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
207 ierr = PetscObjectReference((PetscObject)(lag->subspaces[i]));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),207,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
208 }
209 }
210 ierr = PetscDualSpaceGetDimension_SingleCell_Lagrange(sp, sp->order, &pdimMax);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),210,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
211 pdimMax *= (pStratEnd[depth] - pStratStart[depth]);
212 ierr = PetscMalloc1(pdimMax, &sp->functional)PetscMallocA(1,PETSC_FALSE,212,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(pdimMax)*sizeof(**(&sp->functional)),(&sp
->functional));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),212,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
213 if (!dim) {
214 for (c = 0; c < Nc; ++c) {
215 ierr = PetscQuadratureCreate(PETSC_COMM_SELF((MPI_Comm)0x44000001), &sp->functional[f]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),215,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
216 ierr = PetscCalloc1(Nc, &qweights)PetscMallocA(1,PETSC_TRUE,216,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(Nc)*sizeof(**(&qweights)),(&qweights));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),216,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
217 ierr = PetscQuadratureSetOrder(sp->functional[f], 0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),217,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
218 ierr = PetscQuadratureSetData(sp->functional[f], 0, Nc, 1, NULL((void*)0), qweights);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),218,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
219 qweights[c] = 1.0;
220 ++f;
221 lag->numDof[0]++;
222 }
223 } else {
224 PetscSection section;
225 PetscReal *v0, *hv0, *J, *invJ, detJ, hdetJ;
226 PetscInt *tup;
227
228 ierr = PetscSectionCreate(PETSC_COMM_SELF((MPI_Comm)0x44000001),&section);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),228,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
229 ierr = PetscSectionSetChart(section,pStart,pEnd);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),229,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
230 ierr = PetscCalloc5(dim+1,&tup,dim,&v0,dim,&hv0,dim*dim,&J,dim*dim,&invJ)PetscMallocA(5,PETSC_TRUE,230,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(dim+1)*sizeof(**(&tup)),(&tup),(size_t)(dim
)*sizeof(**(&v0)),(&v0),(size_t)(dim)*sizeof(**(&
hv0)),(&hv0),(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),230,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
231 for (p = pStart; p < pEnd; p++) {
232 PetscInt pointDim, d, nFunc = 0;
233 PetscDualSpace hsp;
234
235 ierr = DMPlexComputeCellGeometryFEM(dm, p, NULL((void*)0), v0, J, invJ, &detJ);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),235,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
236 for (d = 0; d < depth; d++) {if (p >= pStratStart[d] && p < pStratEnd[d]) break;}
237 pointDim = (depth == 1 && d == 1) ? dim : d;
238 hsp = ((pointDim < dim) && lag->subspaces) ? lag->subspaces[dim - pointDim - 1] : NULL((void*)0);
239 if (hsp) {
240 PetscDualSpace_Lag *hlag = (PetscDualSpace_Lag *) hsp->data;
241 DM hdm;
242
243 ierr = PetscDualSpaceGetDM(hsp,&hdm);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),243,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
244 ierr = DMPlexComputeCellGeometryFEM(hdm, 0, NULL((void*)0), hv0, NULL((void*)0), NULL((void*)0), &hdetJ);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),244,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
245 nFunc = lag->numDof[pointDim] = hlag->numDof[pointDim];
246 }
247 if (pointDim == dim) {
248 /* Cells, create for self */
249 PetscInt orderEff = continuous ? (!tensorSpace ? order-1-dim : order-2) : order;
250 PetscReal denom = continuous ? order : (!tensorSpace ? order+1+dim : order+2);
251 PetscReal numer = (!simplex || !tensorSpace) ? 2. : (2./dim);
252 PetscReal dx = numer/denom;
253 PetscInt cdim, d, d2;
254
255 if (orderEff < 0) continue;
256 ierr = PetscDualSpaceGetDimension_SingleCell_Lagrange(sp, orderEff, &cdim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),256,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
257 ierr = PetscArrayzero(tup,dim+1)PetscMemzero(tup,(dim+1)*sizeof(*(tup)));;CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),257,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
258 if (!tensorSpace) {
259 while (!tup[dim]) {
260 for (c = 0; c < Nc; ++c) {
261 ierr = PetscQuadratureCreate(PETSC_COMM_SELF((MPI_Comm)0x44000001), &sp->functional[f]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),261,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
262 ierr = PetscMalloc1(dim, &qpoints)PetscMallocA(1,PETSC_FALSE,262,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(dim)*sizeof(**(&qpoints)),(&qpoints));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),262,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
263 ierr = PetscCalloc1(Nc, &qweights)PetscMallocA(1,PETSC_TRUE,263,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(Nc)*sizeof(**(&qweights)),(&qweights));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),263,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
264 ierr = PetscQuadratureSetOrder(sp->functional[f], 0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),264,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
265 ierr = PetscQuadratureSetData(sp->functional[f], dim, Nc, 1, qpoints, qweights);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),265,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
266 for (d = 0; d < dim; ++d) {
267 qpoints[d] = v0[d];
268 for (d2 = 0; d2 < dim; ++d2) qpoints[d] += J[d*dim+d2]*((tup[d2]+1)*dx);
269 }
270 qweights[c] = 1.0;
271 ++f;
272 }
273 ierr = PetscDualSpaceLatticePointLexicographic_Internal(dim, orderEff, tup);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),273,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
274 }
275 } else {
276 while (!tup[dim]) {
277 for (c = 0; c < Nc; ++c) {
278 ierr = PetscQuadratureCreate(PETSC_COMM_SELF((MPI_Comm)0x44000001), &sp->functional[f]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),278,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
279 ierr = PetscMalloc1(dim, &qpoints)PetscMallocA(1,PETSC_FALSE,279,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(dim)*sizeof(**(&qpoints)),(&qpoints));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),279,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
280 ierr = PetscCalloc1(Nc, &qweights)PetscMallocA(1,PETSC_TRUE,280,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(Nc)*sizeof(**(&qweights)),(&qweights));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),280,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
281 ierr = PetscQuadratureSetOrder(sp->functional[f], 0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),281,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
282 ierr = PetscQuadratureSetData(sp->functional[f], dim, Nc, 1, qpoints, qweights);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),282,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
283 for (d = 0; d < dim; ++d) {
284 qpoints[d] = v0[d];
285 for (d2 = 0; d2 < dim; ++d2) qpoints[d] += J[d*dim+d2]*((tup[d2]+1)*dx);
286 }
287 qweights[c] = 1.0;
288 ++f;
289 }
290 ierr = PetscDualSpaceTensorPointLexicographic_Internal(dim, orderEff, tup);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),290,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
291 }
292 }
293 lag->numDof[dim] = cdim;
294 } else { /* transform functionals from subspaces */
295 PetscInt q;
296
297 for (q = 0; q < nFunc; q++, f++) {
298 PetscQuadrature fn;
299 PetscInt fdim, Nc, c, nPoints, i;
300 const PetscReal *points;
301 const PetscReal *weights;
302 PetscReal *qpoints;
303 PetscReal *qweights;
304
305 ierr = PetscDualSpaceGetFunctional(hsp, q, &fn);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),305,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
306 ierr = PetscQuadratureGetData(fn,&fdim,&Nc,&nPoints,&points,&weights);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),306,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
307 if (fdim != pointDim) SETERRQ2(PETSC_COMM_SELF,PETSC_ERR_PLIB,"Expected height dual space dim %D, got %D",pointDim,fdim)return PetscError(((MPI_Comm)0x44000001),307,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,77,PETSC_ERROR_INITIAL,"Expected height dual space dim %D, got %D"
,pointDim,fdim);
308 ierr = PetscMalloc1(nPoints * dim, &qpoints)PetscMallocA(1,PETSC_FALSE,308,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(nPoints * dim)*sizeof(**(&qpoints)),(&qpoints
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),308,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
309 ierr = PetscCalloc1(nPoints * Nc, &qweights)PetscMallocA(1,PETSC_TRUE,309,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(nPoints * Nc)*sizeof(**(&qweights)),(&qweights
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),309,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
310 for (i = 0; i < nPoints; i++) {
311 PetscInt j, k;
312 PetscReal *qp = &qpoints[i * dim];
313
314 for (c = 0; c < Nc; ++c) qweights[i*Nc+c] = weights[i*Nc+c];
315 for (j = 0; j < dim; ++j) qp[j] = v0[j];
316 for (j = 0; j < dim; ++j) {
317 for (k = 0; k < pointDim; k++) qp[j] += J[dim * j + k] * (points[pointDim * i + k] - hv0[k]);
318 }
319 }
320 ierr = PetscQuadratureCreate(PETSC_COMM_SELF((MPI_Comm)0x44000001), &sp->functional[f]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),320,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
321 ierr = PetscQuadratureSetOrder(sp->functional[f],0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),321,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
322 ierr = PetscQuadratureSetData(sp->functional[f],dim,Nc,nPoints,qpoints,qweights);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),322,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
323 }
324 }
325 ierr = PetscSectionSetDof(section,p,lag->numDof[pointDim]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),325,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
326 }
327 ierr = PetscFree5(tup,v0,hv0,J,invJ)PetscFreeA(5,327,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,&(tup),&(v0),&(hv0),&(J),&(invJ));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),327,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
328 ierr = PetscSectionSetUp(section);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),328,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
329 { /* reorder to closure order */
330 PetscInt *key, count;
331 PetscQuadrature *reorder = NULL((void*)0);
332
333 ierr = PetscCalloc1(f,&key)PetscMallocA(1,PETSC_TRUE,333,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(f)*sizeof(**(&key)),(&key));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),333,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
334 ierr = PetscMalloc1(f*sp->Nc,&reorder)PetscMallocA(1,PETSC_FALSE,334,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(f*sp->Nc)*sizeof(**(&reorder)),(&reorder
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),334,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
335
336 for (p = pStratStart[depth], count = 0; p < pStratEnd[depth]; p++) {
337 PetscInt *closure = NULL((void*)0), closureSize, c;
338
339 ierr = DMPlexGetTransitiveClosure(dm,p,PETSC_TRUE,&closureSize,&closure);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),339,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
340 for (c = 0; c < closureSize; c++) {
341 PetscInt point = closure[2 * c], dof, off, i;
342
343 ierr = PetscSectionGetDof(section,point,&dof);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),343,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
344 ierr = PetscSectionGetOffset(section,point,&off);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),344,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
345 for (i = 0; i < dof; i++) {
346 PetscInt fi = i + off;
347 if (!key[fi]) {
348 key[fi] = 1;
349 reorder[count++] = sp->functional[fi];
350 }
351 }
352 }
353 ierr = DMPlexRestoreTransitiveClosure(dm,p,PETSC_TRUE,&closureSize,&closure);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),353,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
354 }
355 ierr = PetscFree(sp->functional)((*PetscTrFree)((void*)(sp->functional),355,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
) || ((sp->functional) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),355,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
356 sp->functional = reorder;
357 ierr = PetscFree(key)((*PetscTrFree)((void*)(key),357,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
) || ((key) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),357,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
358 }
359 ierr = PetscSectionDestroy(&section);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),359,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
360 }
361 if (pStratEnd[depth] == 1 && f != pdimMax) SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_PLIB, "Number of dual basis vectors %D not equal to dimension %D", f, pdimMax)return PetscError(((MPI_Comm)0x44000001),361,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,77,PETSC_ERROR_INITIAL,"Number of dual basis vectors %D not equal to dimension %D"
,f,pdimMax);
362 if (f > pdimMax) SETERRQ2(PETSC_COMM_SELF, PETSC_ERR_PLIB, "Number of dual basis vectors %D is greater than max size %D", f, pdimMax)return PetscError(((MPI_Comm)0x44000001),362,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,77,PETSC_ERROR_INITIAL,"Number of dual basis vectors %D is greater than max size %D"
,f,pdimMax);
363 ierr = PetscFree2(pStratStart, pStratEnd)PetscFreeA(2,363,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,&(pStratStart),&(pStratEnd));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),363,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.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
367static PetscErrorCode PetscDualSpaceGetDimension_Lagrange(PetscDualSpace sp, PetscInt *dim)
368{
369 DM K;
370 const PetscInt *numDof;
371 PetscInt spatialDim, Nc, size = 0, d;
372 PetscErrorCode ierr;
373
374 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 374; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
375 ierr = PetscDualSpaceGetDM(sp, &K);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),375,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
376 ierr = PetscDualSpaceGetNumDof(sp, &numDof);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),376,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
377 ierr = DMGetDimension(K, &spatialDim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),377,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
378 ierr = DMPlexGetHeightStratum(K, 0, NULL((void*)0), &Nc);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),378,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
379 if (Nc == 1) {ierr = PetscDualSpaceGetDimension_SingleCell_Lagrange(sp, sp->order, dim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),379,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0); 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);}
380 for (d = 0; d <= spatialDim; ++d) {
381 PetscInt pStart, pEnd;
382
383 ierr = DMPlexGetDepthStratum(K, d, &pStart, &pEnd);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),383,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
384 size += (pEnd-pStart)*numDof[d];
385 }
386 *dim = size;
387 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);
388}
389
390static PetscErrorCode PetscDualSpaceGetNumDof_Lagrange(PetscDualSpace sp, const PetscInt **numDof)
391{
392 PetscDualSpace_Lag *lag = (PetscDualSpace_Lag *) sp->data;
393
394 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 394; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
395 *numDof = lag->numDof;
396 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);
397}
398
399static PetscErrorCode PetscDualSpaceGetHeightSubspace_Lagrange(PetscDualSpace sp, PetscInt height, PetscDualSpace *bdsp)
400{
401 PetscDualSpace_Lag *lag = (PetscDualSpace_Lag *) sp->data;
402 PetscErrorCode ierr;
403
404 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 404; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
405 if (height == 0) {
406 *bdsp = sp;
407 } else {
408 DM dm;
409 PetscInt dim;
410
411 ierr = PetscDualSpaceGetDM(sp,&dm);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),411,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
412 ierr = DMGetDimension(dm,&dim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),412,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
413 if (height > dim || height < 0) SETERRQ2(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE,"Asked for dual space at height %D for dimension %D reference element\n",height,dim)return PetscError(((MPI_Comm)0x44000001),413,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,63,PETSC_ERROR_INITIAL,"Asked for dual space at height %D for dimension %D reference element\n"
,height,dim);
414 if (height <= lag->height) {*bdsp = lag->subspaces[height-1];}
415 else {*bdsp = NULL((void*)0);}
416 }
417 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);
418}
419
420#define BaryIndex(perEdge,a,b,c)(((b)*(2*perEdge+1-(b)))/2)+(c) (((b)*(2*perEdge+1-(b)))/2)+(c)
421
422#define CartIndex(perEdge,a,b)(perEdge*(a)+b) (perEdge*(a)+b)
423
424static PetscErrorCode PetscDualSpaceGetSymmetries_Lagrange(PetscDualSpace sp, const PetscInt ****perms, const PetscScalar ****flips)
425{
426
427 PetscDualSpace_Lag *lag = (PetscDualSpace_Lag *) sp->data;
428 PetscInt dim, order, p, Nc;
429 PetscErrorCode ierr;
430
431 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 431; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
432 ierr = PetscDualSpaceGetOrder(sp,&order);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),432,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
433 ierr = PetscDualSpaceGetNumComponents(sp,&Nc);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),433,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
434 ierr = DMGetDimension(sp->dm,&dim);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),434,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
435 if (!dim || !lag->continuous || order < 3) 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);
436 if (dim > 3) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_SUP,"Lagrange symmetries not implemented for dim = %D > 3",dim)return PetscError(((MPI_Comm)0x44000001),436,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,56,PETSC_ERROR_INITIAL,"Lagrange symmetries not implemented for dim = %D > 3"
,dim);
437 if (!lag->symmetries) { /* store symmetries */
438 PetscDualSpace hsp;
439 DM K;
440 PetscInt numPoints = 1, d;
441 PetscInt numFaces;
442 PetscInt ***symmetries;
443 const PetscInt ***hsymmetries;
444
445 if (lag->simplexCell) {
446 numFaces = 1 + dim;
447 for (d = 0; d < dim; d++) numPoints = numPoints * 2 + 1;
448 } else {
449 numPoints = PetscPowInt(3,dim);
450 numFaces = 2 * dim;
451 }
452 ierr = PetscCalloc1(numPoints,&symmetries)PetscMallocA(1,PETSC_TRUE,452,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(numPoints)*sizeof(**(&symmetries)),(&symmetries
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),452,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
453 if (0 < dim && dim < 3) { /* compute self symmetries */
454 PetscInt **cellSymmetries;
455
456 lag->numSelfSym = 2 * numFaces;
457 lag->selfSymOff = numFaces;
458 ierr = PetscCalloc1(2*numFaces,&cellSymmetries)PetscMallocA(1,PETSC_TRUE,458,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(2*numFaces)*sizeof(**(&cellSymmetries)),(&cellSymmetries
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),458,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
459 /* we want to be able to index symmetries directly with the orientations, which range from [-numFaces,numFaces) */
460 symmetries[0] = &cellSymmetries[numFaces];
461 if (dim == 1) {
462 PetscInt dofPerEdge = order - 1;
463
464 if (dofPerEdge > 1) {
465 PetscInt i, j, *reverse;
466
467 ierr = PetscMalloc1(dofPerEdge*Nc,&reverse)PetscMallocA(1,PETSC_FALSE,467,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(dofPerEdge*Nc)*sizeof(**(&reverse)),(&reverse
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),467,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
468 for (i = 0; i < dofPerEdge; i++) {
469 for (j = 0; j < Nc; j++) {
470 reverse[i*Nc + j] = Nc * (dofPerEdge - 1 - i) + j;
471 }
472 }
473 symmetries[0][-2] = reverse;
474
475 /* yes, this is redundant, but it makes it easier to cleanup if I don't have to worry about what not to free */
476 ierr = PetscMalloc1(dofPerEdge*Nc,&reverse)PetscMallocA(1,PETSC_FALSE,476,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(dofPerEdge*Nc)*sizeof(**(&reverse)),(&reverse
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),476,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
477 for (i = 0; i < dofPerEdge; i++) {
478 for (j = 0; j < Nc; j++) {
479 reverse[i*Nc + j] = Nc * (dofPerEdge - 1 - i) + j;
480 }
481 }
482 symmetries[0][1] = reverse;
483 }
484 } else {
485 PetscInt dofPerEdge = lag->simplexCell ? (order - 2) : (order - 1), s;
486 PetscInt dofPerFace;
487
488 if (dofPerEdge > 1) {
489 for (s = -numFaces; s < numFaces; s++) {
490 PetscInt *sym, i, j, k, l;
491
492 if (!s) continue;
493 if (lag->simplexCell) {
494 dofPerFace = (dofPerEdge * (dofPerEdge + 1))/2;
495 ierr = PetscMalloc1(Nc*dofPerFace,&sym)PetscMallocA(1,PETSC_FALSE,495,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(Nc*dofPerFace)*sizeof(**(&sym)),(&sym));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),495,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
496 for (j = 0, l = 0; j < dofPerEdge; j++) {
497 for (k = 0; k < dofPerEdge - j; k++, l++) {
498 i = dofPerEdge - 1 - j - k;
499 switch (s) {
500 case -3:
501 sym[Nc*l] = BaryIndex(dofPerEdge,i,k,j)(((k)*(2*dofPerEdge+1-(k)))/2)+(j);
502 break;
503 case -2:
504 sym[Nc*l] = BaryIndex(dofPerEdge,j,i,k)(((i)*(2*dofPerEdge+1-(i)))/2)+(k);
505 break;
506 case -1:
507 sym[Nc*l] = BaryIndex(dofPerEdge,k,j,i)(((j)*(2*dofPerEdge+1-(j)))/2)+(i);
508 break;
509 case 1:
510 sym[Nc*l] = BaryIndex(dofPerEdge,k,i,j)(((i)*(2*dofPerEdge+1-(i)))/2)+(j);
511 break;
512 case 2:
513 sym[Nc*l] = BaryIndex(dofPerEdge,j,k,i)(((k)*(2*dofPerEdge+1-(k)))/2)+(i);
514 break;
515 }
516 }
517 }
518 } else {
519 dofPerFace = dofPerEdge * dofPerEdge;
520 ierr = PetscMalloc1(Nc*dofPerFace,&sym)PetscMallocA(1,PETSC_FALSE,520,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(Nc*dofPerFace)*sizeof(**(&sym)),(&sym));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),520,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
521 for (j = 0, l = 0; j < dofPerEdge; j++) {
522 for (k = 0; k < dofPerEdge; k++, l++) {
523 switch (s) {
524 case -4:
525 sym[Nc*l] = CartIndex(dofPerEdge,k,j)(dofPerEdge*(k)+j);
526 break;
527 case -3:
528 sym[Nc*l] = CartIndex(dofPerEdge,(dofPerEdge - 1 - j),k)(dofPerEdge*((dofPerEdge - 1 - j))+k);
529 break;
530 case -2:
531 sym[Nc*l] = CartIndex(dofPerEdge,(dofPerEdge - 1 - k),(dofPerEdge - 1 - j))(dofPerEdge*((dofPerEdge - 1 - k))+(dofPerEdge - 1 - j));
532 break;
533 case -1:
534 sym[Nc*l] = CartIndex(dofPerEdge,j,(dofPerEdge - 1 - k))(dofPerEdge*(j)+(dofPerEdge - 1 - k));
535 break;
536 case 1:
537 sym[Nc*l] = CartIndex(dofPerEdge,(dofPerEdge - 1 - k),j)(dofPerEdge*((dofPerEdge - 1 - k))+j);
538 break;
539 case 2:
540 sym[Nc*l] = CartIndex(dofPerEdge,(dofPerEdge - 1 - j),(dofPerEdge - 1 - k))(dofPerEdge*((dofPerEdge - 1 - j))+(dofPerEdge - 1 - k));
541 break;
542 case 3:
543 sym[Nc*l] = CartIndex(dofPerEdge,k,(dofPerEdge - 1 - j))(dofPerEdge*(k)+(dofPerEdge - 1 - j));
544 break;
545 }
546 }
547 }
548 }
549 for (i = 0; i < dofPerFace; i++) {
550 sym[Nc*i] *= Nc;
551 for (j = 1; j < Nc; j++) {
552 sym[Nc*i+j] = sym[Nc*i] + j;
553 }
554 }
555 symmetries[0][s] = sym;
556 }
557 }
558 }
559 }
560 ierr = PetscDualSpaceGetHeightSubspace(sp,1,&hsp);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),560,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
561 ierr = PetscDualSpaceGetSymmetries(hsp,&hsymmetries,NULL((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),561,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
562 if (hsymmetries) {
563 PetscBool *seen;
564 const PetscInt *cone;
565 PetscInt KclosureSize, *Kclosure = NULL((void*)0);
566
567 ierr = PetscDualSpaceGetDM(sp,&K);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),567,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
568 ierr = PetscCalloc1(numPoints,&seen)PetscMallocA(1,PETSC_TRUE,568,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(numPoints)*sizeof(**(&seen)),(&seen));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),568,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
569 ierr = DMPlexGetCone(K,0,&cone);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),569,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
570 ierr = DMPlexGetTransitiveClosure(K,0,PETSC_TRUE,&KclosureSize,&Kclosure);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),570,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
571 for (p = 0; p < numFaces; p++) {
572 PetscInt closureSize, *closure = NULL((void*)0), q;
573
574 ierr = DMPlexGetTransitiveClosure(K,cone[p],PETSC_TRUE,&closureSize,&closure);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),574,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
575 for (q = 0; q < closureSize; q++) {
576 PetscInt point = closure[2*q], r;
577
578 if(!seen[point]) {
579 for (r = 0; r < KclosureSize; r++) {
580 if (Kclosure[2 * r] == point) break;
581 }
582 seen[point] = PETSC_TRUE;
583 symmetries[r] = (PetscInt **) hsymmetries[q];
584 }
585 }
586 ierr = DMPlexRestoreTransitiveClosure(K,cone[p],PETSC_TRUE,&closureSize,&closure);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),586,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
587 }
588 ierr = DMPlexRestoreTransitiveClosure(K,0,PETSC_TRUE,&KclosureSize,&Kclosure);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),588,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
589 ierr = PetscFree(seen)((*PetscTrFree)((void*)(seen),589,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
) || ((seen) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),589,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
590 }
591 lag->symmetries = symmetries;
592 }
593 if (perms) *perms = (const PetscInt ***) lag->symmetries;
594 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);
595}
596
597static PetscErrorCode PetscDualSpaceLagrangeGetContinuity_Lagrange(PetscDualSpace sp, PetscBool *continuous)
598{
599 PetscDualSpace_Lag *lag = (PetscDualSpace_Lag *) sp->data;
600
601 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 601; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
602 PetscValidHeaderSpecific(sp, PETSCDUALSPACE_CLASSID, 1)do { if (!sp) return PetscError(((MPI_Comm)0x44000001),602,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,85,PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",1); if (
!PetscCheckPointer(sp,PETSC_OBJECT)) return PetscError(((MPI_Comm
)0x44000001),602,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,1); if (((PetscObject)(sp))->classid != PETSCDUALSPACE_CLASSID
) { if (((PetscObject)(sp))->classid == -1) return PetscError
(((MPI_Comm)0x44000001),602,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,1); else return PetscError(((MPI_Comm)0x44000001),602,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,62,PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d"
,1); } } while (0);
603 PetscValidPointer(continuous, 2)do { if (!continuous) return PetscError(((MPI_Comm)0x44000001
),603,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,85,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",2); if
(!PetscCheckPointer(continuous,PETSC_CHAR)) return PetscError
(((MPI_Comm)0x44000001),603,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer: Parameter # %d",2);
} while (0);
604 *continuous = lag->continuous;
605 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);
606}
607
608static PetscErrorCode PetscDualSpaceLagrangeSetContinuity_Lagrange(PetscDualSpace sp, PetscBool continuous)
609{
610 PetscDualSpace_Lag *lag = (PetscDualSpace_Lag *) sp->data;
611
612 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 612; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
613 PetscValidHeaderSpecific(sp, PETSCDUALSPACE_CLASSID, 1)do { if (!sp) return PetscError(((MPI_Comm)0x44000001),613,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,85,PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",1); if (
!PetscCheckPointer(sp,PETSC_OBJECT)) return PetscError(((MPI_Comm
)0x44000001),613,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,1); if (((PetscObject)(sp))->classid != PETSCDUALSPACE_CLASSID
) { if (((PetscObject)(sp))->classid == -1) return PetscError
(((MPI_Comm)0x44000001),613,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,1); else return PetscError(((MPI_Comm)0x44000001),613,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,62,PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d"
,1); } } while (0);
614 lag->continuous = continuous;
615 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);
616}
617
618/*@
619 PetscDualSpaceLagrangeGetContinuity - Retrieves the flag for element continuity
620
621 Not Collective
622
623 Input Parameter:
624. sp - the PetscDualSpace
625
626 Output Parameter:
627. continuous - flag for element continuity
628
629 Level: intermediate
630
631.seealso: PetscDualSpaceLagrangeSetContinuity()
632@*/
633PetscErrorCode PetscDualSpaceLagrangeGetContinuity(PetscDualSpace sp, PetscBool *continuous)
634{
635 PetscErrorCode ierr;
636
637 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 637; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
638 PetscValidHeaderSpecific(sp, PETSCDUALSPACE_CLASSID, 1)do { if (!sp) return PetscError(((MPI_Comm)0x44000001),638,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,85,PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",1); if (
!PetscCheckPointer(sp,PETSC_OBJECT)) return PetscError(((MPI_Comm
)0x44000001),638,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,1); if (((PetscObject)(sp))->classid != PETSCDUALSPACE_CLASSID
) { if (((PetscObject)(sp))->classid == -1) return PetscError
(((MPI_Comm)0x44000001),638,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,1); else return PetscError(((MPI_Comm)0x44000001),638,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,62,PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d"
,1); } } while (0);
639 PetscValidPointer(continuous, 2)do { if (!continuous) return PetscError(((MPI_Comm)0x44000001
),639,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,85,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",2); if
(!PetscCheckPointer(continuous,PETSC_CHAR)) return PetscError
(((MPI_Comm)0x44000001),639,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer: Parameter # %d",2);
} while (0);
640 ierr = PetscTryMethod(sp, "PetscDualSpaceLagrangeGetContinuity_C", (PetscDualSpace,PetscBool*),(sp,continuous))0; do { PetscErrorCode (*f)(PetscDualSpace,PetscBool*), __ierr
; __ierr = PetscObjectQueryFunction_Private(((PetscObject)sp)
,("PetscDualSpaceLagrangeGetContinuity_C"),(PetscVoidFunction
*)(&f));do {if (__builtin_expect(!!(__ierr),0)) return PetscError
(((MPI_Comm)0x44000001),640,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0); if (f) {__ierr =
(*f)(sp,continuous);do {if (__builtin_expect(!!(__ierr),0)) return
PetscError(((MPI_Comm)0x44000001),640,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0);} } while(0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),640,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
641 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);
642}
643
644/*@
645 PetscDualSpaceLagrangeSetContinuity - Indicate whether the element is continuous
646
647 Logically Collective on sp
648
649 Input Parameters:
650+ sp - the PetscDualSpace
651- continuous - flag for element continuity
652
653 Options Database:
654. -petscdualspace_lagrange_continuity <bool>
655
656 Level: intermediate
657
658.seealso: PetscDualSpaceLagrangeGetContinuity()
659@*/
660PetscErrorCode PetscDualSpaceLagrangeSetContinuity(PetscDualSpace sp, PetscBool continuous)
661{
662 PetscErrorCode ierr;
663
664 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 664; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
665 PetscValidHeaderSpecific(sp, PETSCDUALSPACE_CLASSID, 1)do { if (!sp) return PetscError(((MPI_Comm)0x44000001),665,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,85,PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",1); if (
!PetscCheckPointer(sp,PETSC_OBJECT)) return PetscError(((MPI_Comm
)0x44000001),665,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,1); if (((PetscObject)(sp))->classid != PETSCDUALSPACE_CLASSID
) { if (((PetscObject)(sp))->classid == -1) return PetscError
(((MPI_Comm)0x44000001),665,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,1); else return PetscError(((MPI_Comm)0x44000001),665,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,62,PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d"
,1); } } while (0);
666 PetscValidLogicalCollectiveBool(sp, continuous, 2)do { PetscErrorCode _7_ierr; PetscMPIInt b1[2],b2[2]; b1[0] =
-(PetscMPIInt)continuous; b1[1] = (PetscMPIInt)continuous; _7_ierr
= (PetscAllreduceBarrierCheck(PetscObjectComm((PetscObject)sp
),2,666,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
) || ((petsc_allreduce_ct += PetscMPIParallelComm((PetscObjectComm
((PetscObject)sp))),0) || MPI_Allreduce((b1),(b2),(2),(((MPI_Datatype
)0x4c000405)),((MPI_Op)(0x58000001)),(PetscObjectComm((PetscObject
)sp)))));do {if (__builtin_expect(!!(_7_ierr),0)) return PetscError
(((MPI_Comm)0x44000001),666,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,_7_ierr,PETSC_ERROR_REPEAT," ");} while (0); if (-b2[0] != b2
[1]) return PetscError(PetscObjectComm((PetscObject)sp),666,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,62,PETSC_ERROR_INITIAL,"Bool value must be same on all processes, argument # %d"
,2); } while (0);
667 ierr = PetscTryMethod(sp, "PetscDualSpaceLagrangeSetContinuity_C", (PetscDualSpace,PetscBool),(sp,continuous))0; do { PetscErrorCode (*f)(PetscDualSpace,PetscBool), __ierr
; __ierr = PetscObjectQueryFunction_Private(((PetscObject)sp)
,("PetscDualSpaceLagrangeSetContinuity_C"),(PetscVoidFunction
*)(&f));do {if (__builtin_expect(!!(__ierr),0)) return PetscError
(((MPI_Comm)0x44000001),667,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0); if (f) {__ierr =
(*f)(sp,continuous);do {if (__builtin_expect(!!(__ierr),0)) return
PetscError(((MPI_Comm)0x44000001),667,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0);} } while(0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),667,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
668 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);
669}
670
671static PetscErrorCode PetscDualSpaceLagrangeGetTensor_Lagrange(PetscDualSpace sp, PetscBool *tensor)
672{
673 PetscDualSpace_Lag *lag = (PetscDualSpace_Lag *)sp->data;
674
675 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 675; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
676 *tensor = lag->tensorSpace;
677 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);
678}
679
680static PetscErrorCode PetscDualSpaceLagrangeSetTensor_Lagrange(PetscDualSpace sp, PetscBool tensor)
681{
682 PetscDualSpace_Lag *lag = (PetscDualSpace_Lag *)sp->data;
683
684 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 684; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
685 lag->tensorSpace = tensor;
686 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);
687}
688
689/*@
690 PetscDualSpaceLagrangeGetTensor - Get the tensor nature of the dual space
691
692 Not collective
693
694 Input Parameter:
695. sp - The PetscDualSpace
696
697 Output Parameter:
698. tensor - Whether the dual space has tensor layout (vs. simplicial)
699
700 Level: intermediate
701
702.seealso: PetscDualSpaceLagrangeSetTensor(), PetscDualSpaceCreate()
703@*/
704PetscErrorCode PetscDualSpaceLagrangeGetTensor(PetscDualSpace sp, PetscBool *tensor)
705{
706 PetscErrorCode ierr;
707
708 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 708; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
3
←
Within the expansion of the macro 'PetscFunctionBegin':
→
709 PetscValidHeaderSpecific(sp, PETSCDUALSPACE_CLASSID, 1)do { if (!sp) return PetscError(((MPI_Comm)0x44000001),709,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,85,PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",1); if (
!PetscCheckPointer(sp,PETSC_OBJECT)) return PetscError(((MPI_Comm
)0x44000001),709,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,1); if (((PetscObject)(sp))->classid != PETSCDUALSPACE_CLASSID
) { if (((PetscObject)(sp))->classid == -1) return PetscError
(((MPI_Comm)0x44000001),709,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,1); else return PetscError(((MPI_Comm)0x44000001),709,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,62,PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d"
,1); } } while (0);
4
←
Within the expansion of the macro 'PetscValidHeaderSpecific':
→
a
Assuming the condition is false
710 PetscValidPointer(tensor, 2)do { if (!tensor) return PetscError(((MPI_Comm)0x44000001),710
,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,85,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",2); if
(!PetscCheckPointer(tensor,PETSC_CHAR)) return PetscError(((
MPI_Comm)0x44000001),710,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer: Parameter # %d",2);
} while (0);
5
←
Within the expansion of the macro 'PetscValidPointer':
→
a
Assuming the condition is false
711 ierr = PetscTryMethod(sp,"PetscDualSpaceLagrangeGetTensor_C",(PetscDualSpace,PetscBool *),(sp,tensor))0; do { PetscErrorCode (*f)(PetscDualSpace,PetscBool *), __ierr
; __ierr = PetscObjectQueryFunction_Private(((PetscObject)sp)
,("PetscDualSpaceLagrangeGetTensor_C"),(PetscVoidFunction*)(&
f));do {if (__builtin_expect(!!(__ierr),0)) return PetscError
(((MPI_Comm)0x44000001),711,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0); if (f) {__ierr =
(*f)(sp,tensor);do {if (__builtin_expect(!!(__ierr),0)) return
PetscError(((MPI_Comm)0x44000001),711,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0);} } while(0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),711,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
6
←
Within the expansion of the macro 'PetscTryMethod':
→
a
Assuming '__ierr' is 0
b
Assuming 'f' is null
7
←
Within the expansion of the macro 'CHKERRQ':
→
712 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);
8
←
Within the expansion of the macro 'PetscFunctionReturn':
→
a
Assuming 'petscstack' is null
713}
714
715/*@
716 PetscDualSpaceLagrangeSetTensor - Set the tensor nature of the dual space
717
718 Not collective
719
720 Input Parameters:
721+ sp - The PetscDualSpace
722- tensor - Whether the dual space has tensor layout (vs. simplicial)
723
724 Level: intermediate
725
726.seealso: PetscDualSpaceLagrangeGetTensor(), PetscDualSpaceCreate()
727@*/
728PetscErrorCode PetscDualSpaceLagrangeSetTensor(PetscDualSpace sp, PetscBool tensor)
729{
730 PetscErrorCode ierr;
731
732 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 732; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
733 PetscValidHeaderSpecific(sp, PETSCDUALSPACE_CLASSID, 1)do { if (!sp) return PetscError(((MPI_Comm)0x44000001),733,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,85,PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",1); if (
!PetscCheckPointer(sp,PETSC_OBJECT)) return PetscError(((MPI_Comm
)0x44000001),733,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,1); if (((PetscObject)(sp))->classid != PETSCDUALSPACE_CLASSID
) { if (((PetscObject)(sp))->classid == -1) return PetscError
(((MPI_Comm)0x44000001),733,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,1); else return PetscError(((MPI_Comm)0x44000001),733,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,62,PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d"
,1); } } while (0);
734 ierr = PetscTryMethod(sp,"PetscDualSpaceLagrangeSetTensor_C",(PetscDualSpace,PetscBool),(sp,tensor))0; do { PetscErrorCode (*f)(PetscDualSpace,PetscBool), __ierr
; __ierr = PetscObjectQueryFunction_Private(((PetscObject)sp)
,("PetscDualSpaceLagrangeSetTensor_C"),(PetscVoidFunction*)(&
f));do {if (__builtin_expect(!!(__ierr),0)) return PetscError
(((MPI_Comm)0x44000001),734,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0); if (f) {__ierr =
(*f)(sp,tensor);do {if (__builtin_expect(!!(__ierr),0)) return
PetscError(((MPI_Comm)0x44000001),734,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0);} } while(0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),734,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
735 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);
736}
737
738static PetscErrorCode PetscDualSpaceInitialize_Lagrange(PetscDualSpace sp)
739{
740 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 740; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
741 sp->ops->destroy = PetscDualSpaceDestroy_Lagrange;
742 sp->ops->view = PetscDualSpaceView_Lagrange;
743 sp->ops->setfromoptions = PetscDualSpaceSetFromOptions_Lagrange;
744 sp->ops->duplicate = PetscDualSpaceDuplicate_Lagrange;
745 sp->ops->setup = PetscDualSpaceSetUp_Lagrange;
746 sp->ops->getdimension = PetscDualSpaceGetDimension_Lagrange;
747 sp->ops->getnumdof = PetscDualSpaceGetNumDof_Lagrange;
748 sp->ops->getheightsubspace = PetscDualSpaceGetHeightSubspace_Lagrange;
749 sp->ops->getsymmetries = PetscDualSpaceGetSymmetries_Lagrange;
750 sp->ops->apply = PetscDualSpaceApplyDefault;
751 sp->ops->applyall = PetscDualSpaceApplyAllDefault;
752 sp->ops->createallpoints = PetscDualSpaceCreateAllPointsDefault;
753 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);
754}
755
756/*MC
757 PETSCDUALSPACELAGRANGE = "lagrange" - A PetscDualSpace object that encapsulates a dual space of pointwise evaluation functionals
758
759 Level: intermediate
760
761.seealso: PetscDualSpaceType, PetscDualSpaceCreate(), PetscDualSpaceSetType()
762M*/
763
764PETSC_EXTERNextern __attribute__((visibility ("default"))) PetscErrorCode PetscDualSpaceCreate_Lagrange(PetscDualSpace sp)
765{
766 PetscDualSpace_Lag *lag;
767 PetscErrorCode ierr;
768
769 PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
; petscstack->line[petscstack->currentsize] = 769; petscstack
->petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);
770 PetscValidHeaderSpecific(sp, PETSCDUALSPACE_CLASSID, 1)do { if (!sp) return PetscError(((MPI_Comm)0x44000001),770,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,85,PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",1); if (
!PetscCheckPointer(sp,PETSC_OBJECT)) return PetscError(((MPI_Comm
)0x44000001),770,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,1); if (((PetscObject)(sp))->classid != PETSCDUALSPACE_CLASSID
) { if (((PetscObject)(sp))->classid == -1) return PetscError
(((MPI_Comm)0x44000001),770,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,1); else return PetscError(((MPI_Comm)0x44000001),770,__func__
,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,62,PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d"
,1); } } while (0);
771 ierr = PetscNewLog(sp,&lag)(PetscMallocA(1,PETSC_TRUE,771,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,(size_t)(1)*sizeof(**(((&lag)))),(((&lag)))) || PetscLogObjectMemory
((PetscObject)sp,sizeof(**(&lag))));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),771,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
772 sp->data = lag;
773
774 lag->numDof = NULL((void*)0);
775 lag->simplexCell = PETSC_TRUE;
776 lag->tensorSpace = PETSC_FALSE;
777 lag->continuous = PETSC_TRUE;
778
779 ierr = PetscDualSpaceInitialize_Lagrange(sp);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),779,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
780 ierr = PetscObjectComposeFunction((PetscObject) sp, "PetscDualSpaceLagrangeGetContinuity_C", PetscDualSpaceLagrangeGetContinuity_Lagrange)PetscObjectComposeFunction_Private((PetscObject) sp,"PetscDualSpaceLagrangeGetContinuity_C"
,(PetscVoidFunction)(PetscDualSpaceLagrangeGetContinuity_Lagrange
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),780,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
781 ierr = PetscObjectComposeFunction((PetscObject) sp, "PetscDualSpaceLagrangeSetContinuity_C", PetscDualSpaceLagrangeSetContinuity_Lagrange)PetscObjectComposeFunction_Private((PetscObject) sp,"PetscDualSpaceLagrangeSetContinuity_C"
,(PetscVoidFunction)(PetscDualSpaceLagrangeSetContinuity_Lagrange
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),781,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
782 ierr = PetscObjectComposeFunction((PetscObject) sp, "PetscDualSpaceLagrangeGetTensor_C", PetscDualSpaceLagrangeGetTensor_Lagrange)PetscObjectComposeFunction_Private((PetscObject) sp,"PetscDualSpaceLagrangeGetTensor_C"
,(PetscVoidFunction)(PetscDualSpaceLagrangeGetTensor_Lagrange
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),782,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
783 ierr = PetscObjectComposeFunction((PetscObject) sp, "PetscDualSpaceLagrangeSetTensor_C", PetscDualSpaceLagrangeSetTensor_Lagrange)PetscObjectComposeFunction_Private((PetscObject) sp,"PetscDualSpaceLagrangeSetTensor_C"
,(PetscVoidFunction)(PetscDualSpaceLagrangeSetTensor_Lagrange
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),783,__func__,"/sandbox/petsc/petsc.master/src/dm/dt/dualspace/impls/lagrange/dspacelagrange.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);
784 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);
785}
786