Home
last modified time | relevance | path

Searched refs:xi (Results 1 – 25 of 63) sorted by relevance

123

/petsc/src/dm/impls/swarm/
H A Dswarmpic_da.c8 PetscReal *xi; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular() local
27 PetscCall(PetscMalloc1(dim * npoints, &xi)); in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
32 xi[dim * cnt + 0] = -1.0 + 0.5 * ds[d] + ii * ds[0]; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
41 xi[dim * cnt + 0] = -1.0 + 0.5 * ds[0] + ii * ds[0]; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
42 xi[dim * cnt + 1] = -1.0 + 0.5 * ds[1] + jj * ds[1]; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
53 xi[dim * cnt + 0] = -1.0 + 0.5 * ds[0] + ii * ds[0]; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
54 xi[dim * cnt + 1] = -1.0 + 0.5 * ds[1] + jj * ds[1]; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
55 xi[dim * cnt + 2] = -1.0 + 0.5 * ds[2] + kk * ds[2]; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
63 *_xi = xi; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
71 PetscReal *xi; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Gauss() local
[all …]
H A Dswarmpic_plex.c17 PetscErrorCode private_DMSwarmSetPointCoordinatesCellwise_PLEX(DM, DM, PetscInt, PetscReal *xi);
161 PetscReal *xi, ds, ds2; in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular() local
181 PetscCall(PetscMalloc1(dim * npoints * npoints, &xi)); in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
187 xi[dim * pcnt + 0] = ii * ds; in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
188 xi[dim * pcnt + 1] = jj * ds; in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
190 xi[dim * pcnt + 0] *= (1.0 - 1.2 * ds2); in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
191 xi[dim * pcnt + 1] *= (1.0 - 1.2 * ds2); in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
193 xi[dim * pcnt + 0] += 0.35 * ds2; in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
194 xi[dim * pcnt + 1] += 0.35 * ds2; in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
205 basis[q][0] = 1.0 - xi[dim * q + 0] - xi[dim * q + 1]; in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
[all …]
/petsc/src/tao/bound/tutorials/
H A Djbearing2.c41 static PetscReal p(PetscReal xi, PetscReal ecc);
166 static PetscReal p(PetscReal xi, PetscReal ecc) in p() argument
168 PetscReal t = 1.0 + ecc * PetscCosScalar(xi); in p()
215 PetscReal xi, v[5]; in FormFunctionGradient() local
247 xi = (i + 1) * hx; in FormFunctionGradient()
248 trule1 = hxhy * (p(xi, ecc) + p(xi + hx, ecc) + p(xi, ecc)) / six; /* L(i,j) */ in FormFunctionGradient()
249 trule2 = hxhy * (p(xi, ecc) + p(xi - hx, ecc) + p(xi, ecc)) / six; /* U(i,j) */ in FormFunctionGradient()
250 trule3 = hxhy * (p(xi, ecc) + p(xi + hx, ecc) + p(xi + hx, ecc)) / six; /* U(i+1,j) */ in FormFunctionGradient()
251 trule4 = hxhy * (p(xi, ecc) + p(xi - hx, ecc) + p(xi - hx, ecc)) / six; /* L(i-1,j) */ in FormFunctionGradient()
330 PetscReal xi, v[5]; in FormHessian() local
[all …]
/petsc/src/ts/tutorials/
H A Dex43.c12 static void Exact(PetscReal t, PetscReal omega, PetscReal xi, PetscReal u0, PetscReal v0, PetscReal… in Exact() argument
15 if (xi < 1) { in Exact()
16 PetscReal a = xi * omega; in Exact()
17 PetscReal w = PetscSqrtReal(1 - xi * xi) * omega; in Exact()
22 } else if (xi > 1) { in Exact()
23 PetscReal w = PetscSqrtReal(xi * xi - 1) * omega; in Exact()
24 PetscReal C1 = (w * u0 + xi * u0 + v0) / (2 * w); in Exact()
25 PetscReal C2 = (w * u0 - xi * u0 - v0) / (2 * w); in Exact()
26 u = C1 * PetscExpReal((-xi + w) * t) + C2 * PetscExpReal((-xi - w) * t); in Exact()
27 …v = C1 * (-xi + w) * PetscExpReal((-xi + w) * t) + C2 * (-xi - w) * PetscExpReal((-xi -… in Exact()
[all …]
H A Dex25.c292 PetscReal xi = i * hx; in FormInitialSolution() local
293 x[i].u = user->uleft * (1. - xi) + user->uright * xi + PetscSinReal(2. * PETSC_PI * xi); in FormInitialSolution()
294 x[i].v = user->vleft * (1. - xi) + user->vright * xi; in FormInitialSolution()
H A Dex35.cxx379 PetscReal xi = vpos[0]; in FormInitialSolution() local
380 …x[dof].u = user->leftbc.u * (1. - xi) + user->rightbc.u * xi + PetscSinReal(2. * PETSC_PI * xi… in FormInitialSolution()
381 x[dof].v = user->leftbc.v * (1. - xi) + user->rightbc.v * xi; in FormInitialSolution()
H A Dex34.c28 PetscReal xi; /* interblock spring constant */ member
55 …f[i].v = -hx * (PetscSqr(user->gammaTilde) * u[i].u + (PetscSqr(user->gamma) / user->xi) * (u[i].… in FormRHSFunction()
231 user.xi = 0.5; in main()
236 PetscCall(PetscOptionsReal("-xi", "Interblock spring constant", "", user.xi, &user.xi, NULL)); in main()
/petsc/src/mat/impls/sbaij/seq/
H A Drelax.h9 #define PetscSparseDensePlusDot_no_function(sum, r, xv, xi, nnz) \ argument
15 sum += *xv++ * r[*xi++]; \
17 sum += *xv++ * r[*xi++]; \
19 sum += *xv++ * r[*xi++]; \
23 sum += xv[0] * r[xi[0]] + xv[1] * r[xi[1]] + xv[2] * r[xi[2]] + xv[3] * r[xi[3]]; \
25 xi += 4; \
29 xi -= nnz; \
34 #define PetscSparseDensePlusDot_no_function(sum, r, xv, xi, nnz) \ argument
38 __i1 = xi[__i]; \
39 __i2 = xi[__i + 1]; \
[all …]
/petsc/src/mat/impls/aij/seq/
H A Daij.h505 #define PetscSparseDenseMinusDot(sum, r, xv, xi, nnz) \ argument
511 sum -= *xv++ * r[*xi++]; \
513 sum -= *xv++ * r[*xi++]; \
515 sum -= *xv++ * r[*xi++]; \
519 sum -= xv[0] * r[xi[0]] + xv[1] * r[xi[1]] + xv[2] * r[xi[2]] + xv[3] * r[xi[3]]; \
521 xi += 4; \
525 xi -= nnz; \
530 #define PetscSparseDenseMinusDot(sum, r, xv, xi, nnz) \ argument
534 __i1 = xi[__i]; \
535 __i2 = xi[__i + 1]; \
[all …]
/petsc/src/dm/tests/
H A Dex36.c87 PetscScalar xi = XX[dim * i]; in DAApplyConformalMapping() local
95 Ni[0] = 0.125 * (1.0 - xi) * (1.0 - eta) * (1.0 - zeta); in DAApplyConformalMapping()
96 Ni[1] = 0.125 * (1.0 + xi) * (1.0 - eta) * (1.0 - zeta); in DAApplyConformalMapping()
97 Ni[2] = 0.125 * (1.0 - xi) * (1.0 + eta) * (1.0 - zeta); in DAApplyConformalMapping()
98 Ni[3] = 0.125 * (1.0 + xi) * (1.0 + eta) * (1.0 - zeta); in DAApplyConformalMapping()
100 Ni[4] = 0.125 * (1.0 - xi) * (1.0 - eta) * (1.0 + zeta); in DAApplyConformalMapping()
101 Ni[5] = 0.125 * (1.0 + xi) * (1.0 - eta) * (1.0 + zeta); in DAApplyConformalMapping()
102 Ni[6] = 0.125 * (1.0 - xi) * (1.0 + eta) * (1.0 + zeta); in DAApplyConformalMapping()
103 Ni[7] = 0.125 * (1.0 + xi) * (1.0 + eta) * (1.0 + zeta); in DAApplyConformalMapping()
171 PetscScalar xi = XX[dim * i]; in DAApplyConformalMapping() local
[all …]
/petsc/src/ksp/ksp/tutorials/
H A Dex70.c46 PetscScalar xi = _xi[0]; in EvaluateBasis_Q1() local
49 N[0] = 0.25 * (1.0 - xi) * (1.0 - eta); in EvaluateBasis_Q1()
50 N[1] = 0.25 * (1.0 + xi) * (1.0 - eta); in EvaluateBasis_Q1()
51 N[2] = 0.25 * (1.0 + xi) * (1.0 + eta); in EvaluateBasis_Q1()
52 N[3] = 0.25 * (1.0 - xi) * (1.0 + eta); in EvaluateBasis_Q1()
57 PetscScalar xi = _xi[0]; in EvaluateBasisDerivatives_Q1() local
65 dN[1][0] = -0.25 * (1.0 - xi); in EvaluateBasisDerivatives_Q1()
66 dN[1][1] = -0.25 * (1.0 + xi); in EvaluateBasisDerivatives_Q1()
67 dN[1][2] = 0.25 * (1.0 + xi); in EvaluateBasisDerivatives_Q1()
68 dN[1][3] = 0.25 * (1.0 - xi); in EvaluateBasisDerivatives_Q1()
[all …]
H A Dex42.c104 PetscReal xi = PetscRealPart(_xi[0]); in ShapeFunctionQ13D_Evaluate() local
108 Ni[0] = 0.125 * (1.0 - xi) * (1.0 - eta) * (1.0 - zeta); in ShapeFunctionQ13D_Evaluate()
109 Ni[1] = 0.125 * (1.0 - xi) * (1.0 + eta) * (1.0 - zeta); in ShapeFunctionQ13D_Evaluate()
110 Ni[2] = 0.125 * (1.0 + xi) * (1.0 + eta) * (1.0 - zeta); in ShapeFunctionQ13D_Evaluate()
111 Ni[3] = 0.125 * (1.0 + xi) * (1.0 - eta) * (1.0 - zeta); in ShapeFunctionQ13D_Evaluate()
113 Ni[4] = 0.125 * (1.0 - xi) * (1.0 - eta) * (1.0 + zeta); in ShapeFunctionQ13D_Evaluate()
114 Ni[5] = 0.125 * (1.0 - xi) * (1.0 + eta) * (1.0 + zeta); in ShapeFunctionQ13D_Evaluate()
115 Ni[6] = 0.125 * (1.0 + xi) * (1.0 + eta) * (1.0 + zeta); in ShapeFunctionQ13D_Evaluate()
116 Ni[7] = 0.125 * (1.0 + xi) * (1.0 - eta) * (1.0 + zeta); in ShapeFunctionQ13D_Evaluate()
121 PetscReal xi = PetscRealPart(_xi[0]); in ShapeFunctionQ13D_Evaluate_dxi() local
[all …]
/petsc/src/snes/utils/dm/
H A Ddminterpolatesnes.c401 PetscReal xi[4]; in DMInterpolate_Triangle_Private() local
412 xi[d] = 0.0; in DMInterpolate_Triangle_Private()
413 …for (PetscInt f = 0; f < ctx->dim; ++f) xi[d] += invJ[d * ctx->dim + f] * 0.5 * PetscRealPart(coor… in DMInterpolate_Triangle_Private()
414 …* ctx->dof + comp] += PetscRealPart(x[(d + 1) * ctx->dof + comp] - x[0 * ctx->dof + comp]) * xi[d]; in DMInterpolate_Triangle_Private()
431 PetscReal xi[4]; in DMInterpolate_Tetrahedron_Private() local
442 xi[d] = 0.0; in DMInterpolate_Tetrahedron_Private()
443 …for (PetscInt f = 0; f < ctx->dim; ++f) xi[d] += invJ[d * ctx->dim + f] * 0.5 * PetscRealPart(coor… in DMInterpolate_Tetrahedron_Private()
444 … ctx->dof + comp] += PetscRealPart(x[order[d] * ctx->dof + comp] - x[0 * ctx->dof + comp]) * xi[d]; in DMInterpolate_Tetrahedron_Private()
543 PetscScalar *xi; in DMInterpolate_Quad_Private() local
586 PetscCall(VecGetArray(real, &xi)); in DMInterpolate_Quad_Private()
[all …]
/petsc/src/dm/dt/fe/interface/
H A Dfegeom.c29 g->xi = p; in PetscFEGeomCreate()
101 (*chunkGeom)->xi = geom->xi; in PetscFEGeomGetChunk()
173 pgeom->xi = geom->xi; in PetscFEGeomGetPoint()
179 …if (pcoords) CoordinatesRefToReal(dE, dim, pgeom->xi, &geom->v[c * Np * dE], pgeom->J, pcoords, pg… in PetscFEGeomGetPoint()
/petsc/src/ts/tests/
H A Dex25.c314 PetscReal xi = i * hx; in FormInitialSolution() local
315 x[i].u = user->uleft * (1. - xi) + user->uright * xi + PetscSinReal(2. * PETSC_PI * xi); in FormInitialSolution()
316 x[i].v = user->vleft * (1. - xi) + user->vright * xi; in FormInitialSolution()
H A Dex27.c316 PetscReal xi[3], xi2, xi3, mag; in QCompute() local
320 DMPlex_WaxpyD_Internal(dim, -1.0, vq, vp, xi); in QCompute()
321 xi2 = DMPlex_DotD_Internal(dim, xi, xi); in QCompute()
325 for (e = 0; e < dim; ++e) Q[d * dim + e] = -xi[d] * xi[e] / xi3; in QCompute()
/petsc/src/snes/impls/ngmres/
H A Dngmresfunc.c45 PetscScalar *xi = ngmres->xi; in SNESNGMRESFormCombinedSolution_Private() local
56 PetscCall(VecMDotBegin(FM, l, Fdot, xi)); in SNESNGMRESFormCombinedSolution_Private()
58 PetscCall(VecMDotEnd(FM, l, Fdot, xi)); in SNESNGMRESFormCombinedSolution_Private()
68 for (i = 0; i < l; i++) beta[i] = nu - xi[i]; in SNESNGMRESFormCombinedSolution_Private()
72 for (i = 0; i < l; i++) H(i, j) = Q(i, j) - xi[i] - xi[j] + nu; in SNESNGMRESFormCombinedSolution_Private()
/petsc/src/ts/tutorials/power_grid/
H A Dex7.c149 PetscScalar xi, yi; in ini_bou() local
175 xi = coors[j][i].x; in ini_bou()
177 if ((xi == mux) && (yi == muy)) p[j][i] = 1.0; in ini_bou()
185 xi = coors[j][i].x; in ini_bou()
187 …/ (1 - rho * rho) * (PetscPowScalar((xi - mux) / sigmax, 2) + PetscPowScalar((yi - muy) / sigmay, … in ini_bou()
307 PetscScalar val[5], xi, yi; in IJacobian() local
321 xi = coors[j][i].x; in IJacobian()
328 …c3 = (user->ws / (2.0 * user->H)) * (user->PM_min - user->Pmax * PetscSinScalar(xi)) / us… in IJacobian()
H A Dex6.c131 PetscScalar xi, yi; in ini_bou() local
149 xi = coors[j][i].x; in ini_bou()
153 …/ (1 - rho * rho) * (PetscPowScalar((xi - mux) / sigmax, 2) + PetscPowScalar((yi - muy) / sigmay, … in ini_bou()
312 PetscScalar val[5], xi, yi; in IJacobian() local
326 xi = coors[j][i].x; in IJacobian()
337 fthetac = user->ws / (2 * user->H) * (user->PM_min - user->Pmax * PetscSinScalar(xi)); in IJacobian()
435 …c3 = (user->ws / (2.0 * user->H)) * (user->PM_min - user->Pmax * PetscSinScalar(xi)) / (2 *… in IJacobian()
/petsc/src/dm/impls/da/
H A Ddainterp.c117 PetscScalar *xi; in DMCreateInterpolation_DA_1D_Q1() local
123 PetscCall(PetscMalloc1(nxi, &xi)); in DMCreateInterpolation_DA_1D_Q1()
124 for (li = 0; li < nxi; li++) xi[li] = -1.0 + (PetscScalar)li * (2.0 / (PetscScalar)(nxi - 1)); in DMCreateInterpolation_DA_1D_Q1()
152 Ni[0] = 0.5 * (1.0 - xi[li]); in DMCreateInterpolation_DA_1D_Q1()
153 Ni[1] = 0.5 * (1.0 + xi[li]); in DMCreateInterpolation_DA_1D_Q1()
159 PetscCall(PetscFree(xi)); in DMCreateInterpolation_DA_1D_Q1()
407 PetscScalar *xi, *eta; in DMCreateInterpolation_DA_2D_Q1() local
413 PetscCall(PetscMalloc1(nxi, &xi)); in DMCreateInterpolation_DA_2D_Q1()
415 for (li = 0; li < nxi; li++) xi[li] = -1.0 + (PetscScalar)li * (2.0 / (PetscScalar)(nxi - 1)); in DMCreateInterpolation_DA_2D_Q1()
455 Ni[0] = 0.25 * (1.0 - xi[li]) * (1.0 - eta[lj]); in DMCreateInterpolation_DA_2D_Q1()
[all …]
/petsc/src/ts/impls/explicit/euler/
H A Deuler.c98 static PetscErrorCode TSComputeLinearStability_Euler(TS ts, PetscReal xr, PetscReal xi, PetscReal *… in TSComputeLinearStability_Euler() argument
102 *yi = xi; in TSComputeLinearStability_Euler()
/petsc/src/dm/field/impls/ds/
H A Ddmfieldds.c1008 for (q = 0; q < Nq; q++) orientPoints[o][q] = -geom->xi[q]; in DMFieldComputeFaceData_DS()
1010 for (q = 0; q < Nq; q++) orientPoints[o][q] = geom->xi[q]; in DMFieldComputeFaceData_DS()
1019 lambda[0] = -(geom->xi[2 * q] + geom->xi[2 * q + 1]) / 2.; in DMFieldComputeFaceData_DS()
1020 lambda[1] = (geom->xi[2 * q] + 1.) / 2.; in DMFieldComputeFaceData_DS()
1021 lambda[2] = (geom->xi[2 * q + 1] + 1.) / 2.; in DMFieldComputeFaceData_DS()
1034 PetscReal xi[2], xio[2]; in DMFieldComputeFaceData_DS() local
1037 xi[0] = geom->xi[2 * q]; in DMFieldComputeFaceData_DS()
1038 xi[1] = geom->xi[2 * q + 1]; in DMFieldComputeFaceData_DS()
1041 xio[0] = xi[1]; in DMFieldComputeFaceData_DS()
1042 xio[1] = -xi[0]; in DMFieldComputeFaceData_DS()
[all …]
/petsc/src/sys/classes/draw/interface/
H A Ddtri.c130 PetscErrorCode PetscDrawTensorContour(PetscDraw draw, int m, int n, const PetscReal xi[], const Pet… in PetscDrawTensorContour() argument
169 if (!xi) { in PetscDrawTensorContour()
175 } else ctx.x = (PetscReal *)xi; in PetscDrawTensorContour()
/petsc/src/dm/impls/plex/
H A Dplexply.c26 PetscInt dim = 2, cdim = 3, Nvp = 0, coordSize, xi = -1, yi = -1, zi = -1, v, c, p; in DMPlexCreatePLYFromFile() local
99 if (matchB) xi = Nvp; in DMPlexCreatePLYFromFile()
161 if (p == xi) coords[v * cdim + 0] = rbuf[0]; in DMPlexCreatePLYFromFile()
/petsc/src/ksp/ksp/utils/dm/
H A Ddmproject.c944 PetscInt xi[3], vi[3]; in InitializeParticles_Regular() local
946 xi[0] = q % n; in InitializeParticles_Regular()
947 xi[1] = (q / n) % n; in InitializeParticles_Regular()
948 xi[2] = (q / PetscSqr(n)) % n; in InitializeParticles_Regular()
949 …for (PetscInt d = 0; d < xdim; ++d) x[p * xdim + d] = lower[d] + (xi[d] + 0.5) * (upper[d] - lower… in InitializeParticles_Regular()
1097 for (PetscInt xi = xs[0] - 1; xi < xs[0] + 3; ++xi) { in DMSwarmRemap_Colella_Internal() local
1099 const PetscInt xip = xi < 0 ? xi + xend[0] : (xi >= xend[0] ? xi - xend[0] : xi); in DMSwarmRemap_Colella_Internal()
1100 …PetscCall(W_3_Interpolation_Private((xmin[0] + (xi + 0.5) * hx[0] - x[p * dim + 0]) / hx[0], &Wx[0… in DMSwarmRemap_Colella_Internal()
1114 …t_FMT ", %" PetscInt_FMT ", %" PetscInt_FMT ") w = %g (%g, %g, %g, %g)\n", xi, xj, vi, vj, (double… in DMSwarmRemap_Colella_Internal()

123