Re: [petsc-users] MatNestGetISs in fortran

2014-07-07 Thread Klaij, Christiaan
Satish, Thanks for your reply! I've been on holiday for some time, do you still want me to test this or has it been fixed with the 3.5 release? Chris From: Satish Balay ba...@mcs.anl.gov Sent: Monday, June 16, 2014 5:41 PM To: Klaij, Christiaan Cc:

[petsc-users] Jacobian construction, DA vs Plex

2014-07-07 Thread Tobin Isaac
Hi, If I have a pointwise Jacobian function f, I know that I can call DMGetDS() and pass f to PetscDSSetJacobian(), and that f will be used by PetscFEIntegrate() and thus by DMPlexSNESComputeJacobianFEM(). It looks like PetscFEIntegrate() is only used by plex and not da. Is there anyway that I

Re: [petsc-users] Jacobian construction, DA vs Plex

2014-07-07 Thread Matthew Knepley
On Mon, Jul 7, 2014 at 5:36 PM, Tobin Isaac tis...@ices.utexas.edu wrote: Hi, If I have a pointwise Jacobian function f, I know that I can call DMGetDS() and pass f to PetscDSSetJacobian(), and that f will be used by PetscFEIntegrate() and thus by DMPlexSNESComputeJacobianFEM(). It looks

Re: [petsc-users] MatNestGetISs in fortran

2014-07-07 Thread Satish Balay
My change is in petsc-3.5. You can upgrade and see if it works for you.. Satish On Mon, 7 Jul 2014, Klaij, Christiaan wrote: Satish, Thanks for your reply! I've been on holiday for some time, do you still want me to test this or has it been fixed with the 3.5 release? Chris

Re: [petsc-users] Jacobian construction, DA vs Plex

2014-07-07 Thread Tobin Isaac
On Mon, Jul 07, 2014 at 05:43:48PM +0200, Matthew Knepley wrote: On Mon, Jul 7, 2014 at 5:36 PM, Tobin Isaac tis...@ices.utexas.edu wrote: Hi, If I have a pointwise Jacobian function f, I know that I can call DMGetDS() and pass f to PetscDSSetJacobian(), and that f will be used by

[petsc-users] [petsc4py] vec.getCUSPArray for cusp vecs

2014-07-07 Thread Ashwin Srinath
Hello, petsc-users, In PETSc, we have VecCUSPGetArray to access the underlying array from a Vec of type 'cusp'. Does petsc4py offer something similar? If not, I'm considering using petsc4py and pyCUDA together in some way, and I'd be extremely grateful for any advice on doing that. Thanks so