<div dir="ltr">Thank you very much, but I need to access this value in my program and not just print it, is there a way to do this? Also, I need to use Cholesky factorization, is there a way to avoid MUMPS so I could access the diagonal?<br>
</div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Feb 25, 2013 at 7:13 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@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"><br>
  Jelena,<br>
<br>
     MUMPS has its own (complicated) data structure for storing the factored matrix, you cannot easily look inside it.<br>
<br>
      However, fortunately MUMPS has a way of accessing the determinate of the matrix directly. Look at src/mat/impls/aij/mpi/mumps.c and search for the word determinant, you will find the command line argument to have it computed and where it will be printed out if you call MatView() on the factored matrix. You mush pass the PetscViewerFormat PETSC_VIEWER_ASCII_INFO to PetscViewerSetFormat() before viewing the matrix.<br>

<br>
   Good luck,<br>
<br>
   Barry<br>
<div class="HOEnZb"><div class="h5"><br>
On Feb 25, 2013, at 6:02 PM, Jelena Slivka <<a href="mailto:slivkaje@gmail.com">slivkaje@gmail.com</a>> wrote:<br>
<br>
> Hello,<br>
> I am trying to find the determinant of the matrix. I have found this thread:<br>
><br>
> <a href="http://lists.mcs.anl.gov/pipermail/petsc-users/2010-July/006716.html" target="_blank">http://lists.mcs.anl.gov/pipermail/petsc-users/2010-July/006716.html</a><br>
><br>
> But, could you please explain how to access the data structure directly?<br>
> So far I have tried the following:<br>
><br>
> #include <petsc-private/matimpl.h><br>
> #include <../src/mat/impls/aij/seq/aij.h><br>
> ...<br>
> // (Solving the system A*x = rhs<br>
> MatSetOption(A, MAT_SPD, PETSC_TRUE);<br>
> Mat F;<br>
> MatGetFactor(A, MATSOLVERMUMPS, MAT_FACTOR_CHOLESKY, &F);<br>
> MatFactorInfo factinfo;<br>
> MatFactorInfoInitialize(&factinfo);<br>
> IS perm, iscol;<br>
> MatGetOrdering(A, MATORDERINGNATURAL, &perm, &iscol);<br>
> MatCholeskyFactorSymbolic(F, A, perm, &factinfo);<br>
> MatCholeskyFactorNumeric(F, A, &factinfo);<br>
> ISDestroy(&iscol);<br>
> ISDestroy(&perm);<br>
><br>
> MatSolve(F, rhs, x);<br>
><br>
> // Now I would like to access the diagonal part of the matrix F, can I do something like:<br>
><br>
> Mat_SeqAIJ *aa = (Mat_SeqAIJ*)F->data;<br>
><br>
> Do you maybe have an example I could look at?<br>
><br>
> Grateful in advance,<br>
> Jelena<br>
<br>
</div></div></blockquote></div><br></div>