Use correct primal values in the sensitivity calls

[jamiebramwell]

As written, the sensitivity calculation calls in the HeatTransfer and SolidMechanics modules will not use the correct primal field values (e.g. temperature) during the backward adjoint pass. See

serac/src/serac/physics/heat_transfer.hpp

Lines 921 to 927 in a2d7924

	std::array<std::function<decltype((*residual_)(DifferentiateWRT<0>{}, temperature_, temperature_rate_,
	shape_displacement_, *parameters_[parameter_indices].state...))()>,
	sizeof...(parameter_indices)>
	d_residual_d_ = {[&]() {
	return (*residual_)(DifferentiateWRT<NUM_STATE_VARS + parameter_indices>{}, temperature_, temperature_rate_,
	shape_displacement_, *parameters_[parameter_indices].state...);
	}...};

.

[btalamini]

We should check that the correct time is used as well. Ideally, we should write tests that fail before patching the bugs.