Parallel Colt 0.7.2

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

java.lang.Object
  extended by cern.colt.matrix.tfloat.algo.solver.AbstractFloatIterativeSolver
      extended by 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
static float sqrteps
           
 
Constructor Summary
FloatMRNSD()
           
 
Method Summary
 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

Parallel Colt 0.7.2

Jump to the Parallel Colt Homepage