implement CPDI2 3D with numerical integration, not done yet