com.numericalmethod.suanshu.matrix.doubles.matrixtype.sparse.solver.iterative.nonstationary
Class SteepestDescentSolver

java.lang.Object
  com.numericalmethod.suanshu.matrix.doubles.matrixtype.sparse.solver.iterative.nonstationary.SteepestDescentSolver

All Implemented Interfaces:

IterativeSolver

public class SteepestDescentSolver
extends java.lang.Object
implements IterativeSolver

The Steepest Descent method (SDM) can solve symmetric n-by-n linear systems. The tolerance is guaranteed if A is symmetric positive definite.

Only left preconditioning is supported in this implementation. It requires in addition that the preconditioner is symmetric and positive definite.

See Also:

Wikipedia: Gradient descent

Nested Class Summary

Nested classes/interfaces inherited from interface com.numericalmethod.suanshu.matrix.doubles.matrixtype.sparse.solver.iterative.IterativeSolver

IterativeSolver.ConvergenceFailure, IterativeSolver.Problem

Field Summary

static int

DEFAULT_RESIDUAL_REFRESH_RATE The algorithm recomputes the residual as b - Axi once per this number of iterations

Constructor Summary

SteepestDescentSolver()

SteepestDescentSolver(int residualRefreshRate)
The solver recomputes the residual as b - Ax_i once per this number of iterations

Method Summary

Vector	solve(IterativeSolver.Problem problem) Solve iteratively Ax = b until the solution is close enough, i.e., the norm of residual (b - Ax) is less than or equal to the specified iteration.
Vector	solve(IterativeSolver.Problem problem, IterationMonitor monitor) Solve iteratively Ax = b until the solution is close enough, i.e., the norm of residual (b - Ax) is less than or equal to the specified iteration.

Methods inherited from class java.lang.Object

clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Field Detail

DEFAULT_RESIDUAL_REFRESH_RATE

public static final int DEFAULT_RESIDUAL_REFRESH_RATE

The algorithm recomputes the residual as b - Ax_i once per this number of iterations

See Also:

Constant Field Values

Constructor Detail

SteepestDescentSolver

public SteepestDescentSolver()

SteepestDescentSolver

public SteepestDescentSolver(int residualRefreshRate)

The solver recomputes the residual as b - Ax_i once per this number of iterations

Parameters:

residualRefreshRate - the number of iterations before the next refresh

Method Detail

solve

public Vector solve(IterativeSolver.Problem problem)
             throws IterativeSolver.ConvergenceFailure

Description copied from interface: IterativeSolver

Solve iteratively

Ax = b

until the solution is close enough, i.e., the norm of residual (b - Ax) is less than or equal to the specified iteration.

Specified by:

solve in interface IterativeSolver

Parameters:

problem - the problem of solving Ax = b

Returns:

the computed solution for the problem

Throws:

IterativeSolver.ConvergenceFailure - if the algorithm fails to converge

solve

public Vector solve(IterativeSolver.Problem problem,
                    IterationMonitor monitor)
             throws IterativeSolver.ConvergenceFailure

Description copied from interface: IterativeSolver

Solve iteratively

Ax = b

until the solution is close enough, i.e., the norm of residual (b - Ax) is less than or equal to the specified iteration.

In each iteration, the newly computed iterate is added to the IterationMonitor for statistics or diagnostic purpose.

Specified by:

solve in interface IterativeSolver

Parameters:

problem - the problem of solving Ax = b

monitor - an IterationMonitor instance

Returns:

the computed solution for the problem

Throws:

IterativeSolver.ConvergenceFailure - if the algorithm fails to converge