<div dir="ltr"><div><div><span style="font-family:courier new,monospace">Hi all,<br><br></span></div><div><span style="font-family:courier new,monospace">I'm trying to solve the Navier-Stokes equations on a staggered grid. I've created three DA vectors U, V and P to store the values of the x- and y-components of velocity and the pressure. Each of these has different global dimensions. I have assembled matrices Gx and Gy to calculate each component of the gradient of the pressure field:<br>


</span></div><div><span style="font-family:courier new,monospace"><br>Gx * P = Rx<br></span></div><div><span style="font-family:courier new,monospace">Gy * P = Ry<br><br></span></div><div><span style="font-family:courier new,monospace">Rx and Ry have the same lengths and distributed structures of U and V. <br>


<br>I would like to create matrix G and vector R such that<br><br></span></div><div><span style="font-family:courier new,monospace">G*P = R<br></span></div><div><span style="font-family:courier new,monospace"></span></div>


<div><span style="font-family:courier new,monospace"><br>where<br><br>G = /Gx\<br></span></div><span style="font-family:courier new,monospace">    \Gy/<br></span></div><span style="font-family:courier new,monospace">and</span><br>


<span style="font-family:courier new,monospace"></span><div><span style="font-family:courier new,monospace"><span style="font-family:courier new,monospace"><br>R = /Rx\<br></span><span style="font-family:courier new,monospace">    \Ry/<br>


<br></span></span></div><div><span style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">and vector R is created from the DMComposite of the DAs of U and V.<br></span></span></div><div>


<span style="font-family:courier new,monospace"><br></span></div><div><span style="font-family:courier new,monospace">What would be the best way to go about doing this?<br><br></span></div><div><span style="font-family:courier new,monospace">Thank you,<br>


Anush<br></span></div></div>