How to utilize reconstruction operator for relaxed HDiv HDG?

I am currently trying to reconstruct solution u for the relaxed HDiv space, as given in https://epubs.siam.org/doi/pdf/10.1137/17M1138078. Can anyone suggest some hints on how to do that or is there any related Example/Operator()? Any input will help.

Wh