Hi guys, <div><br></div><div>I am setting up a matrix (parallel) using DA&#39;s in 3D. </div><div>DA is built upon a width=1, STAR_STENCIL type. </div><div>The nonzero structure of the resulting matrix is hepta-diagonal for the most part, however, since I am imposing boundary conditions using interpolation methods, for some boundary nodes, it is possible that some extra nonzeros are added into new locations. </div>
<div>This obviously reduces the matrix setup and assembly efficiency since the new nonzeros are not preallocated. </div><div>Is there any way that I could do it myself, i.e. preallocate the new nonzeros before matrix setup to get the best performance? </div>
<div><br></div><div>Thanks and have a great weekend, </div><div>Mohamad</div><div><br></div>