SNES: Add flag indicating if final residual is always computedAlthough SNESSolve_XXX always updates the solution variable, it may notbe the case that the residual in vec_func is the residual at th
SNES: Add flag indicating if final residual is always computedAlthough SNESSolve_XXX always updates the solution variable, it may notbe the case that the residual in vec_func is the residual at the newsolution. Track this behaviour explicitly for all concreteimplementations. We will use this in SNESSolve_FAS to determine if theresidual must be computed before transferring it between levels.
show more ...
removed a bunch of unneeded inclusions of xxximpl.h which exposed private data where it did not have to be exposedalso made more systematic the inclusion of the base xxximpl.h in the derived yyyimp
removed a bunch of unneeded inclusions of xxximpl.h which exposed private data where it did not have to be exposedalso made more systematic the inclusion of the base xxximpl.h in the derived yyyimpl.hso the inclusion is done in the include file rather than in each source code that includes it
Updated script for detecting visible symbolsRemoved unneeded visibility of some symbolsFixed some namespacingRemoved the use of unneeded implementation includes in the TAO code
Merge branch 'pr401/tferma/virs-pcfieldsplit-fix/master'
When SNES VI fails with a divered linear solve destroy all the objects before breaking from the iteration loop.
A fix for fieldsplit pcs when using snes virs.When reset pc fieldsplit saves a copy of the original IS. SNES VIRS then calls a new function, PCFieldSplitRestrictIS, and suppliesthe inactive set t
A fix for fieldsplit pcs when using snes virs.When reset pc fieldsplit saves a copy of the original IS. SNES VIRS then calls a new function, PCFieldSplitRestrictIS, and suppliesthe inactive set to create new ISs for the fieldsplit. This is then reindexed and passed down to any nested fieldsplits.
fixed up formatting of references in manualpages to consistent styleReferences:+ 1. - ref1. 2. - ref2- 3. - ref3previously the formatting was varied and ugly
broken nonlinear solver table links fixedReported-by: Andrew McRae <A.T.T.McRae@bath.ac.uk>Also made References: formatting consistent by always listing authors first
begun adding more monitoring functionality to SNESLineSearch
snesvi should only use DM to compute bounds if it can compute boundsfixes for shashi example, converges to some degree with rs method
Begin renaming the PetscOptions objects for the future where PetscOptions will not be a singleton
Introduce MPIU_Allreduce() to detect cases where all processes that share an MPI comm do not call the same sequence of MPI_Allreduce()This is a bug that can be difficult to track down.Suggested-by
Introduce MPIU_Allreduce() to detect cases where all processes that share an MPI comm do not call the same sequence of MPI_Allreduce()This is a bug that can be difficult to track down.Suggested-by: Eric Chamberland <Eric.Chamberland@giref.ulaval.ca>Suggested-by: Patrick Lacasse <placasse@giref.ulaval.ca>
added error checking for when PetscViewer is not passed to the standard monitor routinesReported-by: Garth N. Wells <gnw20@cam.ac.uk>Removed outdated output in example test file
changed SNESVIMonitor to only display residual norm with %g
initial locks with concurrencykit and OpenMPcleanup of many monitors that used nonthread safe constructors when no viewer passed incleanup of examples that generated unused variable messages when
initial locks with concurrencykit and OpenMPcleanup of many monitors that used nonthread safe constructors when no viewer passed incleanup of examples that generated unused variable messages when compiled with --with-log=0
Have SNESVIGetInactiveSet return the *actual* current inactive set.Previously it returned the previous inactive set.
fixed definition of active constraints to only include f > 0 or f < 0 to match citationremoved dead vi examplesremoved dead vi routineupdated vi examples outputadded command line option to zero e
fixed definition of active constraints to only include f > 0 or f < 0 to match citationremoved dead vi examplesremoved dead vi routineupdated vi examples outputadded command line option to zero equality of x[] value with its constraint instead of hardwired 1.e-8Reported-by: Asbjørn Nilsen Riseth <riseth@maths.ox.ac.uk>
The default bounds for SNESVISetVariableBounds() in the manual page were reversed from correct valuesReported-by: Mark Lohry <mlohry@gmail.com>
Infrastructure that allows failures in PCSetUp(), PCApply(), MatMult() etc due to, for example a zero pivot or function evaluation outside its domain, to propagate up and become a KSP_DIVERGED inste
Infrastructure that allows failures in PCSetUp(), PCApply(), MatMult() etc due to, for example a zero pivot or function evaluation outside its domain, to propagate up and become a KSP_DIVERGED instead of generating an error that stops the program. In response to Issue 96. This includes failures in MatCreateSNESMF() applications due to domain errors. The mechanism to propagate some errors is by setting Info or Nan into the output vector and using the norm or inner product reductions in SNES or KSP to propagate the error condition to all processes and then handling them immediately after the norm or inner product.This allows, for example, ODE integrators to try again with a smaller time-step if the PCSetUp failed instead of requiring a complete restart of the run with other options.Currently some error conditions, such as function domain error in a line search may not get progated up using the correct SNESConvergedReasonSee src/snes/examples/tests/ex69.c for the handling of several conditions
Support the case where the user specifies a function to set the bounds.
If the user passes in variable bounds, just set them on snes->xl and snes->xu;don't set the VI type to SNESVINEWTONRSLS.
updated PETSc directory layout to match standard packaging strategiesinclude/petsc finclude,private,mpiunilib/petsc confbin/petsc*
Fix obviously-wrong aspects of documentation on SNESVINEWTONRSLS and SNESVINEWTONSSLS.
fixed memory leak introduced in SNESVIMonitorResidual()fixed automatic use of -snes_vi_monitor_residual when -snes_vi_monitor was set
if SNESSolve_VINEWTONRSLS was broken out early then make sure the DM is cleaned up and VI related information removed
12345678910>>...16