cern.colt.matrix.tfloat.algo.solver

Class FloatMRNSD

java.lang.Object

cern.colt.matrix.tfloat.algo.solver.AbstractFloatIterativeSolver

cern.colt.matrix.tfloat.algo.solver.FloatMRNSD

All Implemented Interfaces:

FloatIterativeSolver

public class FloatMRNSD
extends AbstractFloatIterativeSolver

MRNSD is Modified Residual Norm Steepest Descent method used for solving large-scale, ill-posed inverse problems of the form: b = A*x + noise. This algorithm is nonnegatively constrained.

References:

[1] J. Nagy, Z. Strakos, "Enforcing nonnegativity in image reconstruction algorithms" in Mathematical Modeling, Estimation, and Imaging, David C. Wilson, et.al., Eds., 4121 (2000), pg. 182--190.

[2] L. Kaufman, "Maximum likelihood, least squares and penalized least squares for PET", IEEE Trans. Med. Imag. 12 (1993) pp. 200--214.

Author:

Piotr Wendykier (piotr.wendykier@gmail.com)

Field Summary

Fields

Modifier and Type	Field and Description
static float	sqrteps

Constructor Summary

Constructors

Constructor and Description
FloatMRNSD()

Method Summary

Methods

Modifier and Type	Method and Description
FloatMatrix1D	solve(FloatMatrix2D A, FloatMatrix1D b, FloatMatrix1D x) Solves the given problem, writing result into the vector.

Methods inherited from class cern.colt.matrix.tfloat.algo.solver.AbstractFloatIterativeSolver

getIterationMonitor, getPreconditioner, setIterationMonitor, setPreconditioner

Methods inherited from class java.lang.Object

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

Field Detail

sqrteps

public static final float sqrteps

Constructor Detail

FloatMRNSD

public FloatMRNSD()

Method Detail

solve

public FloatMatrix1D solve(FloatMatrix2D A,
                  FloatMatrix1D b,
                  FloatMatrix1D x)
                    throws IterativeSolverFloatNotConvergedException

Description copied from interface: FloatIterativeSolver

Solves the given problem, writing result into the vector.

Parameters:

A - Matrix of the problem

b - Right hand side

x - Solution is stored here. Also used as initial guess

Returns:

The solution vector x

Throws:

IterativeSolverFloatNotConvergedException