Lines Matching refs:fieldA
139 for (PetscInt fieldA = 0; fieldA < ctx->num_species; fieldA++) { in LandauFormJacobian_Internal() local
140 …Eq_m[fieldA] = ctx->Ez * ctx->t_0 * ctx->charges[fieldA] / (ctx->v_0 * ctx->masses[fieldA]); /* no… in LandauFormJacobian_Internal()
141 …if (dim == 2) Eq_m[fieldA] *= 2 * PETSC_PI; /* ad… in LandauFormJacobian_Internal()
407 … for (PetscInt fieldA = 0, f_off = ctx->species_offset[grid]; fieldA < loc_Nf; ++fieldA) { in LandauFormJacobian_Internal() local
409 gg2[fieldA][d2] = gg2_temp[d2] * nu_alpha[fieldA + f_off]; in LandauFormJacobian_Internal()
410 … (d3 = 0; d3 < LANDAU_DIM; d3++) gg3[fieldA][d2][d3] = -gg3_temp[d2][d3] * nu_alpha[fieldA + f_off… in LandauFormJacobian_Internal()
414 …for (PetscInt fieldA = 0, f_off = ctx->species_offset[grid]; fieldA < loc_Nf; ++fieldA) gg2[fieldA… in LandauFormJacobian_Internal() local
416 for (PetscInt fieldA = 0; fieldA < loc_Nf; ++fieldA) { in LandauFormJacobian_Internal() local
418 g2[fieldA][d] = 0.0; in LandauFormJacobian_Internal()
420 g2[fieldA][d] += invJj[d * dim + d2] * gg2[fieldA][d2]; in LandauFormJacobian_Internal()
421 g3[fieldA][d][d2] = 0.0; in LandauFormJacobian_Internal()
423 …for (dp = 0; dp < dim; ++dp) g3[fieldA][d][d2] += invJj[d * dim + d3] * gg3[fieldA][d3][dp] * invJ… in LandauFormJacobian_Internal()
425 g3[fieldA][d][d2] *= wj; in LandauFormJacobian_Internal()
427 g2[fieldA][d] *= wj; in LandauFormJacobian_Internal()
433 for (PetscInt fieldA = 0; fieldA < loc_Nf; ++fieldA) { in LandauFormJacobian_Internal() local
435 g0[fieldA] = wj * shift * 2. * PETSC_PI; // move this to below and remove g0 in LandauFormJacobian_Internal()
437 g0[fieldA] = wj * shift; // move this to below and remove g0 in LandauFormJacobian_Internal()
443 PetscInt fieldA, d, f, d2, g; in LandauFormJacobian_Internal() local
446 for (fieldA = 0; fieldA < loc_Nf; fieldA++) { in LandauFormJacobian_Internal()
448 const PetscInt i = fieldA * Nb + f; /* Element matrix row */ in LandauFormJacobian_Internal()
450 const PetscInt j = fieldA * Nb + g; /* Element matrix column */ in LandauFormJacobian_Internal()
454 elemMat[fOff] += DIq[f * dim + d] * g2[fieldA][d] * BJq[g]; in LandauFormJacobian_Internal()
455 …for (d2 = 0; d2 < dim; ++d2) elemMat[fOff] += DIq[f * dim + d] * g3[fieldA][d][d2] * DIq[g * dim +… in LandauFormJacobian_Internal()
458 elemMat[fOff] += BJq[f] * g0[fieldA] * BJq[g]; in LandauFormJacobian_Internal()
482 …PetscInt fieldA, q, f, g, d, nr, nc, rows0[LANDAU_MAX_Q_FACE] = {0}, cols0[LANDAU_MAX_Q_FACE] =… in LandauFormJacobian_Internal() local
486 for (fieldA = 0; fieldA < loc_Nf; fieldA++) { in LandauFormJacobian_Internal()
487 LandauIdx *const Idxs = &maps[grid].gIdx[loc_elem][fieldA][0]; in LandauFormJacobian_Internal()
517 const PetscInt i = fieldA * Nb + f; /* Element matrix row */ in LandauFormJacobian_Internal()
518 const PetscInt j = fieldA * Nb + g; /* Element matrix column */ in LandauFormJacobian_Internal()
522 …_elem_offsets[elem_offset[num_grids]] + coo_elem_offsets[glb_elem_idx] + fieldA * fullNb2 + fullNb… in LandauFormJacobian_Internal()
1518 for (PetscInt fieldA = 0; fieldA < Nf[grid]; fieldA++) { in CreateStaticData() local
1524 elMat[(fieldA * Nb + q) * totDim + fieldA * Nb + q] = 1; in CreateStaticData()
1528 … for (PetscInt q = 0; q < numindices; ++q) maps[grid].gIdx[eidx][fieldA][q] = indices[q]; in CreateStaticData()
1536 maps[grid].gIdx[eidx][fieldA][q] = plex_batch[indices[f]]; in CreateStaticData()
1538 maps[grid].gIdx[eidx][fieldA][q] = indices[f]; in CreateStaticData()
1545 …maps[grid].gIdx[eidx][fieldA][q] = -maps[grid].num_reduced - 1; // store (-)index: id = -(idx+1): … in CreateStaticData()
1575 …I = %22.16e (LANDAU_MAX_Q_FACE=%d, #face=%" PetscInt_FMT ")\n", eidx, q, fieldA, (double)sum, LAND… in CreateStaticData()
1595 …if (fieldA == 0) { // cache full Nb for this element, on this … in CreateStaticData()
1656 for (PetscInt fieldA = 0; fieldA < Nf[grid]; fieldA++) { in CreateStaticData() local
1657 const LandauIdx *const Idxs = &maps[grid].gIdx[ej][fieldA][0]; in CreateStaticData()
1673 … = b_id * coo_elem_offsets[ncellsTot] + coo_elem_offsets[glb_elem_idx] + fieldA * fullNb2 + fullNb… in CreateStaticData()