Is there a function which returns one of the properties accessible through
cubeprop() (e.g. electron density), but at any arbitrary 3D coordinate, given a wavefunction? I assume the information itself could be easily computed, given that
cubeprop() can return this property if you choose the grid correctly, but this of course comes along with the entire cubefile. I’m thinking of something like this:
```E, wfn = psi4.energy(level_of_theory, return_wfn=True) result = pointprop(wfn, [x, y, z], tasks="density")```
oeprop() seems to provide some of this functionality, but not for instance electron density.
Thanks and best regards,