FEniCS/basix

Check array order consistency - RowMajor vs ColumnMajor

Closed this issue · 2 comments

IgorBaratta commented

Example:
basix::tabulate receives a pointer to the memory of a set of points [npoints, tdim], which is expected to be in row-major order, maps it into a row-major eigen array, but in the subsequent function calls it's treated as a column-major array.

garth-wells commented

This has caused some issues. Plan is to strip Eigen from the interface and make everything row-major.

  • 👍1
garth-wells commented

Eigen has been removed.