TAOLinearBoundSolver.h

Note

The documentation on this page was automatically extracted from the DOLFIN C++ code and may need to be edited or expanded.

class TAOLinearBoundSolver

Parent class(es)

• Variable
• PETScObject

This class provides a bound constrained solver for a linear variational inequality defined by a matrix A and a vector b. It solves the problem:

Find \(x_l\leq x\leq x_u\) such that \((Ax-b)\cdot (y-x)\geq 0,\; \forall x_l\leq y\leq x_u\)

It is a wrapper for the TAO bound constrained solver.

Example

# Assemble the linear system
A, b = assemble_system(a, L, bc)
# Define the constraints
constraint_u = Constant(1.)
constraint_l = Constant(0.)
u_min = interpolate(constraint_l, V)
u_max = interpolate(constraint_u, V)
# Define the function to store the solution
usol=Function(V)
# Create the TAOLinearBoundSolver
solver=TAOLinearBoundSolver("tao_gpcg","gmres")
# Set some parameters
solver.parameters["monitor_convergence"]=True
solver.parameters["report"]=True
# Solve the problem
solver.solve(A, usol.vector(), b , u_min.vector(), u_max.vector())
info(solver.parameters,True)

TAOLinearBoundSolver(const std::string method="default", const std::string ksp_type="default", const std::string pc_type="default")

Create TAO bound constrained solver

std::size_t solve(const GenericMatrix& A, GenericVector& x, const GenericVector& b, const GenericVector& xl, const GenericVector& xu)

Solve the linear variational inequality defined by A and b with xl =< x <= xu

std::size_t solve(const PETScMatrix& A, PETScVector& x, const PETScVector& b, const PETScVector& xl, const PETScVector& xu)

Solve the linear variational inequality defined by A and b with xl =< x <= xu

void set_ksp(const std::string ksp_type="default")

Set PETSC Krylov Solver (ksp) used by TAO

static std::map<std::string, std::string> methods()

Return a list of available Tao solver methods

static std::map<std::string, std::string> krylov_solvers()

Return a list of available krylov solvers

static std::map<std::string, std::string> preconditioners()

Return a list of available preconditioners

static Parameters default_parameters()

Default parameter values