Can anyone tell me whether the concept I understood is correct or not.<div><br></div><div>If I am having a MATLAB function func(Matrix A, Vector x, Vector b, Matrix Q, Matrix M, int max_it, int tol, Vec* sol, int* err, int* iter, int* flag)   to be converted into a KSP solver, then I have to do these following steps.</div>
<div><br></div><div>1. In KSPSetUp_func(KSP ksp), initialize  KSPDefaultGetWork(ksp,11);  // This tells the KSP that it is going to take 11 arguments</div><div>2. In KSPSolve_func(KSP ksp), declare </div><div>    Mat A, Q, M;</div>
<div>    Vec x, b, *sol;</div><div>    PetscInt max_it, *flag, *iter;</div><div>    PetscScalar tol, *err;</div><div><br></div><div>    A  = ksp-&gt;work[0];</div><div>    x   = ksp-&gt;work[1];</div><div>    b   = ksp-&gt;work[2];</div>
<div><div>    Q  = ksp-&gt;work[3];</div><div>    M   = ksp-&gt;work[4];</div><div>    max_it   = ksp-&gt;work[5];</div><div>    tol  = ksp-&gt;work[6];</div><div>    *sol   = ksp-&gt;work[7];</div><div>    *err   = ksp-&gt;work[8];</div>
<div>    *iter   = ksp-&gt;work[9];</div><div>    *flag   = ksp-&gt;work[10];</div><div><br></div><div>------------------------------</div><div>Is the above said is correct. Is that the way to get the arguments of the KSPSolve.</div>
<div><br></div><div>Moreover, I have a doubt. What are these variables that are been used in cg.c</div><div><div>  </div><div>  stored_max_it = ksp-&gt;max_it;</div><div>  X             = ksp-&gt;vec_sol;</div><div>  B             = ksp-&gt;vec_rhs;</div>
</div><div><br></div><div>Where are these variables (ksp-&gt;max_it, ksp-&gt;vec_sol and ksp-&gt;vec_rhs) getting initialized.</div><div><br></div><div>-Vedaprakash</div><div><br></div><div class="gmail_quote">On Wed, Oct 6, 2010 at 7:56 PM, vedaprakash subramanian <span dir="ltr">&lt;<a href="mailto:vedaprakashsubramanian@gmail.com">vedaprakashsubramanian@gmail.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">I need my KSP to take in Tolerence, flag, err, 3 matrices, right hand side vector. How can I make it to take all the parameters.<br>
<br><div class="gmail_quote">On Mon, Sep 20, 2010 at 11:00 AM,  <span dir="ltr">&lt;<a href="mailto:petsc-users-request@mcs.anl.gov" target="_blank">petsc-users-request@mcs.anl.gov</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send petsc-users mailing list submissions to<br>
        <a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="https://lists.mcs.anl.gov/mailman/listinfo/petsc-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/petsc-users</a><br>
or, via email, send a message with subject or body &#39;help&#39; to<br>
        <a href="mailto:petsc-users-request@mcs.anl.gov" target="_blank">petsc-users-request@mcs.anl.gov</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:petsc-users-owner@mcs.anl.gov" target="_blank">petsc-users-owner@mcs.anl.gov</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than &quot;Re: Contents of petsc-users digest...&quot;<br>
<br>
<br>
Today&#39;s Topics:<br>
<br>
   1.  How to pass the parameters for KSP (vedaprakash subramanian)<br>
   2. Re:  How to pass the parameters for KSP (Barry Smith)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Sun, 19 Sep 2010 19:05:48 -0600<br>
From: vedaprakash subramanian &lt;<a href="mailto:vedaprakashsubramanian@gmail.com" target="_blank">vedaprakashsubramanian@gmail.com</a>&gt;<br>
Subject: [petsc-users] How to pass the parameters for KSP<br>
To: <a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a><br>
Message-ID:<br>
        &lt;AANLkTinjiGMTMVB1yVx6JEvPhjGi_=<a href="mailto:T%2BQY6HK7wnHgOa@mail.gmail.com" target="_blank">T+QY6HK7wnHgOa@mail.gmail.com</a>&gt;<br>
Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
<br>
I am converting a MATLAB function into a KSP solver. I am doing it similar<br>
to BiCGStab. But I wanted to know how to pass the arguments of the function<br>
into KSP solver.<br>
<br>
Thanks,<br>
Vedaprakash<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: &lt;<a href="http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100919/9d902e74/attachment-0001.htm" target="_blank">http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100919/9d902e74/attachment-0001.htm</a>&gt;<br>


<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Sun, 19 Sep 2010 21:49:47 -0500<br>
From: Barry Smith &lt;<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>&gt;<br>
Subject: Re: [petsc-users] How to pass the parameters for KSP<br>
To: PETSc users list &lt;<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>&gt;<br>
Message-ID: &lt;<a href="mailto:27B820EF-3A6A-424B-8A43-2A43E82089B3@mcs.anl.gov" target="_blank">27B820EF-3A6A-424B-8A43-2A43E82089B3@mcs.anl.gov</a>&gt;<br>
Content-Type: text/plain; charset=us-ascii<br>
<br>
<br>
 What arguments? Do you mean the right hand side x and the matrix? Or do you mean parameters like the relative tolerance in convergence?<br>
<br>
  See src/ksp/ksp/examples/tutorials/ex1.c for a simple example.<br>
<br>
   Barry<br>
<br>
On Sep 19, 2010, at 8:05 PM, vedaprakash subramanian wrote:<br>
<br>
&gt; I am converting a MATLAB function into a KSP solver. I am doing it similar to BiCGStab. But I wanted to know how to pass the arguments of the function into KSP solver.<br>
&gt;<br>
&gt; Thanks,<br>
&gt; Vedaprakash<br>
<br>
<br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
petsc-users mailing list<br>
<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/petsc-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/petsc-users</a><br>
<br>
<br>
End of petsc-users Digest, Vol 21, Issue 32<br>
*******************************************<br>
</blockquote></div><br>
</blockquote></div><br></div>