[petsc-users] petscfe_type opencl, variable_coefficient field

Hi, Building off of snes/ex12, I'm trying to use -variable_coefficient field with -petscfe_type opencl -mat_petscfe_type opencl constant coefficients "-variable_coefficient none" works fine, but I am getting NaN in my residual vector for spatially varying coefficients. I tried changing N_t m

Re: [petsc-users] Elasticity tensor in ex52

Thank you. On Mon, Apr 21, 2014 at 10:24 AM, Matthew Knepley wrote: > On Mon, Apr 21, 2014 at 10:23 AM, Miguel Angel Salazar de Troya < > salazardetr...@gmail.com> wrote: > >> Thank you. Now it makes sense. Just to confirm though, back to the weak >> form I wrote above in index notation: >> >>

Re: [petsc-users] Elasticity tensor in ex52

On Mon, Apr 21, 2014 at 10:23 AM, Miguel Angel Salazar de Troya < salazardetr...@gmail.com> wrote: > Thank you. Now it makes sense. Just to confirm though, back to the weak > form I wrote above in index notation: > > phi_{i,j} C_{ijkl} u_{k,l} > > The g3 indices are equivalent to the indices of th

Re: [petsc-users] Elasticity tensor in ex52

Thank you. Now it makes sense. Just to confirm though, back to the weak form I wrote above in index notation: phi_{i,j} C_{ijkl} u_{k,l} The g3 indices are equivalent to the indices of this weak form in this manner: ic = i jc = k id = j jd = l Is this correct? This is what I understand when you

Re: [petsc-users] an ambiguity on a Lanczos solver for a symmetric system -- advice needed

Sorry, ignore my previous email. You are not solving eigenvalue problem. Hong On Mon, Apr 21, 2014 at 9:59 AM, Hong Zhang wrote: > Umut, > Have you tried slepc for eigenvalue problems? > Why do you need mumps in your eigensolver? Shift-and-invert? > > Hong > > On Mon, Apr 21, 2014 at 9:39 AM, Mat

Re: [petsc-users] an ambiguity on a Lanczos solver for a symmetric system -- advice needed

Umut, Have you tried slepc for eigenvalue problems? Why do you need mumps in your eigensolver? Shift-and-invert? Hong On Mon, Apr 21, 2014 at 9:39 AM, Matthew Knepley wrote: > On Sat, Apr 19, 2014 at 2:13 PM, Umut Tabak wrote: >> >> Dear all, > > > For any timing question, we need to see the ou

Re: [petsc-users] Elasticity tensor in ex52

On Sun, Apr 20, 2014 at 4:17 PM, Miguel Angel Salazar de Troya < salazardetr...@gmail.com> wrote: > I understand. So if we had a linear elastic material, the weak form would > be something like > > phi_{i,j} C_{ijkl} u_{k,l} > > where the term C_{ijkl} u_{k,l} would correspond to the term f_1(U, >

Re: [petsc-users] an ambiguity on a Lanczos solver for a symmetric system -- advice needed

On Sat, Apr 19, 2014 at 2:13 PM, Umut Tabak wrote: > Dear all, > For any timing question, we need to see the output f -log_summary. Also, if you have significant time in routines you wrote, we need you to create PETSc events for these. Matt > I am experiencing lately some issues with a symm

Re: [petsc-users] libpetsc.so: undefined reference to `_gfortran_transfer_character_write@GFORTRAN_1.4'

Thanks Satish! I am trying to fix the issue following your suggestions. I will let you know. Thanks On Fri, Apr 18, 2014 at 12:21 PM, Satish Balay wrote: > detailed logs would have been useful. > > Mostlikely you have multiple version of libgfortran.so - and the wrong > version got linked in.