<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
On 01.11.2011 16:14, Jed Brown wrote:
<blockquote
cite="mid:CAM9tzSmqTvzn1+wPrq0-2Ap8fyUHtL40zivQ-8UQY9-zAW9hxw@mail.gmail.com"
type="cite">
<div class="gmail_quote">On Tue, Nov 1, 2011 at 09:08, Alexander
Grayver <span dir="ltr"><<a moz-do-not-send="true"
href="mailto:agrayver@gfz-potsdam.de">agrayver@gfz-potsdam.de</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt
0.8ex; border-left: 1px solid rgb(204, 204, 204);
padding-left: 1ex;">
<div id=":wa">I solve 3D vector PDE using FD, so my solution
vector b has size 3N, where N stands for total number of
grid cells.<br>
Each part of b contains solution for corresponding
dimension, i.e. fx = b(1:N), fy = b(N+1:2N), fz =
b(2N+1:3*N).<br>
</div>
</blockquote>
<div><br>
</div>
<div>You set the number of degrees of freedom per node to 3.
Then these vectors are interlaced in b.</div>
</div>
</blockquote>
<br>
Thanks Jed, looks easy. I've read about this parameter in
documentation, but it wasn't obvious how to use it. <br>
May be it's worth to mention it shortly in "2.4.4 Accessing the
Vector Entries for DMDA Vectors" or somewhere else.<br>
<br>
How then can I ensure correct ordering between my system matrix and
vectors in case of dof > 1? <br>
Can I say that dof_i always corresponds to g((i-1)*N+1:i*N) in
global DM vector?<br>
<br>
Regards,<br>
Alexander<br>
</body>
</html>