date:20130227

[petsc-dev] Use KSPSetUpNorms_Private() in KSPSetFromOptions()

2013-02-27 Thread Dmitry Karpeev

move them to KSPSetUp() then the user cannot > reset them. > >Barry > > > > KSPSetUpNorms_Private() will be just about the first thing called in > KSPSetUp() and will compute the appropriate normtype and pcside, > > but now based on complete KSP opt

[petsc-dev] Deflated Krylov solvers for PETSc

2013-02-27 Thread Jed Brown

t; Email: jiechen at mcs.anl.gov > Homepage: http://www.mcs.anl.gov/~jiechen > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130227/5cad9baf/attachment.html>

[petsc-dev] Deflated Krylov solvers for PETSc

2013-02-27 Thread Matthew Knepley

t, IL 60439 >> Phone: (630) 252-3313 >> Email: jiechen at mcs.anl.gov >> Homepage: http://www.mcs.anl.gov/~jiechen >> > > -- 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: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130227/1a57bcc0/attachment.html>

[petsc-dev] Deflated Krylov solvers for PETSc

2013-02-27 Thread Jie Chen

My personal experience is to use approximate eigenvectors to deflate conjugate gradient. In that regard, it is kinda like DGMRES, except that the iterative method is different. A lot of situations I have seen use eigenvectors, but just as what Vuik does, you never know when someone will propose

[petsc-dev] Use KSPSetUpNorms_Private() in KSPSetFromOptions()

2013-02-27 Thread Dmitry Karpeev

te() is called > after KSP type has been set to KSPPREONLY, > > > but BEFORE there was a chance to set the PC side! > > > > > > Is there any reason to keep these calls to KSPSetUpNorms_Private() in > KSPSetFromOptions()? > > > >I think they are there so that users can overwrite them in the rest > of KSPSetFromOptions(). If you move them to KSPSetUp() then the user cannot > reset them. > > > >Barry > > > > > > > KSPSetUpNorms_Private() will be just about the first thing called in > KSPSetUp() and will compute the appropriate normtype and pcside, > > > but now based on complete KSP options. > > > > > > Not to mention that calling XXXSetUp_YYY() XXXSetFromOptions() seems > generally contrary to PETSc's design > > > philosophy. > > > > > > Dmitry. > > > > > > > > > > > > > > > > > > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130227/84d4f6a7/attachment-0001.html>

[petsc-dev] Deflated Krylov solvers for PETSc

2013-02-27 Thread Barry Smith

The PETSc Krylov methods are supposed to be implemented using KSP_MatMult(), KSP_PCApply() etc instead of directly using MatMult() and PCApply() etc. This means, in theory, one can add support for techniques that modify the application of the multiply or preconditioner in a single place and

[petsc-dev] Vector operations on ghosted vectors

2013-02-27 Thread John Mousel

I'd think something like VecCopyGhosted would be fairly useful. John -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130227/476d9aba/attachment.html>

[petsc-dev] How to view scatter plot of petsc matrix

2013-02-27 Thread Nystrom, William D

Is there a way to view a scatter plot of the non-zero elements of a petsc matrix after it and the right hand side have been created in petsc - so that the matrix structure can be visually inspected? Thanks, Dave -- Dave Nystrom LANL HPC-5 Phone: 505-667-7913 Email: wdn at lanl.gov Smail: Mail St