<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Aug 5, 2015 at 4:15 AM, Nicolas Pozin <span dir="ltr"><<a href="mailto:nicolas.pozin@inria.fr" target="_blank">nicolas.pozin@inria.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="font-family:times new roman,new york,times,serif;font-size:12pt;color:#000000"><div>Hello,<br></div><div><br></div><div>I'm trying to solve a system with a matrix free operator and through conjugate gradient method. <br></div><div>To make ideas clear, I set up the following simple example (I am using petsc-3.6) and I get this error message :</div></div></div></blockquote><div><br></div><div>Yes, you are passing a C++ function userMult, so the compiler sticks "this" in as the first argument. We do not</div><div>recommend this kind of wrapping.</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="font-family:times new roman,new york,times,serif;font-size:12pt;color:#000000"><div><div>"</div><div>[0]PETSC ERROR: --------------------- Error Message ------------------------------------<br>[0]PETSC ERROR: Invalid argument!<br>[0]PETSC ERROR: Wrong type of object: Parameter # 1!<br>[0]PETSC ERROR: ------------------------------------------------------------------------<br>[0]PETSC ERROR: Petsc Release Version 3.4.3, Oct, 15, 2013 <br>[0]PETSC ERROR: See docs/changes/index.html for recent updates.<br>[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>[0]PETSC ERROR: See docs/index.html for manual pages.<br>[0]PETSC ERROR: ------------------------------------------------------------------------<br>[0]PETSC ERROR: ./test on a ubuntu_release named pl-59080 by npozin Wed Aug 5 10:55:26 2015<br>[0]PETSC ERROR: Libraries linked from /home/npozin/Felisce_libraries/petsc_3.4.3/ubuntu_release/lib<br>[0]PETSC ERROR: Configure run at Wed Jul 22 16:18:36 2015<br>[0]PETSC ERROR: Configure options PETSC_ARCH=ubuntu_release --with-cxx=g++ --with-fc=gfortran --with-cc=gcc --with-x=0 --download-openmpi --download-f-blas-lapack --download-superlu --download-superlu_dist --with-superlu_dist=1 --download-metis --download-mumps --download-parmetis --with-superlu_dist=1 --download-boost --with-boost=1 --download-scalapack with-external-packages-dir=/home/npozin/Felisce_libraries/petsc_3.4.3/packages<br>[0]PETSC ERROR: ------------------------------------------------------------------------<br>[0]PETSC ERROR: MatShellGetContext() line 202 in /home/npozin/Felisce_libraries/petsc_3.4.3/src/mat/impls/shell/shell.c<br>End userMult<br>[0]PETSC ERROR: MatMult() line 2179 in /home/npozin/Felisce_libraries/petsc_3.4.3/src/mat/interface/matrix.c<br>[0]PETSC ERROR: KSP_MatMult() line 204 in /home/npozin/Felisce_libraries/petsc_3.4.3/include/petsc-private/kspimpl.h<br>[0]PETSC ERROR: KSPSolve_CG() line 219 in /home/npozin/Felisce_libraries/petsc_3.4.3/src/ksp/ksp/impls/cg/cg.c<br>[0]PETSC ERROR: KSPSolve() line 441 in /home/npozin/Felisce_libraries/petsc_3.4.3/src/ksp/ksp/interface/itfunc.c<br>"</div><div><br></div><div>I don't understand where the problem comes from with the matrix argument of MatShellGetContext.</div><div>Any idea on what I do wrong?</div><div><br></div><div>Thanks a lot,</div>Nicolas</div><div><br></div><div><br></div><div><br></div><div>#include <iostream><br>#include <petscksp.h><br><br>using namespace std;<br><br><br>typedef struct {<br> int val;<br>} MyCtx;<br><br><br>class ShellClass {<br> Mat matShell;<br> KSP ksp;<br> PC pc;<br> Vec x;<br> Vec b;<br> <br>public:<br> void userMult(Mat Amat, Vec x, Vec y) {<br> cout << "Inside userMult" << endl;<br> <br> MyCtx *ctx;<br> MatShellGetContext(Amat, (void *) ctx);<br> <br> cout << "End userMult" << endl;<br> }<br><br> void solveShell() {<br> // context<br> MyCtx *ctx = new MyCtx;<br> ctx->val = 42;<br><br> // pc<br> PCCreate(PETSC_COMM_WORLD, &pc);<br> PCSetType(pc, PCNONE);<br><br> // ksp<br> KSPCreate(PETSC_COMM_WORLD, &ksp);<br> KSPSetType(ksp, KSPCG);<br> KSPSetPC(ksp, pc);<br> KSPSetFromOptions(ksp);<br><br> // matshell<br> int m = 10;<br> int n = 10;<br> MatCreateShell(PETSC_COMM_WORLD, m, n, PETSC_DETERMINE, PETSC_DETERMINE, ctx, &matShell);<br> MatShellSetOperation(matShell, MATOP_MULT, (void(*)(void))&ShellClass::userMult);<br><br><br> // create vectors<br> MatCreateVecs(matShell, &x, 0);<br> VecDuplicate(x, &b);<br> VecSet(b, 1.);<br><br> // set operators<br> KSPSetOperators(ksp, matShell, matShell);<br><br> // solve (call to userMult)<br> KSPSolve(ksp, b, x);<br> }<br>};<br><br><br><br>int main(int argc, char** argv) {<br> PetscInitialize(&argc, &argv, NULL, NULL);<br> <br> ShellClass foo;<br> foo.solveShell();<br> <br> PetscFinalize();<br> return 0;<br>}<br><br></div></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">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</div>
</div></div>