C and C++ ABI

Solvers live in Rust. Geometry hosts that call one outer iteration at a time hold an opaque rgmin_solver_t; scripts that want a one-shot solve still call rgmin_minimize. Every vector is a dlpk DLManagedTensorVersioned tensor.

Session (eOn / ASE step)

#include "xts.h"

rgmin_control_t ctrl = { .maxiter = 200, .gtol = 1e-8, .istep = 1.0, .memory = 8 };
rgmin_solver_t *s = rgmin_solver_create(RGMIN_LBFGS, &ctrl, n);
rgmin_report_t report;
while (!done) {
    rgmin_solver_set_maxmove(s, max_move);
    rgmin_solver_set_qn_step(s, RGMIN_QN_NEWTON); /* or RGMIN_QN_LBFGS + P as H0 */
    rgmin_solver_set_accept(s, RGMIN_ACCEPT_NONE);
    rgmin_solver_step_hess_fg(s, evalgrad, hess, user, x, &report);
}
rgmin_solver_free(s);

rgmin_solver_forget drops pairs / conjugacy / moments. Newton and RFO use rgmin_solver_step_hess_fg (fused energy and gradient, one host potential call). Split eval / grad callbacks remain. accept is none (take the clipped step), energy, or nonmonotone. Callbacks are arguments of each step, not stored on the handle.

One-shot

#include "xts.h"

rgmin_control_t ctrl = { .maxiter = 200, .gtol = 1e-8, .istep = 1.0, .memory = 8 };
rgmin_report_t report;
rgmin_status_t st = rgmin_minimize(eval, grad, user, x, &ctrl, RGMIN_LBFGS, &report);

A non-CPU tensor returns RGMIN_UNSUPPORTED_DEVICE. The ABI does not change when a CUDA path lands.

C++

#include <xts/optimize.hpp>

xts::optimize::Control ctrl;
auto report = xts::optimize::minimize(eval, grad, user, x, ctrl,
                                      xts::optimize::Method::Lbfgs);

xtensor callers include xts/xtensor.hpp, which adapts xt::xarray to the same rgmin_minimize entry. There is no solver in the header.

Manifold

rgmin_solver_set_manifold(s, RGMIN_MANIFOLD_SPHERE);

Tokens: RGMIN_MANIFOLD_EUCLIDEAN (default), RIGID_QUOTIENT (R^{3N}/SE(3), Sella Cartesian), MW_RIGID (Page–McIver / Sella IRC Eckart), SPHERE, SO3 (length 9), STIEFEL (St(n,1)), SE3 (length 12). See retract onto an embedded manifold. rgmin_solver_set_masses supplies N atomic masses for MW_RIGID.