Why physicists still use Fortran (2015)(moreisdifferent.com)

    real, dimension(:,:), allocatable :: name_of_array
    allocate(name_of_array(xdim, ydim))

    auto *my_matrix = new MathArray<double, 10, 10>();

    typedef struct {
        size_t rows, cols;
        void values[];
    } mat;

    inline mat *matrix_make(size_t rows, size_t cols, size_t value_size) {
        mat *m = malloc(sizeof(size_t) + sizeof(size_t) + (rows * cols * value_size));
        m->rows = rows;
        m->cols  = cols;
        memset(m->values , 0, rows * cols * value_size);
        return m;
    }

    mat *matrix = matrix_make(10, 10, sizeof(double));

    #define MATRIX_MAKE(r, c, type) matrix_make(r, c, sizeof(type))

   1. We already know it
   2. We're not going to get grant money to rewrite a library
   3. We've built a bunch of computer clusters and have to justify what we spent (rather than buying 2 GPUs for your workstation)
   4. We've all spent a lot of time learning how to use MPI that we're never getting back.