org.lsmp.djep.jama.Inverse Maven / Gradle / Ivy

Go to download

Show more of this group Show more artifacts with this name
Show all versions of jep Show documentation

JEP is a Java library for parsing and evaluating mathematical expressions.

The newest version!

/* @author rich
 * Created on 15-Feb-2005
 *
 * See LICENSE.txt for license information.
 */
package org.lsmp.djep.jama;

import java.util.Stack;

import org.nfunk.jep.ParseException;
import org.nfunk.jep.function.PostfixMathCommand;
import org.lsmp.djep.vectorJep.Dimensions;
import org.lsmp.djep.vectorJep.values.*;
import org.lsmp.djep.vectorJep.function.*;

/**
 * Find the inverses of a matrix.
 * Serves a wrapper around the Jama linear algebra function.
 * @see http://math.nist.gov/javanumerics/jama/
 * 
 * @author Rich Morris
 * Created on 15-Feb-2005
 */
public class Inverse extends PostfixMathCommand implements UnaryOperatorI 
{
	public Inverse()
	{
		this.numberOfParameters = 1;
	}


	public void run(Stack s) throws ParseException
	{
		Object o = s.pop();
		if(!(o instanceof Matrix))
			throw new ParseException("inverse: can only be applied to a matrix");
		Jama.Matrix m = JamaUtil.toJama((Matrix) o);
		Jama.Matrix inv = m.inverse();
		Matrix res = JamaUtil.fromJama(inv);
		s.push(res);
	}

	public Dimensions calcDim(Dimensions ldim)
	{
		int rows = ldim.getFirstDim();
		int cols = ldim.getLastDim();
		if(rows < cols) return Dimensions.valueOf(rows,rows);
		else if(rows > cols) return Dimensions.valueOf(cols,cols);
		else return ldim;
	}

	public MatrixValueI calcValue(MatrixValueI res, MatrixValueI lhs)
		throws ParseException
	{
		if(!(lhs instanceof Matrix))
			throw new ParseException("inverse: can only be applied to a matrix");
		if(!(res instanceof Matrix))
			throw new ParseException("inverse: result should be a matrix");
		Jama.Matrix m = JamaUtil.toJama((Matrix) lhs);
		Jama.Matrix inv = m.inverse();
		JamaUtil.fromJama(inv,(Matrix) res);
		return res;
	}

}