Указатель на функцию с Eigen

Я довольно хорошо разбираюсь в Python, но я новичок в C++ и таких вещах, как указатели. Я пытаюсь написать несколько кодов для решения ODE с помощью пакета Eigen для линейной алгебры (позже мне нужно будет иметь дело с множеством матриц, поэтому я планирую начать с него). У меня есть следующий код для RK4, и они работают:

#include "../eigen-eigen-b3f3d4950030/Eigen/Dense"
using namespace Eigen;

VectorXd Func(const VectorXd& a)
{ // equations for solving simple harmonic oscillator
    Vector2d ans;
    ans(0) = a(1);    //  dy/dt
    ans(1) = -a(0);   //  d2y/dt2
    return ans;
}

MatrixXd RK4(VectorXd Func(const VectorXd& y), const Ref<const VectorXd>& y0, double h, int step_num)
{
    MatrixXd y(step_num, y0.rows());
    y.row(0) = y0;

    for (int i=1; i<step_num; i++){
        VectorXd y_old = y.row(i-1).transpose();
        VectorXd k1 = h*Func(y_old);
        VectorXd k2 = h*Func(y_old+k1/2);
        VectorXd k3 = h*Func(y_old+k2/2);
        VectorXd k4 = h*Func(y_old+k3);
        VectorXd dy = (k1 + 2*k2 + 2*k3 + k4)/6;
        y.row(i) = y.row(i-1) + dy.transpose();
    }
    return y;
}

int main()
{
    Vector2d v1;
    v1(0) = 1.4;    v1(1) = -0.1;
    double h = 0.1;
    int step_num = 50;

    MatrixXd sol = RK4(Func,v1,h,step_num);
    return 0;
}

У меня следующие вопросы:

1. Что означает & в аргументе функции? Пройти по ссылке? Я просто скопировал код из официальная документация, но я не уверен, понимаю ли я каждый бит в аргументах функций RK4, таких как VectorXd Func(const VectorXd& y). Существуют ли альтернативные способы принятия Eigen :: MatrixXd и функций, которые принимают Eigen :: MatrixXd в качестве аргументов функции?

2. Насколько я понимаю, мы не можем вернуть весь 2D-массив из функции, и мы возвращаем только первый элемент массива (поправьте меня, если я ошибаюсь). А как насчет Eigen::MatrixX? Что мы на самом деле проходим / возвращаем? Первый элемент матрицы или совершенно новый объект, определенный библиотекой Eigen?

3. Я не уверен, что эти коды написаны эффективно. Могу ли я что-нибудь сделать, чтобы оптимизировать эту часть? (Просто интересно, сделал ли я что-нибудь, что может значительно замедлить скорость).

Спасибо