改章節筆者在廣東游玩的時候突然想到的...明天就有想寫幾篇關于函數說明的文章，所以回家到之后就奮筆疾書的寫出來發布了

????

/*=======================================================================

????

* 第二天了，算打弄完這個去吃劈柴院的灌湯包，海上起霧了，要是不冷了，去趟只聽過沒去過的棧橋

????

* 金剛狼3上映有一段時間了吧，晚上去試試。。

????

=======================================================================*/

????我看過cublas的說明文檔，它在解講的時候先從總體上說明cublas的基本存儲則規，所以這里，我說幾點，這些容內將在實際中體具習練。

????

????

2.1 實際之前的cublas幾點說明

????

• cublas實現三個別級的BLAS（Basic Linear Algebra Subprograms）的函數能功，關于BLAS的說明，大家可以在 這里 找到。
• 本節中述講blas1中的某個函數的應用。
• cublas接受的參數，按照存儲點地分，有三種：設備上，主機上，設備或者主機上。按照數學類型分：量標（對應單變量，常量），量向（對應數組），矩陣（對應數組）。
• cublas中量向，矩陣 從1開始引索 。這點要注意的是，cublas自己的函數這樣實現，我們自己的數據還是按照自己的語言標準來做的，就是說，你用C語言的數組，還是從0開始，但是用cublas庫的時候，庫本身自己從1開始。

????

????

2.2 境環外額配置

????上一篇中紹介了如何配置cublas的編譯境環，不過要應用cublas，還需要cuda的runtime支撐，所以你需要在鏈接器的依附項中添加cudart.lib的外額依附。體具添加法方，與上一篇中法方相似，不再贅述。

????另外，源碼中應用cuda runtime函數的文件中要添加#include “cuda_runtime.h”

????

????

2.3 驗實函數說明

????明天我們要驗實的庫函數是cublasIsamax(...)，它要實現的能功是得獲某個量向中第一個最大模的引索號，注意這個引索號是從1開始的。

????函數型原是：cublasStatus_t cublasIsamax(cublasHandle_t handle, int n,?const float *x, int incx, int *result)

????參數說明：

????

• 回返值，含包函數的誤錯碼代，這是個枚舉類型，值為CUBLAS_STATUS_SUCCESS時，時常正，否則說明發生了誤錯。
• handle，表現是cublas的上下文句柄。由cublasStatus_t?cublasCreate(cublasHandle_t *handle)函數創立，由cublasStatus_t?cublasDestroy(cublasHandle_t handle)釋放。
• n，表現量向x中要引索的元素的總個數。
• x，量向x。他的存儲 必須存儲在GPU上 。
• incx，要引索的x中元素的增量。
• result，寄存結果。他的存 儲可以在CPU上也可以在GPU上 。

????數學上解釋下算法：

????此函數回返第一個i，使得 |Im(x[j])| + |Re(x[j])| 為最大值，其中i=1,...,n，并且，j=1 + (i-1)*incx.

????

????

2.4 基本驗實

????碼代我貼在了上面

    #include <stdio.h>

#include <stdlib.h>

#include <time.h>



#include "cuda_runtime.h"

#include "cublas_v2.h"

#define DATA_SIZE 1024



#define CUBLAS_ERROR_CHECK(sdata) if(CUBLAS_STATUS_SUCCESS!=sdata){printf("ERROR at:%s:%d\n",__FILE__,__LINE__);exit(-1);}



void InitVector(float *vec, int n, int down, int up)

{

	int i;

	srand((unsigned int)time(NULL));

	for (i = 0; i < n; i++)

	{

		vec[i] = (float)(rand()%(up-down+1)+down);

	}

	return ;

}

void myIsamax(int n, float *x, int incx, int *result)

{

	int i = 0;

	*result = 0;

	for (i = 1; i < n; i+=incx)

	{

		if(x[i]>x[*result]) *result=i; 

	}

}

int main(int argc, char **argv)

{

	cublasHandle_t cubhandle;

	cublasStatus_t cubStatus = CUBLAS_STATUS_SUCCESS;



	float *h_datax;

	float *d_datax;

	int    h_result;

	int    d_result;

	int    data_size = DATA_SIZE;

	

	// initialize cublas resource

	cubStatus = cublasCreate(&cubhandle); CUBLAS_ERROR_CHECK(cubStatus)



	// allocate memory CPU + GPU

	h_datax = (float*)malloc(data_size*sizeof(float));

	cudaMalloc((void**)&d_datax, data_size*sizeof(float));



	// initialize cpu memory

	InitVector(h_datax, data_size, 0, 100);



	// copy cpu value to gpu memory

	cubStatus = cublasSetVector(data_size, sizeof(float), h_datax, 1, d_datax, 1);

	CUBLAS_ERROR_CHECK(cubStatus)

	

	// execute GPU version

	cubStatus = cublasIsamax(cubhandle, data_size, d_datax, 1, &d_result); CUBLAS_ERROR_CHECK(cubStatus)

	// execute CPU version



	myIsamax(data_size, h_datax, 1, &h_result);





	printf("host   : h_datax[%d]=%.2f\n", h_result, h_datax[h_result]);

	printf("device : d_datax[%d]=%.2f\n", d_result-1, h_datax[d_result-1]);





	// Free memory

	free(h_datax);h_datax=0x0;

	cudaFree(d_datax);d_datax=0x0;



	cublasDestroy(cubhandle);

	getchar();

	return 0;

}
  

????總結

????為盡的地方，歡送大家批駁針指。

????

????

????

????

????

????

文章結束給大家分享下程序員的一些笑話語錄： 這年頭的互聯網真是娛樂了中國，網民們從各種各樣的“門”里鉆來鉆去，又有好多“哥”好多“帝”，值得大家品味不已……網絡經典語錄，關于IT與互聯網，經典與您分享！

函數說明玩玩CUBLAS(2)——level1函數