How to compute the JVP for multiple-valued outputs?

[RohanGautam]

I have a c++ function with a std::vec as input and a (larger) std::vec for it's output. I couldn't find a clear way to compute the JVP of the output in this multi-input, multi-output case. Seems like it should be possible (the JVP is computed in the Enzyme-JAX interop library example), but I couldn't find any way to do this from the documentation.

A sample fuction signature:

std::vector<double> my_function(std::vector<double> input)
{//code returning another std::vector}

Any pointers to relavant resources would be appreciated!

[ZuseZ4]

I am not aware of support for it, so you'll likely need to write a wrapper around it which also accepts the output vector by reference.