gfortran

1열

1답변

BLAS 서브 루틴 dgemm, dgemv 및 ddot가 스칼라에서 작동하지 않습니까?

행렬 * 행렬, 행렬 * 벡터 및 벡터 * 벡터를 계산하는 BLAS의 서브 루틴 dgemm, dgemv 및 ddot를 사용하는 포트란 서브 루틴이 있습니다. 나는 m * m 행렬과 m * 1 벡터를가집니다. 어떤 경우에는 m = 1이다. 이러한 서브 루틴은 이러한 경우 제대로 작동하지 않는 것으로 보입니다. 그들은 오류를주지 않지만 결과에 수치 적 불안정성