fixing bug

NCAR · Oct 20, 2023 · ab01538 · ab01538
1 parent 782de06
commit ab01538
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/solver/linear_solver.cu b/src/solver/linear_solver.cu
@@ -119,7 +119,7 @@ __global__ void SolveKernel(SolveDevice* device)
     cudaMemcpy(d_lower_matrix, sparseMatrix.lower_matrix_, sizeof(double)*sparseMatrix.lower_matrix_size_, cudaMemcpyHostToDevice);
     cudaMemcpy(d_upper_matrix, sparseMatrix.upper_matrix_, sizeof(double)* sparseMatrix.upper_matrix_size_, cudaMemcpyHostToDevice);
     cudaMemcpy(d_b, denseMatrix.b_, sizeof(double)* denseMatrix.n_grids_* denseMatrix.b_column_counts_, cudaMemcpyHostToDevice);
-    cudaMemcpy(d_x, denseMatrix.x_, sizeof(double)* denseMatrix.n_grids * denseMatrix.x_column_counts_, cudaMemcpyHostToDevice);
+    cudaMemcpy(d_x, denseMatrix.x_, sizeof(double)* denseMatrix.n_grids_ * denseMatrix.x_column_counts_, cudaMemcpyHostToDevice);
 
     cudaMemcpy(&(device->nLij_Lii_), &d_nLij_Lii, sizeof(std::pair<size_t, size_t>*),cudaMemcpyHostToDevice);
     cudaMemcpy(&(device->Lij_yj_), &d_Lij_yj, sizeof(std::pair<size_t, size_t>*), cudaMemcpyHostToDevice);