<div dir="ltr">Hello!<div><br></div><div>I have an MPI code in which a linear system is created and solved with PETSc. It works in sequential run but when I use multiple cores the VecAssemblyBegin/End give segmentation fault. Here's a sample of my code:</div><div><br></div><div><div>call PetscInitialize(PETSC_NULL_CHARACTER,perr)</div><div>      </div><div>      ind(1) = 3*nnod_loc*max_time_deg</div><div>      call VecCreate(PETSC_COMM_WORLD,feP,perr)</div><div>      call VecSetSizes(feP,PETSC_DECIDE,ind,perr)</div><div>      call VecSetFromOptions(feP,perr)</div><div>      </div><div>      do in = nnod_loc</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">           </span>do jt = 1,mm</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">                   </span>ind(1) = 3*((in -1)*max_time_deg + (jt-1))</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">                     </span>fval(1) = fe(3*((in -1)*max_time_deg + (jt-1)) +1)</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">                     </span>call VecSetValues(feP,1,ind,fval(1),INSERT_VALUES,perr)</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">                        </span>ind(1) = 3*((in -1)*max_time_deg + (jt-1)) +1</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">                  </span>fval(1) = fe(3*((in -1)*max_time_deg + (jt-1)) +2)</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">                     </span>call VecSetValues(feP,1,ind,fval(1),INSERT_VALUES,perr)</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">                        </span>ind(1) = 3*((in -1)*max_time_deg + (jt-1)) +2</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">                  </span>fval(1) = fe(3*((in -1)*max_time_deg + (jt-1)) +3)</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">                     </span>call VecSetValues(feP,1,ind,fval(1),INSERT_VALUES,perr)</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">                </span> enddo</div><div>         enddo</div><div>      enddo</div><div>      call VecAssemblyBegin(feP,perr)</div><div>      call VecAssemblyEnd(feP,perr)<br></div></div><div><br></div><div>The vector has 640.000 elements more or less but I am running on a high performing computer so there shouldn't be memory issues. Does anyone know where is the problem and how can I fix it?</div><div><br></div><div>Thank you,</div><div>Francesco Migliorini</div></div>