Commit ba9d7aab authored by Giovanni La Mura's avatar Giovanni La Mura
Browse files

Reset cumulative cross-section at 0 in scale iteration

parent 62e4b64f

src/include/clu_subs.h

+3 −0
Original line number Diff line number Diff line
@@ -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];