<div dir="ltr"><div class="gmail_quote">On Thu, Apr 7, 2011 at 09:42, Thomas Witkowski <span dir="ltr"><<a href="mailto:thomas.witkowski@tu-dresden.de">thomas.witkowski@tu-dresden.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":4kw">The problem is that the data field (vec->data) of the nested vectors is empty. I created the nested vectors exactly in the same way as in /src/ksp/ksp/examples/tests/ex22.c. When I would make a VecGetArray(myvec, &vecptr) on a nested vector myvec, than vecptr is not a vailde pointer. </div>
</blockquote><div><br></div><div>Indeed, this is why this test uses no preconditioner. I do not recommend VecNest for general use, there are many places where it will be problematic. So go ahead and use the MatNest, but just create a normal Vec.</div>
<div><br></div><div>Eventually, I want to remove the explicit scatters in the PCFieldSplit implementation in favor of VecGetSubVector(). When that is done, VecNest will provide a way to eliminate some copies (at the expense of memory layout so it will still not always be faster), but its use should still remain hidden from user code. Don't bother with it now.</div>
</div></div>