#include <Eigen/Dense>
int main(){
auto c = Eigen::ArrayXi::LinSpaced(11, 0, 10);
auto r = Eigen::pow(3.7, c);
}
but if I switch the ArrayXi for ArrayXd it compiles (because the scalar types match), but that is not what I want.
Really, what I need is to call my own implementation of pow(double, int), because it is MUCH faster than pow(double, double).
In the end, I would like some guidance on how I can craft an Eigen _expression_ for this use case that will not introduce any intermediate allocation.
Thanks,
Ian