gistfile1.txt

   double* derivs_p = get_derivatives();
    for (int i = 0; i < outdim; i++)
    {
        for (int j = 0; j < indim; j++)
        {
            derivs_p[i * indim + j] += source_p[i] * input_p[j];
        }
    }