package org.matheclipse.core.reflection.system;

import org.matheclipse.core.basic.Config;
import org.matheclipse.core.eval.exception.Validate;
import org.matheclipse.core.eval.interfaces.AbstractFunctionEvaluator;
import org.matheclipse.core.expression.F;
import org.matheclipse.core.generic.Functors;
import org.matheclipse.core.interfaces.IAST;
import org.matheclipse.core.interfaces.IComplexNum;
import org.matheclipse.core.interfaces.IExpr;
import org.matheclipse.core.interfaces.INum;
import org.matheclipse.core.interfaces.ISymbol;

public class Chop extends AbstractFunctionEvaluator {
	public final double DEFAULT_CHOP_DELTA = 10E-10;

	public Chop() {
	}

	@Override
	public IExpr evaluate(final IAST ast) {
		Validate.checkRange(ast, 2, 3);

		IExpr arg1 = ast.arg1();
		double delta = DEFAULT_CHOP_DELTA;
		if (ast.size() == 3 && ast.arg2() instanceof INum) {
			delta = ((INum) ast.arg2()).getRealPart();
		}
		try {
			arg1 = F.eval(arg1);
			if (arg1 instanceof INum) {
				if (F.isZero(((INum) arg1).getRealPart(), delta)) {
					return F.C0;
				}
				return arg1;
			} else if (arg1 instanceof IComplexNum) {
				if (F.isZero(((IComplexNum) arg1).getRealPart(), delta)) {
					if (F.isZero(((IComplexNum) arg1).getImaginaryPart(), delta)) {
						return F.C0;
					}
					return F.complexNum(0.0, ((IComplexNum) arg1).getImaginaryPart());
				}
				if (F.isZero(((IComplexNum) arg1).getImaginaryPart(), delta)) {
					return F.complexNum(((IComplexNum) arg1).getRealPart(), 0.0);
				}
				return arg1;
			} else if (arg1.isAST()) {
				IAST list = (IAST)arg1;
				// Chop[{a,b,c}] -> {Chop[a],Chop[b],Chop[c]}
				return list.map(Functors.replace1st(F.Chop(F.Null)));
			}
			return arg1;
		} catch (Exception e) {
			if (Config.SHOW_STACKTRACE) {
				e.printStackTrace();
			}
		}

		return null;
	}

	public void setUp(final ISymbol symbol) {
		symbol.setAttributes(ISymbol.NHOLDALL | ISymbol.LISTABLE);
	}
}

 

Ads help maintain this website.