- weight_rt_1 *= tav_viol_Rtav7*pow(dd->rm3tav[pair], seven_three)+
- instant_viol_Rtav7*pow(dd->rt[pair], -7);
+ weight_rt_1 *= tav_viol_Rtav7*std::pow(dd->rm3tav[pair], seven_three)+
+ instant_viol_Rtav7*std::pow(dd->rt[pair], static_cast<real>(-7));