Use correct primal values in the sensitivity calls

[Jamie Bramwell]

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...); 
     }...}; 

.

[Brandon Talamini]

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