diff --git a/src/libnptm/clu_subs.cpp b/src/libnptm/clu_subs.cpp index d83384a7441d61e0e3d23ab02506985be1009709..cf643b1a9bebcb00e1ed020642b325d52526c9fb 100644 --- a/src/libnptm/clu_subs.cpp +++ b/src/libnptm/clu_subs.cpp @@ -2156,7 +2156,7 @@ void scr0(double vk, double exri, C1 *c1, C1_AddOns *c1ao, C3 *c3, C4 * c4) { nvtxRangePush("scr0 outer loop"); #endif -#pragma omp parallel for reduction(+:scs, ecs, acs, tfsas) + //#pragma omp parallel for reduction(+:scs, ecs, acs, tfsas) for (int i14 = 1; i14 <= c4->nsph; i14++) { int iogi = c1->iog[i14 - 1]; if (iogi >= i14) { @@ -2196,12 +2196,12 @@ void scr0(double vk, double exri, C1 *c1, C1_AddOns *c1ao, C3 *c3, C4 * c4) { acs += c1->sabs[iogi - 1]; tfsas += c1->fsas[iogi - 1]; } // i14 loop - c3->scs = scs; - c3->ecs = ecs; - c3->acs = acs; - c3->tfsas = tfsas; + c3->scs = scs; + c3->ecs = ecs; + c3->acs = acs; + c3->tfsas = tfsas; #ifdef USE_NVTX - nvtxRangePop(); + nvtxRangePop(); #endif } @@ -2233,7 +2233,7 @@ void scr2( #ifdef USE_NVTX nvtxRangePush("scr2 outer loop 1"); #endif -#pragma omp parallel for reduction(+:tsas00, tsas10, tsas01, tsas11) + //#pragma omp parallel for reduction(+:tsas00, tsas10, tsas01, tsas11) for (int i14 = 1; i14 <= c4->nsph; i14++) { int i = i14 - 1; int iogi = c1->iog[i14 - 1]; @@ -2299,6 +2299,11 @@ void scr2( #ifdef USE_NVTX nvtxRangePop(); #endif + + + + + c3->tsas[0][0] = tsas00; c3->tsas[1][0] = tsas10; c3->tsas[0][1] = tsas01;