计算矩阵各元素指数的最有效方法

我正在从Matlab迁移到C+GSL,我想知道

B[i][j] = exp(A[i][j])

请注意,这与矩阵指数不同:

B = exp(A)

我刚刚找到了蛮力解决方案(几个for循环),但是 有没有更聪明的方法?

编辑

Drew Hall解决方案的结果

所有结果都来自1024x1024 for(for) 在每次迭代中 double 赋值(复数)。 .

• 当考虑{Row,Column}-主要模式来存储矩阵时的结果:
  • 223.22 ms,当以行主模式在内环中的列上循环时(情况2)。
  • gsl_matrix_complex_set

案例1的源代码 :

for(i=0; i<Nx; i++)
{
    for(j=0; j<Ny; j++)
    {
        /* Operations to obtain c_value (including exponentiation) */
        matrix[2*(i*s_tda + j)] = GSL_REAL(c_value);
        matrix[2*(i*s_tda + j)+1] = GSL_IMAG(c_value);
    }
}

案例2的源代码

for(i=0; i<Nx; i++)
{
    for(j=0; j<Ny; j++)
    {
        /* Operations to obtain c_value (including exponentiation) */
        matrix->data[2*(j*s_tda + i)] = GSL_REAL(c_value);
        matrix->data[2*(j*s_tda + i)+1] = GSL_IMAG(c_value);
    }
}

案例3的源代码 :

for(i=0; i<Nx; i++)
{
    for(j=0; j<Ny; j++)
    {
        /* Operations to obtain c_value (including exponentiation) */
        gsl_matrix_complex_set(matrix, i, j, c_value);
    }
}