expose vectorial RZ from c code

would be nice to get directly at the result of `VectorialPSF::integral_` ... there are slight differences between the pure python result and the c++ result, and I suspect it is only related to interpolation at the end (differences decrease as the dxy pixel size gets smaller