/*
* CombinedOperatorSchedule.java
*
* Copyright (c) 2002-2015 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.inference.operators;
import dr.math.MathUtils;
import java.util.List;
import java.util.Vector;
/**
* Package: CombinedOperatorSchedule
* Description:
* <p/>
* <p/>
* Created by
*
* @author Alexander V. Alekseyenko (alexander.alekseyenko@gmail.com)
* Date: Dec 23, 2009
* Time: 10:17:28 AM
*/
public class CombinedOperatorSchedule implements OperatorSchedule {
List<OperatorSchedule> operatorSchedules = null;
public CombinedOperatorSchedule() {
operatorSchedules = new Vector<OperatorSchedule>();
}
public void addOperatorSchedule(OperatorSchedule os) {
operatorSchedules.add(os);
}
public int getScheduleCount() {
return operatorSchedules.size();
}
public int getNextOperatorIndex() {
final int v = MathUtils.nextInt(operatorSchedules.size());
int offset = 0;
for (int i = 0; i < v; ++i) {
offset += operatorSchedules.get(i).getOperatorCount();
}
return offset + operatorSchedules.get(v).getNextOperatorIndex();
}
public void reset() {
for (OperatorSchedule os : operatorSchedules) {
for (int i = 0; i < os.getOperatorCount(); ++i) {
os.getOperator(i).reset();
}
}
}
public int getOperatorCount() {
int operatorCount = 0;
for (OperatorSchedule os : operatorSchedules) {
operatorCount += os.getOperatorCount();
}
return operatorCount;
}
public MCMCOperator getOperator(int index) {
for (OperatorSchedule os : operatorSchedules) {
int opCount = os.getOperatorCount();
if (index < opCount) {
return os.getOperator(index);
} else {
index -= opCount;
}
}
// if we reach here the index must be out of bounds return null
return null;
}
public void addOperator(MCMCOperator op) {
System.err.println("ERROR: addOperator() should not be called on CombinedOperatorSchedule, " +
"use individual operator schedules to add operators. This call will have no affect!");
}
public void addOperators(List<MCMCOperator> v) {
System.err.println("ERROR: addOperators() should not be called on CombinedOperatorSchedule, " +
"use individual operator schedules to add operators. This call will have no affect!");
}
public void operatorsHasBeenUpdated() {
for (OperatorSchedule os : operatorSchedules) {
os.operatorsHasBeenUpdated();
}
}
public double getOptimizationTransform(double d) {
if (operatorSchedules.size() > 0)
return operatorSchedules.get(0).getOptimizationTransform(d);
else
return 0;
}
public long getMinimumAcceptAndRejectCount() {
long minCount = Long.MAX_VALUE;
for (OperatorSchedule os : operatorSchedules) {
if (os.getMinimumAcceptAndRejectCount() < minCount) {
minCount = os.getMinimumAcceptAndRejectCount();
}
}
return minCount;
}
}