Feature/custom bsg gridoperator

Extract the C++ implementation of the optimized blockstructured gridoperator from my personal pdelab fork. I think the implementation is more useful in this repository. Now the optimized GO can be used without any changes to pdelab.