Not sure what you are asking. Why not just think of the grid as being transposed?<br>Something has to be numbered first. If you had a more sophisticated change, I<br>would recommend using an AO to renumber the matrix. This is what happens<br>
when you go from PETSc ordering to Natural ordering.<br><br> Matt<br><br><div class="gmail_quote">On Tue, Feb 2, 2010 at 9:18 AM, ilyas yilmaz <span dir="ltr"><<a href="mailto:ilyasy@chalmers.se">ilyasy@chalmers.se</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><br>
Hi,<br>
<br>
As far as I understand, DAs create matrix structure which is based on<br>
natural ordering of mesh, ie.,<br>
<br>
(i,j)<br>
cell(0,0) --> 0th row of DA matrix<br>
cell(1,0) --> 1st row of DA matrix<br>
cell(2,0) --> 2nd row of DA matrix and so on<br>
<br>
and for coefficients : (i,j-1) (i-1,j) (i,j) (i+1,j) (i,j+1)<br>
<br>
My question is :<br>
----------------<br>
I need such a matrix stucture for my application that<br>
<br>
(i,j)<br>
cell(0,0) --> 0th row of DA matrix<br>
cell(0,1) --> 1st row of DA matrix<br>
cell(0,2) --> 2nd row of DA matrix and so on<br>
<br>
and for coefficients : (i-1,j) (i,j-1) (i,j) (i,j+1) (i+1,j)<br>
<br>
How can I tell this DAs?<br>
<br>
Note:<br>
To simplify and to stay away from processor ordering issue,<br>
I assume that I am running DAs code on 1 cpu.<br>
<br>
<br>
Thank you for providing PETSc library and for your help.<br>
<br>
Regards,<br>
<font color="#888888">Ilyas<br>
<br>
</font></blockquote></div><br><br clear="all"><br>-- <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<br>