60 …const PetscInt MBS = M * BS;             /* MBS=bs. We turn MBS into a compile-time const when EQ=…  in d_Pack()  local
63     PetscInt t = (opt ? MapTidToIndex(opt, tid) : (idx ? idx[tid] : start + tid)) * MBS;  in d_Pack()
64     PetscInt s = tid * MBS;  in d_Pack()
65     for (PetscInt i = 0; i < MBS; i++) buf[s + i] = data[t + i];  in d_Pack()
72   const PetscInt M = (EQ) ? 1 : bs / BS, MBS = M * BS;  in d_UnpackAndOp()  local
76     PetscInt t = (opt ? MapTidToIndex(opt, tid) : (idx ? idx[tid] : start + tid)) * MBS;  in d_UnpackAndOp()
77     PetscInt s = tid * MBS;  in d_UnpackAndOp()
78     for (PetscInt i = 0; i < MBS; i++) op(data[t + i], buf[s + i]);  in d_UnpackAndOp()
85   const PetscInt M = (EQ) ? 1 : bs / BS, MBS = M * BS;  in d_FetchAndOp()  local
89 …Int r = (rootopt ? MapTidToIndex(rootopt, tid) : (rootidx ? rootidx[tid] : rootstart + tid)) * MBS;  in d_FetchAndOp()
90     PetscInt l = tid * MBS;  in d_FetchAndOp()
91     for (PetscInt i = 0; i < MBS; i++) leafbuf[l + i] = op(rootdata[r + i], leafbuf[l + i]);  in d_FetchAndOp()
98   const PetscInt M = (EQ) ? 1 : bs / BS, MBS = M * BS;  in d_ScatterAndOp()  local
124     s *= MBS;  in d_ScatterAndOp()
125     t *= MBS;  in d_ScatterAndOp()
126     for (PetscInt i = 0; i < MBS; i++) op(dst[t + i], src[s + i]);  in d_ScatterAndOp()
133   const PetscInt M = (EQ) ? 1 : bs / BS, MBS = M * BS;  in d_FetchAndOpLocal()  local
137 …Int r = (rootopt ? MapTidToIndex(rootopt, tid) : (rootidx ? rootidx[tid] : rootstart + tid)) * MBS;  in d_FetchAndOpLocal()
138 …Int l = (leafopt ? MapTidToIndex(leafopt, tid) : (leafidx ? leafidx[tid] : leafstart + tid)) * MBS;  in d_FetchAndOpLocal()
139     for (PetscInt i = 0; i < MBS; i++) leafupdate[l + i] = op(rootdata[r + i], leafdata[l + i]);  in d_FetchAndOpLocal()