<div dir="ltr">Hello, <div>In the MPI context, after assembling the distributed matrix A (matmpiaij) and the right-hand-side b, I am trying to apply the 1st kind boundary condition using MatZeroRows() and VecSetValues(), for A and b respectively. </div><div>The pseudo-code is:</div><div>=========</div><div><font face="arial, sans-serif">for (<span style="background-color:rgb(255,255,255)">int key = 0; key < BCNodes_Length; key++){</span></font></div><div><span style="background-color:rgb(255,255,255)">      // retrieving the global row position</span></div>      pos = BCNodes[key];  <div><span style="background-color:rgb(255,255,255)">      // Set all elements in that row 0 except the one on the diagonal to be 1.0</span></div><div>      MatZeroRows(A, 1, &pos, 1.0, NULL, NULL);      </div><div><span style="background-color:rgb(255,255,255)">}</span></div><div><div style="line-height:19px">MatAssemblyBegin(A, MAT_FINAL_ASSEMBLY);</div><div style="line-height:19px">MatAssemblyEnd(A, MAT_FINAL_ASSEMBLY);</div><div style="line-height:19px">=========</div><div><br></div>For BCNodes_Length = 10^4, the FOR loop timing is 8 seconds.<br>For BCNodes_Length = 15*10^4, the FOR loop timing is 3000 seconds.<br>I am using two computational nodes and each having 12 cores.<br><br>My questions are:<br>1) Is the timing plausible? Is the MatZeroRows() function so costly?<br>2) Any suggestions to apply the 1st kind boundary conditions for a better performance?<br><br>Thanks,<br>Qiyue Lu</div></div>