h******6 发帖数: 5 | 1 本来矩阵向量运算的code都是自己乱写的, 现在想学用一下 blas lapack 这样的标准
库。 所以先找了一段 ATLAS 的源代码, 不是很懂, 有牛人能解释下下面的这段 代
码么? 特别是为什么 Mjoin 函数名后面能跟两个参数列表? 多谢!
#include "atlas_misc.h"
#include "atlas_level1.h"
#ifdef TREAL
void Mjoin(PATL,rot)(const int N, TYPE *X, const int incX,
TYPE *Y, const int incY, const TYPE c, const TYPE s)
{
int i;
TYPE tmp;
if (c != ATL_rone || s != ATL_rzero)
{
if (incX == 1 && incY == 1)
{
for (i=0; i != N; i++)
{
tmp = c * |
|