Fix E2PFromCSR

Author: isuruf

sumpy/e2p.py

@@ -330,10 +330,15 @@ def get_kernel(self, max_ntargets_in_one_box):
                               if=run_itgt}
                         end
                     end
-                    result[iknl, itgt] = result[iknl, itgt] + result_temp[ \
+                    for itgt_offset
+                        <> itgt2 = itgt_start + itgt_offset {id=init_itgt_for_write}
+                        <> run_itgt2 = itgt_start + itgt_offset < itgt_end  \
+                                {id=init_cond_for_write}
+                        result[iknl, itgt2] = result[iknl, itgt2] + result_temp[ \
                             itgt_offset, iknl] * kernel_scaling \
-                        {dep=update_result:init_result,id=write_result, \
-                         dup=iknl,if=run_itgt}
+                            {dep=update_result:init_result,id=write_result, \
+                             dup=iknl,if=run_itgt2}
+                    end
                 end
                 """],
                 [
@@ -386,7 +391,8 @@ def get_optimized_kernel(self, max_ntargets_in_one_box):
         knl = lp.privatize_temporaries_with_inames(knl,
             "itgt_offset_outer", "result_temp")
         knl = lp.duplicate_inames(knl, "itgt_offset_outer", "id:init_result")
-        knl = lp.duplicate_inames(knl, "itgt_offset_outer", "id:write_result")
+        knl = lp.duplicate_inames(knl, "itgt_offset_outer",
+            "id:write_result or id:init_itgt_for_write or id:init_cond_for_write")
         knl = lp.add_inames_to_insn(knl, "dummy",
             "id:init_box* or id:fetch_src_box or id:fetch_center "
             "or id:kernel_scaling")