Home
last modified time | relevance | path

Searched refs:dd (Results 1 – 25 of 45) sorted by relevance

12

/petsc/src/ksp/ksp/tutorials/
H A Dex59.c66 static PetscErrorCode BuildCSRGraph(DomainData dd, PetscInt **xadj, PetscInt **adjncy) in BuildCSRGraph() argument
76 for (k = 0; k < dd.zm_l; k++) { in BuildCSRGraph()
79 kend_csr = k < dd.zm_l - 1 ? k + 2 : k + 1; in BuildCSRGraph()
80 if (k == 0 || k == dd.zm_l - 1) { in BuildCSRGraph()
85 for (j = 0; j < dd.ym_l; j++) { in BuildCSRGraph()
87 jend_csr = j < dd.ym_l - 1 ? j + 2 : j + 1; in BuildCSRGraph()
88 if (j == 0 || j == dd.ym_l - 1) { in BuildCSRGraph()
93 for (i = 0; i < dd.xm_l; i++) { in BuildCSRGraph()
95 iend_csr = i < dd.xm_l - 1 ? i + 2 : i + 1; in BuildCSRGraph()
96 if (i == 0 || i == dd.xm_l - 1) { in BuildCSRGraph()
[all …]
/petsc/src/dm/impls/da/
H A Ddaltol.c19 DM_DA *dd = (DM_DA *)da->data; in DMLocalToLocalCreate_DA() local
24 if (dd->ltol) PetscFunctionReturn(PETSC_SUCCESS); in DMLocalToLocalCreate_DA()
30 PetscCall(VecScatterCopy(dd->gtol, &dd->ltol)); in DMLocalToLocalCreate_DA()
32 left = dd->xs - dd->Xs; in DMLocalToLocalCreate_DA()
33 PetscCall(PetscMalloc1(dd->xe - dd->xs, &idx)); in DMLocalToLocalCreate_DA()
34 for (j = 0; j < dd->xe - dd->xs; j++) idx[j] = left + j; in DMLocalToLocalCreate_DA()
36 left = dd->xs - dd->Xs; in DMLocalToLocalCreate_DA()
37 down = dd->ys - dd->Ys; in DMLocalToLocalCreate_DA()
38 up = down + dd->ye - dd->ys; in DMLocalToLocalCreate_DA()
39 PetscCall(PetscMalloc1((dd->xe - dd->xs) * (up - down), &idx)); in DMLocalToLocalCreate_DA()
[all …]
H A Dda.c23 DM_DA *dd = (DM_DA *)da->data; in DMDASetSizes() local
35 dd->M = M; in DMDASetSizes()
36 dd->N = N; in DMDASetSizes()
37 dd->P = P; in DMDASetSizes()
58 DM_DA *dd = (DM_DA *)da->data; in DMDASetNumProcs() local
66 dd->m = m; in DMDASetNumProcs()
67 dd->n = n; in DMDASetNumProcs()
68 dd->p = p; in DMDASetNumProcs()
72 if ((dd->m > 0) && (dd->n < 0)) { in DMDASetNumProcs()
73 dd->n = size / dd->m; in DMDASetNumProcs()
[all …]
H A Ddacreate.c5 DM_DA *dd = (DM_DA *)da->data; in DMSetFromOptions_DA() local
12 …PetscCheck(dd->M >= 0, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "Dimension must… in DMSetFromOptions_DA()
13 …PetscCheck(dd->N >= 0, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "Dimension must… in DMSetFromOptions_DA()
14 …PetscCheck(dd->P >= 0, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "Dimension must… in DMSetFromOptions_DA()
17 …Int("-da_grid_x", "Number of grid points in x direction", "DMDASetSizes", dd->M, &dd->M, NULL, 1)); in DMSetFromOptions_DA()
18 …Int("-da_grid_y", "Number of grid points in y direction", "DMDASetSizes", dd->N, &dd->N, NULL, 1)); in DMSetFromOptions_DA()
19 …Int("-da_grid_z", "Number of grid points in z direction", "DMDASetSizes", dd->P, &dd->P, NULL, 1)); in DMSetFromOptions_DA()
21 …verlap", "Decomposition overlap in all directions", "DMDASetOverlap", dd->xol, &dd->xol, &flg, 0)); in DMSetFromOptions_DA()
22 if (flg) PetscCall(DMDASetOverlap(da, dd->xol, dd->xol, dd->xol)); in DMSetFromOptions_DA()
23 …overlap_x", "Decomposition overlap in x direction", "DMDASetOverlap", dd->xol, &dd->xol, NULL, 0)); in DMSetFromOptions_DA()
[all …]
H A Ddadestroy.c9 DM_DA *dd = (DM_DA *)da->data; in DMDestroy_DA() local
14 PetscCall(PetscFree(dd->startghostedout[i])); in DMDestroy_DA()
15 PetscCall(PetscFree(dd->startghostedin[i])); in DMDestroy_DA()
16 PetscCall(PetscFree(dd->startout[i])); in DMDestroy_DA()
17 PetscCall(PetscFree(dd->startin[i])); in DMDestroy_DA()
20 PetscCall(VecScatterDestroy(&dd->gtol)); in DMDestroy_DA()
21 PetscCall(VecScatterDestroy(&dd->ltol)); in DMDestroy_DA()
22 PetscCall(VecDestroy(&dd->natural)); in DMDestroy_DA()
23 PetscCall(VecScatterDestroy(&dd->gton)); in DMDestroy_DA()
24 PetscCall(AODestroy(&dd->ao)); in DMDestroy_DA()
[all …]
H A Ddaindex.c15 DM_DA *dd = (DM_DA *)da->data; in DMDAGetNatural_Private() local
18 Nlocal = (dd->xe - dd->xs); in DMDAGetNatural_Private()
19 if (dim > 1) Nlocal *= (dd->ye - dd->ys); in DMDAGetNatural_Private()
20 if (dim > 2) Nlocal *= (dd->ze - dd->zs); in DMDAGetNatural_Private()
25 for (i = dd->xs; i < dd->xe; i++) { in DMDAGetNatural_Private()
30 for (j = dd->ys; j < dd->ye; j++) { in DMDAGetNatural_Private()
31 for (i = dd->xs; i < dd->xe; i++) { in DMDAGetNatural_Private()
33 lidx[lict++] = i + j * dd->M * dd->w; in DMDAGetNatural_Private()
37 for (k = dd->zs; k < dd->ze; k++) { in DMDAGetNatural_Private()
38 for (j = dd->ys; j < dd->ye; j++) { in DMDAGetNatural_Private()
[all …]
H A Ddaview.c110 DM_DA *dd = (DM_DA *)da->data; in DMDAGetInfo() local
116 if (dd->Mo < 0) *M = dd->M; in DMDAGetInfo()
117 else *M = dd->Mo; in DMDAGetInfo()
120 if (dd->No < 0) *N = dd->N; in DMDAGetInfo()
121 else *N = dd->No; in DMDAGetInfo()
124 if (dd->Po < 0) *P = dd->P; in DMDAGetInfo()
125 else *P = dd->Po; in DMDAGetInfo()
127 if (m) *m = dd->m; in DMDAGetInfo()
128 if (n) *n = dd->n; in DMDAGetInfo()
129 if (p) *p = dd->p; in DMDAGetInfo()
[all …]
H A Ddagtol.c9 DM_DA *dd = (DM_DA *)da->data; in DMGlobalToLocalBegin_DA() local
15 PetscCall(VecScatterBegin(dd->gtol, g, l, mode, SCATTER_FORWARD)); in DMGlobalToLocalBegin_DA()
21 DM_DA *dd = (DM_DA *)da->data; in DMGlobalToLocalEnd_DA() local
27 PetscCall(VecScatterEnd(dd->gtol, g, l, mode, SCATTER_FORWARD)); in DMGlobalToLocalEnd_DA()
33 DM_DA *dd = (DM_DA *)da->data; in DMLocalToGlobalBegin_DA() local
40 PetscCall(VecScatterBegin(dd->gtol, l, g, ADD_VALUES, SCATTER_REVERSE)); in DMLocalToGlobalBegin_DA()
42 …PetscCheck(dd->bx == DM_BOUNDARY_GHOSTED || dd->bx == DM_BOUNDARY_NONE || dd->s <= 0 || dd->m != 1… in DMLocalToGlobalBegin_DA()
43 …PetscCheck(dd->bx == DM_BOUNDARY_GHOSTED || dd->by == DM_BOUNDARY_NONE || dd->s <= 0 || dd->n != 1… in DMLocalToGlobalBegin_DA()
44 …PetscCheck(dd->bx == DM_BOUNDARY_GHOSTED || dd->bz == DM_BOUNDARY_NONE || dd->s <= 0 || dd->p != 1… in DMLocalToGlobalBegin_DA()
45 PetscCall(VecScatterBegin(dd->gtol, l, g, INSERT_VALUES, SCATTER_REVERSE_LOCAL)); in DMLocalToGlobalBegin_DA()
[all …]
H A Ddacorn.c41 DM_DA *dd = (DM_DA *)da->data; in DMDASetFieldName() local
45 …PetscCheck(nf >= 0 && nf < dd->w, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid field number… in DMDASetFieldName()
46 …PetscCheck(dd->fieldname, PetscObjectComm((PetscObject)da), PETSC_ERR_ORDER, "You should call DMSe… in DMDASetFieldName()
47 PetscCall(PetscFree(dd->fieldname[nf])); in DMDASetFieldName()
48 PetscCall(PetscStrallocpy(name, &dd->fieldname[nf])); in DMDASetFieldName()
72 DM_DA *dd = (DM_DA *)da->data; in DMDAGetFieldNames() local
75 *names = (const char *const *)dd->fieldname; in DMDAGetFieldNames()
100 DM_DA *dd = (DM_DA *)da->data; in DMDASetFieldNames() local
105 …PetscCheck(dd->fieldname, PetscObjectComm((PetscObject)da), PETSC_ERR_ORDER, "You should call DMSe… in DMDASetFieldNames()
107 …PetscCheck(nf == dd->w + 1, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid number of fields %… in DMDASetFieldNames()
[all …]
H A Ddaghost.c40 DM_DA *dd = (DM_DA *)da->data; in DMDAGetGhostCorners() local
46 w = dd->w; in DMDAGetGhostCorners()
47 if (x) *x = dd->Xs / w + dd->xo; in DMDAGetGhostCorners()
48 if (y) *y = dd->Ys + dd->yo; in DMDAGetGhostCorners()
49 if (z) *z = dd->Zs + dd->zo; in DMDAGetGhostCorners()
50 if (m) *m = (dd->Xe - dd->Xs) / w; in DMDAGetGhostCorners()
51 if (n) *n = (dd->Ye - dd->Ys); in DMDAGetGhostCorners()
52 if (p) *p = (dd->Ze - dd->Zs); in DMDAGetGhostCorners()
H A Ddasub.c106 DM_DA *dd = (DM_DA *)da->data; in DMDAGetRay() local
119 PetscCall(PetscMalloc1(dd->w, &indices)); in DMDAGetRay()
120 indices[0] = dd->w * gp; in DMDAGetRay()
121 for (i = 1; i < dd->w; ++i) indices[i] = indices[i - 1] + 1; in DMDAGetRay()
122 PetscCall(AOApplicationToPetsc(ao, dd->w, indices)); in DMDAGetRay()
124 PetscCall(VecSetBlockSize(*newvec, dd->w)); in DMDAGetRay()
125 PetscCall(VecSetSizes(*newvec, dd->w, PETSC_DETERMINE)); in DMDAGetRay()
127 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, dd->w, indices, PETSC_OWN_POINTER, &is)); in DMDAGetRay()
134 PetscCall(PetscMalloc1(dd->w * dd->M, &indices)); in DMDAGetRay()
135 indices[0] = gp * dd->M * dd->w; in DMDAGetRay()
[all …]
H A Dda1.c13 DM_DA *dd = (DM_DA *)da->data; in DMView_DA_1d() local
55 …m %" PetscInt_FMT " w %" PetscInt_FMT " s %" PetscInt_FMT "\n", rank, dd->M, dd->m, dd->w, dd->s)); in DMView_DA_1d()
62 double ymin = -1, ymax = 1, xmin = -1, xmax = dd->M, x; in DMView_DA_1d()
81 …for (xmin_tmp = 0; xmin_tmp < dd->M; xmin_tmp++) PetscCall(PetscDrawLine(draw, (double)xmin_tmp, y… in DMView_DA_1d()
83 xmax = dd->M - 1; in DMView_DA_1d()
95 xmin = dd->xs / dd->w; in DMView_DA_1d()
96 xmax = (dd->xe / dd->w) - 1; in DMView_DA_1d()
102 base = dd->base / dd->w; in DMView_DA_1d()
125 DM_DA *dd = (DM_DA *)da->data; in DMSetUp_DA_1D() local
126 const PetscInt M = dd->M; in DMSetUp_DA_1D()
[all …]
H A Ddadist.c9 DM_DA *dd = (DM_DA *)da->data; in DMCreateGlobalVector_DA() local
15 PetscCall(VecSetSizes(*g, dd->Nlocal, PETSC_DETERMINE)); in DMCreateGlobalVector_DA()
16 PetscCall(VecSetBlockSize(*g, dd->w)); in DMCreateGlobalVector_DA()
18 if (dd->Nlocal < da->bind_below) { in DMCreateGlobalVector_DA()
61 DM_DA *dd = (DM_DA *)da->data; in DMDACreateNaturalVector() local
66 if (dd->natural) { in DMDACreateNaturalVector()
67 PetscCall(PetscObjectGetReference((PetscObject)dd->natural, &cnt)); in DMDACreateNaturalVector()
69 PetscCall(PetscObjectReference((PetscObject)dd->natural)); in DMDACreateNaturalVector()
70 *g = dd->natural; in DMDACreateNaturalVector()
71 } else PetscCall(VecDuplicate(dd->natural, g)); in DMDACreateNaturalVector()
[all …]
H A Ddalocal.c47 DM_DA *dd = (DM_DA *)da->data; in DMCreateLocalVector_DA() local
53 PetscCall(VecSetSizes(*g, dd->nlocal, PETSC_DETERMINE)); in DMCreateLocalVector_DA()
54 PetscCall(VecSetBlockSize(*g, dd->w)); in DMCreateLocalVector_DA()
56 if (dd->nlocal < da->bind_below) { in DMCreateLocalVector_DA()
62 …if (dd->w == 1 && da->dim == 2) PetscCall(PetscObjectComposeFunction((PetscObject)*g, "PetscMatlab… in DMCreateLocalVector_DA()
413 DM_DA *dd = (DM_DA *)da->data; in DMDAGetArray() local
419 if (dd->arrayghostedin[i]) { in DMDAGetArray()
420 *iptr = dd->arrayghostedin[i]; in DMDAGetArray()
421 iarray_start = (char *)dd->startghostedin[i]; in DMDAGetArray()
422 dd->arrayghostedin[i] = NULL; in DMDAGetArray()
[all …]
H A Ddareg.c5 DM_DA *dd = (DM_DA *)da->data; in DMSetUp_DA() local
9 …Check(dd->w >= 1, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "Must have 1 or more… in DMSetUp_DA()
10 …PetscCheck(dd->s >= 0, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "Stencil width … in DMSetUp_DA()
12 PetscCall(PetscCalloc1(dd->w + 1, &dd->fieldname)); in DMSetUp_DA()
13 PetscCall(PetscCalloc1(da->dim, &dd->coordinatename)); in DMSetUp_DA()
21 da->Nf = dd->w; in DMSetUp_DA()
H A Dda2.c8 DM_DA *dd = (DM_DA *)da->data; in DMView_DA_2d() local
50 …t_FMT " w %" PetscInt_FMT " s %" PetscInt_FMT "\n", rank, dd->M, dd->N, dd->m, dd->n, dd->w, dd->s… in DMView_DA_2d()
57 double ymin = -1 * dd->s - 1, ymax = dd->N + dd->s; in DMView_DA_2d()
58 double xmin = -1 * dd->s - 1, xmax = dd->M + dd->s; in DMView_DA_2d()
77 ymax = dd->N - 1; in DMView_DA_2d()
78 …for (xmin = 0; xmin < dd->M; xmin++) PetscCall(PetscDrawLine(draw, xmin, ymin, xmin, ymax, PETSC_D… in DMView_DA_2d()
80 xmax = dd->M - 1; in DMView_DA_2d()
81 …for (ymin = 0; ymin < dd->N; ymin++) PetscCall(PetscDrawLine(draw, xmin, ymin, xmax, ymin, PETSC_D… in DMView_DA_2d()
89 xmin = dd->xs / dd->w; in DMView_DA_2d()
90 xmax = (dd->xe - 1) / dd->w; in DMView_DA_2d()
[all …]
H A Dda3.c13 DM_DA *dd = (DM_DA *)da->data; in DMView_DA_3d() local
56 … w %" PetscInt_FMT " s %" PetscInt_FMT "\n", rank, dd->M, dd->N, dd->P, dd->m, dd->n, dd->p, dd->w… in DMView_DA_3d()
76 PetscReal ymin = -1.0, ymax = (PetscReal)dd->N; in DMView_DA_3d()
77 PetscReal xmin = -1.0, xmax = (PetscReal)((dd->M + 2) * dd->P), x, y, ycoord, xcoord; in DMView_DA_3d()
94 for (k = 0; k < dd->P; k++) { in DMView_DA_3d()
96 ymax = (PetscReal)(dd->N - 1); in DMView_DA_3d()
97 …for (xmin = (PetscReal)(k * (dd->M + 1)); xmin < (PetscReal)(dd->M + (k * (dd->M + 1))); xmin++) P… in DMView_DA_3d()
98 xmin = (PetscReal)(k * (dd->M + 1)); in DMView_DA_3d()
99 xmax = xmin + (PetscReal)(dd->M - 1); in DMView_DA_3d()
100 …for (ymin = 0; ymin < (PetscReal)dd->N; ymin++) PetscCall(PetscDrawLine(draw, xmin, ymin, xmax, ym… in DMView_DA_3d()
[all …]
H A Ddagtona.c34 DM_DA *dd = (DM_DA *)da->data; in DMDAGlobalToNaturalAllCreate() local
42 …PetscCall(VecCreateMPIWithArray(PetscObjectComm((PetscObject)da), dd->w, dd->Nlocal, PETSC_DETERMI… in DMDAGlobalToNaturalAllCreate()
47 PetscCall(VecCreateSeqWithArray(PETSC_COMM_SELF, dd->w, N, NULL, &tmplocal)); in DMDAGlobalToNaturalAllCreate()
75 DM_DA *dd = (DM_DA *)da->data; in DMDANaturalAllToGlobalCreate() local
76 PetscInt M, m = dd->Nlocal, start; in DMDANaturalAllToGlobalCreate()
88 …PetscCall(VecCreateMPIWithArray(PetscObjectComm((PetscObject)da), dd->w, m, PETSC_DETERMINE, NULL,… in DMDANaturalAllToGlobalCreate()
93 PetscCall(VecCreateSeqWithArray(PETSC_COMM_SELF, dd->w, M, NULL, &tmplocal)); in DMDANaturalAllToGlobalCreate()
H A Ddascatter.c30 DM_DA *dd = (DM_DA *)da->data; in DMDAGetScatter() local
34 if (gtol) *gtol = dd->gtol; in DMDAGetScatter()
36 if (!dd->ltol) PetscCall(DMLocalToLocalCreate_DA(da)); in DMDAGetScatter()
37 *ltol = dd->ltol; in DMDAGetScatter()
H A Ddagetelem.c302 DM_DA *dd = (DM_DA *)da->data; in DMDASetElementType() local
310 if (dd->elementtype != etype) { in DMDASetElementType()
311 PetscCall(PetscFree(dd->e)); in DMDASetElementType()
312 PetscCall(ISDestroy(&dd->ecorners)); in DMDASetElementType()
314 dd->elementtype = etype; in DMDASetElementType()
315 dd->ne = 0; in DMDASetElementType()
316 dd->nen = 0; in DMDASetElementType()
317 dd->e = NULL; in DMDASetElementType()
340 DM_DA *dd = (DM_DA *)da->data; in DMDAGetElementType() local
348 *etype = dd->elementtype; in DMDAGetElementType()
[all …]
H A Dgr2.c587 DM_DA *dd = (DM_DA *)da->data; in DMDAArrayMPIIO() local
608 PetscCall(PetscMPIIntCast(dd->w, &dof)); in DMDAArrayMPIIO()
610 PetscCall(PetscMPIIntCast(dd->M, gsizes + 1)); in DMDAArrayMPIIO()
611 PetscCall(PetscMPIIntCast(dd->N, gsizes + 2)); in DMDAArrayMPIIO()
612 PetscCall(PetscMPIIntCast(dd->P, gsizes + 3)); in DMDAArrayMPIIO()
614 PetscCall(PetscMPIIntCast((dd->xe - dd->xs) / dof, lsizes + 1)); in DMDAArrayMPIIO()
615 PetscCall(PetscMPIIntCast(dd->ye - dd->ys, lsizes + 2)); in DMDAArrayMPIIO()
616 PetscCall(PetscMPIIntCast(dd->ze - dd->zs, lsizes + 3)); in DMDAArrayMPIIO()
618 PetscCall(PetscMPIIntCast(dd->xs / dof, lstarts + 1)); in DMDAArrayMPIIO()
619 PetscCall(PetscMPIIntCast(dd->ys, lstarts + 2)); in DMDAArrayMPIIO()
[all …]
/petsc/src/dm/impls/stag/tests/
H A Dex7.c79 PetscInt dd, extra[3]; in main() local
89 …for (dd = 0; dd < 3; ++dd) expected *= (bnd[dd] ? stencilWidth + 1 + extra[dd] : 2 * stencilWidth … in main()
98 …for (dd = 1; dd < 3; ++dd) expected *= (bnd[dd] ? stencilWidth + 1 + extra[dd] : 2 * stencilWidth … in main()
107 …for (dd = 0; dd < 3; dd += 2) expected *= (bnd[dd] ? stencilWidth + 1 + extra[dd] : 2 * stencilWid… in main()
116 for (dd = 0; dd < 2; ++dd) expected *= ((bnd[dd] ? 1 : 2) * stencilWidth + 1); in main()
125 …for (dd = 0; dd < 2; ++dd) expected *= (bnd[dd] ? stencilWidth + 1 + extra[dd] : 2 * stencilWidth … in main()
134 for (dd = 0; dd < 3; dd += 2) expected *= ((bnd[dd] ? 1 : 2) * stencilWidth + 1); in main()
143 for (dd = 1; dd < 3; ++dd) expected *= ((bnd[dd] ? 1 : 2) * stencilWidth + 1); in main()
152 for (dd = 0; dd < 3; ++dd) expected *= ((bnd[dd] ? 1 : 2) * stencilWidth + 1); in main()
H A Dex10.c73 PetscInt dd, extra[2]; in main() local
81 …for (dd = 0; dd < 2; ++dd) expected *= (bnd[dd] ? stencilWidth + 1 + extra[dd] : 2 * stencilWidth … in main()
108 for (dd = 0; dd < 2; ++dd) expected *= ((bnd[dd] ? 1 : 2) * stencilWidth + 1); in main()
/petsc/src/ksp/ksp/impls/cg/
H A Dcgeig.c46 PetscReal *dd, *ee; in KSPComputeExtremeSingularValues_CG() local
57 dd = cgP->dd; in KSPComputeExtremeSingularValues_CG()
62 dd[j] = PetscRealPart(d[j]); in KSPComputeExtremeSingularValues_CG()
68 PetscCallBLAS("LAPACKREALstev", LAPACKREALstev_("N", &bn, dd, ee, NULL, &ldz, NULL, &lierr)); in KSPComputeExtremeSingularValues_CG()
71 *emin = dd[0]; in KSPComputeExtremeSingularValues_CG()
72 *emax = dd[n - 1]; in KSPComputeExtremeSingularValues_CG()
/petsc/src/snes/tests/
H A Dex20.c139 PetscScalar tleft, tright, beta, td, ad, dd, fd = 0.0, tu, au, du = 0.0, fu = 0.0; in FormFunction() local
196 dd = PetscPowScalar(ad, beta); in FormFunction()
197 fd = dd * (t0 - td); in FormFunction()
237 dd = PetscPowScalar(ad, beta); in FormFunction()
238 fd = dd * (t0 - td); in FormFunction()
285 dd = PetscPowScalar(ad, beta); in FormFunction()
286 fd = dd * (t0 - td); in FormFunction()
322 dd = PetscPowScalar(ad, beta); in FormFunction()
323 fd = dd * (t0 - td); in FormFunction()
359 dd = PetscPowScalar(ad, beta); in FormFunction()
[all …]

12