<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Feb 17, 2014 at 5:42 AM, Torquil Macdonald Sørensen <span dir="ltr"><<a href="mailto:torquil@gmail.com" target="_blank">torquil@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div></div>How well is PETSC suited for solving elliptic PDEs on four-dimensional domains? Looking at the list of commands, I see e.g. DMDACreate{1d,2d,3d} with no 4d version. On the other hand, I see you have VecGetArray{1d,2d,3d,4d}.<br>

<br>Is there e.g. a convenient way of using the PETSCs data structures for solving the Poisson equation on a 4d "cube"?<br></div></div></blockquote><div><br></div><div>There is no 4d DA. The problem is that the routines for communication are not generic in dimension. I see at least</div>
<div>two avenues:</div><div><br></div><div>  1) Use DMPlex</div><div><br></div><div>       This would have the communication done, but none of the nice array extraction routines</div><div><br></div><div>  2) Write the 4d routines for communication setup</div>
<div><br></div><div>       This sounds like a fair amount of work (there have already been some failures), but its not hard to understand. We</div><div>       could help if this is what you want to do.</div><div><br></div>
<div>Are you looking to use 4d arrays to interact with it?</div><div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">
<div></div><div>Best regards<span class="HOEnZb"><font color="#888888"><br>Torquil Sørensen<br><br></font></span></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>