All Downloads are FREE. Search and download functionalities are using the official Maven repository.

org.chocosolver.solver.variables.impl.RealVarImpl Maven / Gradle / Ivy

There is a newer version: 4.10.17
Show newest version
/*
 * This file is part of choco-solver, http://choco-solver.org/
 *
 * Copyright (c) 2020, IMT Atlantique. All rights reserved.
 *
 * Licensed under the BSD 4-clause license.
 *
 * See LICENSE file in the project root for full license information.
 */
package org.chocosolver.solver.variables.impl;

import org.chocosolver.memory.IStateDouble;
import org.chocosolver.solver.ICause;
import org.chocosolver.solver.Model;
import org.chocosolver.solver.exception.ContradictionException;
import org.chocosolver.solver.exception.SolverException;
import org.chocosolver.solver.expression.continuous.arithmetic.CArExpression;
import org.chocosolver.solver.variables.RealVar;
import org.chocosolver.solver.variables.delta.NoDelta;
import org.chocosolver.solver.variables.events.RealEventType;
import org.chocosolver.solver.variables.impl.scheduler.RealEvtScheduler;
import org.chocosolver.util.iterators.EvtScheduler;
import org.chocosolver.util.objects.RealInterval;

import java.util.List;
import java.util.TreeSet;

/**
 * An implementation of RealVar, variable for continuous constraints (solved using IBEX).
 * 
* * @author Charles Prud'homme * @since 18/07/12 */ public class RealVarImpl extends AbstractVariable implements RealVar { private IStateDouble LB, UB; private double precision; public RealVarImpl(String name, double lb, double ub, double precision, Model model) { super(name, model); this.LB = model.getEnvironment().makeFloat(lb); this.UB = model.getEnvironment().makeFloat(ub); this.precision = precision; } @Override public double getPrecision() { return precision; } @Override public void silentlyAssign(RealInterval bounds) { silentlyAssign(bounds.getLB(), bounds.getUB()); } @Override public void silentlyAssign(double lb, double ub) { this.LB.set(lb); this.UB.set(ub); } @Override public double getLB() { return LB.get(); } @Override public double getUB() { return UB.get(); } @Override public void intersect(double l, double u, ICause cause) throws ContradictionException { updateBounds(l, u, cause); } @Override public boolean updateLowerBound(double value, ICause cause) throws ContradictionException { assert cause != null; // TODO ICause antipromo = cause; double old = this.getLB(); if (old < value) { if (this.getUB() < value) { // TODO solver.getExplainer().updateLowerBound(this, old, value, antipromo); this.contradiction(cause, MSG_LOW); } else { LB.set(value); this.notifyPropagators(RealEventType.INCLOW, cause); // TODO solver.getExplainer().updateLowerBound(this, old, value, antipromo); return true; } } return false; } @Override public boolean updateUpperBound(double value, ICause cause) throws ContradictionException { assert cause != null; // TODO ICause antipromo = cause; double old = this.getUB(); if (old > value) { if (this.getLB() > value) { // TODO solver.getExplainer().updateUpperBound(this, old, value, antipromo); this.contradiction(cause, MSG_UPP); } else { UB.set(value); this.notifyPropagators(RealEventType.DECUPP, cause); // TODO solver.getExplainer().updateUpperBound(this, old, value, antipromo); return true; } } return false; } @Override public boolean updateBounds(double lowerbound, double upperbound, ICause cause) throws ContradictionException { assert cause != null; // TODO ICause antipromo = cause; double oldlb = this.getLB(); double oldub = this.getUB(); if (oldlb < lowerbound || oldub > upperbound) { if (oldub < lowerbound || oldlb > upperbound) { // TODO solver.getExplainer()... this.contradiction(cause, MSG_BOUND); } else { RealEventType e = RealEventType.VOID; if (oldlb < lowerbound) { LB.set(lowerbound); e = RealEventType.INCLOW; } if (oldub > upperbound) { UB.set(upperbound); e = (e == RealEventType.INCLOW) ? RealEventType.BOUND : RealEventType.DECUPP; } this.notifyPropagators(e, cause); // TODO solver.getExplainer()... return true; } } return false; } @Override public boolean isInstantiated() { double lb = LB.get(); double ub = UB.get(); return ub - lb < precision || nextValue(lb) >= ub; } private double nextValue(double x) { if (x < 0) { return Double.longBitsToDouble(Double.doubleToLongBits(x) - 1); } else if (x == 0) { return Double.longBitsToDouble(1L); } else if (x < Double.POSITIVE_INFINITY) { return Double.longBitsToDouble(Double.doubleToLongBits(x) + 1); } else { return x; // nextValue(infty) = infty } } @Override public NoDelta getDelta() { return NoDelta.singleton; } @Override public void createDelta() { throw new SolverException("Unable to create delta for RealVar!"); } @Override public int getTypeAndKind() { return VAR | REAL; } @Override protected EvtScheduler createScheduler() { return new RealEvtScheduler(); } @Override public String toString() { return String.format("%s = [%.16f .. %.16f]", name, getLB(), getUB()); } @Override public void tighten() { } @Override public void project(ICause cause) throws ContradictionException { } @Override public void collectVariables(TreeSet set) { set.add(this); } @Override public void subExps(List list) { list.add(this); } @Override public boolean isolate(RealVar var, List wx, List wox) { return var == this; } @Override public void init() { // void } }




© 2015 - 2024 Weber Informatics LLC | Privacy Policy