dalon-work / fortcpp

FortC++ provides header-only arrays that act like Fortran arrays

FortC++

FortC++ is a templated header library designed to add Fortran-style arrays into C++. This is not the first or the best such offering. Here is a list of others:

• ndarray
• Blitz++
• Boost::MultiArray
• Eigen

This library is not intended to be a linear algebra package. It is intended to create user-friendly arrays which provide whole-array operations, slicing, views, multi-d indexing, and bounds checking. Whole-array operations do not generate temporaries (with appropriate inlining), keeping the overhead of using the arrays down.

FortC++ supports any rank of array using variadic templates. Column-Major ordering of the data is the default, but Row-Major ordering is fully supported. There are two main components to FortCpp:

FortCpp::Alloc< Type, Rank, Options >
FortCpp::Fixed< Type, Dims... >

The Alloc class specifies an allocatable array, whose dimensions are only known at runtime. This is the main workhorse of FortC++. The Fixed class specifies a static array whose size must be determined at compile time. These two classes are fully compatible with each other.

---

How to use FortC++

Because FortC++ is a header-file only library, using the library is as straightforward as placing the library in your project's directory and including FortCpp.h.

For a more detailed introduction to FortC++, see the tutorial.

fortcpp's Issues

where

This will have to act similar to NumPy's where routine.

MultiD Iterator

Implement a multiD Iterator, compare it's speed with normal usage, etc.

Lbound shift

Create a wrapper class to specify the index of the lower bounds of an array

Vector Indexing

Create new class that allows for Vector Indexing, similar to Fortran:

real,dimension(10) :: a
a = 1.0
a( [1,3,4,7] ) = 0.0

I'm thinking of using Array1d as a linear indexer, which can then be used in functions like where

Row-Major Slicing

This might require a lot or rewriting.

Align data

Need to default align dynamically allocated arrays

Fixed Size Arrays

Order-specific offset

Specialize the offsets to account for Col vs. Row Major arrays.

struct member/function access

By overloading the ->* operator, it might be possible to create arrays that reference a particular member in a struct, or to apply a member function to all elements in an array.

Contiguous Slice Produces Contiguous Array

Create LHS resize

Make a function that allows the LHS of an operation to be resized to the size of the RHS.
Something like this:

Array1d<float> a,b(10),c(10);
a.resize() = b+c;

Unit Tests

Create a small testing framework to make sure everything is working correctly, at all times.