Commit e6d94516 authored by Alain O' Miniussi's avatar Alain O' Miniussi
Browse files

use c++17 pow

parent 30790af3
......@@ -161,9 +161,9 @@ namespace fargOCA {
arr3d<real>::HostMirror heta2 = kk::create_mirror_view(eta2);
kk::parallel_for("init_stell_heating", grid.fullRange(),
LBD(int i, int h, int j) {
real denom = 1+16*dt*diffusion.PlankOpacity(i,h,j)*SIGMARAD*pow(temperature(i,h,j),3.0)/specificHeat;
real denom = 1+16*dt*diffusion.PlankOpacity(i,h,j)*SIGMARAD*pow(temperature(i,h,j),3)/specificHeat;
heta1(i,h,j) = (temperature(i,h,j)
+ dt*(12*diffusion.PlankOpacity(i,h,j)/specificHeat*SIGMARAD*pow(temperature(i,h,j),4.0)
+ dt*(12*diffusion.PlankOpacity(i,h,j)/specificHeat*SIGMARAD*pow(temperature(i,h,j),4)
+ (diffusion.stellarRadiation(i,h,j)+qplus(i,h,j)+planetRadiation(i,h,j))/(density(i,h,j)*specificHeat))) / denom;
heta2(i,h,j) = dt*C0*diffusion.PlankOpacity(i,h,j)/(specificHeat*denom);
});
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment