/*******************************************************************************
*   Copyright 2012 Analog Devices, Inc.
*
*   Licensed under the Apache License, Version 2.0 (the "License");
*   you may not use this file except in compliance with the License.
*   You may obtain a copy of the License at
*
*       http://www.apache.org/licenses/LICENSE-2.0
*
*   Unless required by applicable law or agreed to in writing, software
*   distributed under the License is distributed on an "AS IS" BASIS,
*   WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
*   See the License for the specific language governing permissions and
*   limitations under the License.
********************************************************************************/

package com.analog.lyric.math;

import org.apache.commons.math3.linear.RealMatrix;
import org.apache.commons.math3.linear.SingularValueDecomposition;

import com.analog.lyric.dimple.exceptions.DimpleException;

public class LyricSingularValueDecomposition extends SingularValueDecomposition
{
	public LyricSingularValueDecomposition(RealMatrix arg0)
	{
		super(checkMatrix(arg0));
	}
	
	private static RealMatrix checkMatrix(RealMatrix m)
	{
		for (int i = 0; i < m.getRowDimension(); i++)
		{
			for (int j = 0; j < m.getColumnDimension(); j++)
			{
				if (Double.isNaN(m.getEntry(i,j)) || Double.isInfinite(m.getEntry(i,j)))
				{
					throw new DimpleException("cannot do SVD on matrix that contains NaN or infinite");
				}
			}
		}
		return m;
	}
}
