Clarabel.cpp/examples/cpp/example_expcone.cpp at a42179aa67f6f8ba42168668b8ab3f3ab2cb4c43 · oxfordcontrol/Clarabel.cpp · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
#include "utils.h"

#include <clarabel.hpp>
#include <Eigen/Eigen>
#include <cmath>
#include <cstdint>
#include <vector>

using namespace clarabel;
using namespace std;
using namespace Eigen;

int main(void)
{
    SparseMatrix<double> P = MatrixXd::Zero(3, 3).sparseView();
    P.makeCompressed();

    Vector<double, 3> q{ -1.0, 0.0, 0.0 };

    MatrixXd A_dense(5, 3);
    A_dense <<
        -1.,  0.,  0.,
         0., -1.,  0.,
         0.,  0., -1.,
         0.,  1.,  0.,
         0.,  0.,  1.;

    SparseMatrix<double> A = A_dense.sparseView();
    A.makeCompressed();

    Vector<double, 5> b { 0., 0., 0., 1., exp(5.0) };

    vector<SupportedConeT<double>> cones
    {
        ExponentialConeT<double>(),
        ZeroConeT<double>(2),
    };

    // Settings
    DefaultSettings<double> settings = DefaultSettingsBuilder<double>::default_settings()
                                           .verbose(true)
                                           .build();

    // Build solver
    DefaultSolver<double> solver(P, q, A, b, cones, settings);

    // Solve
    solver.solve();

    // Get solution
    DefaultSolution<double> solution = solver.solution();
    utils::print_solution(solution);

    return 0;
}