import clarabel
import numpy as np
from scipy import sparse
n = 20000
P = sparse.identity(n).tocsc()
A1 = sparse.identity(n).tocsc()
A2 = -sparse.identity(n).tocsc()
A = sparse.vstack([A1, A2]).tocsc()
q = np.ones(n)
b = np.ones(2*n)
cones = [clarabel.NonnegativeConeT(b.size)]
settings = clarabel.DefaultSettings()
solver = clarabel.DefaultSolver(P, q, A, b, cones, settings)
solution = solver.solve()