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->work[0];</div><div> x = ksp->work[1];</div><div> b = ksp->work[2];</div>
<div><div> Q = ksp->work[3];</div><div> M = ksp->work[4];</div><div> max_it = ksp->work[5];</div><div> tol = ksp->work[6];</div><div> *sol = ksp->work[7];</div><div> *err = ksp->work[8];</div>
<div> *iter = ksp->work[9];</div><div> *flag = ksp->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->max_it;</div><div> X = ksp->vec_sol;</div><div> B = ksp->vec_rhs;</div>
</div><div><br></div><div>Where are these variables (ksp->max_it, ksp->vec_sol and ksp->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"><<a href="mailto:vedaprakashsubramanian@gmail.com">vedaprakashsubramanian@gmail.com</a>></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"><<a href="mailto:petsc-users-request@mcs.anl.gov" target="_blank">petsc-users-request@mcs.anl.gov</a>></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 'help' 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 "Re: Contents of petsc-users digest..."<br>
<br>
<br>
Today'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 <<a href="mailto:vedaprakashsubramanian@gmail.com" target="_blank">vedaprakashsubramanian@gmail.com</a>><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>
<AANLkTinjiGMTMVB1yVx6JEvPhjGi_=<a href="mailto:T%2BQY6HK7wnHgOa@mail.gmail.com" target="_blank">T+QY6HK7wnHgOa@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<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: <<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>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Sun, 19 Sep 2010 21:49:47 -0500<br>
From: Barry Smith <<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>><br>
Subject: Re: [petsc-users] How to pass the parameters for KSP<br>
To: PETSc users list <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
Message-ID: <<a href="mailto:27B820EF-3A6A-424B-8A43-2A43E82089B3@mcs.anl.gov" target="_blank">27B820EF-3A6A-424B-8A43-2A43E82089B3@mcs.anl.gov</a>><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>
> 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>
><br>
> Thanks,<br>
> 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>