"Zhang, Junchao" <[email protected]> writes: > On Tue, Jan 22, 2019 at 4:08 PM Jed Brown > <[email protected]<mailto:[email protected]>> wrote: >> It is not supported at this time. What does your use case look like? >> Do roots have degree greater than 1? > > Yes. Imagine a vecscatter x[0]->y[0], x[1]->y[0], x[1]->y[1],x[2]->y[1]. I > build an SF for SCATTER_FORWARD. Now I wan to do SCATTER_REVERSE with > ADD_VALUES. > To solve this problem without creating another SF and without breaking > current SF API, I propose to add PetscSFBcastAndOp(sf, unit, rootdata, > leafdata, op)
Is this pattern needed for some algorithm or application or is this about implementing the full VecScatter interface (modulo ill-defined operations) in terms of SF? If the former, I'd like to understand it. If the latter, I'm okay with you extending the interface as proposed.
