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>