solver with test-case

src/cvx/ball/solver.py

@@ -0,0 +1,25 @@
+import cvxpy as cp
+import numpy as np
+
+from .utils.circle import Circle
+
+
+def min_circle_cvx(points, **kwargs):
+    # Use con_1 if no constraint construction is defined
+    # cvxpy variable for the radius
+    r = cp.Variable(name="Radius")
+    # cvxpy variable for the midpoint
+    x = cp.Variable(points.shape[1], name="Midpoint")
+    objective = cp.Minimize(r)
+    constraints = [
+        cp.SOC(
+            r * np.ones(points.shape[0]),
+            points - cp.outer(np.ones(points.shape[0]), x),
+            axis=1,
+        )
+    ]
+
+    problem = cp.Problem(objective=objective, constraints=constraints)
+    problem.solve(**kwargs)
+
+    return Circle(radius=float(r.value), center=x.value)