Searched refs:qj (Results 1 – 4 of 4) sorted by relevance
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex54f.F90 | 28 PetscInt qj, qi, ne, M, Istart, Iend, geq
112 qj = geq/(ne + 1); qi = mod(geq, (ne + 1))
113 x = h*qi - 1.0; y = h*qj - 1.0 ! lower left corner (-1,-1)
114 if (qi < ne .and. qj < ne) then
143 if (qj > 0) then
160 if (qj > 0) then ! set rhs
|
| /petsc/src/snes/tutorials/ |
| H A D | ex16.c | 282 void DeformationGradient(Field *ex, PetscInt qi, PetscInt qj, PetscInt qk, PetscScalar *invJ, Petsc… in DeformationGradient() argument
294 PetscInt bidx = NEB * idx + qi + NQ * qj + NQ * NQ * qk; in DeformationGradient()
311 void DeformationGradientJacobian(PetscInt qi, PetscInt qj, PetscInt qk, PetscInt ii, PetscInt jj, P… in DeformationGradientJacobian() argument
316 PetscInt bidx = NEB * idx + qi + NQ * qj + NQ * NQ * qk; in DeformationGradientJacobian()
437 void QuadraturePointGeometricJacobian(CoordField *ec, PetscInt qi, PetscInt qj, PetscInt qk, PetscS… in QuadraturePointGeometricJacobian() argument
446 PetscInt bidx = NEB * idx + qi + NQ * qj + NQ * NQ * qk; in QuadraturePointGeometricJacobian()
468 PetscInt i, j, k, l, ii, jj, kk, ll, qi, qj, qk, m; in FormElementJacobian() local
486 for (qj = 0; qj < NQ; qj++) { in FormElementJacobian()
488 QuadraturePointGeometricJacobian(ec, qi, qj, qk, J); in FormElementJacobian()
490 scl = vol * wts[qi] * wts[qj] * wts[qk]; in FormElementJacobian()
[all …]
|
| /petsc/src/ts/utils/dmplexlandau/kokkos/ |
| H A D | landau.kokkos.cxx | 699 for (int qj = 0; qj < Nq; qj++) { // look at others integration points in LandauKokkosJacobian() local
700 const PetscReal *BJq = &d_BB[qj * Nb], *DIq = &d_DD[qj * Nb * dim]; in LandauKokkosJacobian()
702 t += DIq[f * dim + d] * g2(d, fieldA, qj) * BJq[g]; in LandauKokkosJacobian()
703 …for (int d2 = 0; d2 < dim; ++d2) t += DIq[f * dim + d] * g3(d, d2, fieldA, qj) * DIq[g * dim + d2]; in LandauKokkosJacobian()
754 for (int qj = 0; qj < Nq; qj++) { // look at others integration points in LandauKokkosJacobian() local
755 const PetscReal *BJq = &d_BB[qj * Nb]; in LandauKokkosJacobian()
756 const PetscInt jpidx_glb = jpidx_0 + qj; in LandauKokkosJacobian()
|
| /petsc/src/ts/utils/dmplexlandau/ |
| H A D | plexland.c | 334 for (PetscInt qj = 0; qj < Nq; ++qj) { in LandauFormJacobian_Internal() local
335 const PetscInt jpidx_glb = ip_offset[grid] + qj + loc_elem * Nq; in LandauFormJacobian_Internal()
339 const PetscReal *const invJj = &invJe[qj * dim * dim]; in LandauFormJacobian_Internal()
444 const PetscReal *BJq = &BB[qj * Nb], *DIq = &DD[qj * Nb * dim]; in LandauFormJacobian_Internal()
913 PetscInt Nq, Nb, *Nb2, cStart, cEnd, c, dim, qj, k; in adaptToleranceFEM() local
944 for (qj = 0; qj < Nq; ++qj) { in adaptToleranceFEM()
945 …tt = PetscSqr(v0[dim * qj + 0]) + PetscSqr(v0[dim * qj + 1]) + PetscSqr((dim == 3) ? v0[dim * qj +… in adaptToleranceFEM()
951 …ell %" PetscInt_FMT ", qp %" PetscInt_FMT "/%" PetscInt_FMT "\n", grid, (double)r, c, qj + 1, Nq)); in adaptToleranceFEM()
957 …MT ", qp %" PetscInt_FMT "/%" PetscInt_FMT ", d=%e\n", grid, (double)r, c, qj + 1, Nq, (double)(r … in adaptToleranceFEM()
1760 for (PetscInt qj = 0; qj < Nq; qj++, outer_ipidx++) { in CreateStaticData() local
[all …]
|