Eh bien, le titre dit, je fais suite genre d'opérations en Fortran:BLAS daxpy et dcopy quand inc = 1, sont-ils plus rapides que l'utilisation d'un (:, t) = b (:, t)?
a(:,t) = b(:,t)
c(:,t) = x(i,t)*d(:,t)
Y at-il des avantages de l'utilisation daxpy et dcopy de BLAS dans les sous-routines ce cas où inc = 1?
Aucune différence que vous auriez jamais remarqué. BLAS doit être compatible avec Fortran 77, que je suis sûr que ne possédait pas ces fonctionnalités de fantaisie.
Ces sous-routines sont là pour faire des tableaux ou des matrices en copiant une ligne de code, parce que c'est fait beaucoup. Les cycles ont tendance à s'embourber dans d'autres routines, comme la matrice inverse, donc la copie n'est généralement pas un problème de performance.
Si vous êtes préoccupé par les performances, il vous suffit de les coder de manière raisonnable. Alors ce que je ferais est interrupt it a few times. Cela vous montrera où va le temps. S'il passe beaucoup de temps à copier, il vous le dira. Sinon, il vous le dira.