improved GL integration.

This commit is contained in:
suter_a 2025-05-27 15:59:08 +02:00
parent e525abe614
commit d1f025a8c2

View File

@ -78,8 +78,10 @@ PLGKT_LF::PLGKT_LF(std::vector<double> &param, const double tmax) : fParam(param
*/
int PLGKT_LF::DynamicLGKTLF()
{
std::vector<double> rr={0.2, 0.4, 0.6, 0.8, 1.0, 2.0, 3.0, 4.0, 5.0, 6.25, 7.5, 8.75, 10.0, 32.5, 55.0, 77.5, 100.0};
std::vector<double> ww={0.4, 1.6, 0.8, 1.6, 2.4, 8.0, 4.0, 8.0, 4.5, 10.0, 5.0, 10.0, 47.5, 180.0, 90.0, 180.0, 45.0};
// std::vector<double> rr={0.2, 0.4, 0.6, 0.8, 1.0, 2.0, 3.0, 4.0, 5.0, 6.25, 7.5, 8.75, 10.0, 32.5, 55.0, 77.5, 100.0};
// std::vector<double> ww={0.4, 1.6, 0.8, 1.6, 2.4, 8.0, 4.0, 8.0, 4.5, 10.0, 5.0, 10.0, 47.5, 180.0, 90.0, 180.0, 45.0};
std::vector<double> rr={0.2, 0.4, 0.6, 0.8, 1.0, 1.25, 1.5, 1.75, 2.0, 2.5, 3.0, 4.0, 5.0, 7.5, 10.0, 12.8125, 15.625, 18.4375, 21.25, 26.875, 32.5, 43.75, 55.0, 77.5, 100.0};
std::vector<double> ww={0.4, 1.6, 0.8, 1.6, 0.9, 2.0, 1.0, 2.0, 1.5, 4.0, 3.0, 8.0, 7.0, 20.0, 10.625, 22.5, 11.25, 22.5, 16.875, 45.0, 33.75, 90.0, 67.5, 180.0, 45.0};
std::vector<double> pp={fParam[0], 0.0, fParam[2]};
std::vector<double> pol;