Матрица GSL gsl_blas_sgemm умножается очень неточно при n> 140 - PullRequest
Новая
       28

Матрица GSL gsl_blas_sgemm умножается очень неточно при n> 140

1 голос
/ 08 мая 2011

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

1 Ответ

4 голосов
/ 09 мая 2011

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

#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 n\n\tWhere 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 = %f\n", 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
...