[petsc-dev] [petsc-users] MATNEST with shell matrices

2013-02-21 Thread Boyce Griffith
idding doesn't preserve things like > norms, so you end up with tangled semantics.) At present all solvers operate only on fixed grid configurations. I do not expect this to change any time soon. -- Boyce -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130221/ec9220a7/attachment-0001.html>

[petsc-dev] [petsc-users] MATNEST with shell matrices

2013-02-21 Thread Jed Brown
next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130221/666fa6b4/attachment.html>

[petsc-dev] [petsc-users] MATNEST with shell matrices

2013-02-21 Thread Barry Smith
The way I view changing sizes inside a (non)linear solver is that one does an approximate (non)linear solver, decides to adapt the mesh, then does another approximate (non)linear solver etc. So even though formally they may not be stopping the SNESSolve() and restarting one can think of it th

[petsc-dev] [petsc-users] MATNEST with shell matrices

2013-02-21 Thread Jed Brown
they regrid all vectors. But regridding doesn't preserve things like norms, so you end up with tangled semantics.) -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130221/32c1fb17/attachment.html>

[petsc-dev] Alternatives to cholmod

2013-02-21 Thread Jed Brown
ld be a good package for this scenario. Are you having problems with it or just looking for other (likely slower) packages to spend time installing? -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130221/24039b76/attachment.html>

[petsc-dev] Alternatives to cholmod

2013-02-21 Thread Dave Nystrom
Jed Brown writes: > On Thu, Feb 21, 2013 at 4:47 PM, Dave Nystrom > wrote: > > > I've been using cholmod from SuiteSparse for over a year now to solve one > > of > > my linear systems that is difficult to solve with an iterative method. Are > > there other sparse direct solvers, such as

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Richard Tran Mills
I will add DMShellSetGlobalToLocal (and DMShellSetLocalToGlobal). One other thing: We also need to be do LocalToLocal operations. For the structured grid case, we call DMDALocalToLocalBegin/End, and for the unstructured grid case we construct our own VecScatter context essentially just how it i

[petsc-dev] Alternatives to cholmod

2013-02-21 Thread Dave Nystrom
I've been using cholmod from SuiteSparse for over a year now to solve one of my linear systems that is difficult to solve with an iterative method. Are there other sparse direct solvers, such as clique, that I should try sometime? My linear system is SPD and I am solving it on a single node. My

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Barry Smith
On Feb 21, 2013, at 3:31 PM, Richard Tran Mills wrote: > I will add DMShellSetGlobalToLocal (and DMShellSetLocalToGlobal). > > One other thing: We also need to be do LocalToLocal operations. For the > structured grid case, we call DMDALocalToLocalBegin/End, and for the > unstructured grid cas

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Matthew Knepley
ell. >> > > Yes > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -- next part -- An HTML attachment was scrubbed... URL: <

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Jed Brown
ng PetscSF, it may be a fine choice. -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130221/82529888/attachment.html>

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Barry Smith
Jed, The issue is not allowing or even often using SF. The issue is that DMGlobalToLocal() already has the flexibility to use a globaltolocal function pointer. And I object to hiding that or pretending it doesn't exist and saying one HAS to use SF especially if one already has code that

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Richard Tran Mills
Mills, Ph.D. Computational Earth Scientist | Joint Assistant Professor Hydrogeochemical Dynamics Team | EECS and Earth & Planetary Sciences Oak Ridge National Laboratory | University of Tennessee, Knoxville E-mail: rmills at ornl.gov V: 865-241-3198 http://climate.ornl.gov/~rmills -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130221/b1f0f352/attachment.html>

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Jed Brown
e alternative to use VecScatter is also fine, but sometimes it is easier to use PetscSF. The Point SF, in particular, is very useful for handling integer data in the setup phase. Its API is richer than VecScatter, in addition to supporting more data types. -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130221/3f3a4af6/attachment.html>

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Barry Smith
On Feb 21, 2013, at 1:39 PM, Jed Brown wrote: > > On Thu, Feb 21, 2013 at 1:31 PM, Barry Smith wrote: > I don't understand. There is a DMGlobalToLocalBegin/End() which includes > > } else { > ierr = (*dm->ops->globaltolocalbegin)(dm,g,mode == INSERT_ALL_VALUES ? > INSERT_VALUES : (mod

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Jed Brown
,PetscErrorCode (*end)(DM,Vec,InsertMode,Vec)); -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130221/6fe675cb/attachment.html>

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Richard Tran Mills
Hi Folks, I want to set a LocalToGlobal mapping in a DMShell object. Should we add a DMShell routine to do this, or just add a DMSetLocalToGlobalMapping() routine? Thanks, Richard -- Richard Tran Mills, Ph.D. Computational Earth Scientist | Joint Assistant Professor Hydrogeochemical Dyn

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Barry Smith
On Feb 21, 2013, at 1:23 PM, Jed Brown wrote: > > On Thu, Feb 21, 2013 at 1:10 PM, Matthew Knepley wrote: > Yes, he needs the local to global scatters. Rich, I think the easiest way to > do this now may be > to define your data layout as a PetscSection (these are all local) using > DMSetDefa

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Jed Brown
s, then we need to add the DMShell functions to get your functions called. -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130221/cb7a3be7/attachment-0001.html>

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Jed Brown
DMShell. > Yes -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130221/0b83367b/attachment.html>

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Peter Brune
/climate.ornl.gov/~rmills >>> >>> >> >> >> -- >> Richard Tran Mills, Ph.D. >> Computational Earth Scientist | Joint Assistant Professor >> Hydrogeochemical Dynamics Team | EECS and Earth & Planetary Sciences >> Oak Ridge National Laboratory | University of Tennessee, Knoxville >> E-mail: rmills at ornl.gov V: 865-241-3198 http://climate.ornl.gov/~rmills >> >> > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130221/5ef7cfc6/attachment.html>

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Jed Brown
like Rich wants to run his scatters using DMGlobalToLocalBegin/End, and doesn't really care about about the ISLocalToGlobalMapping. -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130221/8812a85d/attachment.html>

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Jed Brown
-- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130221/b27b401e/attachment.html>

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Dmitry Karpeev
| EECS and Earth & Planetary Sciences > Oak Ridge National Laboratory | University of Tennessee, Knoxville > E-mail: rmills at ornl.gov V: 865-241-3198 http://climate.ornl.gov/~rmills > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130221/e2c56679/attachment-0001.html>

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Peter Brune
mical Dynamics Team | EECS and Earth & Planetary Sciences >> Oak Ridge National Laboratory | University of Tennessee, Knoxville >> E-mail: rmills at ornl.gov V: 865-241-3198 http://climate.ornl.gov/~** >> rmills <http://climate.ornl.gov/~rmills> >> >> >

[petsc-dev] Need mechanism to set LocalToGlobal mapping in a DMShell

2013-02-21 Thread Jed Brown
limate.ornl.gov/~**rmills<http://climate.ornl.gov/~rmills> > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130221/87e69ab0/attachment.html>

[petsc-dev] [SIAM-SC] DOE Workshop on Applied Mathematics Research for Exascale Computing

2013-02-21 Thread Barry Smith
On Feb 21, 2013, at 7:00 AM, Karl Rupp wrote: > Hi guys, > > it might be interesting to write a position paper from PETSc's point of view. > It's just that I'm afraid that nobody will actually care about the issues of > solving really interesting problems? You are correct, no one involved

[petsc-dev] Fwd: [SIAM-SC] DOE Workshop on Applied Mathematics Research for Exascale Computing

2013-02-21 Thread Karl Rupp
Hi guys, it might be interesting to write a position paper from PETSc's point of view. It's just that I'm afraid that nobody will actually care about the issues of solving really interesting problems... Best regards, Karli Original Message Subject: [SIAM-SC] DOE Workshop on