Ecosyste.ms: Awesome
An open API service indexing awesome lists of open source software.
https://github.com/rurlus/carma
Converters between Armadillo matrices (C++) and Numpy arrays using Pybind11
https://github.com/rurlus/carma
armadillo c-plus-plus carma cpp numpy pybind11 python3 rcpparmadillo
Last synced: about 6 hours ago
JSON representation
Converters between Armadillo matrices (C++) and Numpy arrays using Pybind11
- Host: GitHub
- URL: https://github.com/rurlus/carma
- Owner: RUrlus
- License: apache-2.0
- Created: 2020-01-23T19:17:19.000Z (almost 5 years ago)
- Default Branch: stable
- Last Pushed: 2024-11-09T17:11:41.000Z (about 1 month ago)
- Last Synced: 2024-12-15T17:12:06.358Z (7 days ago)
- Topics: armadillo, c-plus-plus, carma, cpp, numpy, pybind11, python3, rcpparmadillo
- Language: C++
- Homepage:
- Size: 1.93 MB
- Stars: 88
- Watchers: 3
- Forks: 24
- Open Issues: 4
-
Metadata Files:
- Readme: README.md
- Changelog: CHANGELOG.md
- License: LICENSE
- Citation: CITATION.cff
Awesome Lists containing this project
README
A C++ header only library providing conversions between Numpy arrays and Armadillo matrices.
|
Documentation
|
[](https://github.com/RUrlus/carma/actions/workflows/linux.yml)
[](https://github.com/RUrlus/carma/actions/workflows/macos.yml)
[](https://github.com/RUrlus/carma/actions/workflows/windows.yml)
[](https://coveralls.io/github/RUrlus/carma?branch=stable)
[](https://carma.readthedocs.io/en/latest/?badge=latest)
[](https://github.com/RUrlus/carma/blob/stable/LICENSE)
[](https://github.com/RUrlus/carma/releases)
[](https://zenodo.org/badge/latestdoi/235867993)
## Introduction
CARMA provides fast bidirectional conversions between [Numpy](https://numpy.org) arrays and [Armadillo](http://arma.sourceforge.net/docs.html) matrices, vectors and cubes, much like [RcppArmadillo](https://github.com/RcppCore/RcppArmadillo) does for R and Armadillo.
The library extends the impressive [pybind11](https://pybind11.readthedocs.io/en/stable/intro.html) library with support for Armadillo.
For details on Pybind11 and Armadillo refer to their respective documentation [1](https://pybind11.readthedocs.io/en/stable/intro.html), [2](http://arma.sourceforge.net/docs.html).
## Installation
CARMA is a header only library that relies on two other header only libraries, Armadillo and Pybind11.
It can be integrated in a CMake build using `ADD_SUBDIRECTORY()` or installation which provides an interface library target `carma::carma` that has been linked with Python, Numpy, Pybind11 and Armadillo. See [build configuration](https://carma.readthedocs.io/en/stable/building.html) for details.
It can be installed using:
```bash
mkdir build
cd build
# optionally with -DCMAKE_INSTALL_PREFIX:PATH=
cmake -DCARMA_INSTALL_LIB=ON ..
cmake --build . --config Release --target install
```
You can than include it in a project using:
```cmake
FIND_PACKAGE(carma CONFIG REQUIRED)
TARGET_LINK_LIBRARIES( PRIVATE carma::carma)
```
### CMake subdirectory
Alternatively you can forgo installing CARMA and directly use it as CMake subdirectory.
For Pybind11 and or Armadillo we create target(s) based on user settable version, see [build configuration](https://carma.readthedocs.io/en/stable/building.html), when they are not defined.
To link with CARMA:
```cmake
ADD_SUBDIRECTORY(extern/carma)
TARGET_LINK_LIBRARIES( PRIVATE carma::carma)
```
CARMA and Armadillo can then be included using:
```C++
#include
#include
```
CARMA provides a number of configurations that can be set in the `carma_config.cmake` file at the root of the directory or passed to CMake, see [Configuration](https://carma.readthedocs.io/en/stable/configuration.html) and [Build configuration](https://carma.readthedocs.io/en/stable/building.html) documentation sections for details.
## Requirements
Numpy v2.* is supported by CARMA >= v0.8.0 which requires a compiler with support for C++14 and supports:
* Python 3.8 -- 3.12
* Numpy >= 1.14
* **Pybind11 >= v2.12.0**
* Armadillo >= 10.5.2
CARMA <= v0.7 requires a compiler with support for C++14 and supports:
* Python 3.8 -- 3.12
* **Numpy >= 1.14 < 2.0**
* **Pybind11 >= v2.6.0 < v2.12.0**
* Armadillo >= 10.5.2
CARMA makes use of Armadillo's `ARMA_ALIEN_MEM_ALLOC` and `ARMA_ALIEN_MEM_FREE` functionality introduced in version 10.5.2 to use Numpy's (de)allocator.
### Considerations
In order to achieve fast conversions the default behaviour is avoid copying both from and to Numpy whenever possible and reasonable.
This allows very low overhead conversions but it impacts memory safety and requires user vigilance.
If you intend to return the memory of the input array back as another array, you must make sure to either copy or steal the memory on the conversion in or copy the memory out.
If you don't the memory will be aliased by the two Numpy arrays and bad things will happen.
A second consideration is memory layout. Armadillo is optimised for column-major (Fortran order) memory whereas Numpy defaults to row-major (C order).
The default behaviour is to automatically convert, read copy, C-order arrays to F-order arrays upon conversion to Armadillo. Users should note that the library will not convert back to C-order when returning.
For details see the documentation section [Memory Management](https://carma.readthedocs.io/en/latest/memory_management.html).
### Example
On a high level CARMA provides two ways to work Numpy arrays in Armadillo:
Automatic conversion saves a bit on code but provides less flexibility with
regards to when to copy and when not.
Manual conversion should be used when you need fine grained control.
Combining the two; we use automatic conversion on the conversion in and manual when
creating the tuple for the way out.
```cpp
#include
#include
#include
#include
#include
namespace py = pybind11;
py::tuple ols(arma::mat& X, arma::colvec& y) {
// We borrow the data underlying the numpy arrays
int n = X.n_rows, k = X.n_cols;
arma::colvec coeffs = arma::solve(X, y);
arma::colvec resid = y - X * coeffs;
double sig2 = arma::as_scalar(arma::trans(resid) * resid / (n-k));
arma::colvec std_errs = arma::sqrt(sig2 * arma::diagvec( arma::inv(arma::trans(X)*X)) );
// We take ownership of the memory from the armadillo objects and
// return to python as a tuple containing two Numpy arrays.
return py::make_tuple(
carma::col_to_arr(coeffs),
carma::col_to_arr(std_errs)
);
}
// adapted from https://gallery.rcpp.org/articles/fast-linear-model-with-armadillo/
```
Which can be called using:
```python
y = np.linspace(1, 100, num=100) + np.random.normal(0, 0.5, 100)
X = np.hstack((np.ones(100)[:, None], np.arange(1, 101)[:, None]))
coeff, std_err = carma.ols(X, y)
```
The repository contains tests, examples and CMake build instructions that can be used as an reference.
### About
This project was created by Ralph Urlus. Significant improvements to the project have been contributed by [Pascal H.](https://github.com/hpwxf)
### License
CARMA is provided under a Apache 2.0 license that can be found in the LICENSE file. By using, distributing, or contributing to this project, you agree to the terms and conditions of this license.