Блог

Просто чтобы проверить, a умножил матрицу 150×150 всех единиц на себя, используя функцию gsl_blas_sgemm. Я получаю ошибки до 38. Почему это так неточно?

Вы уверены, что правильно это называете? У меня нет проблем:

 #include <stdio.h>
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_blas.h>

int main(int argc, char **argv) {
    gsl_matrix_float *a, *b;
    float err, maxerr;
    int n;

    if (argc != 2) {
        fprintf(stderr,"Usage: %s nntWhere n is size of matrix.n", argv[0]);
        exit(-1);
    }

    n = atoi(argv[1]);
    if (n<2 || n > 1000) {
        n = 100;
        fprintf(stderr,"%s: invalid n = %s; using n = %d.n",
                argv[0], argv[1], n );
    }

    a = gsl_matrix_float_alloc(n,n);
    b = gsl_matrix_float_alloc(n,n);

    for (int i=0; i<n; i  )
        for (int j=0; j<n; j  )
            gsl_matrix_float_set(a,i,j,1.);

    gsl_blas_sgemm (CblasNoTrans, CblasNoTrans, 1., a, a, 0., b);

    maxerr = 0.;
    for (int i=0; i<n; i  ) {
        for (int j=0; j<n; j  ) {
            err = gsl_matrix_float_get(b,i,j)-(float)n;
            if (err > maxerr) maxerr = err;
        }
    }

    printf("Maxerr = %fn", maxerr);

    gsl_matrix_float_free(a);
    gsl_matrix_float_free(b);
    return 0;
}
  

и затем

 $ gcc -o gslblas gslblas.c -lgsl -lgslcblas -std=c99
$ ./gslblas 150
Maxerr = 0.000000