Commit 356ea70a authored by Alain O' Miniussi's avatar Alain O' Miniussi
Browse files

in Bell Lin opacity, replace pow(rho,0.3333..) with more accurate cbrt(rho).

A regold was needed.
parent cc138d7a
1.885256000000e+02 5.090310207946e-02 3.075250735101e-02 5.936078698619e-05 1.538611058814e-06 2.140170351974e-07
1.885884000000e+02 5.090310333225e-02 3.075250735101e-02 5.935841492450e-05 1.533140170143e-06 2.138868398047e-07
1.886512000000e+02 5.090310458427e-02 3.075250735101e-02 5.935553416218e-05 1.527696902510e-06 2.137689891434e-07
1.887140000000e+02 5.090310583656e-02 3.075250735101e-02 5.935291431282e-05 1.522284121568e-06 2.136610733579e-07
1.887768000000e+02 5.090310708755e-02 3.075250735101e-02 5.935069446804e-05 1.516905226332e-06 2.135650924277e-07
1.888396000000e+02 5.090310833830e-02 3.075250735101e-02 5.934825375185e-05 1.511565199670e-06 2.134826990804e-07
1.889024000000e+02 5.090310958883e-02 3.075250735101e-02 5.934583132623e-05 1.506272111326e-06 2.134130700832e-07
1.889652000000e+02 5.090311083855e-02 3.075250735101e-02 5.934338505050e-05 1.501037945662e-06 2.133556533222e-07
1.890280000000e+02 5.090311208815e-02 3.075250735101e-02 5.934096371280e-05 1.495874304519e-06 2.133109098342e-07
1.890908000000e+02 5.090311333535e-02 3.075250735101e-02 5.933853050061e-05 1.490790309111e-06 2.132783026807e-07
1.891536000000e+02 5.090311457709e-02 3.075250735101e-02 5.933611843238e-05 1.485789238487e-06 2.132563302117e-07
1.892164000000e+02 5.090311581535e-02 3.075250735101e-02 5.933364100743e-05 1.480869173416e-06 2.132448006004e-07
1.892792000000e+02 5.090311705259e-02 3.075250735101e-02 5.933111982003e-05 1.476027149791e-06 2.132445884475e-07
1.886512000000e+02 5.090310458427e-02 3.075250735101e-02 5.935553416219e-05 1.527696902510e-06 2.137689891434e-07
1.887140000000e+02 5.090310583656e-02 3.075250735101e-02 5.935291431283e-05 1.522284121568e-06 2.136610733580e-07
1.887768000000e+02 5.090310708755e-02 3.075250735101e-02 5.935069446805e-05 1.516905226332e-06 2.135650924277e-07
1.888396000000e+02 5.090310833830e-02 3.075250735101e-02 5.934825375186e-05 1.511565199670e-06 2.134826990804e-07
1.889024000000e+02 5.090310958883e-02 3.075250735101e-02 5.934583132624e-05 1.506272111326e-06 2.134130700833e-07
1.889652000000e+02 5.090311083855e-02 3.075250735101e-02 5.934338505051e-05 1.501037945662e-06 2.133556533223e-07
1.890280000000e+02 5.090311208815e-02 3.075250735101e-02 5.934096371281e-05 1.495874304519e-06 2.133109098343e-07
1.890908000000e+02 5.090311333535e-02 3.075250735101e-02 5.933853050063e-05 1.490790309112e-06 2.132783026808e-07
1.891536000000e+02 5.090311457709e-02 3.075250735101e-02 5.933611843240e-05 1.485789238487e-06 2.132563302118e-07
1.892164000000e+02 5.090311581535e-02 3.075250735101e-02 5.933364100745e-05 1.480869173416e-06 2.132448006004e-07
1.892792000000e+02 5.090311705259e-02 3.075250735101e-02 5.933111982005e-05 1.476027149791e-06 2.132445884475e-07
301 9.94394087846767e-01 -1.07582432469537e-01 2.33115184434714e-04 1.07296827848511e-01 9.93906473518895e-01 1.00082045447214e-03 0.00000000000000e+00 1.88588400000006e+02 1.00036420875420e+00
302 9.99165616512461e-01 -4.49951231093621e-02 2.96062086413061e-04 4.46170151150747e-02 9.98669995821453e-01 1.00228412822682e-03 0.00000000000000e+00 1.88651200000006e+02 1.00046514733816e+00
303 9.99993973785480e-01 1.77685345644442e-02 3.58963328607858e-04 -1.82412938041172e-02 9.99519699685176e-01 9.99455574572164e-04 0.00000000000000e+00 1.88714000000006e+02 1.00048438632118e+00
304 9.96875605487044e-01 8.04625415326816e-02 4.21565004407284e-04 -8.10339413134328e-02 9.96451461025401e-01 9.92774957027987e-04 0.00000000000000e+00 1.88776800000006e+02 1.00042253623984e+00
304 9.96875605487044e-01 8.04625415326816e-02 4.21565004407284e-04 -8.10339413134327e-02 9.96451461025401e-01 9.92774957027987e-04 0.00000000000000e+00 1.88776800000006e+02 1.00042253623984e+00
305 9.89822302663278e-01 1.42841137424119e-01 4.83639035438303e-04 -1.43516964113233e-01 9.89477002542290e-01 9.82708705388290e-04 0.00000000000000e+00 1.88839600000007e+02 1.00027778101015e+00
306 9.78861137068390e-01 2.04659795146272e-01 5.44990148065365e-04 -2.05447638830394e-01 9.78623866437111e-01 9.69803483512203e-04 0.00000000000000e+00 1.88902400000007e+02 1.00004590495116e+00
307 9.64034355508449e-01 2.65676218228488e-01 6.05461964106819e-04 -2.66585574050126e-01 9.63935530837187e-01 9.54759529862352e-04 0.00000000000000e+00 1.88965200000007e+02 9.99719597953966e-01
308 9.45399196207236e-01 3.25651347400553e-01 6.64950248195395e-04 -3.26693902154045e-01 9.45471603356290e-01 9.38550616991676e-04 0.00000000000000e+00 1.89028000000007e+02 9.99287285012958e-01
307 9.64034355508449e-01 2.65676218228487e-01 6.05461964106819e-04 -2.66585574050126e-01 9.63935530837187e-01 9.54759529862352e-04 0.00000000000000e+00 1.88965200000007e+02 9.99719597953966e-01
308 9.45399196207236e-01 3.25651347400553e-01 6.64950248195395e-04 -3.26693902154045e-01 9.45471603356290e-01 9.38550616991675e-04 0.00000000000000e+00 1.89028000000007e+02 9.99287285012958e-01
309 9.23027624415959e-01 3.84350385312382e-01 7.23426643560922e-04 -3.85540673063419e-01 9.23308218140609e-01 9.22646769560822e-04 0.00000000000000e+00 1.89090800000008e+02 9.98730968466328e-01
310 8.97005968283864e-01 4.41543858774352e-01 7.80983167234996e-04 -4.42900669500748e-01 8.97538866893497e-01 9.09475799732708e-04 0.00000000000000e+00 1.89153600000008e+02 9.98022133760266e-01
311 8.67434413632281e-01 4.97008761340118e-01 8.37922761666330e-04 -4.98558173540879e-01 8.68276216890349e-01 9.03500216471810e-04 0.00000000000000e+00 1.89216400000008e+02 9.97113421966821e-01
312 8.34426257951122e-01 5.50529879282855e-01 8.94967049892886e-04 -5.52312143873279e-01 8.35656411450723e-01 9.14138339456619e-04 0.00000000000000e+00 1.89279200000009e+02 9.95919793406288e-01
311 8.67434413632281e-01 4.97008761340118e-01 8.37922761666330e-04 -4.98558173540879e-01 8.68276216890349e-01 9.03500216471809e-04 0.00000000000000e+00 1.89216400000008e+02 9.97113421966822e-01
312 8.34426257951122e-01 5.50529879282855e-01 8.94967049892886e-04 -5.52312143873279e-01 8.35656411450723e-01 9.14138339456618e-04 0.00000000000000e+00 1.89279200000009e+02 9.95919793406288e-01
......@@ -811,7 +811,7 @@ namespace fargOCA {
} else if (temp < t678*std::pow(rho,2.267e-1) || rho <= 1.0e-10){
// to avoid overflow
real ts4 = 1.e-4*temp;
real rho13 = std::pow(rho,0.333333333);
real rho13 = std::cbrt(rho);
  • @elena, I replaced pow(rho,0.333..) with cbrt(rho) which is a better way to extract \sqrt[3]{rho}, assuming that was the value we meant to compute.

    I had to regold 3D_stellar_radiation and the radiative step test.

    Can you confirm that's ok ?

Please register or sign in to reply
real rho23 = rho13*rho13;
// disjoint opacity laws for 5, 6, and 7.
......
No preview for this file type
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