Commits

Rio Yokota committed 74d5028

Fixed a few bugs, but convergence is still off for parallel run.

  • Participants
  • Parent commits 4d787ce

Comments (0)

Files changed (2)

   //nbodyGn(ni,xi,yi,ui,nb,xm,ym,um,dxe,dye,re,eps);
 
   ic = 0;
-  for (i=info.xs+1; i<info.xs+info.xm-1; i++) {
-    for (j=info.ys+1; j<info.ys+info.ym-1; j++,ic++) {
-      Y[j][i] = ui[ic];
+  for (i=info.xs; i<info.xs+info.xm; i++) { // Loop over interior points
+    for (j=info.ys; j<info.ys+info.ym; j++) {
+      if (i != 0 && i != info.mx-1 && j != 0 && j != info.my-1) {
+        Y[j][i] = ui[ic];
+        ic++;
+      }
     }
   }
 
 Poisson2D: Poisson2D.o fmm.o
 	${CLINKER} $^ -o $@ ${PETIGA_LIB} -L../../exafmm2d/wrappers -lfmm -ltbb
 	${RM} -f $<
-	mpirun -np 1 ./Poisson2D -ksp_monitor -ksp_norm_type UNPRECONDITIONED -draw_pause 1
+	mpirun -np 2 ./Poisson2D -ksp_monitor -ksp_norm_type UNPRECONDITIONED -ksp_max_it 100 -draw_pause 1
 
 clean::
 	-@${RM} ${TARGETS}