using namespace std;

int ymin, int ymax, float dx, float dy, float dxxinv, float dyyinv) {//MODIFY to suit algorithm

}

// Set grid size and number of iterations

const int save_iters = 20;

dim3 dBlock(m);

for (k = 1; k <= total_iters; k++) {

if (k % 2 == 0) {

}

}

cudaMemcpy(a, d_a, n* m * sizeof(float), cudaMemcpyDeviceToHost);

outputImage("jacobi out", a);