1. update CPDI2, the gradient can now be with respect to initial configuration
2. implemented comparison of particle based integration and quadrature points based integration,
   particle-based is 2x faster