<div><div class="gmail_msg">You need to swap the order of your function calls.</div><div class="gmail_msg">Call VecSetSizes() before VecSetType()</div><div class="gmail_msg"><br class="gmail_msg"></div><div class="gmail_msg">Thanks,</div><div class="gmail_msg">  Dave</div></div><div><div class="gmail_msg"><br class="gmail_msg"></div><div class="gmail_msg"><br class="gmail_msg"><div class="gmail_quote gmail_msg"><div class="gmail_msg">On Wed, 4 Jan 2017 at 23:21, Manuel Valera <<a href="mailto:mvalera@mail.sdsu.edu" class="gmail_msg" target="_blank">mvalera@mail.sdsu.edu</a>> wrote:<br class="gmail_msg"></div><blockquote class="gmail_quote gmail_msg" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="gmail_msg">Hello all, happy new year,<div class="gmail_msg"><br class="gmail_msg"></div><div class="gmail_msg">I'm working on parallelizing my code, it worked and provided some results when i just called more than one processor, but created artifacts because i didn't need one image of the whole program in each processor, conflicting with each other. </div><div class="gmail_msg"><br class="gmail_msg"></div><div class="gmail_msg">Since the pressure solver is the main part i need in parallel im chosing mpi to run everything in root processor until its time to solve for pressure, at this point im trying to create a distributed vector using either </div><div class="gmail_msg"><br class="gmail_msg"></div><div class="gmail_msg"><p style="margin:0px;font-size:11px;line-height:normal;font-family:menlo" class="gmail_msg"><span style="font-variant-ligatures:no-common-ligatures" class="gmail_msg">     </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(52,187,199)" class="gmail_msg">call</span><span style="font-variant-ligatures:no-common-ligatures" class="gmail_msg"> VecCreateMPI(PETSC_COMM_WORLD,PETSC_DECIDE,nbdp,xp,ierr)</span></p></div><div class="gmail_msg"><span style="font-variant-ligatures:no-common-ligatures" class="gmail_msg">or</span></div><div class="gmail_msg"><span style="font-variant-ligatures:no-common-ligatures" class="gmail_msg"><p style="margin:0px;font-size:11px;line-height:normal;font-family:menlo" class="gmail_msg"><span style="font-variant-ligatures:no-common-ligatures;color:rgb(52,187,199)" class="gmail_msg">     call</span><span style="font-variant-ligatures:no-common-ligatures" class="gmail_msg"> VecCreate(PETSC_COMM_WORLD,xp,ierr); CHKERRQ(ierr)</span></p><p style="margin:0px;font-size:11px;line-height:normal;font-family:menlo;color:rgb(83,48,225)" class="gmail_msg"><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)" class="gmail_msg">     </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(52,187,199)" class="gmail_msg">call</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(0,0,0)" class="gmail_msg"> VecSetType(xp,VECMPI,ierr)  </span></p><p style="margin:0px;font-size:11px;line-height:normal;font-family:menlo" class="gmail_msg"><span style="font-variant-ligatures:no-common-ligatures" class="gmail_msg">     </span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(52,187,199)" class="gmail_msg">call</span><span style="font-variant-ligatures:no-common-ligatures" class="gmail_msg"> VecSetSizes(xp,PETSC_DECIDE,nbdp,ierr); CHKERRQ(ierr)</span><br class="gmail_msg"></p><p style="margin:0px;font-size:11px;line-height:normal;font-family:menlo" class="gmail_msg"><span style="font-variant-ligatures:no-common-ligatures" class="gmail_msg"><br class="gmail_msg"></span></p><p style="margin:0px;font-size:11px;line-height:normal;font-family:menlo" class="gmail_msg"><span style="font-variant-ligatures:no-common-ligatures" class="gmail_msg"><br class="gmail_msg"></span></p><p style="margin:0px;line-height:normal" class="gmail_msg"><span style="font-variant-ligatures:no-common-ligatures" class="gmail_msg"><font face="arial, helvetica, sans-serif" class="gmail_msg">In both cases program hangs at this point, something it never happened on the naive way i described before. I've made sure the global size, nbdp, is the same in every processor. What can be wrong?</font></span></p><p style="margin:0px;line-height:normal" class="gmail_msg"><span style="font-variant-ligatures:no-common-ligatures" class="gmail_msg"><font face="arial, helvetica, sans-serif" class="gmail_msg"><br class="gmail_msg"></font></span></p><p style="margin:0px;line-height:normal" class="gmail_msg"><span style="font-variant-ligatures:no-common-ligatures" class="gmail_msg"><font face="arial, helvetica, sans-serif" class="gmail_msg">Thanks for your kind help,</font></span></p><p style="margin:0px;line-height:normal" class="gmail_msg"><span style="font-variant-ligatures:no-common-ligatures" class="gmail_msg"><font face="arial, helvetica, sans-serif" class="gmail_msg"><br class="gmail_msg"></font></span></p><p style="margin:0px;line-height:normal" class="gmail_msg"><span style="font-variant-ligatures:no-common-ligatures" class="gmail_msg"><font face="arial, helvetica, sans-serif" class="gmail_msg">Manuel.</font></span></p></span></div></div><br class="gmail_msg"><br class="gmail_msg"></blockquote></div></div></div>