[petsc-dev] SNESSetSolution and SNESSolve

[Barry Smith]

    I suspect Matt added it so he could have 

ierr = SNESGetDM(snes, &dm);CHKERRQ(ierr);
ierr = DMShellSetGlobalVector(dm, u);CHKERRQ(ierr);

and not get complaints from DMShell.  I absolutely hate this ad hoc mess; surely there is a better way to handle it. We need to remove this thing, reproduce the problem Matt was having and find a better fix.

  Barry




> On May 26, 2015, at 5:47 AM, Lisandro Dalcin <dalcinl at gmail.com> wrote:
> 
> Matt, you added SNESSetSolution(), however SNESSolve(snes,NULL,NULL)
> does not reuse vec_sol if it is already set. IMHO, it should.
> 
> What other people think about it?
> 
> 
> -- 
> Lisandro Dalcin
> ============
> Research Scientist
> Computer, Electrical and Mathematical Sciences & Engineering (CEMSE)
> Numerical Porous Media Center (NumPor)
> King Abdullah University of Science and Technology (KAUST)
> http://numpor.kaust.edu.sa/
> 
> 4700 King Abdullah University of Science and Technology
> al-Khawarizmi Bldg (Bldg 1), Office # 4332
> Thuwal 23955-6900, Kingdom of Saudi Arabia
> http://www.kaust.edu.sa
> 
> Office Phone: +966 12 808-0459