gsl 复矩阵 * 复向量

Question

谁能帮我弄清楚如何做这个操作。我正在查看文档，但自从我不得不做任何线性代数类型的东西以来，我有点迷茫了。

我有一个名为 Y 的 14x14 复数矩阵和一个名为 I 的复数向量。我需要将它们相乘并将结果设置为一个名为 IL 的复数向量。

到目前为止，我已经弄清楚我需要使用：

gsl_blas_zgemv (CBLAS_TRANSPOSE_t TransA, const gsl_complex alpha, const 
                 gsl_matrix_complex * A, const gsl_vector_complex * x, const  
                 gsl_complex beta, gsl_vector_complex * y)

但我不知道去哪里了。不知道去哪里。像这样的东西？但是什么是 alpha 和 beta？

gsl_blas_zgemv(CblasNoTrans, ???, &Y, &I, ???, IL); 

Answer 1

没有看到你实际在做什么，看起来你没有包含正确的标题。因为以下对我来说没有错误：

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#include <gsl/gsl_blas.h>
#include <gsl/gsl_complex_math.h>

int
main(int argc, char **argv)
{

    int size = 14;
    gsl_matrix_complex *A = NULL;
    gsl_vector_complex *x = NULL;
    gsl_vector_complex *y = NULL;

    A = gsl_matrix_complex_alloc(size, size);
    x = gsl_vector_complex_alloc(size);
    y = gsl_vector_complex_alloc(size);

    gsl_matrix_complex_set_all(A, GSL_COMPLEX_ONE);
    gsl_vector_complex_set_all(x, GSL_COMPLEX_ONE);
    gsl_vector_complex_set_all(y, GSL_COMPLEX_ZERO);

    gsl_blas_zgemv(CblasNoTrans, GSL_COMPLEX_ONE, A, x,
            GSL_COMPLEX_ZERO, y);

    return(EXIT_SUCCESS);

}

并编译

gcc -o test -I/opt/local/include/gsl -L/opt/local/lib -lgsl -lgslcblas test.c

（是的，我在使用 MacPorts 的 Mac 上）。