class SimpleSpace : public Gecode::Space
{
public:
	Gecode::FloatVarArray mFloatVariables;
        // The float channel variables are channeled to the int variables and used in float expressions.
	Gecode::FloatVarArray mIntVariableChannels;
        // Int variables serve as discrete option selectors.
	Gecode::IntVarArray mIntVariables;

public:
	SimpleSpace(
		const size_t floatVariableCount,
		const size_t intVariableCount) noexcept
		:
		mFloatVariables(*this, floatVariableCount, 0.f, 10000.f),
		mIntVariableChannels(*this, intVariableCount, 0.f, 10.f),
		mIntVariables(*this, intVariableCount, 0, 10)
	{}

	SimpleSpace(SimpleSpace &other)
		:
		Space(other)
	{
		mFloatVariables.update(*this, other.mFloatVariables);
		mIntVariableChannels.update(*this, other.mIntVariableChannels);
		mIntVariables.update(*this, other.mIntVariables);
	}

	// Inherited via Space
	virtual Space *copy(void) override
	{
		return new SimpleSpace(*this);
	}
};

// Initialize the space with 2 float variables and 1 int variable.
SimpleSpace home(2, 1);

// Bind the int variable to the float channel variable.
// The int variable is used like a bool variable in this context, but can be used to express more options.
Gecode::channel(home, home.mIntVariables[0], home.mIntVariableChannels[0]);
Gecode::dom(home, home.mIntVariableChannels[0], 0.f, 1.f);
Gecode::dom(home, home.mIntVariables[0], 0, 1);

// Declare temporary float variables.
// The limits and potentially other relations for the float variables are given,
// but I want the option to switch their values and constraints.
Gecode::FloatVarArray floatVariables(home, 2, 0.f, 10000.f);
Gecode::dom(home, floatVariables[0], 2.f, 2.f);
Gecode::dom(home, floatVariables[1], 3.f, 4.f);

// Limit domain of actual variables.
Gecode::dom(home, home.mFloatVariables[0], 2.f, 4.f);
Gecode::dom(home, home.mFloatVariables[1], 2.f, 4.f);

// Flip accordingly.
// I'm using two multivariate quadratic expressions here to determine the value of the actual variables.
Gecode::FloatVar flip = home.mIntVariableChannels[0];
const auto flipExpr = flip * (floatVariables[1] - floatVariables[0]);
// This relation goes through correctly. If I make the domain of floatVariables[0] not tight in this case,
// it will also throw the exception.
Gecode::rel(home, home.mFloatVariables[0] == floatVariables[0] + flipExpr);
// This one errors out as described earlier.
Gecode::rel(home, home.mFloatVariables[1] == floatVariables[1] - flipExpr); // EXCEPTION

// Flip to be determined later...

set(CMAKE_DEBUG_POSTFIX "d" CACHE PATH "")

project(GECODE)