Source code for ocean.maxsat._variables._feature

from collections.abc import Mapping

from ...feature import Feature
from ...feature._keeper import FeatureKeeper
from ...typing import Key
from .._base import BaseModel, Var


[docs] class FeatureVar(Var, FeatureKeeper): """MaxSAT variable bundle associated with a single parsed feature.""" X_VAR_NAME_FMT: str = "x[{name}]" _x: int _u: Mapping[Key, int] _mu: Mapping[Key, int] def __init__(self, feature: Feature, name: str) -> None: Var.__init__(self, name=name) FeatureKeeper.__init__(self, feature=feature)
[docs] def build(self, model: BaseModel) -> None: if self.is_binary: self._x = self._add_x(model) if self.is_numeric: self._mu = self._add_mu(model) # Add exactly-one constraint for mu variables (exactly one interval) model.add_exactly_one(list(self._mu.values())) if self.is_one_hot_encoded: self._u = self._add_u(model)
[docs] def xget(self, code: Key | None = None, mu: Key | None = None) -> int: if mu is not None and code is not None: msg = "Cannot get both 'mu' and 'code' at the same time" raise ValueError(msg) if self.is_one_hot_encoded: return self._xget_one_hot_encoded(code) if code is not None: msg = "Get by code is only supported for one-hot encoded features" raise ValueError(msg) if self.is_numeric: return self._xget_numeric(mu) if mu is not None: msg = "Get by 'mu' is only supported for numeric features" raise ValueError(msg) return self._x
def _add_x(self, model: BaseModel) -> int: if not self.is_binary: msg = "The '_add_x' method is only supported for binary features" raise ValueError(msg) name = self.X_VAR_NAME_FMT.format(name=self._name) return self._add_binary(model, name) def _add_u(self, model: BaseModel) -> Mapping[Key, int]: name = self._name.format(name=self._name) u = self._add_one_hot_encoded(model=model, name=name) model.add_exactly_one(list(u.values())) return u def _add_one_hot_encoded( self, model: BaseModel, name: str, ) -> Mapping[Key, int]: return { code: model.add_var(name=f"{name}[{code}]") for code in self.codes } def _add_mu(self, model: BaseModel) -> Mapping[Key, int]: name = self._name.format(name=self._name) if self.is_discrete: # For discrete features: one mu variable per level (value) # mu[i] means value == levels[i] n_values = len(self.levels) return { lv: model.add_var(name=f"{name}[{lv}]") for lv in range(n_values) } # For continuous features: n-1 mu variables for n levels (intervals) # mu[i] means value in interval (levels[i], levels[i+1]] n_intervals = len(self.levels) - 1 return { lv: model.add_var(name=f"{name}[{lv}]") for lv in range(n_intervals) } @staticmethod def _add_binary(model: BaseModel, name: str) -> int: return model.add_var(name=name) def _xget_one_hot_encoded(self, code: Key | None) -> int: if code is None: msg = "Code is required for one-hot encoded features get" raise ValueError(msg) if code not in self.codes: msg = f"Code '{code}' not found in the feature codes" raise ValueError(msg) return self._u[code] def _xget_numeric(self, mu: Key | None) -> int: if mu is None: msg = "mu is required to get numeric features" raise ValueError(msg) if self.is_discrete: # For discrete: mu[i] represents value levels[i] n_values = len(self.levels) if mu not in range(n_values): msg = f"mu '{mu}' not in values (0 to {n_values - 1})" raise ValueError(msg) else: # For continuous: mu[i] represents interval (levels[i], levels[i+1]] n_intervals = len(self.levels) - 1 if mu not in range(n_intervals): msg = f"mu '{mu}' not in intervals (0 to {n_intervals - 1})" raise ValueError(msg) return self._mu[mu]