由买买提看人间百态

topics

全部话题 - 话题: mjoin
(共0页)
h******6
发帖数: 5
1
来自主题: Computation版 - 一段 ATLAS 的源代码
本来矩阵向量运算的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 *
a**a
发帖数: 416
2
来自主题: Computation版 - 一段 ATLAS 的源代码
想通过看代码来学习复杂算法,就像通过拆发动机来学习发动机原理一样,都是很吃力
的事情。建议你找本专门讲线性代数算法的书,比如著名的Matrix Computation. 这样
可以省事很多。
这段代码嘛, Mjoin不是函数,而是一个宏。它根据给出的参数来生成函数名字。
(共0页)