<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
To see how this can be done with DM field for a PDE with DOF=2, you can refer to RHSJacobian() in src/ts/tutorials/advection-diffusion-reaction/ex5.c<br class="">
<div><br class="">
</div>
<div>Hong (Mr.)</div>
<div><br class="">
<blockquote type="cite" class="">
<div class="">On May 29, 2020, at 11:15 AM, Lucas Banting <<a href="mailto:bantingl@myumanitoba.ca" class="">bantingl@myumanitoba.ca</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class="">
Hello,</div>
<div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class="">
<br class="">
</div>
<div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class="">
I have a structured grid problem with 5 unknowns per cell: U,V,P,T,F, this is a CFD code but I don't think these details are necessary to answer my question.<span class="Apple-converted-space"> </span><br class="">
</div>
<div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class="">
I was wondering how exactly I should use MatSetValuesStencil() to fill my matrix.</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Verdana; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
I<span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""><span class="Apple-converted-space"> </span>have 5 degrees of freedom, so I need to set up to 25 coefficients per grid cell. I don't
understand how i,j, and c map to all coefficients.</span></div>
<div style="caret-color: rgb(0, 0, 0); font-family: Verdana; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""><br class="">
</span></div>
<div style="caret-color: rgb(0, 0, 0); font-family: Verdana; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">For a 1 degree of freedom system to fill my 9 point stencil I did:</span></div>
<div style="caret-color: rgb(0, 0, 0); font-family: Verdana; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""><br class="">
</span></div>
<div style="caret-color: rgb(0, 0, 0); font-family: Verdana; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""></span><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">
values(1)= ASW(II,JJ)</span><span class=""><br class="">
</span>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> values(2)= AS (II,JJ)</span><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> values(3)= ASE(II,JJ)</span><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> values(4)= AW (II,JJ)</span><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> values(5)= AP (II,JJ)</span><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> values(6)= AE (II,JJ)</span><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> values(7)= ANW(II,JJ)</span><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> values(8)= AN (II,JJ)</span><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> values(9)= ANE(II,JJ)</span><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> <span class="Apple-converted-space"> </span></span><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> idxm(MatStencil_i,1) = II;idxm(MatStencil_j,1) = JJ</span></div>
<div class=""><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> idxn(MatStencil_i,1) = II-1;idxn(MatStencil_j,1) = JJ-1</span><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> idxn(MatStencil_i,2) = II ;idxn(MatStencil_j,2) = JJ-1</span><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> idxn(MatStencil_i,3) = II+1;idxn(MatStencil_j,3) = JJ-1</span><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> idxn(MatStencil_i,4) = II-1 ;idxn(MatStencil_j,4) = JJ</span><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> idxn(MatStencil_i,5) = II ;idxn(MatStencil_j,5) = JJ</span><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> idxn(MatStencil_i,6) = II+1;idxn(MatStencil_j,6) = JJ</span><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> idxn(MatStencil_i,7) = II-1 ;idxn(MatStencil_j,7) = JJ+1</span><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> idxn(MatStencil_i,8) = II ;idxn(MatStencil_j,8) = JJ+1</span><br class="">
</div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> idxn(MatStencil_i,9) = II+1;idxn(MatStencil_j,9) = JJ+1</span><br class="">
</div>
<span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""> call MatSetValuesStencil(A,1,idxm,9,idxn,values,INSERT_VALUES,ierr)</span></div>
<div style="caret-color: rgb(0, 0, 0); font-family: Verdana; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""><br class="">
</span></div>
<span style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">Which
seemed to work just fine. For 5 degrees of freedom, instead of just having a single AP coefficient for example, I have an AP matrix:</span><span style="caret-color: rgb(0, 0, 0); font-family: Verdana; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class=""></span>
<div style="caret-color: rgb(0, 0, 0); font-family: Verdana; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""><br class="">
</span></div>
<div style="caret-color: rgb(0, 0, 0); font-family: Verdana; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">
<table cellspacing="0" cellpadding="1" style="border-collapse: collapse;" class="">
<tbody class="">
<tr style="background-color: rgb(255, 255, 255);" class="">
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_uu</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_uv</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_up</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_ut</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_uf</span><br class="">
</td>
</tr>
<tr style="background-color: rgb(255, 255, 255);" class="">
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_vu</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_vv</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_vp</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_vt</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_vf</span><br class="">
</td>
</tr>
<tr style="background-color: rgb(255, 255, 255);" class="">
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_pu</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_pv</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_pp</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_pt</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_pf</span><br class="">
</td>
</tr>
<tr style="background-color: rgb(255, 255, 255);" class="">
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_tu</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_tv</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_tp</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_tt</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_tf</span><br class="">
</td>
</tr>
<tr style="background-color: rgb(255, 255, 255);" class="">
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_fu</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_fv</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_fp</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_ft</span><br class="">
</td>
<td style="width: 100px; border: 1px solid rgb(171, 171, 171);" class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">ap_ff</span><br class="">
</td>
</tr>
</tbody>
</table>
</span></div>
<div style="caret-color: rgb(0, 0, 0); font-family: Verdana; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<br class="">
</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Verdana; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">In 1 degree of freedom, i and j corresponded to the solution variable in the grid. For multi degree of freedom, I don't
understand how values of i, j, and c could distinguish from an ap_uv and ap_vv coefficent for example, wouldn't they both be at i, j ,c=2?<br class="">
</span></div>
<div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class=""></span></div>
<span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">Is there anyway I can use MatSetValuesStencil() to fill in my 9 point stencil with my 5x5 matrix coefficients? To clarify I have
nine 5x5 matrices for each cell which correspond to the 5 unknowns per cell.<br class="">
</span></div>
<div style="caret-color: rgb(0, 0, 0); font-family: Verdana; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<br class="">
</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Verdana; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">Thanks,</span></div>
<div style="caret-color: rgb(0, 0, 0); font-family: Verdana; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<br class="">
</div>
<div style="caret-color: rgb(0, 0, 0); font-family: Verdana; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; background-color: rgba(0, 0, 0, 0);" class="">Lucas</span></div>
</div>
</blockquote>
</div>
<br class="">
</body>
</html>