     if (y >= height-1)
       p[threadIdx.x+1][threadIdx.y+1].y = 0.0f;
 
-//        syncthreads();
+    //        syncthreads();
 
     float divergence = p[threadIdx.x+1][threadIdx.y+1].x - p[threadIdx.x][threadIdx.y+1].x +
-                       p[threadIdx.x+1][threadIdx.y+1].y - p[threadIdx.x+1][threadIdx.y].y;
+        p[threadIdx.x+1][threadIdx.y+1].y - p[threadIdx.x+1][threadIdx.y].y;
 
     // update primal variable
     u = (u + tau*(divergence + lambda*f))/(1.0f+tau*lambda);

    else
      theta = 1.0f;

    update_shared_single<<<dimGrid, dimBlock>>>(device_u->data(), device_u_->data(), device_p->data(),
                                                tau, theta, sigma, lambda, width, height,
                                                device_u->stride(), device_p->stride());
    sigma /= theta;

    if (threadIdx.y==BSY-1)
      u_[threadIdx.x][threadIdx.y+1] = tex2D(rof_tex_u_, x+0.5f, y+1.5f);

    //    syncthreads();

    float2 grad_u = make_float2(u_[threadIdx.x+1][threadIdx.y] - u_[threadIdx.x][threadIdx.y],
                                u_[threadIdx.x][threadIdx.y+1] - u_[threadIdx.x][threadIdx.y] );
