<div dir="ltr"><div>Hi,<br></div><div><br></div><div>Building off of snes/ex12,</div><div><br></div><div>I'm trying to use</div><div><br></div><div>-variable_coefficient field  with -petscfe_type opencl -mat_petscfe_type  opencl</div>
<div><br></div><div>constant coefficients "-variable_coefficient none" works fine, but I am getting NaN in my residual vector for spatially varying coefficients.</div><div><br></div><div>I tried changing N_t manually for N_comp_aux =1, also tried setting my variable coefficient field to a constant at each quadrature point to see if any difference. </div>
<div><br></div><div>Is there any other place I may look ? </div><div><br></div><div><a href="https://bitbucket.org/petsc/petsc/src/fb3a4ffdea0449562fb0f58dad4b44552b6d3589/src/dm/dt/interface/dtfe.c?at=master#cl-4924">https://bitbucket.org/petsc/petsc/src/fb3a4ffdea0449562fb0f58dad4b44552b6d3589/src/dm/dt/interface/dtfe.c?at=master#cl-4924</a><br>
</div><div><br></div><pre style="font-family:'Bitstream Vera Sans Mono','DejaVu Sans Mono',Monaco,monospace;font-size:12px;line-height:1.4000000000000001;margin-top:0px;margin-bottom:0px;padding:0px;color:rgb(51,51,51)">
  if (useAux) {
<a name="cl-4922" style="color:rgb(59,115,175)"></a>  ierr = PetscSNPrintfCount(string_tail, end_of_buffer - string_tail,
<a name="cl-4923" style="color:rgb(59,115,175)"></a>"    /* Load coefficients a_i for this cell */\n"
<a name="cl-4924" style="color:rgb(59,115,175);display:block;background-color:rgb(235,242,249);height:16px;width:980px"></a>"    /* TODO: This should not be N_t here, it should be N_bc*N_comp_aux */\n"
<a name="cl-4925" style="color:rgb(59,115,175)"></a>"    a_i[tidx] = coefficientsAux[Goffset+batch*N_t+tidx];\n",
<a name="cl-4926" style="color:rgb(59,115,175)"></a>                            &count);STRING_ERROR_CHECK("Message to short");
<a name="cl-4927" style="color:rgb(59,115,175)"></a>  }</pre><div><br></div><div><br></div><div><div>$ ./ireSolver -run_type full -dim 3 -petscspace_order 1 -mat_petscspace_order 0 -variable_coefficient field  -snes_type ksponly  -snes_monitor -snes_converged_reason -ksp_converged_reason -ksp_rtol 1.e-12 -pc_type bjacobi -vtk geserviceStudy0025axt1.vtk -f meshIRElores.e -sourcelandmark sourcelandmarks.vtk -targetlandmark targetlandmarks.vtk -petscfe_type opencl -mat_petscfe_type  opencl -info -info_exclude null,vec,mat -petscfe_num_blocks 2 -petscfe_num_batches 2 </div>
</div><div>.</div><div>.</div><div>.</div><div>[0] PetscFEIntegrateResidual_OpenCL(): GPU layout grid(32,59,1) block(8,1,1) with 2 batches<br></div><div><div>[0] PetscFEIntegrateResidual_OpenCL():  N_t: 8, N_cb: 2</div><div>
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div><div>[0]PETSC ERROR: Floating point exception</div><div>[0]PETSC ERROR: Vec entry at local location 120 is not-a-number or infinite at end of function: Parameter number 3</div>
<div>[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.</div><div>[0]PETSC ERROR: Petsc Development GIT revision: v3.4.4-4111-g7562b2e  GIT Date: 2014-04-18 18:37:01 -0500</div>
<div>[0]PETSC ERROR: ./ireSolver on a arch-precise-gcc-4.6.3-dbg named maeda by fuentes Mon Apr 21 22:23:25 2014</div><div>[0]PETSC ERROR: Configure options --with-shared-libraries --with-clanguage=c++ --CFLAGS=-O0 --CXXFLAGS=-O0 --download-ctetgen --download-triangle --with-debugging=yes --with-blas-lapack-dir=/opt/apps/MKL/12.1/lib/intel64/ --with-exodusii-lib="[/usr/lib/x86_64-linux-gnu/libexoIIv2.so]" --with-netcdf-dir=/usr/lib --with-hdf5-dir=/usr/lib --with-c2html=0 --with-exodusii-include=/usr/include --with-opencl-include=/usr/include/ --with-opencl-lib=/usr/lib/libOpenCL.so --download-viennacl=yes</div>
<div>[0]PETSC ERROR: #1 VecValidValues() line 34 in /opt/apps/PETSc/petsc-usr/src/vec/vec/interface/rvector.c</div><div>[0]PETSC ERROR: #2 SNESComputeFunction() line 2096 in /opt/apps/PETSc/petsc-usr/src/snes/interface/snes.c</div>
<div>[0]PETSC ERROR: #3 SNESSolve_KSPONLY() line 24 in /opt/apps/PETSc/petsc-usr/src/snes/impls/ksponly/ksponly.c</div><div>[0]PETSC ERROR: #4 SNESSolve() line 3794 in /opt/apps/PETSc/petsc-usr/src/snes/interface/snes.c</div>
<div>[0]PETSC ERROR: #5 main() line 808 in /home/fuentes/github/DMPlexApplications/ireSolver.c</div><div>[0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint@mcs.anl.gov----------</div>
<div>application called MPI_Abort(MPI_COMM_WORLD, 72) - process 0</div><div><br></div></div><div><br></div><div><br></div><div><br></div><div>thanks,</div><div>David</div></div>