Add a blocksparse multiplication implementation

This is a first step, it may still need more work.
I don't understand "bptr" mean and do. Make a small docstring for that fct.