325   PetscInt i, j, Nk, Mk;  in PetscDTAltVPullback()  local
331     PetscCall(PetscDTBinomialInt(M, PetscAbsInt(k), &Mk));  in PetscDTAltVPullback()
339         for (j = 0; j < Mk; j++) sum += L[j * Nk + i] * w[j];  in PetscDTAltVPullback()
348         for (j = 0; j < Mk; j++) sum += L[(Mk - 1 - j) * Nk + (Nk - 1 - i)] * w[j] * mult[j];  in PetscDTAltVPullback()
360 …for (j = 0; j < Mk; j++) sum += (L[pairs[j][0] * N + pairs[i][0]] * L[pairs[j][1] * N + pairs[i][1…  in PetscDTAltVPullback()
375 …for (j = 0; j < Mk; j++) sum += (L[pairs[offj + j][0] * N + pairs[offi + i][0]] * L[pairs[offj + j…  in PetscDTAltVPullback()
392     PetscCall(PetscDTBinomialInt(M, PetscAbsInt(k), &Mk));  in PetscDTAltVPullback()
398       PetscCall(PetscMalloc1(Mk, &walloc));  in PetscDTAltVPullback()
406     for (i = 0; i < Mk; i++) {  in PetscDTAltVPullback()
454   PetscInt   Nk, Mk, Nf, i, j, l, p;  in PetscDTAltVPullbackMatrix()  local
466     PetscCall(PetscDTBinomialInt(M, PetscAbsInt(k), &Mk));  in PetscDTAltVPullbackMatrix()
472         for (j = 0; j < Mk; j++) Lstar[i * Mk + j] = L[j * Nk + i];  in PetscDTAltVPullbackMatrix()
476 …for (j = 0; j < Mk; j++) Lstar[i * Mk + j] = L[(Mk - 1 - j) * Nk + (Nk - 1 - i)] * mult[i] * mult[…  in PetscDTAltVPullbackMatrix()
486 …for (j = 0; j < Mk; j++) Lstar[i * Mk + j] = L[pairs[j][0] * N + pairs[i][0]] * L[pairs[j][1] * N …  in PetscDTAltVPullbackMatrix()
498 …for (j = 0; j < Mk; j++) Lstar[i * Mk + j] = L[pairs[offj + j][0] * N + pairs[offi + i][0]] * L[pa…  in PetscDTAltVPullbackMatrix()
510     PetscCall(PetscDTBinomialInt(M, PetscAbsInt(k), &Mk));  in PetscDTAltVPullbackMatrix()
514     for (i = 0; i < Nk * Mk; i++) Lstar[i] = 0.;  in PetscDTAltVPullbackMatrix()
515     for (i = 0; i < Mk; i++) {  in PetscDTAltVPullbackMatrix()
520       iidx = negative ? Mk - 1 - i : i;  in PetscDTAltVPullbackMatrix()
536           Lstar[jidx * Mk + iidx] += prod;  in PetscDTAltVPullbackMatrix()