gurobipy 12.0.2__cp310-cp310-win_amd64.whl → 13.0.0b1__cp310-cp310-win_amd64.whl

gurobipy/__init__.py

@@ -1,4 +1,55 @@
-__version__ = "12.0.2"
+"""
+Python interface to the Gurobi Optimizer
+
+The Gurobi Optimizer is a mathematical optimization software library for
+solving mixed-integer linear, quadratic, and nonlinear optimization problems.
+
+This package comes with a trial license that allows you to solve problems of
+limited size. As a student or staff member of an academic institution, you
+qualify for a free, full product license. For more information visit
+https://www.gurobi.com/academia/academic-program-and-licenses/
+
+For a commercial evaluation, you can request an evaluation license by
+visiting https://www.gurobi.com/free-trial/.
+
+Other useful resources to get started:
+
+- Gurobi Documentation: https://docs.gurobi.com/
+- Gurobi Community: https://support.gurobi.com/hc/en-us/community/topics/
+
+A simple example, formulating and solving a mixed-integer linear program:
+
+#  maximize
+#        x +   y + 2 z
+#  subject to
+#        x + 2 y + 3 z <= 4
+#        x +   y       >= 1
+#        x, y, z binary
+
+import gurobipy as gp
+
+with gp.Env() as env, gp.Model(env=env) as model:
+
+    # Create variables
+    x = model.addVar(vtype='B', name="x")
+    y = model.addVar(vtype='B', name="y")
+    z = model.addVar(vtype='B', name="z")
+
+    # Set objective function
+    model.setObjective(x + y + 2 * z, gp.GRB.MAXIMIZE)
+
+    # Add constraints
+    model.addConstr(x + 2 * y + 3 * z <= 4)
+    model.addConstr(x + y >= 1)
+
+    # Solve it!
+    model.optimize()
+
+    print(f"Optimal objective value: {model.ObjVal}")
+    print(f"Solution values: x={x.X}, y={y.X}, z={z.X}")
+"""
+
+__version__ = "13.0.0b1"
 
 from gurobipy._batch import Batch
 
@@ -25,14 +76,11 @@ from gurobipy._core import (
     disposeDefaultEnv,
     getParamInfo,
     gurobi,
-    help,
-    models,
     paramHelp,
     read,
     readParams,
     resetParams,
     setParam,
-    system,
     writeParams,
 )
 
@@ -57,4 +105,35 @@ from gurobipy._matrixapi import (
 
 from gurobipy._model import Model
 
+from gurobipy._load_model import loadModel
+
 import gurobipy.nlfunc as nlfunc
+
+
+# fmt: off
+__all__ = [
+    # _batch
+    "Batch",
+    # _core
+    "Column", "Constr", "Env", "GenConstr", "GenExpr", "LinExpr", "QConstr",
+    "QuadExpr", "SOS", "Var", "NLExpr", "TempConstr", "tuplelist", "tupledict",
+    "abs_", "all_", "and_", "any_", "max_", "min_", "norm", "or_",
+    "disposeDefaultEnv", "getParamInfo", "gurobi", "paramHelp", "read",
+    "readParams", "resetParams", "setParam", "writeParams",
+    # _exception
+    "GurobiError",
+    # _grb
+    "GRB",
+    # _helpers
+    "multidict", "quicksum",
+    # _matrixapi
+    "MConstr", "MGenConstr", "MLinExpr", "MQConstr", "MQuadExpr", "MVar",
+    "MNLExpr", "concatenate", "hstack", "vstack",
+    # _model
+    "Model",
+    # _load_model
+    "loadModel",
+    # _nlfunc
+    "nlfunc",
+]
+# fmt: on

gurobipy/__init__.pyi

@@ -85,7 +85,6 @@ def getParamInfo(
     __paramname: str
 ) -> Tuple[str, Type[Any], Any, Any, Any, Any]: ...
 def help(argument: Any = None) -> None: ...
-def models() -> None: ...
 @overload
 def multidict(
     __data: Mapping[_T, float]
@@ -151,11 +150,12 @@ def resetParams() -> None: ...
 def setParam(paramname: str, newvalue: float) -> None: ...
 @overload
 def setParam(paramname: str, newvalue: str) -> None: ...
-def system(__command: str) -> None: ...
 def writeParams(__filename: str) -> None: ...
 
 class AttrConstClass:
+    BarStatus: str = ...
     BarIterCount: str = ...
+    PDHGIterCount: str = ...
     BarX: str = ...
     BarPi: str = ...
     BatchErrorCode: str = ...
@@ -260,6 +260,7 @@ class AttrConstClass:
     MinRHS: str = ...
     ModelName: str = ...
     ModelSense: str = ...
+    NLBarIterCount: str = ...
     NodeCount: str = ...
     NumBinVars: str = ...
     NumConstrs: str = ...
@@ -267,6 +268,7 @@ class AttrConstClass:
     NumIntVars: str = ...
     NumNZs: str = ...
     NumObj: str = ...
+    NumObjPasses: str = ...
     NumPWLObjVars: str = ...
     NumQCNZs: str = ...
     NumQConstrs: str = ...
@@ -282,18 +284,40 @@ class AttrConstClass:
     ObjN: str = ...
     ObjNAbsTol: str = ...
     ObjNCon: str = ...
+    ObjNPass: str = ...
+    ObjNPassIterCount: str = ...
+    ObjNPassMIPGap: str = ...
     ObjNName: str = ...
+    ObjNPassNodeCount: str = ...
+    ObjNPassObjBound: str = ...
+    ObjNPassObjVal: str = ...
+    ObjNPassOpenNodeCount: str = ...
     ObjNPriority: str = ...
     ObjNRelTol: str = ...
+    ObjNPassRuntime: str = ...
+    ObjNPassStatus: str = ...
     ObjNVal: str = ...
     ObjNWeight: str = ...
+    ObjNPassWork: str = ...
     ObjVal: str = ...
     OpenNodeCount: str = ...
+    PoolNBoundVio: str = ...
+    PoolNBoundVioIndex: str = ...
+    PoolNBoundVioSum: str = ...
+    PoolNConstrVio: str = ...
+    PoolNConstrVioIndex: str = ...
+    PoolNConstrVioSum: str = ...
+    PoolNIntVio: str = ...
+    PoolNIntVioIndex: str = ...
+    PoolNIntVioSum: str = ...
+    PoolNMaxVio: str = ...
     PStart: str = ...
     PWLObjCvx: str = ...
     Partition: str = ...
     Pi: str = ...
     PoolIgnore: str = ...
+    PoolNObjVal: str = ...
+    PoolNX: str = ...
     PoolObjBound: str = ...
     PoolObjVal: str = ...
     PreFixVal: str = ...
@@ -425,6 +449,13 @@ class CallbackClass:
     MULTIOBJ_NODLFT: int = ...
     MULTIOBJ_RUNTIME: int = ...
     MULTIOBJ_WORK: int = ...
+    PDHG: int = ...
+    PDHG_COMPL: int = ...
+    PDHG_DUALINF: int = ...
+    PDHG_DUALOBJ: int = ...
+    PDHG_ITRCNT: int = ...
+    PDHG_PRIMINF: int = ...
+    PDHG_PRIMOBJ: int = ...
     POLLING: int = ...
     PRESOLVE: int = ...
     PRE_BNDCHG: int = ...
@@ -658,6 +689,7 @@ class GRB:
     METHOD_DETERMINISTIC_CONCURRENT_SIMPLEX: int = ...
     METHOD_DUAL: int = ...
     METHOD_NONE: int = ...
+    METHOD_PDHG: int = ...
     METHOD_PRIMAL: int = ...
     MINIMIZE: int = ...
     NODE_LIMIT: int = ...
@@ -1291,6 +1323,7 @@ class MVar:
     PWLObjCvx: int = ...
     Partition: int = ...
     PoolIgnore: int = ...
+    PoolNX: float = ...
     RC: float = ...
     SALBLow: float = ...
     SALBUp: float = ...
@@ -1455,6 +1488,7 @@ class MVar:
 
 class Model:
     BarIterCount: int = ...
+    PDHGIterCount: int = ...
     BoundSVio: float = ...
     BoundSVioIndex: int = ...
     BoundSVioSum: float = ...
@@ -1526,6 +1560,7 @@ class Model:
     MinRHS: float = ...
     ModelName: str = ...
     ModelSense: int = ...
+    NLBarIterCount: int = ...
     NodeCount: float = ...
     NumBinVars: int = ...
     NumConstrs: int = ...
@@ -1533,6 +1568,7 @@ class Model:
     NumIntVars: int = ...
     NumNZs: int = ...
     NumObj: int = ...
+    NumObjPasses: int = ...
     NumPWLObjVars: int = ...
     NumQCNZs: int = ...
     NumQConstrs: int = ...
@@ -1546,14 +1582,35 @@ class Model:
     ObjCon: float = ...
     ObjNAbsTol: float = ...
     ObjNCon: float = ...
+    ObjNPass: int = ...
+    ObjNPassIterCount: float = ...
+    ObjNPassMIPGap: float = ...
+    ObjNPassNodeCount: float = ...
     ObjNName: str = ...
+    ObjNPassObjBound: float = ...
+    ObjNPassObjVal: float = ...
+    ObjNPassOpenNodeCount: float = ...
     ObjNPriority: int = ...
     ObjNRelTol: float = ...
+    ObjNPassRuntime: float = ...
+    ObjNPassStatus: int = ...
     ObjNVal: float = ...
     ObjNWeight: float = ...
+    ObjNPassWork: float = ...
     ObjVal: float = ...
     OpenNodeCount: float = ...
     Params: ParamClass = ...
+    PoolNBoundVio: float = ...
+    PoolNBoundVioIndex: int = ...
+    PoolNBoundVioSum: float = ...
+    PoolNConstrVio: float = ...
+    PoolNConstrVioIndex: int = ...
+    PoolNConstrVioSum: float = ...
+    PoolNIntVio: float = ...
+    PoolNIntVioIndex: int = ...
+    PoolNIntVioSum: float = ...
+    PoolNMaxVio: float = ...
+    PoolNObjVal: float = ...
     PoolObjBound: float = ...
     PoolObjVal: float = ...
     Runtime: float = ...
@@ -2136,6 +2193,11 @@ class Model:
     def fixed(self) -> Model: ...
     # no type hinting for scipy.sparse
     def getA(self) -> Any: ...
+    def getQ(self) -> Any: ...
+    @overload
+    def getQCMatrices(self, qc: QConstr) -> Tuple[Any, Any]: ...
+    @overload
+    def getQCMatrices(self, qc: MQConstr) -> Tuple[Any, Any]: ...
     @overload
     def getAttr(self, attrname: str) -> Any: ...
     @overload
@@ -2329,6 +2391,8 @@ class Model:
     @overload
     def setAttr(self, attrname: str, arg1: str) -> None: ...
     @overload
+    def setAttr(self, attrname: str, arg1: Sequence[_Scalar]) -> None: ...
+    @overload
     def setAttr(
         self,
         attrname: str,
@@ -2555,6 +2619,7 @@ class ParamClass:
     DualReductions: int = ...
     FeasRelaxBigM: float = ...
     FeasibilityTol: float = ...
+    FixVarsInIndicators: int = ...
     FlowCoverCuts: int = ...
     FlowPathCuts: int = ...
     FuncMaxVal: float = ...
@@ -2572,8 +2637,10 @@ class ParamClass:
     ImproveStartGap: float = ...
     ImproveStartNodes: float = ...
     ImproveStartTime: float = ...
+    ImproveStartWork: float = ...
     InfProofCuts: int = ...
     InfUnbdInfo: int = ...
+    InheritParams: int = ...
     IntFeasTol: float = ...
     IntegralityFocus: int = ...
     IterationLimit: float = ...
@@ -2601,6 +2668,11 @@ class ParamClass:
     NLPHeur: int = ...
     NetworkAlg: int = ...
     NetworkCuts: int = ...
+    NLBarCFeasTol: float = ...
+    NLBarDFeasTol: float = ...
+    NLBarIterLimit: int = ...
+    NLBarPFeasTol: float = ...
+    NoRelHeurSolutions: int = ...
     NoRelHeurTime: float = ...
     NoRelHeurWork: float = ...
     NodeLimit: float = ...
@@ -2613,8 +2685,14 @@ class ParamClass:
     OBBT: int = ...
     ObjNumber: int = ...
     ObjScale: float = ...
+    OptimalityTarget: int = ...
     OptimalityTol: float = ...
     OutputFlag: int = ...
+    PDHGIterLimit: float = ...
+    PDHGRelTol: float = ...
+    PDHGAbsTol: float = ...
+    PDHGConvTol: float = ...
+    PDHGGPU: int = ...
     PSDCuts: int = ...
     PSDTol: float = ...
     PartitionPlace: int = ...
@@ -2658,6 +2736,8 @@ class ParamClass:
     SolutionNumber: int = ...
     StartNodeLimit: int = ...
     StartNumber: int = ...
+    StartTimeLimit: float = ...
+    StartWorkLimit: float = ...
     StrongCGCuts: int = ...
     SubMIPCuts: int = ...
     SubMIPNodes: int = ...
@@ -2740,6 +2820,7 @@ class ParamConstClass:
     DualReductions: str = ...
     FeasRelaxBigM: str = ...
     FeasibilityTol: str = ...
+    FixVarsInIndicators: str = ...
     FlowCoverCuts: str = ...
     FlowPathCuts: str = ...
     FuncMaxVal: str = ...
@@ -2757,8 +2838,10 @@ class ParamConstClass:
     ImproveStartGap: str = ...
     ImproveStartNodes: str = ...
     ImproveStartTime: str = ...
+    ImproveStartWork: str = ...
     InfProofCuts: str = ...
     InfUnbdInfo: str = ...
+    InheritParams: str = ...
     IntFeasTol: str = ...
     IntegralityFocus: str = ...
     IterationLimit: str = ...
@@ -2786,6 +2869,11 @@ class ParamConstClass:
     NLPHeur: str = ...
     NetworkAlg: str = ...
     NetworkCuts: str = ...
+    NLBarCFeasTol: str = ...
+    NLBarDFeasTol: str = ...
+    NLBarIterLimit: str = ...
+    NLBarPFeasTol: str = ...
+    NoRelHeurSolutions: str = ...
     NoRelHeurTime: str = ...
     NoRelHeurWork: str = ...
     NodeLimit: str = ...
@@ -2799,7 +2887,13 @@ class ParamConstClass:
     ObjNumber: str = ...
     ObjScale: str = ...
     OptimalityTol: str = ...
+    OptimalityTarget: str = ...
     OutputFlag: str = ...
+    PDHGIterLimit: str = ...
+    PDHGRelTol: str = ...
+    PDHGAbsTol: str = ...
+    PDHGConvTol: str = ...
+    PDHGGPU: str = ...
     PSDCuts: str = ...
     PSDTol: str = ...
     PartitionPlace: str = ...
@@ -2843,6 +2937,8 @@ class ParamConstClass:
     SolutionNumber: str = ...
     StartNodeLimit: str = ...
     StartNumber: str = ...
+    StartTimeLimit: str = ...
+    StartWorkLimit: str = ...
     StrongCGCuts: str = ...
     SubMIPCuts: str = ...
     SubMIPNodes: str = ...
@@ -3106,6 +3202,7 @@ class Var:
     PWLObjCvx: int = ...
     Partition: int = ...
     PoolIgnore: int = ...
+    PoolNX: float = ...
     RC: float = ...
     SALBLow: float = ...
     SALBUp: float = ...

gurobipy/_attrconst.py

@@ -34,7 +34,8 @@ class AttrConstClass(object):
     objVal: Objective value for current solution
     objBound: Best available lower bound on solution objective
     poolObjBound: Bound on objective for undiscovered solutions
-    poolObjVal: Retrieve the objective value of an alternate MIP solution
+    poolNObjVal: Retrieve the objective value of an alternate MIP solution
+    poolObjVal: Deprecated since v13 - use poolNObjVal instead
     MIPGap: Current relative MIP optimality gap
     runtime: Runtime (in seconds) for most recent optimization
     work: Work (in work units) for most recent optimization
@@ -44,6 +45,8 @@ class AttrConstClass(object):
     iterCount: Number of simplex iterations performed
     nodeCount: Number of branch-and-cut nodes explored
     barIterCount: Number of barrier iterations performed
+    nlbarIterCount: Number of NL barrier iterations performed
+    PDHGIterCount: Number of PDHG iterations performed
     isMIP: Indicates whether the model is MIP (1) or not (0)
     isQP: Indicates whether the model has a quadratic objective
     isQCP: Indicates whether the model has quadratic constraints
@@ -82,7 +85,8 @@ class AttrConstClass(object):
     varName: Variable name
     x: Solution value
     rc: Reduced cost
-    xn: Solution value in alternate MIP solution
+    poolNX: Solution value in alternate MIP solution
+    xn: Deprecated since v13 - use poolNX instead
     start: Start vector (use GRB.UNDEFINED to leave a value unspecified)
     vBasis: Basis status
     varHintVal: Variable hint value
@@ -177,6 +181,18 @@ class AttrConstClass(object):
     ObjNName: names of multi-objectives
     NumObj: number of multi-objectives
 
+    NumObjPasses: number of optimization passes during the multi-objective solve
+    ObjNPass: optimization pass in which the selected objective function was solved
+    ObjNPassIterCount: simplex iteration count for a pass during the multi-objective solve
+    ObjNPassMIPGap: MIP gap for a pass during the multi-objective solve
+    ObjNPassNodeCount: number of explored nodes for a pass during the multi-objective solve
+    ObjNPassObjBound: objective bound for a pass during the multi-objective solve
+    ObjNPassObjVal: objective value for a pass during the multi-objective solve
+    ObjNPassOpenNodeCount: number of unexplored nodes for a pass during the multi-objective solve
+    ObjNPassRuntime: run time for a pass during the multi-objective solve
+    ObjNPassStatus: status for a pass during the multi-objective solve
+    ObjNPassWork: work done for a pass during the multi-objective solve
+
   Multi-scenarios
 
     ScenNLB: modification to lower bound vector in current scenario
@@ -245,6 +261,7 @@ class AttrConstClass(object):
   ObjBound        = "ObjBound"
   ObjBoundC       = "ObjBoundC"
   PoolObjBound    = "PoolObjBound"
+  PoolNObjVal     = "PoolNObjVal"
   PoolObjVal      = "PoolObjVal"
   MIPGap          = "MIPGap"
   Runtime         = "Runtime"
@@ -253,7 +270,10 @@ class AttrConstClass(object):
   SolCount        = "SolCount"
   IterCount       = "IterCount"
   NodeCount       = "NodeCount"
+  BarStatus       = "BarStatus"
   BarIterCount    = "BarIterCount"
+  NLBarIterCount  = "NLBarIterCount"
+  PDHGIterCount   = "PDHGIterCount"
   FarkasProof     = "FarkasProof"
   NumStart        = "NumStart"
   TuneResultCount = "TuneResultCount"
@@ -272,6 +292,7 @@ class AttrConstClass(object):
   VarName        = "VarName"
   X              = "X"
   RC             = "RC"
+  PoolNX         = "PoolNX"
   Xn             = "Xn"
   BarX           = "BarX"
   BarPi          = "BarPi"
@@ -399,6 +420,19 @@ class AttrConstClass(object):
   IntVioIndex          = "IntVioIndex"
   IntVioSum            = "IntVioSum"
 
+# Solution pool quality attributes
+
+  PoolNMaxVio          = "PoolNMaxVio"
+  PoolNBoundVio        = "PoolNBoundVio"
+  PoolNBoundVioIndex   = "PoolNBoundVioIndex"
+  PoolNBoundVioSum     = "PoolNBoundVioSum"
+  PoolNConstrVio       = "PoolNConstrVio"
+  PoolNConstrVioIndex  = "PoolNConstrVioIndex"
+  PoolNConstrVioSum    = "PoolNConstrVioSum"
+  PoolNIntVio          = "PoolNIntVio"
+  PoolNIntVioIndex     = "PoolNIntVioIndex"
+  PoolNIntVioSum       = "PoolNIntVioSum"
+
 # Multi-objective attributes
 
   ObjN                 = "ObjN"
@@ -411,6 +445,18 @@ class AttrConstClass(object):
   ObjNName             = "ObjNName"
   NumObj               = "NumObj"
 
+  NumObjPasses          = "NumObjPasses"
+  ObjNPass              = "ObjNPass"
+  ObjNPassIterCount     = "ObjNPassIterCount"
+  ObjNPassMIPGap        = "ObjNPassMIPGap"
+  ObjNPassNodeCount     = "ObjNPassNodeCount"
+  ObjNPassObjBound      = "ObjNPassObjBound"
+  ObjNPassObjVal        = "ObjNPassObjVal"
+  ObjNPassOpenNodeCount = "ObjNPassOpenNodeCount"
+  ObjNPassRuntime       = "ObjNPassRuntime"
+  ObjNPassStatus        = "ObjNPassStatus"
+  ObjNPassWork          = "ObjNPassWork"
+
 # Multi-scenario attributes
 
   ScenNLB              = "ScenNLB"

gurobipy/_callbackconst.py

@@ -19,6 +19,7 @@ class CallbackConstClass(object):
     MIPNODE:  MIP node information available
     MULTIOBJ: In multi-objective optimization
     IIS:      In IIS computation
+    NLBAR:    In NL barrier
     MESSAGE:  Optimizer output a message
 
   Your method can call Model.cbGet() to obtain detailed information
@@ -92,6 +93,11 @@ class CallbackConstClass(object):
     IIS_BOUNDMIN: Minimum number of bounds in IIS
     IIS_BOUNDMAX: Maximum number of bounds in IIS
     IIS_BOUNDGUESS: Estimated number of bounds in IIS
+    NLBAR_ITRCNT: NL barrier iteration count
+    NLBAR_PRIMOBJ: NL barrier iterate primal objective
+    NLBAR_PRIMINF: NL barrier iterate primal infeasibility
+    NLBAR_DUALINF: NL barrier iterate dual infeasibility
+    NLBAR_COMPL: NL barrier iterate complementarity violation
     MSG_STRING: Output message
 
   Your callback method can call other methods on the model object:
@@ -113,6 +119,8 @@ class CallbackConstClass(object):
   BARRIER       =     7
   MULTIOBJ      =     8
   IIS           =     9
+  PDHG          =    10
+  NLBAR         =    11
 
   PRE_COLDEL     =  1000
   PRE_ROWDEL     =  1001
@@ -183,3 +191,16 @@ class CallbackConstClass(object):
   IIS_BOUNDMIN    = 9004
   IIS_BOUNDMAX    = 9005
   IIS_BOUNDGUESS  = 9006
+
+  PDHG_ITRCNT     = 10001
+  PDHG_PRIMOBJ    = 10002
+  PDHG_DUALOBJ    = 10003
+  PDHG_PRIMINF    = 10004
+  PDHG_DUALINF    = 10005
+  PDHG_COMPL      = 10006
+
+  NLBAR_ITRCNT  =  11001
+  NLBAR_PRIMOBJ =  11002
+  NLBAR_PRIMINF =  11003
+  NLBAR_DUALINF =  11004
+  NLBAR_COMPL   =  11005

gurobipy/_grb.py

@@ -51,23 +51,25 @@ class GRB(object):
 
 # Status codes
 
-  LOADED          = 1
-  OPTIMAL         = 2
-  INFEASIBLE      = 3
-  INF_OR_UNBD     = 4
-  UNBOUNDED       = 5
-  CUTOFF          = 6
-  ITERATION_LIMIT = 7
-  NODE_LIMIT      = 8
-  TIME_LIMIT      = 9
-  SOLUTION_LIMIT  = 10
-  INTERRUPTED     = 11
-  NUMERIC         = 12
-  SUBOPTIMAL      = 13
-  INPROGRESS      = 14
-  USER_OBJ_LIMIT  = 15
-  WORK_LIMIT      = 16
-  MEM_LIMIT       = 17
+  LOADED             = 1
+  OPTIMAL            = 2
+  INFEASIBLE         = 3
+  INF_OR_UNBD        = 4
+  UNBOUNDED          = 5
+  CUTOFF             = 6
+  ITERATION_LIMIT    = 7
+  NODE_LIMIT         = 8
+  TIME_LIMIT         = 9
+  SOLUTION_LIMIT     = 10
+  INTERRUPTED        = 11
+  NUMERIC            = 12
+  SUBOPTIMAL         = 13
+  INPROGRESS         = 14
+  USER_OBJ_LIMIT     = 15
+  WORK_LIMIT         = 16
+  MEM_LIMIT          = 17
+  LOCALLY_OPTIMAL    = 18
+  LOCALLY_INFEASIBLE = 19
 
 # Batch status codes
 
@@ -142,6 +144,8 @@ class GRB(object):
   OPCODE_LOG2         = 15
   OPCODE_LOG10        = 16
   OPCODE_LOGISTIC     = 17
+  OPCODE_TANH         = 18
+  OPCODE_SIGNPOW      = 19
 
 # Basis status
 
@@ -169,9 +173,9 @@ class GRB(object):
 
 # Version number
 
-  VERSION_MAJOR     = 12
+  VERSION_MAJOR     = 13
   VERSION_MINOR     = 0
-  VERSION_TECHNICAL = 2
+  VERSION_TECHNICAL = 0
 
 # Errors
 
@@ -236,6 +240,7 @@ class GRB(object):
   METHOD_CONCURRENT                       = 3
   METHOD_DETERMINISTIC_CONCURRENT         = 4
   METHOD_DETERMINISTIC_CONCURRENT_SIMPLEX = 5
+  METHOD_PDHG                             = 6
 
 # BarHomogeneous parameter values