tried to directly compute deformation gradient in CPDI2, no evident improvement