Re: [petsc-users] MatNestGetISs in fortran

2014-06-13 Thread Klaij, Christiaan
Perhaps this message from May 27 slipped through the email cracks as Matt puts it? Chris dr. ir. Christiaan Klaij CFD Researcher Research Development E mailto:c.kl...@marin.nl T +31 317 49 33 44 MARIN 2, Haagsteeg, P.O. Box 28, 6700 AA Wageningen, The Netherlands T +31 317 49 39 11, F +31

[petsc-users] Interpreting -log_summary, amount of communication

2014-06-13 Thread Åsmund Ervik
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Dear PETSc, First of all, bug report for the manual (for petsc-current): in Fig. 20 and 21, something has not gone well with \href and listings, so I can't understand those figures properly. I read the chapter in the manual, 12.1, but it din't

Re: [petsc-users] Ensuring non-negative solutions

2014-06-13 Thread Jed Brown
Gautam Bisht gbi...@lbl.gov writes: PETSc, I had few questions regarding adaptive time stepping inTS. - Can temporal adaptivity in TS ensure non-negative solution of a PDE? If short steps are sufficient to prevent negative solutions, then it can. The simplest way is to set your own step

[petsc-users] Multi-DOF DMDA Vec

2014-06-13 Thread Anush Krishnan
Hello petsc-users, I created a vector using DMDACreate with 3 degrees of freedom. Is it possible for me to access each vector corresponding to a degree of freedom? Seeing that I need to access the array as [k][j][i][dof], does it mean that the values of each component are not contiguous? Also,

Re: [petsc-users] Multi-DOF DMDA Vec

2014-06-13 Thread Jed Brown
Anush Krishnan an...@bu.edu writes: Hello petsc-users, I created a vector using DMDACreate with 3 degrees of freedom. Is it possible for me to access each vector corresponding to a degree of freedom? Seeing that I need to access the array as [k][j][i][dof], does it mean that the values of

Re: [petsc-users] Multi-DOF DMDA Vec

2014-06-13 Thread Anush Krishnan
Thanks, Jed! On 13 June 2014 14:43, Jed Brown j...@jedbrown.org wrote: Anush Krishnan an...@bu.edu writes: Hello petsc-users, I created a vector using DMDACreate with 3 degrees of freedom. Is it possible for me to access each vector corresponding to a degree of freedom? Seeing that

Re: [petsc-users] Multi-DOF DMDA Vec

2014-06-13 Thread Barry Smith
On Jun 13, 2014, at 1:43 PM, Jed Brown j...@jedbrown.org wrote: Anush Krishnan an...@bu.edu writes: Hello petsc-users, I created a vector using DMDACreate with 3 degrees of freedom. Is it possible for me to access each vector corresponding to a degree of freedom? Seeing that I need to

Re: [petsc-users] Interpreting -log_summary, amount of communication

2014-06-13 Thread Barry Smith
On Jun 13, 2014, at 6:57 AM, Åsmund Ervik asmund.er...@ntnu.no wrote: -BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Dear PETSc, First of all, bug report for the manual (for petsc-current): in Fig. 20 and 21, something has not gone well with \href and listings, so I can't understand

Re: [petsc-users] Multi-DOF DMDA Vec

2014-06-13 Thread Anush Krishnan
On 13 June 2014 14:59, Barry Smith bsm...@mcs.anl.gov wrote: On Jun 13, 2014, at 1:43 PM, Jed Brown j...@jedbrown.org wrote: Anush Krishnan an...@bu.edu writes: Hello petsc-users, I created a vector using DMDACreate with 3 degrees of freedom. Is it possible for me to access each

Re: [petsc-users] Multi-DOF DMDA Vec

2014-06-13 Thread Jed Brown
Anush Krishnan an...@bu.edu writes: With regard to the interlaced memory performing better: If I used three vectors created from the same DMDA for each degree of freedom, how different would that be in performance compared to a fully interlaced vector? Wouldn't cache reuse be about the same

Re: [petsc-users] Multi-DOF DMDA Vec

2014-06-13 Thread Barry Smith
The main reason to “pull out” a single component is, for example, to solve a linear system for that single component; that is, to work on that single component a great deal. You wouldn’t pull out the individual components to iterate on them all together. Barry On Jun 13, 2014, at 9:09

Re: [petsc-users] Multi-DOF DMDA Vec

2014-06-13 Thread Anush Krishnan
On 13 June 2014 22:22, Barry Smith bsm...@mcs.anl.gov wrote: The main reason to “pull out” a single component is, for example, to solve a linear system for that single component; that is, to work on that single component a great deal. You wouldn’t pull out the individual components to