anonymous_namespace{PETScNonlinearSolver.cpp} Namespace Reference

Classes
struct	PetscContext

Functions
void	printConvergenceReason (SNES const &snes_solver, SNESConvergedReason const &reason)
PetscErrorCode	updateResidual (SNES, Vec x, Vec petsc_r, void *petsc_context)
PetscErrorCode	updateJacobian (SNES, Vec, Mat J, Mat, void *petsc_context)

Function Documentation

printConvergenceReason()

void anonymous_namespace{PETScNonlinearSolver.cpp}::printConvergenceReason	(	SNES const &	snes_solver,
		SNESConvergedReason const &	reason )

Definition at line 33 of file PETScNonlinearSolver.cpp.

{
    const char* strreason;
    SNESGetConvergedReasonString(snes_solver, &strreason);
    if (reason < 0)
    {
        INFO("PETScSNES diverged reason: {:s}.", strreason);
    }
    else
    {
        INFO("PETScSNES converged reason: {:s}.", strreason);
    }
}

References INFO().

updateJacobian()

PetscErrorCode anonymous_namespace{PETScNonlinearSolver.cpp}::updateJacobian	(	SNES	,
		Vec	,
		Mat	J,
		Mat	,
		void *	petsc_context )

Definition at line 78 of file PETScNonlinearSolver.cpp.

{
    DBUG("PETScNonlinearSolver: jacobian callback called.");
    // Assume the system is already assembled.
    // Copy the results into petsc vectors.
 
    auto context = static_cast<PetscContext*>(petsc_context);
    MatAssemblyBegin(J, MAT_FINAL_ASSEMBLY);
    MatAssemblyEnd(J, MAT_FINAL_ASSEMBLY);
    MatCopy(context->J->getRawMatrix(), J, DIFFERENT_NONZERO_PATTERN);
    return 0;
}

References DBUG().

updateResidual()

PetscErrorCode anonymous_namespace{PETScNonlinearSolver.cpp}::updateResidual	(	SNES	,
		Vec	x,
		Vec	petsc_r,
		void *	petsc_context )

Definition at line 48 of file PETScNonlinearSolver.cpp.

{
    auto context = static_cast<PetscContext*>(petsc_context);
 
    DBUG("PETScNonlinearSolver: residual callback called.");
 
    VecCopy(x, context->x[context->process_id]->getRawVector());
 
    // Assemble in ogs context.
    BaseLib::RunTime time_assembly;
    time_assembly.start();
    context->system->assemble(context->x, context->x_prev, context->process_id);
 
    INFO("[time] Assembly took {} s.", time_assembly.elapsed());
    context->system->getResidual(*context->x[context->process_id],
                                 *context->x_prev[context->process_id],
                                 *context->r);
    context->r->finalizeAssembly();
    context->J->finalizeAssembly();
 
    context->system->getJacobian(*context->J);
    context->system->applyKnownSolutionsNewton(
        *context->J, *context->r, *context->x[context->process_id]);
 
    VecCopy(context->r->getRawVector(), petsc_r);
 
    return 0;
}

References DBUG(), BaseLib::RunTime::elapsed(), INFO(), and BaseLib::RunTime::start().