added pressure derivatives to compute_pressure()