Commit b6a10e5b authored by Mulas, Giacomo's avatar Mulas, Giacomo
Browse files

comment out two parallel omp loops in scr0 and scr2 that cause problems

parent 8db80d37

src/libnptm/clu_subs.cpp

+12 −7
Original line number Diff line number Diff line
@@ -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) {
@@ -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;