From ba9d7aab15b7e547f7f7696b25f36259b84ae61c Mon Sep 17 00:00:00 2001
From: Giovanni La Mura <giovanni.lamura@inaf.it>
Date: Thu, 7 Dec 2023 18:09:52 +0100
Subject: [PATCH] Reset cumulative cross-section at 0 in scale iteration
---
src/include/clu_subs.h | 3 +++
1 file changed, 3 insertions(+)
diff --git a/src/include/clu_subs.h b/src/include/clu_subs.h
index c025610a..df1d1955 100644
--- a/src/include/clu_subs.h
+++ b/src/include/clu_subs.h
@@ -1918,6 +1918,9 @@ void scr0(double vk, double exri, C1 *c1, C1_AddOns *c1ao, C3 *c3, C4 * c4) {
std::complex<double> sum21, rm, re, csam;
csam = -(ccs / (exri * vk)) * std::complex<double>(0.0, 0.5);
//double scs = 0.0, ecs = 0.0, acs = 0.0;
+ c3->scs = 0.0;
+ c3->ecs = 0.0;
+ c3->acs = 0.0;
c3->tfsas = cc0;
for (int i14 = 1; i14 <= c4->nsph; i14++) {
int iogi = c1->iog[i14 - 1];
--
GitLab