Thank you very much, Satish! Ill try it<div><br></div><div>Alexey<br><br><div class="gmail_quote">8 апреля 2011 г. 3:32 пользователь Satish Balay <span dir="ltr">&lt;<a href="mailto:balay@mcs.anl.gov">balay@mcs.anl.gov</a>&gt;</span> написал:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div><div></div><div class="h5">On Fri, 8 Apr 2011, Алексей Рязанов wrote:<br>
<br>
&gt; Hello.<br>
&gt;<br>
&gt; When I create vectors using<br>
&gt;<br>
&gt; VecCreate(PETSC_COMM_WORLD,&amp;u);<br>
&gt; VecSetSizes(u,PETSC_DECIDE, VecSize);<br>
&gt; VecSetFromOptions(u);<br>
&gt; VecDuplicate(u,&amp;b);<br>
&gt;<br>
&gt; and matrix using<br>
&gt;<br>
&gt; MatCreate(PETSC_COMM_WORLD,&amp;A);<br>
&gt; MatSetSizes(A,PETSC_DECIDE,PETSC_DECIDE,VecSize,VecSize);<br>
&gt; MatSetFromOptions(A);<br>
&gt;<br>
&gt; PETSc distributes their elements in a proper identical way among processors,<br>
&gt; so I can use procedures like<br>
&gt;<br>
&gt; MatMult(A,u,b);<br>
&gt;<br>
&gt; and<br>
&gt;      KSPSolve(ksp,b,x);<br>
&gt; Ofcourse after matrix assembling and initialization of KSP and PC<br>
&gt;<br>
&gt; KSPCreate(PETSC_COMM_WORLD,&amp;ksp);<br>
&gt; KSPSetOperators(ksp,A,A,DIFFERENT_NONZERO_PATTERN);<br>
&gt;<br>
&gt; And thats great and works amazingly!<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; But now I&#39;ve created DA vectors &quot;u&quot; and &quot;b&quot; and assembled them through the<br>
&gt; natural grid indexing.<br>
&gt;<br>
&gt; And I need to solve the same SLE Au=b, where A is a Laplacian.<br>
&gt;<br>
&gt; How should I create and assemble the A matrix according to my DA vector to<br>
&gt; use the same functionality?<br>
<br>
</div></div>Create u,b with DAGetGlobalVector() and A with DAGetMatrix() and they<br>
will match the DA. For eg: check: src/snes/examples/tutorials/ex5.c<br>
[or some of the examples in src/dm/da/examples]<br>
<br>
Satish<br>
<div class="im"><br>
&gt;<br>
&gt; Thank you!<br>
&gt;<br>
&gt; Alexey Ryazanov<br>
&gt; ______________________________________<br>
&gt; Nuclear Safety Institute of Russian Academy of Sciences<br>
</div>&gt; &lt;<a href="http://www.ibrae.ac.ru/" target="_blank">http://www.ibrae.ac.ru/</a>&gt;<br>
&gt; </blockquote></div><br></div>