Keep rank.m

Truncates the singular value decomposition at the specified rank and reassembles the matrix.

Syntax

    A=keep_rank(A,rank)

Arguments

        A     -  real or complex matrix, will be 
                 converted to full if a sparse
                 matrix is received

        rank  -  number of singular values to keep

Outputs

        A     -  filtered matrix, returned as full

See also

Kernel_utilities#Numerical_infrastructure

Version 2.5, authors: Ilya Kuprov