<html><body><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; color: #000000"><div><div>Hi all,</div><div><br></div><div>First, I have to apologize for I lack some experience on this subject, I have tried looking around to solve my issue without much success.</div><div><br data-mce-bogus="1"></div><div>I am assembling a 3D Laplacian matrix on a DMDA using MatSetValuesStencil. The following code, which inserts one row at a time, works as expected:</div><div><br data-mce-bogus="1"></div><div>for (PetscInt iz = info.zs; iz < info.zs + info.zm; ++iz) {</div><div> for (PetscInt iy = info.ys; iy < info.ys + info.ym; ++iy) {</div><div> for (PetscInt ix = info.xs; ix < info.xs + info.xm; ++ix) {</div><div> const std::array<MatStencil, 1> rows{MatStencil{iz, iy, ix, 0}};</div><div> const std::array<MatStencil, 7> cols{{</div><div> {iz, iy, ix, 0},</div><div> {iz, iy, ix - 1, 0},</div><div> {iz, iy, ix + 1, 0},</div><div> {iz, iy - 1, ix, 0},</div><div> {iz, iy + 1, ix, 0},</div><div> {iz - 1, iy, ix, 0},</div><div> {iz + 1, iy, ix, 0}</div><div> }};</div><div><br></div><div> CHKERRQ(MatSetValuesStencil(matrix_A,</div><div> rows.size(), rows.data(),</div><div> cols.size(), cols.data(),</div><div> laplacian.data(), INSERT_VALUES));</div><div> }</div><div> }</div><div>}</div></div><div><br></div><div>How would you batch the whole x dimension as a single <!--StartFragment--><span style="color: #000000; font-family: arial, helvetica, sans-serif; font-size: 16px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: 2; text-align: start; text-indent: 0px; text-transform: none; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; white-space: normal; background-color: #ffffff; text-decoration-thickness: initial; text-decoration-style: initial; text-decoration-color: initial; display: inline !important; float: none;" data-mce-style="color: #000000; font-family: arial, helvetica, sans-serif; font-size: 16px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: 2; text-align: start; text-indent: 0px; text-transform: none; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; white-space: normal; background-color: #ffffff; text-decoration-thickness: initial; text-decoration-style: initial; text-decoration-color: initial; display: inline !important; float: none;">MatSetValuesStencil call ?</span></div><div><span style="color: #000000; font-family: arial, helvetica, sans-serif; font-size: 16px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: 2; text-align: start; text-indent: 0px; text-transform: none; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; white-space: normal; background-color: #ffffff; text-decoration-thickness: initial; text-decoration-style: initial; text-decoration-color: initial; display: inline !important; float: none;" data-mce-style="color: #000000; font-family: arial, helvetica, sans-serif; font-size: 16px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: 2; text-align: start; text-indent: 0px; text-transform: none; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; white-space: normal; background-color: #ffffff; text-decoration-thickness: initial; text-decoration-style: initial; text-decoration-color: initial; display: inline !important; float: none;"><br data-mce-bogus="1"></span></div><div><span style="color: #000000; font-family: arial, helvetica, sans-serif; font-size: 16px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: 2; text-align: start; text-indent: 0px; text-transform: none; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; white-space: normal; background-color: #ffffff; text-decoration-thickness: initial; text-decoration-style: initial; text-decoration-color: initial; display: inline !important; float: none;" data-mce-style="color: #000000; font-family: arial, helvetica, sans-serif; font-size: 16px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: 2; text-align: start; text-indent: 0px; text-transform: none; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; white-space: normal; background-color: #ffffff; text-decoration-thickness: initial; text-decoration-style: initial; text-decoration-color: initial; display: inline !important; float: none;">Each time I try, I endup getting errors such as:</span><!--EndFragment--><div style="clear: both;" data-mce-style="clear: both;"><div style="clear: both;" data-mce-style="clear: both;"> Inserting a new nonzero at global row/column (..., ...) into matrix</div><div><br></div></div></div><div><span style="color: #000000; font-family: arial, helvetica, sans-serif; font-size: 16px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: 2; text-align: start; text-indent: 0px; text-transform: none; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; white-space: normal; background-color: #ffffff; text-decoration-thickness: initial; text-decoration-style: initial; text-decoration-color: initial; display: inline !important; float: none;" data-mce-style="color: #000000; font-family: arial, helvetica, sans-serif; font-size: 16px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: 2; text-align: start; text-indent: 0px; text-transform: none; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; white-space: normal; background-color: #ffffff; text-decoration-thickness: initial; text-decoration-style: initial; text-decoration-color: initial; display: inline !important; float: none;"><br data-mce-bogus="1"></span></div><div><span style="color: #000000; font-family: arial, helvetica, sans-serif; font-size: 16px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: 2; text-align: start; text-indent: 0px; text-transform: none; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; white-space: normal; background-color: #ffffff; text-decoration-thickness: initial; text-decoration-style: initial; text-decoration-color: initial; display: inline !important; float: none;" data-mce-style="color: #000000; font-family: arial, helvetica, sans-serif; font-size: 16px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: 2; text-align: start; text-indent: 0px; text-transform: none; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; white-space: normal; background-color: #ffffff; text-decoration-thickness: initial; text-decoration-style: initial; text-decoration-color: initial; display: inline !important; float: none;">Regards,</span></div><div>-- <br></div><div><div><strong>Etienne Malaboeuf</strong><br>Ingénieur de recherche HPC, Département Calcul Intensif (DCI)<br><hr>Centre Informatique National de l'Enseignement Supérieur (CINES)<br>950 rue de Saint Priest, 34097 Montpellier<br>tel : (334) 67 14 14 02<br>web : <a href="https://urldefense.us/v3/__https://www.cines.fr__;!!G_uCfscf7eWS!dEs8i8p3qAHap5UIsgKxq7Nfj5KU_Mvp9mYj2fNWI5QwZEW8qEuMqPCDekMeZYZHxKcIdmfpIl-BxaWBE-6j-BIpOA$">https://www.cines.fr</a> | <a href="https://urldefense.us/v3/__https://dci.dci-gitlab.cines.fr/webextranet__;!!G_uCfscf7eWS!dEs8i8p3qAHap5UIsgKxq7Nfj5KU_Mvp9mYj2fNWI5QwZEW8qEuMqPCDekMeZYZHxKcIdmfpIl-BxaWBE-6NNu-h_Q$">https://dci.dci-gitlab.cines.fr/webextranet</a><br></div></div></div></body></html>