dep-selector-libgecode 1.0.0.alpha.0

data/ext/libgecode3/vendor/gecode-3.7.3/contribs/qecode/QCOPPlus.hh

@@ -0,0 +1,166 @@
+/****   , [ QCOPPlus.hh ], 
+ Copyright (c) 2009 Universite d'Orleans - Jeremie Vautard 
+ 
+ Permission is hereby granted, free of charge, to any person obtaining a copy
+ of this software and associated documentation files (the "Software"), to deal
+ in the Software without restriction, including without limitation the rights
+ to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+ copies of the Software, and to permit persons to whom the Software is
+ furnished to do so, subject to the following conditions:
+ 
+ The above copyright notice and this permission notice shall be included in
+ all copies or substantial portions of the Software.
+ 
+ THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+ IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+ AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+ LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+ OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
+ THE SOFTWARE.
+ *************************************************************************/
+
+#ifndef __QECODE_QCOP__
+#define __QECODE_QCOP__
+
+#include "myspace.hh"
+#include "vartype.hh"
+#include "OptVar.hh"
+#include <vector>
+
+
+using namespace std;
+using namespace Gecode;
+class QECODE_VTABLE_EXPORT Qcop {
+private:
+	struct Opts {
+		vector<OptVar*> vars; // Several Optvars in universal scopes, one in existential ones
+		int opt_type;        // Used in existential scopes : 0 for ANY, 1 for MIN, 2 for MAX.
+	};
+	int* nvar;
+	int n;
+	int nbSpaces;
+//	void** v;
+	VarType* type_of_v;
+	
+	bool* Quantifiers;
+	MySpace** rules;
+	MySpace* goal;
+	Opts* optim;
+	int* nbVarBySpace;
+	int* whichSpaceOwns;
+	bool* varInitialised;
+	
+	int currentDeclareSpace;
+	Qcop() {}
+public:
+	forceinline QECODE_EXPORT int nv() {return n;}
+	
+	/** \brief  Default constructor for a restricted-quantified space
+	 *
+	 *  This is the first step for building a QCSP+/QCOP+ problem. The prefix is defined here.
+	 *  @param ns The number of scopes in the prefix.
+	 *  @param quant Array of booleans which contains the quantifier of every scope (true for universal, false for existential).
+	 *  @param nv Array of integer which contains the number of variables by scope.
+	 */    
+	QECODE_EXPORT Qcop(int ns, bool* quant,int* nv);
+	QECODE_EXPORT ~Qcop();
+	
+	/** \brief Defines an integer variable in the quantified space
+	 *
+	 *  Defines an integer variable in the quantifies space using a fully declared domain.
+	 *  @param var Number of the variable to be  defined.
+	 *  @param dom The initial domain of the variable.
+	 */    
+	QECODE_EXPORT void QIntVar(int var,int min,int max);
+	
+	/** \brief Defines an integer variable in the quantified space
+	 *
+	 *  Defines an integer variable in the quantifies space using a fully declared domain.
+	 *  @param var Number of the variable to be  defined.
+	 *  @param dom The initial domain of the variable.
+	 */
+	QECODE_EXPORT void QIntVar(int var,IntSet dom);
+	
+	/** \brief Defines a boolean variable in the quantified space
+	 *
+	 *  Defines a boolean variable with a full initial domain in the quantified space.
+	 *  @param var Number of the variable to be defined.
+	 */
+	QECODE_EXPORT void QBoolVar(int var);
+	
+	/** \brief Returns the current declarating space
+	 *
+	 *  Return the space we are currently declarating constraints in. nextScope() is used to set the nextspace as the current one.
+	 */    
+	QECODE_EXPORT MySpace* space();
+	
+	/** \brief Returns an integer variable from the space we are currently declaring
+	 *
+	 *  Returns an integer variable from the cpace we are currently declaring. Will abort if the variable is not integer.
+	 *  @param n The number of the variable to return.
+	 */
+	QECODE_EXPORT IntVar var(int n);
+	
+	/** \brief Returns a boolean variable from the space we are currently declaring
+	 *
+	 * Returns a boolean variable from the space we are currently declaring. Will abort if the variable is not boolean.
+	 *  @param n The number of the variable to return.
+	 */
+	QECODE_EXPORT BoolVar bvar(int n);
+	
+	/** \brief Switch to the next space for constraint declaration
+	 *
+	 *   Switches to the next space for constraint declaration. var, bvar and space methods will now refer to the next rule space (or to the goal space if the current space was the last rule one).
+	 *  Returns the new space number, or -1 if it was called while there was no next space to declare constraints in.
+	 */
+	QECODE_EXPORT int nextScope();
+	
+	/** \brief Finalizes the object construction
+	 *
+	 * Finalizes the restricted-quantified space construction. Must be invoked once all the variable and all the constraints have been declared, and before the search is performed on this space.
+	 * calling this method is not mandatory, although it is recommanded to besure that the problem have been well built.
+	 * For the existential scopes on which an optimization condition have not been defined yet, this method will post a "Any" optimization condition (i.e. do not optimize).
+	 */
+	QECODE_EXPORT void makeStructure();
+	
+	/** \brief returns an aggregate of the problem
+	 *
+	 * Creates an aggregate at a given universal scope. This aggregate is an optimization variable that an existential scope can use for optimizing.
+	 * @param scope the scope where this aggregate will be defined. Must be an unversal scope
+	 * @param opt The optimization variable we want to aggregate. There must not be any universal scope between an agregate and thisoptimization variable.
+	 * @param agg The aggregator function that will be used for this aggregate.
+	 */
+	QECODE_EXPORT OptVar* getAggregate(int scope, OptVar* opt, Aggregator* agg);
+	
+	/** returns an existential optimization variable
+	 *
+	 * Creates an optimization variable that will have the value of an existential variable defined in the problem. 
+	 * @param var the number of the existential variable that must be considered
+	 */
+	QECODE_EXPORT OptVar* getExistential(int var);
+	
+	/** set an optimization condition on an existential scope
+	 *
+	 * set an optimization condition on a given existential scope of the problem. An optimizaiton condition is composed of an optimization variable (aggregate or existential variable), and of an
+	 * optimization type. 
+	 * @param scope the scope on which the optimization condition is posted. Must be existential.
+	 * @param optType the optimization type of the condision we post. 0 is for ANY, 1 is for MIN, 2 is for MAX.
+	 * @param var the optimization variable to be minimized/maximized
+	 */
+	QECODE_EXPORT void optimize(int scope, int optType, OptVar* var);
+	//    QECODE_EXPORT void print();
+	
+	QECODE_EXPORT forceinline bool qt_of_var(int v) {	return Quantifiers[whichSpaceOwns[v]];} ///< returns uantifier of variable 'v'
+	QECODE_EXPORT forceinline bool quantification(int scope) {return Quantifiers[scope];} ///< returns quantifier of scope 'scope'
+	QECODE_EXPORT int spaces(); ///< returns the number of scopes of the problem
+	QECODE_EXPORT forceinline int nbVarInScope(int scope) {return nbVarBySpace[scope];}///< returns the number of variables present in scope 'scope'
+	QECODE_EXPORT MySpace* getSpace(int scope);///< returns a copy of the Gecode::Space corresponding to the 'scope'-th restricted quantifier of the problem
+	QECODE_EXPORT MySpace* getGoal();
+	QECODE_EXPORT int getOptType(int scope); ///< returns the optimization type of the 'scope'-th scope of the problem
+	QECODE_EXPORT OptVar* getOptVar(int scope);///< returns the optimization variable of the 'scope'-th scope of the problem
+	QECODE_EXPORT Qcop* clone(); ///< makes a copy of the quantified problem
+	QECODE_EXPORT forceinline int whoOwns(int var) {return whichSpaceOwns[var];}
+};
+
+#endif

data/ext/libgecode3/vendor/gecode-3.7.3/contribs/qecode/QCSPPlusUnblockable.cc

@@ -0,0 +1,239 @@
+/****   , [ QCSPPlusUnblockable.cc ], 
+Copyright (c) 2009 Universite d'Orleans - Jeremie Vautard 
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in
+all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
+THE SOFTWARE.
+ *************************************************************************/
+
+#include "QCSPPlusUnblockable.hh"
+
+QcspUnblockable::QcspUnblockable(int ns,bool* quant,int* nv) {
+    n=0;
+    for (int i=0;i<ns;i++) {
+        n += nv[i];
+    }
+    nbSpaces=ns;
+    //    cout<<"QCSPU construction : "<<n<<" vriables"<<endl;
+    v=new void*[n];
+    type_of_v = new VarType[n];
+    Quantifiers = quant;
+    whichSpaceOwns=new int[n];
+    nbVarBySpace=new int[nbSpaces];
+    //    cout<<"Before rule creation : "<<arul<<endl;
+    arul=new MySpace(n);
+    //    cout<<"after rule creation : "<<arul<<endl;
+    nbVarBySpace[0]=nv[0];
+    for (int i=1;i<n;i++) {
+        nbVarBySpace[i]=nbVarBySpace[i-1]+nv[i];
+    }
+
+        arul=new MySpace(n);
+    goal = new MySpace(n);
+
+    
+    for(unsigned int i=0;i<n;i++) {
+        int lespace=0;
+        while (nbVarBySpace[lespace]<=i) lespace++;
+        whichSpaceOwns[i]=lespace;
+    }
+
+    varInitialised=new bool[n];
+    for (unsigned int i=0;i<n;i++) varInitialised[i]=false;
+    currentDeclareSpace=0;
+    
+    vars = new vector<int>[nbSpaces];
+    bvars = new vector<int>[nbSpaces];
+}
+
+
+QcspUnblockable::~QcspUnblockable() {
+    delete arul;
+    delete goal;
+}
+
+int QcspUnblockable::spaces() {
+    return nbSpaces;
+}
+
+void QcspUnblockable::QIntVar(int var,int min,int max) {
+    if (varInitialised[var]) {
+        cout<<"Variable "<<var<<"  Already created !!"<<endl;
+        abort();
+    }
+    //    cout<<"Qintvar : arul = "<<arul<<endl;
+    
+    arul->v[var] = new IntVar(*arul,min,max);
+    arul->type_of_v[var] = VTYPE_INT;
+
+    goal->v[var] = new IntVar(*goal,min,max);
+    goal->type_of_v[var] = VTYPE_INT;
+    varInitialised[var]=true;
+    type_of_v[var]=VTYPE_INT;
+}
+
+void QcspUnblockable::QIntVar(int var,IntSet dom) {
+    if (varInitialised[var]) {
+        cout<<"Variable "<<var<<"  Already created !!"<<endl;
+        abort();
+    }
+    
+    arul->v[var] = new IntVar(*arul,dom);
+    arul->type_of_v[var] = VTYPE_INT;
+    goal->v[var] = new IntVar(*goal,dom);
+    goal->type_of_v[var] = VTYPE_INT;
+    varInitialised[var]=true;
+    type_of_v[var]=VTYPE_INT;
+}
+
+
+void QcspUnblockable::QBoolVar(int var) {
+    if (varInitialised[var]) {
+        cout<<"Variable "<<var<<" Already created !!"<<endl;
+        abort();
+    }
+    
+    arul->v[var] = new BoolVar(*arul,0,1);
+    arul->type_of_v[var]=VTYPE_BOOL;
+    goal->v[var] = new BoolVar(*goal,0,1);
+    goal->type_of_v[var]=VTYPE_BOOL;
+    varInitialised[var]=true;
+    type_of_v[var]=VTYPE_BOOL;
+}
+
+MySpace* QcspUnblockable::space() {
+    if (currentDeclareSpace<nbSpaces) {
+      //cout<<"Return space arul"<<endl; cout.flush();
+        return arul;
+}
+    if (currentDeclareSpace==nbSpaces) {
+      //cout<<"Return space Goal"<<endl; cout.flush();
+        return goal;
+}
+cout<<"Return null in space()"<<endl;
+    return NULL;
+}
+
+
+IntVar QcspUnblockable::var(int n) {
+    if (!varInitialised[n]) {
+        cout<<"Variable "<<n<<" not initialized !"<<endl;
+        abort();
+    }
+    if (type_of_v[n] != VTYPE_INT) {
+        cout<<"Variable "<<n<<" is not INT"<<endl;
+        abort();
+    }
+    return *(static_cast<IntVar*>(space()->v[n]));
+}
+
+BoolVar QcspUnblockable::bvar(int n) {
+    if (!varInitialised[n]) {
+        cout<<"Variable "<<n<<" not initialized !"<<endl;
+        abort();
+    }
+    if (type_of_v[n] != VTYPE_BOOL){
+        cout<<"Variable "<<n<<" is not BOOL"<<endl;
+        abort();
+    }
+    return *(static_cast<BoolVar*>(space()->v[n]));
+}
+
+int QcspUnblockable::nextScope() {
+    if (currentDeclareSpace == -1) return -1;
+    currentDeclareSpace++;
+    if (currentDeclareSpace>nbSpaces) return -1;
+    return currentDeclareSpace;
+}
+
+void QcspUnblockable::makeStructure() {
+    for (unsigned int i=0;i<n;i++) {
+        if (varInitialised[i] == false) {
+            cout<<"Can't make structure : variable "<<i<<" not initialised"<<endl;
+            abort();
+        }
+    }
+    for (unsigned int i=0;i<nbSpaces;i++) {
+        unsigned int nbint=0;
+        unsigned int nbbool=0;
+        for (unsigned int j=0;j<nbVarBySpace[i];j++) {
+            if (type_of_v[j] == VTYPE_INT) 
+                nbint++;
+            else
+                nbbool++;
+        }
+        nbint=0;nbbool=0;
+        for (unsigned int j=0;j<nbVarBySpace[i];j++) {
+            if (type_of_v[j] == VTYPE_INT) {
+                (vars[i]).push_back(j);
+            }
+            else {
+                (bvars[i]).push_back(j);
+            }
+        }
+    }
+}
+
+forceinline bool QcspUnblockable::qt_of_var(int v) {
+    return Quantifiers[whichSpaceOwns[v]];
+}
+
+MySpace* QcspUnblockable::getSpace(int scope) {
+    if (scope<0 || scope>nbSpaces) {
+        cout<<"I return NULL coz of bad scope value (<0)"<<endl;
+        return NULL;
+    }
+    if (scope==nbSpaces) {
+        if (goal->status() == SS_FAILED) {
+            cout<<"I return NULL coz goal is failed"<<endl;
+            return NULL;
+        }
+        //        cout<<"I return the goal"<<endl;
+        return static_cast<MySpace*>(goal->clone());
+    }
+    if (arul->status() == SS_FAILED) {
+        cout<<"I return NULL coz scope "<<scope<<" is failed"<<endl;
+        return NULL;
+    }
+    //    cout<<"I return the rule "<<scope<<endl;
+    MySpace* ret= (static_cast<MySpace*>(arul->clone()));
+    IntVarArgs iva(vars[scope].size());
+    BoolVarArgs bva(bvars[scope].size());
+    //    cout << "sizes : " <<iva.size() << " " << bva.size()<<endl;
+    for (int i=0;i<iva.size();i++) {
+        int idx = (vars[scope])[i];
+        iva[i]=*( static_cast<IntVar*>(ret->v[idx]) );
+    }
+    for (int i=0;i<bva.size();i++) {
+        int idx = (bvars[scope])[i];
+        bva[i]=*( static_cast<BoolVar*>(ret->v[idx]) );
+    }
+    br->branch(ret,iva,bva);
+    
+    return ret;
+}
+
+
+MySpace* QcspUnblockable::getGoal() {
+    if (goal->status() == SS_FAILED) return NULL;
+    return static_cast<MySpace*>(goal->clone());
+}
+
+
+void QcspUnblockable::branch(UnblockableBranching* b) {
+    br=b;
+}

data/ext/libgecode3/vendor/gecode-3.7.3/contribs/qecode/QCSPPlusUnblockable.hh

@@ -0,0 +1,140 @@
+/****   , [ QCSPPlusUnblockable.hh ], 
+Copyright (c) 2009 Universite d'Orleans - Jeremie Vautard 
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in
+all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
+THE SOFTWARE.
+ *************************************************************************/
+
+#ifndef __QECODE_UNBLOCKABLE__
+#define __QECODE_UNBLOCKABLE__
+
+#include <iostream>
+#include<vector>
+#include "qecode.hh"
+#include "vartype.hh"
+#include "UnblockableBranching.hh"
+#include "myspace.hh"
+
+
+using namespace std;
+using namespace Gecode;
+
+class QECODE_VTABLE_EXPORT QcspUnblockable {
+private : 
+    
+    int n;
+    int nbSpaces;
+    void** v;
+    VarType* type_of_v;
+    int* whichSpaceOwns;
+    unsigned int currentDeclareSpace;
+    bool* Quantifiers;
+    MySpace* goal;
+    int* nbVarBySpace;
+    bool* varInitialised;
+    UnblockableBranching* br;
+    vector<int>* vars;
+    vector<int>* bvars;
+    MySpace* arul;
+    
+public:
+        QECODE_EXPORT int nv() {return n;}
+    
+    /** \brief  Default constructor for a restricted-quantified space
+        *
+        *  This is the first step for building a QCSP+/QCOP+ problem. The prefix is defined here.
+        *  @param ns The number of scopes in the prefix.
+        *  @param quant Array of booleans which contains the quantifier of every scope (true for universal, false for existential).
+        *  @param nv Array of integer which contains the number of variables by scope.
+        */    
+    QECODE_EXPORT QcspUnblockable(int ns, bool* quant,int* nv);
+    QECODE_EXPORT ~QcspUnblockable();
+    
+    /** \brief Defines an integer variable in the quantified space
+        *
+        *  Defines an integer variable in the quantifies space using a fully declared domain.
+        *  @param var Number of the variable to be  defined.
+        *  @param dom The initial domain of the variable.
+        */    
+    QECODE_EXPORT void QIntVar(int var,int min,int max);
+    
+    /** \brief Defines an integer variable in the quantified space
+        *
+        *  Defines an integer variable in the quantifies space using a fully declared domain.
+        *  @param var Number of the variable to be  defined.
+        *  @param dom The initial domain of the variable.
+        */
+    QECODE_EXPORT void QIntVar(int var,IntSet dom);
+    
+    /** \brief Defines a boolean variable in the quantified space
+        *
+        *  Defines a boolean variable with a full initial domain in the quantified space.
+        *  @param var Number of the variable to be defined.
+        */
+    QECODE_EXPORT void QBoolVar(int var);
+    
+    /** \brief Returns the current declarating space
+        *
+        *  Return the space we are currently declarating constraints in. nextScope() is used to set the nextspace as the current one.
+        */    
+    QECODE_EXPORT MySpace* space();
+    
+    /** \brief Returns an integer variable from the space we are currently declaring
+        *
+        *  Returns an integer variable from the cpace we are currently declaring. Will abort if the variable is not integer.
+        *  @param n The number of the variable to return.
+        */
+    QECODE_EXPORT IntVar var(int n);
+    
+    /** \brief Returns a boolean variable from the space we are currently declaring
+        *
+        * Returns a boolean variable from the space we are currently declaring. Will abort if the variable is not boolean.
+        *  @param n The number of the variable to return.
+        */
+    QECODE_EXPORT BoolVar bvar(int n);
+    
+    /** \brief Switch to the next space for constraint declaration
+        *
+        *   Switches to the next space for constraint declaration. var, bvar and space methods will now refer to the next rule space (or to the goal space if the current space was the last rule one).
+        *  Returns the new space number, or -1 if it was called while there was no next space to declare constraints in.
+        */
+    QECODE_EXPORT int nextScope();
+    
+    /** \brief Finalizes the object construction
+        *
+        * Finalizes the restricted-quantified space construction. Must be invoked once all the variable and all the constraints have been declared, and before the search is performed on this space.
+        * calling this method is not mandatory, although it is recommanded to besure that the problem have been well built.
+        * For the existential scopes on which an optimization condition have not been defined yet, this method will post a "Any" optimization condition (i.e. do not optimize).
+        */
+    QECODE_EXPORT void makeStructure();
+    
+    QECODE_EXPORT bool qt_of_var(int v); ///< returns uantifier of variable 'v'
+    QECODE_EXPORT bool quantification(int scope) {return Quantifiers[scope];} ///< returns quantifier of scope 'scope'
+    QECODE_EXPORT int spaces(); ///< returns the number of scopes of the problem
+    QECODE_EXPORT int nbVarInScope(int scope) {return nbVarBySpace[scope];}///< returns the number of variables present in scope 'scope'
+    QECODE_EXPORT MySpace* getSpace(int scope);///< returns a copy of the Gecode::Space corresponding to the scope-th restricted quantifier of the problem
+    QECODE_EXPORT MySpace* getGoal();
+        
+        /** \brief sets the branching heuristic
+            * 
+            * Declares the branching heuristic that will be used to solve this problem. 
+        */
+    QECODE_EXPORT void branch(UnblockableBranching* b); 
+};
+
+#endif