up #460:
 - introduce domain and image maps in operator in order to handle image and
 domain space which don't have the same size
 - change interface of OperatorPCD in order to retrieve also the velocity space
 and extract the proper information