1 files changed, 0 insertions, 133 deletions
diff --git a/polymatrix/expression/mixins/determinantexprmixin.py b/polymatrix/expression/mixins/determinantexprmixin.py
deleted file mode 100644
index 4659ded..0000000
--- a/polymatrix/expression/mixins/determinantexprmixin.py
+++ /dev/null
@@ -1,133 +0,0 @@
-from __future__ import annotations
-
-import abc
-import collections
-import dataclasses
-import typing
-
-if typing.TYPE_CHECKING:
-    from polymatrix.expressionstate.abc import ExpressionState
-
-from polymatrix.polymatrix.init import init_poly_matrix
-from polymatrix.expression.mixins.expressionbasemixin import ExpressionBaseMixin
-from polymatrix.polymatrix.abc import PolyMatrix
-
-
-class DeterminantExprMixin(ExpressionBaseMixin):
-    @property
-    @abc.abstractmethod
-    def underlying(self) -> ExpressionBaseMixin: ...
-
-    # NP: delete this code?
-    # # overwrites the abstract method of `ExpressionBaseMixin`
-    # @property
-    # def shape(self) -> tuple[int, int]:
-    #     return self.underlying.shape[0], 1
-
-    # overwrites the abstract method of `ExpressionBaseMixin`
-    def apply(
-        self,
-        state: ExpressionState,
-    ) -> tuple[ExpressionState, PolyMatrix]:
-        # raise Exception('not implemented')
-
-        if self in state.cache:
-            return state, state.cache[self]
-
-        state, underlying = self.underlying.apply(state=state)
-
-        assert underlying.shape[0] == underlying.shape[1]
-
-        inequality_data = {}
-        auxillary_equations = {}
-
-        index_start = state.n_param
-        rel_index = 0
-
-        # NP: consider providing a reference eg algorithm name if it has one
-        # NP: or at least a minimal overview / explaination.
-
-        for row in range(underlying.shape[0]):
-            # FIXME: type annotations of dictionaries inside method
-            polynomial = collections.defaultdict(float)
-
-            # f in f-v^T@x-r^2
-            # terms = underlying.get_poly(row, row)
-            try:
-                underlying_poly = underlying.get_poly(row, row)
-
-            # NP: this is polymatrix leaking the abstraction, you should not
-            # NP: need to care about these problems here. => Fix inside PolyMatrix
-            except KeyError:
-                pass
-            else:
-                for monomial, value in underlying_poly.items():
-                    polynomial[monomial] += value
-
-            for inner_row in range(row):
-                # -v^T@x in f-v^T@x-r^2
-                # terms = underlying.get_poly(row, inner_row)
-                try:
-                    underlying_poly = underlying.get_poly(row, inner_row)
-                except KeyError:
-                    pass
-                else:
-                    for monomial, value in underlying_poly.items():
-                        new_monomial = monomial + (index_start + rel_index + inner_row,)
-                        polynomial[new_monomial] -= value
-
-                # auxillary terms
-                # ---------------
-
-                auxillary_polynomial = collections.defaultdict(float)
-
-                for inner_col in range(row):
-                    # P@x in P@x-v
-                    key = tuple(reversed(sorted((inner_row, inner_col))))
-                    try:
-                        underlying_poly = underlying.get_poly(*key)
-                    except KeyError:
-                        pass
-                    else:
-                        for monomial, value in underlying_poly.items():
-                            new_monomial = monomial + (
-                                index_start + rel_index + inner_col,
-                            )
-                            auxillary_polynomial[new_monomial] += value
-
-                # -v in P@x-v
-                try:
-                    underlying_poly = underlying.get_poly(row, inner_row)
-                except KeyError:
-                    pass
-                else:
-                    for monomial, value in underlying_poly.items():
-                        auxillary_polynomial[monomial] -= value
-
-                x_variable = index_start + rel_index + inner_row
-                assert x_variable not in state.auxillary_equations
-                auxillary_equations[x_variable] = dict(auxillary_polynomial)
-
-            rel_index += row
-            inequality_data[row, 0] = dict(polynomial)
-
-        # NP: Type checker does not like the fact that state is of type
-        # NP: ExpressionState, while register returns ExpressionStateMixinis this
-        # NP: was a common problem I had while coding mdpoly, that eventually led
-        # NP: me to abandoning mixins because they are a pain to get right with typing 
-        # NP: such that the type checker (mypy) is happy. This problem happens all over the
-        # NP: place, I don't have a quick solution.
-        state = state.register(rel_index)
-
-        poly_matrix = init_poly_matrix(
-            data=inequality_data,
-            shape=(underlying.shape[0], 1),
-        )
-
-        state = dataclasses.replace(
-            state,
-            auxillary_equations=state.auxillary_equations | auxillary_equations,
-            cache=state.cache | {self: poly_matrix},
-        )
-
-        return state, poly_matrix