PetscErrorCode trapz_4nynznt(const Vec &q, PetscScalar &I, SPE &data)
trapezoidal rule on Vec with 4,ny,nz,nt values
Definition: trapz.cpp:110
int trapz(const Vec &q, PetscScalar &I, SPE &data)
trapezoidal rule on Vec with ny,nz values
Definition: trapz.cpp:10