Experiment with calculations
Create granulate function with output arguments
Use optimized solution for all operations in region_divide_granularity.hpp
Add some asserts for safety
Removing un-used file
Updating comments/doxygen