The memory layout in Lcqpow is row major (not documented) - python bindings use the wrong order for Eigen types - matrices have to be transposed to avoid pybind11 temporary copy (which would have wrong layout), see any Python example
The memory layout in Lcqpow is row major (not documented)