About rotational correlation function

Dear lammps users,

How can one use lammps to calculate rotational autocorrelation function? I know fix ave/correlate can do stress ACF using the green-kubo scripts. Does something like this exist for rotational ACF? I could not find anything similar in the examples. Thank you for your help.

Dear lammps users,

How can one use lammps to calculate rotational autocorrelation function?

by programming a compute/fix that computes such an auto-correlation function.

axel.