The signature of the computeRho
method you want to call is:
T olb::ConstCell<T,DESCRIPTOR>::computeRho() const
i.e. it returns the value of rho as a value, not by writing it to a given pointer.
This means that you need to write e.g.:
T rho = cell.neighbor({n0,n1}).computeRho();