Parallel Colt 0.7.2

cern.colt.matrix.tdouble.algo.solver
Class DoubleMRNSD

java.lang.Object
  extended by cern.colt.matrix.tdouble.algo.solver.AbstractDoubleIterativeSolver
      extended by cern.colt.matrix.tdouble.algo.solver.DoubleMRNSD
All Implemented Interfaces:
DoubleIterativeSolver

public class DoubleMRNSD
extends AbstractDoubleIterativeSolver

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 double sqrteps
           
 
Constructor Summary
DoubleMRNSD()
           
 
Method Summary
 DoubleMatrix1D solve(DoubleMatrix2D A, DoubleMatrix1D b, DoubleMatrix1D x)
          Solves the given problem, writing result into the vector.
 
Methods inherited from class cern.colt.matrix.tdouble.algo.solver.AbstractDoubleIterativeSolver
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 double sqrteps
Constructor Detail

DoubleMRNSD

public DoubleMRNSD()
Method Detail

solve

public DoubleMatrix1D solve(DoubleMatrix2D A,
                            DoubleMatrix1D b,
                            DoubleMatrix1D x)
                     throws IterativeSolverDoubleNotConvergedException
Description copied from interface: DoubleIterativeSolver
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:
IterativeSolverDoubleNotConvergedException

Parallel Colt 0.7.2

Jump to the Parallel Colt Homepage