<div dir="ltr"><div><div>Hi all,<br><br></div>For the linear solver, if I use the standard 7-point differencing scheme and use 2-point boundary conditions, it works fine. But when I change the boundary condition, involving 3-point differencing form at the boundary, such as x[i+2]-2*x[i+1]+x[i]=0, assuming i is boundary point number, I always get the errors like <br>[0]PETSC ERROR: Argument out of range<br>[0]PETSC ERROR: Inserting a new nonzero into matrix<br><br>My boundary code is as follows.<br>if(i==0) {<br> v[0] =  1.0;    col[0].i = i;   <br> v[1] = -2.0;    col[1].i = i+1; <br> v[2] =  1.0;    col[2].i = i+2; <br> MatSetValuesStencil(appctx->A,1,&row,3,col,v,INSERT_VALUES); <br>}<br><br></div>Does anyone have some suggestions? Thank you very much in advance.<br clear="all"><div><div><div><br>-- <br><div class="gmail_signature"><div dir="ltr"><div>Best regards,</div><div><br></div>Feng</div></div>
</div></div></div></div>