/* * MaskFlipOperatorParser.java * * Copyright (c) 2002-2016 Alexei Drummond, Andrew Rambaut and Marc Suchard * * This file is part of BEAST. * See the NOTICE file distributed with this work for additional * information regarding copyright ownership and licensing. * * BEAST is free software; you can redistribute it and/or modify * it under the terms of the GNU Lesser General Public License as * published by the Free Software Foundation; either version 2 * of the License, or (at your option) any later version. * * BEAST is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * GNU Lesser General Public License for more details. * * You should have received a copy of the GNU Lesser General Public * License along with BEAST; if not, write to the * Free Software Foundation, Inc., 51 Franklin St, Fifth Floor, * Boston, MA 02110-1301 USA */ package dr.inferencexml.operators; import dr.inference.model.Parameter; import dr.inference.operators.MCMCOperator; import dr.inference.operators.MaskMoveOperator; import dr.xml.*; import java.util.ArrayList; import java.util.List; /** * */ public class MaskMoveOperatorParser extends AbstractXMLObjectParser { public static final String MASK_FLIP_OPERATOR = "maskMoveOperator"; public static final String CUT_POINT = "cutPoint"; public static final String CUT_MASK = "cutMask"; public static final String SELECT_BEFORE = "selectBefore"; public static final String SELECT_AFTER = "selectAfter"; // public static final String BEFORE_VALUE = "before"; // public static final String AFTER_VALUE = "after"; public String getParserName() { return MASK_FLIP_OPERATOR; } public Object parseXMLObject(XMLObject xo) throws XMLParseException { double weight = xo.getDoubleAttribute(MCMCOperator.WEIGHT); // Parameter parameter = (Parameter) xo.getChild(Parameter.class); List<Parameter> masks = new ArrayList<Parameter>(); for (int i = 0; i < xo.getChildCount(); ++i) { if (xo.getChild(i) instanceof Parameter) { masks.add((Parameter) xo.getChild(i)); } } Parameter cutPoint = (Parameter) xo.getElementFirstChild(CUT_POINT); // Parameter before = (Parameter) xo.getElementFirstChild(BEFORE_VALUE); // Parameter after = (Parameter) xo.getElementFirstChild(AFTER_VALUE); double[] before = xo.getChild(CUT_MASK).getDoubleArrayAttribute(SELECT_BEFORE); double[] after = xo.getChild(CUT_MASK).getDoubleArrayAttribute(SELECT_AFTER); int[] beforeList = new int[before.length]; for (int i = 0; i < before.length; ++i) { beforeList[i] = ((int) (before[i] -1.0 + 0.5)); // Switch to 0-index } int[] afterList = new int[after.length]; for (int i = 0; i < before.length; ++i) { afterList[i] = ((int) (after[i] -1.0 + 0.5)); } if (beforeList.length != afterList.length) { throw new XMLParseException("selectBefore length != selectAfter length"); } if (!MaskMoveOperator.checkMaskValues(masks, cutPoint, beforeList, afterList)) { String name = xo.hasId() ? xo.getId() : "null"; throw new XMLParseException("Bad initialization parameter values in " + name); } return new MaskMoveOperator(masks, cutPoint, beforeList, afterList, weight); } //************************************************************************ // AbstractXMLObjectParser implementation //************************************************************************ public String getParserDescription() { return "This element returns a mask-flip operator on a set of given parameters."; } public Class getReturnType() { return MaskMoveOperator.class; } public XMLSyntaxRule[] getSyntaxRules() { return rules; } private final XMLSyntaxRule[] rules = { AttributeRule.newDoubleRule(MCMCOperator.WEIGHT), new ElementRule(Parameter.class, 1, Integer.MAX_VALUE), new ElementRule(CUT_POINT, new XMLSyntaxRule[] { new ElementRule(Parameter.class), }), new ElementRule(CUT_MASK, new XMLSyntaxRule[] { AttributeRule.newDoubleArrayRule(SELECT_BEFORE), AttributeRule.newDoubleArrayRule(SELECT_AFTER), }) // new ElementRule(BEFORE_VALUE, new XMLSyntaxRule[] { // new ElementRule(Parameter.class), // }), // new ElementRule(AFTER_VALUE, new XMLSyntaxRule[] { // new ElementRule(Parameter.class), // }), }; }