337           PetscInt       d, d2, dp, d3, IPf_idx;  in LandauFormJacobian_Internal()  local
345               for (d2 = 0; d2 < LANDAU_DIM; d2++) gg3_temp[d][d2] = 0;  in LandauFormJacobian_Internal()
383                   for (d2 = 0; d2 < 2; d2++) {  in LandauFormJacobian_Internal()
386                       gg2_temp[d2] += Uk[d2][d3] * temp1[d3];  in LandauFormJacobian_Internal()
388                       gg3_temp[d2][d3] += Ud[d2][d3] * temp2;  in LandauFormJacobian_Internal()
392                   for (d2 = 0; d2 < 3; ++d2) {  in LandauFormJacobian_Internal()
395                       gg2_temp[d2] += U[d2][d3] * temp1[d3];  in LandauFormJacobian_Internal()
397                       gg3_temp[d2][d3] += U[d2][d3] * temp2;  in LandauFormJacobian_Internal()
408               for (d2 = 0; d2 < LANDAU_DIM; d2++) {  in LandauFormJacobian_Internal()
409                 gg2[fieldA][d2] = gg2_temp[d2] * nu_alpha[fieldA + f_off];  in LandauFormJacobian_Internal()
410 …for (d3 = 0; d3 < LANDAU_DIM; d3++) gg3[fieldA][d2][d3] = -gg3_temp[d2][d3] * nu_alpha[fieldA + f_…  in LandauFormJacobian_Internal()
419                 for (d2 = 0; d2 < dim; ++d2) {  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()
443             PetscInt         fieldA, d, f, d2, g;  in LandauFormJacobian_Internal()  local
455 …for (d2 = 0; d2 < dim; ++d2) elemMat[fOff] += DIq[f * dim + d] * g3[fieldA][d][d2] * DIq[g * dim +…  in LandauFormJacobian_Internal()