Home
last modified time | relevance | path

Searched refs:q (Results 1 – 25 of 164) sorted by relevance

1234567

/petsc/src/dm/impls/swarm/
H A Dswarmpic_da.c72 PetscInt d, q, npoints_q; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Gauss() local
78 for (q = 0; q < npoints_q; q++) { in private_DMSwarmCreateCellLocalCoords_DA_Q1_Gauss()
79 for (d = 0; d < dim; d++) xi[dim * q + d] = quadrature_xi[dim * q + d]; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Gauss()
91 PetscInt nel, npe, e, q, k, d; in private_DMSwarmInsertPointsUsingCellDM_DA_Q1() local
132 for (q = 0; q < npoints_q; q++) { in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
133 PetscCall(PetscMalloc1(npe, &basis[q])); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
137 basis[q][0] = 0.5 * (1.0 - xi[dim * q + 0]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
138 basis[q][1] = 0.5 * (1.0 + xi[dim * q + 0]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
141 basis[q][0] = 0.25 * (1.0 - xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
142 basis[q][1] = 0.25 * (1.0 + xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
[all …]
H A Dswarmpic_plex.c22 PetscQuadrature q, fq; in private_PetscFECreateDefault_scalar_pk1() local
70 PetscCall(PetscDTStroudConicalQuadrature(dim, 1, quadPointsPerEdge, -1.0, 1.0, &q)); in private_PetscFECreateDefault_scalar_pk1()
73 PetscCall(PetscDTGaussTensorQuadrature(dim, 1, quadPointsPerEdge, -1.0, 1.0, &q)); in private_PetscFECreateDefault_scalar_pk1()
76 PetscCall(PetscFESetQuadrature(*fem, q)); in private_PetscFECreateDefault_scalar_pk1()
78 PetscCall(PetscQuadratureDestroy(&q)); in private_PetscFECreateDefault_scalar_pk1()
86 PetscInt q, npoints_q, e, nel, pcnt, ps, pe, d, k, r, Nfc; in private_DMSwarmInsertPointsUsingCellDM_PLEX_SubDivide() local
140 for (q = 0; q < npoints_q; q++) { in private_DMSwarmInsertPointsUsingCellDM_PLEX_SubDivide()
143 …for (k = 0; k < nbasis; k++) swarm_coor[dim * pcnt + d] += T->T[0][q * nbasis + k] * PetscRealPart… in private_DMSwarmInsertPointsUsingCellDM_PLEX_SubDivide()
160 PetscInt ii, jj, q, npoints_q, e, nel, npe, pcnt, ps, pe, d, k, nfaces, Nfc; in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular() local
202 for (q = 0; q < npoints_q; q++) { in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
[all …]
/petsc/src/ts/tutorials/network/
H A DpipeImpls.c26 x[i].q = Q0; in PipeComputeSteadyState()
41 return (x[i + 1].q - x[i].q) / dx; in dqdx()
43 return (x[i].q - x[i - 1].q) / dx; in dqdx()
45 return (x[i + 1].q - x[i - 1].q) / (2 * dx); in dqdx()
74 qavg = (xold[i + 1].q + xold[i - 1].q) / 2.0; in PipeIFunctionLocal_Lax()
75 qa = xold[i - 1].q; in PipeIFunctionLocal_Lax()
76 qb = xold[i + 1].q; in PipeIFunctionLocal_Lax()
82 …f[2 * (i - 1) + 2] = (x[i].q - xold_i) + dt * (GRAV * pipe->A * dhdx(xold, i, ilast, dx) + pipe->R… in PipeIFunctionLocal_Lax()
90 …start].q - xold[start + 1].q - ((GRAV * A) / a) * (x[start].h - xold[start + 1].h) + dt * R * xold… in PipeIFunctionLocal_Lax()
91 …ilast].q - xold[ilast - 1].q + ((GRAV * A) / a) * (x[ilast].h - xold[ilast - 1].h) + dt * R * xold… in PipeIFunctionLocal_Lax()
/petsc/src/dm/impls/plex/
H A Dpointqueue.c5 DMPlexPointQueue q; in DMPlexPointQueueCreate() local
9 PetscCall(PetscCalloc1(1, &q)); in DMPlexPointQueueCreate()
10 q->size = size; in DMPlexPointQueueCreate()
11 PetscCall(PetscMalloc1(q->size, &q->points)); in DMPlexPointQueueCreate()
12 q->num = 0; in DMPlexPointQueueCreate()
13 q->front = 0; in DMPlexPointQueueCreate()
14 q->back = q->size - 1; in DMPlexPointQueueCreate()
15 *queue = q; in DMPlexPointQueueCreate()
21 DMPlexPointQueue q = *queue; in DMPlexPointQueueDestroy() local
24 PetscCall(PetscFree(q->points)); in DMPlexPointQueueDestroy()
[all …]
H A Dplexpreallocate.c23 PetscInt p, q, a, aSize, *offsets, aStart, aEnd, *inverse, iSize, *adj, adjSize; in DMPlexComputeAnchorAdjacencies() local
49 for (q = 0; q < dof; q++) { in DMPlexComputeAnchorAdjacencies()
52 a = anchors[off + q]; in DMPlexComputeAnchorAdjacencies()
87 q = inverse[iOff + i]; in DMPlexComputeAnchorAdjacencies()
88 …PetscCall(DMPlexGetAdjacency_Internal(dm, q, useCone, useClosure, PETSC_TRUE, &numAdjQ, &tmpAdjQ)); in DMPlexComputeAnchorAdjacencies()
121 q = inverse[iOff + i]; in DMPlexComputeAnchorAdjacencies()
122 …PetscCall(DMPlexGetAdjacency_Internal(dm, q, useCone, useClosure, PETSC_TRUE, &numAdjQ, &tmpAdjQ)); in DMPlexComputeAnchorAdjacencies()
213 for (PetscInt q = 0; q < numAdj; q++) { in AdjancencyContainsLeafRootPair() local
214 const PetscInt padj = tmpAdj[q]; in AdjancencyContainsLeafRootPair()
224 for (PetscInt q = 0; q < numAdj; q++) { in AdjancencyContainsLeafRootPair() local
[all …]
H A Dplexindices.c51 PetscInt *points = NULL, numPoints, p, q, dof, cldof, cloff; in DMPlexCreateClosureIndex() local
56 for (p = 0, q = 0; p < numPoints * 2; p += 2) { in DMPlexCreateClosureIndex()
60 clPoints[cloff + q * 2] = points[p]; in DMPlexCreateClosureIndex()
61 clPoints[cloff + q * 2 + 1] = points[p + 1]; in DMPlexCreateClosureIndex()
62 ++q; in DMPlexCreateClosureIndex()
67 …cCheck(q * 2 == cldof, PetscObjectComm((PetscObject)dm), PETSC_ERR_PLIB, "Invalid size for closure… in DMPlexCreateClosureIndex()
/petsc/src/dm/field/impls/shell/
H A Ddmfieldshell.c55 PetscInt q; in DMFieldShellEvaluateFEDefault() local
58 for (q = 0; q < Nc; q++) { in DMFieldShellEvaluateFEDefault()
63 for (i = 0; i < dimC; i++) d[i] = sD[(p * Nc + q) * dimC + i]; in DMFieldShellEvaluateFEDefault()
64 for (i = 0; i < dimC; i++) sD[(p * Nc + q) * dimC + i] = 0.; in DMFieldShellEvaluateFEDefault()
66 …for (j = 0; j < dimC; j++) sD[(p * Nc + q) * dimC + i] += geom->J[(p * dimC + j) * dimC + i] * d[j… in DMFieldShellEvaluateFEDefault()
72 PetscInt q; in DMFieldShellEvaluateFEDefault() local
75 for (q = 0; q < Nc; q++) { in DMFieldShellEvaluateFEDefault()
80 for (i = 0; i < dimC; i++) d[i] = rD[(p * Nc + q) * dimC + i]; in DMFieldShellEvaluateFEDefault()
81 for (i = 0; i < dimC; i++) rD[(p * Nc + q) * dimC + i] = 0.; in DMFieldShellEvaluateFEDefault()
83 …for (j = 0; j < dimC; j++) rD[(p * Nc + q) * dimC + i] += geom->J[(p * dimC + j) * dimC + i] * d[j… in DMFieldShellEvaluateFEDefault()
[all …]
/petsc/src/mat/tests/
H A Dex177.c12 PetscInt m, n, M, N, p = 1, q = 1, i, j; in main() local
31 PetscCall(PetscOptionsGetInt(NULL, NULL, "-q", &q, NULL)); in main()
32 PetscCall(PetscMalloc2(p * q, &S, p * q, &T)); in main()
33 for (i = 0; i < p * q; i++) S[i] = 0; in main()
36 for (j = 0; j < q; j++) { in main()
38 S[i + p * j] = ((PetscReal)((i + 1) * (j + 1))) / ((PetscReal)(p + q)); in main()
39 T[i + p * j] = ((PetscReal)((p - i) + j)) / ((PetscReal)(p * q)); in main()
46 PetscCall(MatCreateKAIJ(A, p, q, S, T, &TA)); in main()
68 PetscCall(MatCreateKAIJ(A, p, q, NULL, T, &TA)); in main()
90 PetscCall(MatCreateKAIJ(A, p, q, S, NULL, &TA)); in main()
[all …]
H A Dex248.c9 PetscInt n = 10, m = 20, p = 7, q = 17; in main() local
15 PetscCall(MatCreateDense(PETSC_COMM_SELF, p, q, p, q, NULL, &Bd)); in main()
27 PetscCall(MatCreateKAIJ(A, p, q, NULL, Bv, &K)); in main()
36 PetscCall(MatKAIJSetT(K, p, q, Bv)); in main()
/petsc/src/mat/impls/kaij/
H A Dkaij.c92 if (n) *n = b->q; in MatKAIJGetS()
123 if (n) *n = b->q; in MatKAIJGetSRead()
203 if (n) *n = b->q; in MatKAIJGetT()
234 if (n) *n = b->q; in MatKAIJGetTRead()
346 PetscErrorCode MatKAIJSetS(Mat A, PetscInt p, PetscInt q, const PetscScalar S[]) in MatKAIJSetS() argument
353 PetscCall(PetscMalloc1(p * q, &a->S)); in MatKAIJSetS()
354 PetscCall(PetscArraycpy(a->S, S, p * q)); in MatKAIJSetS()
358 a->q = q; in MatKAIJSetS()
383 if (a->p != a->q) { in MatKAIJGetScaledIdentity()
423 PetscErrorCode MatKAIJSetT(Mat A, PetscInt p, PetscInt q, const PetscScalar T[]) in MatKAIJSetT() argument
[all …]
/petsc/src/dm/dt/fe/impls/basic/
H A Dfebasic.c60 PetscInt Nc, Nq, q, k, c; in PetscFESetUp_Basic() local
70 for (q = 0; q < Nq; ++q) { in PetscFESetUp_Basic()
71 …for (c = 0; c < Nc; ++c) fem->invV[j * pdim + k] += Bf[(q * pdim + k) * Nc + c] * weights[q * Nc +… in PetscFESetUp_Basic()
176 PetscInt qNc, Nq, q; in PetscFEIntegrate_Basic() local
221 for (q = 0; q < Nq; ++q) { in PetscFEIntegrate_Basic()
226 …CoordinatesRefToReal(dE, dim, fegeom.xi, &cgeom->v[e * Np * dE], fegeom.J, &quadPoints[q * dim], x… in PetscFEIntegrate_Basic()
228 fegeom.v = &cgeom->v[(e * Np + q) * dE]; in PetscFEIntegrate_Basic()
229 fegeom.J = &cgeom->J[(e * Np + q) * dE * dE]; in PetscFEIntegrate_Basic()
230 fegeom.invJ = &cgeom->invJ[(e * Np + q) * dE * dE]; in PetscFEIntegrate_Basic()
231 fegeom.detJ = &cgeom->detJ[e * Np + q]; in PetscFEIntegrate_Basic()
[all …]
/petsc/src/ksp/ksp/tutorials/
H A Dex70.c526 PetscInt dim, nel, npe, q, k, d, ncurr, Nfc; in DMSwarmPICInsertPointsCellwise() local
542 for (q = 0; q < npoints; q++) { in DMSwarmPICInsertPointsCellwise()
543 PetscCall(PetscMalloc1(npe, &basis[q])); in DMSwarmPICInsertPointsCellwise()
547 basis[q][0] = 0.5 * (1.0 - xi[dim * q + 0]); in DMSwarmPICInsertPointsCellwise()
548 basis[q][1] = 0.5 * (1.0 + xi[dim * q + 0]); in DMSwarmPICInsertPointsCellwise()
551 basis[q][0] = 0.25 * (1.0 - xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]); in DMSwarmPICInsertPointsCellwise()
552 basis[q][1] = 0.25 * (1.0 + xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]); in DMSwarmPICInsertPointsCellwise()
553 basis[q][2] = 0.25 * (1.0 + xi[dim * q + 0]) * (1.0 + xi[dim * q + 1]); in DMSwarmPICInsertPointsCellwise()
554 basis[q][3] = 0.25 * (1.0 - xi[dim * q + 0]) * (1.0 + xi[dim * q + 1]); in DMSwarmPICInsertPointsCellwise()
558 … basis[q][0] = 0.125 * (1.0 - xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]) * (1.0 - xi[dim * q + 2]); in DMSwarmPICInsertPointsCellwise()
[all …]
/petsc/src/vec/vec/tests/
H A Dex35.c7 PetscInt cnt, i, j, k, l, m = 2, n = 3, p = 4, q = 5; in main() local
13 PetscCall(VecCreateSeq(PETSC_COMM_WORLD, m * n * p * q, &x)); in main()
14 PetscCall(VecGetArray4d(x, m, n, p, q, 0, 0, 0, 0, &xx)); in main()
19 for (l = 0; l < q; l++) xx[i][j][k][l] = cnt++; in main()
21 PetscCall(VecRestoreArray4d(x, m, n, p, q, 0, 0, 0, 0, &xx)); in main()
/petsc/share/petsc/datafiles/meshes/
H A Dgmsh-3d-binary-32.msh44 ���n���o���p���q���r���s���t���u���v���w���x���y���z���������?@������?�������?������?������?����…
60 …q���a���b���e���r���b���a���f���s���b���d���e���t���d���b���h���u���c���b���f���v���b���c���g���w�…
61 ��<���K���E���}��� ��q���x���K���}��� ��-���X���]�������
H A Dmesh-3d-box-innersphere_bin.msh7 ���\3&�q<H]���w3��������?���\3&�q<H]���w3�������п������������������������?�����������������…
17 q����rJe�ĿE�����2$%&����}w�ҹ�}m�wLq�?F�����k*�����Z���q���ߗʿG����"E���?=��Yи?/Đ��@˿H���#���,�…
56 …q���������������g������ ���q���p���������������h������ ���m���q���������������i����…
H A Dgmsh-3d-binary-64.msh39 �������n�������o�������p�������q�������r�������s�������t�������u�������v�������w�������x�������y���…
54 …q�������b�������e�������f�������r�������e�������b�������g�������s�������d�������b�������f�������t�…
56 …q�������o�������|�������3������O�������J�������n�������~�������4������V�������a�������i���������…
/petsc/src/dm/impls/plex/tests/
H A Ddmplexcomputecellgeometryfem.F9034 PetscQuadrature :: q
40 PetscCallA(PetscDTSimplexQuadrature(numDim, tria_qorder, PETSCDTSIMPLEXQUAD_DEFAULT, q, ierr))
42 PetscCallA(DMPlexComputeCellGeometryFEM(dm, cStart, q, v, J, invJ, detJ, ierr))
45 PetscCallA(PetscQuadratureDestroy(q, ierr))
60 PetscQuadrature :: q
74 PetscCallA(PetscDTGaussTensorQuadrature(numDim, nc, npoints, a, b, q, ierr))
76 PetscCallA(DMPlexComputeCellGeometryFEM(dm, cStart, q, v, J, invJ, detJ, ierr))
81 PetscCallA(PetscQuadratureDestroy(q, ierr))
H A Dex75.c8 PetscInt dim = 1, d, cStart, cEnd, c, q, degree = 2, coordSize, offset; in main() local
52 for (q = 0; q < degree + 1; ++q) { in main()
54 …; d < dim; ++d) coords[offset + dim * q + d] = 100. + (PetscScalar)c + (1.0 / (PetscScalar)degree)… in main()
/petsc/src/dm/field/impls/ds/
H A Ddmfieldds.c576 PetscInt i, j, k, l, m, q, r; in DMFieldEvaluateFV_DS() local
600 for (q = 0; q < dim; q++) { in DMFieldEvaluateFV_DS()
601 … for (r = 0; r < dim; r++) pH[l][m] += invJ[q * dimC + l] * invJ[r * dimC + m] * spH[q * dim + r]; in DMFieldEvaluateFV_DS()
634 for (q = 0; q < dim; q++) { in DMFieldEvaluateFV_DS()
635 … for (r = 0; r < dim; r++) pH[l][m] += invJ[q * dimC + l] * invJ[r * dimC + m] * rpH[q * dim + r]; in DMFieldEvaluateFV_DS()
829 PetscInt Np = geom->numPoints, q, dE = geom->dimEmbed, d; in DMFieldComputeFaceData_DS() local
831 for (q = 0; q < Np; ++q) in DMFieldComputeFaceData_DS()
832 for (d = 0; d < dE; ++d) geom->n[(p * Np + q) * dE + d] = -geom->n[(p * Np + q) * dE + d]; in DMFieldComputeFaceData_DS()
886 PetscInt numSupp, s, q; in DMFieldComputeFaceData_DS() local
892 for (q = 0; q < Nq * dE * dE; q++) { in DMFieldComputeFaceData_DS()
[all …]
/petsc/src/dm/dt/tests/
H A Dex15.c31 …for (PetscInt q = 0; q < num_points; q++) I_quad += weights[q] * eval[i * num_points + q] * eval[j… in testQuadrature() local
/petsc/src/binding/petsc4py/demo/legacy/kspsolve/
H A Dpetsc-cg.py10 q = b.duplicate()
21 A.mult(d, q)
22 alpha = delta / d.dot(q)
24 r.axpy(-alpha, q)
/petsc/src/dm/dt/interface/ftn-custom/
H A Dzdtf90.c14 PETSC_EXTERN void petscquadraturegetdata_(PetscQuadrature *q, PetscInt *dim, PetscInt *Nc, PetscInt… in petscquadraturegetdata_() argument
18 *ierr = PetscQuadratureGetData(*q, dim, Nc, npoints, &points, &weights); in petscquadraturegetdata_()
25 PETSC_EXTERN void petscquadraturerestoredata_(PetscQuadrature *q, PetscInt *dim, PetscInt *Nc, Pets… in petscquadraturerestoredata_() argument
32 PETSC_EXTERN void petscquadraturesetdata_(PetscQuadrature *q, PetscInt *dim, PetscInt *Nc, PetscInt… in petscquadraturesetdata_() argument
40 *ierr = PetscQuadratureSetData(*q, *dim, *Nc, *npoints, points, weights); in petscquadraturesetdata_()
/petsc/src/dm/dt/interface/
H A Ddt.c60 PetscErrorCode PetscQuadratureCreate(MPI_Comm comm, PetscQuadrature *q) in PetscQuadratureCreate() argument
63 PetscAssertPointer(q, 2); in PetscQuadratureCreate()
65 …PetscCall(PetscHeaderCreate(*q, PETSCQUADRATURE_CLASSID, "PetscQuadrature", "Quadrature", "DT", co… in PetscQuadratureCreate()
66 (*q)->ct = DM_POLYTOPE_UNKNOWN; in PetscQuadratureCreate()
67 (*q)->dim = -1; in PetscQuadratureCreate()
68 (*q)->Nc = 1; in PetscQuadratureCreate()
69 (*q)->order = -1; in PetscQuadratureCreate()
70 (*q)->numPoints = 0; in PetscQuadratureCreate()
71 (*q)->points = NULL; in PetscQuadratureCreate()
72 (*q)->weights = NULL; in PetscQuadratureCreate()
[all …]
/petsc/src/binding/petsc4py/test/
H A Dtest_ksp_py.py64 q = d.duplicate()
65 self.work += [d, q]
70 r, z, d, q = self.work
78 A.mult(d, q)
79 alpha = delta / d.dot(q)
81 r.axpy(-alpha, q)
/petsc/src/dm/impls/patch/
H A Dpatch.c50 …P, rM, rN, rP, halo = 1, sxb, syb, szb, sxr, syr, szr, exr, eyr, ezr, mxb, myb, mzb, i, j, k, l, q; in DMPatchZoom() local
114 q = 0; in DMPatchZoom()
126 localPoints[q] = lp; in DMPatchZoom()
127 PetscCall(PetscFindInt(indices[q], size + 1, ranges, &ir)); in DMPatchZoom()
130 remotePoints[q].rank = r < 0 ? -(r + 1) - 1 : r; in DMPatchZoom()
131 remotePoints[q].index = indices[q] - ranges[remotePoints[q].rank]; in DMPatchZoom()
132 ++q; in DMPatchZoom()
141 …PetscCall(PetscSFSetGraph(*sfzr, dof * M * N * P, q, localPoints, PETSC_COPY_VALUES, remotePoints,… in DMPatchZoom()
155 q = 0; in DMPatchZoom()
158 for (i = sxb; i < sxb + mxb; ++i, ++q) { in DMPatchZoom()
[all …]

1234567