diff --git a/src/solver/linear_solver.cu b/src/solver/linear_solver.cu
index e658a4d0c..71995cf75 100644
--- a/src/solver/linear_solver.cu
+++ b/src/solver/linear_solver.cu
@@ -133,7 +133,7 @@ __global__ void SolveKernel(SolveDevice* device)
     
     //kernel call 
     size_t num_block = (denseMatrix.n_grids_ + BLOCK_SIZE - 1) / BLOCK_SIZE;
-    AddForcingTermsKernel<<<num_block, BLOCK_SIZE>>>(device);
+    Solve<<<num_block, BLOCK_SIZE>>>(device);
     cudaDeviceSynchronize();
 
     //clean up