 	double a = 1.0;
 	double b = 0.0;
 
-	culasStat = cublasDgemv(handle, CUBLAS_OP_N, m, n, 
+	cublasStat = cublasDgemv(handle, CUBLAS_OP_N, m, n, 
 							&a, mat, n, onesVec, 1, 
 							&b, sum, 1);
 

	double *partW1Cost, *partW2Cost;

	// allocate memeory space
	partW1Cost = malloc(sizeof(double));
	partW2Cost = malloc(sizeof(double));
	
	// copy valu for dvice to host
	cudaStat = cudaMemcpy(partW1Cost, totSumW1, sizeof(double), cudaMemcpyDeviceToHost);

/* --- SUBJECT TO CHANGE; TESTING STAGE --- */	
void squareMatrix(double *mat, int m, int n);
__global__ void squareElement(double *mat, int size);
void rowSum(cublasHandle_t handle, double *mat, int m, int n, double *sum);
void colSum(cublasHandle_t handle, double *mat, int m, int n, double *sum);
/* --- END SUBJECT TO CHANGE AREA --- */


