check number of elements in oogs callback function

Nek5000 · Aug 9, 2020 · 9a1a70b · 9a1a70b
1 parent cde0f8c
commit 9a1a70b
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 0 deletions.
diff --git a/src/libP/solvers/elliptic/src/ellipticBuildMultigridLevel.c b/src/libP/solvers/elliptic/src/ellipticBuildMultigridLevel.c
@@ -994,6 +994,7 @@ elliptic_t* ellipticBuildMultigridLevel(elliptic_t* baseElliptic, int Nc, int Nf
   if(options.compareArgs("THREAD MODEL", "OPENMP")) oogsMode = OOGS_DEFAULT;
   auto callback = [&]() // hardwaire to FP64 const coeff
     {
+      if(mesh->NlocalGatherElements == 0) return;
       occa::kernel &partialAxKernel = elliptic->partialAxKernel;
       partialAxKernel(mesh->NlocalGatherElements,
                       mesh->o_localGatherElementList,

diff --git a/src/libP/solvers/elliptic/src/ellipticSolveSetup.c b/src/libP/solvers/elliptic/src/ellipticSolveSetup.c
@@ -1009,6 +1009,7 @@ void ellipticSolveSetup(elliptic_t* elliptic, occa::properties &kernelInfo)
   if(options.compareArgs("THREAD MODEL", "OPENMP")) oogsMode = OOGS_DEFAULT;
   auto callback = [&]() // hardwired to FP64 variable coeff
     {
+      if(mesh->NlocalGatherElements == 0) return;
       occa::kernel &partialAxKernel = elliptic->partialAxKernel;
       if(elliptic->blockSolver)
           partialAxKernel(mesh->NlocalGatherElements,