我完全被难住了。我有一个用c编写的相当大的递归程序,它调用cblas_dgemm()。结果通过正常工作的程序独立验证。cblas_dgemm - 只有当(β)是2的幂时才有效
C = alpha*A*B + beta*C
在使用随机矩阵反复试验和的参数的所有可能的组合程序给出正确答案仅当ABS(测试版)= 2^N(1,2,4,8 ...)。任何值都适用于alpha。对于测试版,任何其他正面/负面,奇数/偶数值都会在10-30%的时间内给出正确的答案。
我使用的是Ubuntu 10.04,GCC 4.4.x,我试过系统安装了blas/cblas/atlas以及手动编译的地图集。
任何提示或建议将不胜感激。我惊讶于潜藏在这个地方的美妙慷慨(聪明)的人们。
感谢你的所有提前,
拉斯
你能提供你称为cblas_dgemm的确切代码吗?而且,您是否尝试直接打入Fortran常规DGEMM?我不记得在dgemm之前... – 2010-09-19 22:28:01