<div dir="ltr"><div><div><div>Hi, I managed to plot it in matlab with the PetscBinaryRead with options: vc=PetscBinaryRead('outvec','complex',true);<br><br></div>I missed the option 'complex' previously. <br><br></div>It is working fine now. I will ask another separate query so that this can be closed. <br><br></div>Venkatesh<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, May 16, 2015 at 1:16 AM, venkatesh g <span dir="ltr"><<a href="mailto:venkateshgk.j@gmail.com" target="_blank">venkateshgk.j@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi, I am now running the solver with the mumps <div><br></div><div>/ex7 -f1 adt -f2 bdt -st_type sinvert -st_ksp_type preonly -st_pc_type lu -st_pc_factor_mat_solver_package mumps -evecs outvec<br></div><div><br></div><div>However, my eigenvalue is not right by 0.5% and the eigenvector is not correct when I plotted it. </div><div>It has to match the matlab result as it is a benchmark from a book.</div><div><br></div><div>Kindly advice. </div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Venkatesh</div><div><br></div><div><br></div></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Fri, May 15, 2015 at 8:24 PM, venkatesh g <span dir="ltr"><<a href="mailto:venkateshgk.j@gmail.com" target="_blank">venkateshgk.j@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Ah. I made mistake by looking at the old manual. I am now reading the latest version of the manual. <br><br></div>I will do this and get back. <br><span><font color="#888888"><br></font></span></div><span><font color="#888888">Venkatesh<br></font></span></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, May 15, 2015 at 6:14 PM, Jose E. Roman <span dir="ltr"><<a href="mailto:jroman@dsic.upv.es" target="_blank">jroman@dsic.upv.es</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">You have to choose MUMPS as the solver. ONCE AGAIN, read section 3.4.1 of the users manual.<br>
<span><font color="#888888">Jose<br>
</font></span><div><div><br>
<br>
El 15/05/2015, a las 13:48, venkatesh g escribió:<br>
<br>
> Ok, I followed this process and compiled ex7 successfully.<br>
><br>
> Now I give the command: './ex7 -f1 A11 -f2 B11 -eps_nev 4 -eps_target 0.5 -st_type sinvert'<br>
><br>
> The error is the zero pivot problem:<br>
><br>
>  Reading COMPLEX matrices from binary files...<br>
> [0]PETSC ERROR: --------------------- Error Message ------------------------------------<br>
> [0]PETSC ERROR: Detected zero pivot in LU factorization:<br>
> see <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html#ZeroPivot" target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html#ZeroPivot</a>!<br>
> [0]PETSC ERROR: Zero pivot row 420 value 0 tolerance 2.22045e-14!<br>
> [0]PETSC ERROR: ------------------------------------------------------------------------<br>
> [0]PETSC ERROR: Petsc Release Version 3.3.0, Patch 7, Sat May 11 22:15:24 CDT 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: ./ex7 on a linux-gnu named <a href="http://earth.ceas.iisc.ernet.in" target="_blank">earth.ceas.iisc.ernet.in</a> by venkatesh Fri May 15 17:07:02 2015<br>
> [0]PETSC ERROR: Libraries linked from /cluster/share/venkatesh/petsc-3.3-p7/linux-gnu-intel-cmplx-batch/lib<br>
> [0]PETSC ERROR: Configure run at Fri May 15 15:10:57 2015<br>
> [0]PETSC ERROR: Configure options --with-cc=gcc --with-fc=ifort --with-blas-lapack-dir=/opt/intel/composer_xe_2013.5.192/mkl/lib/intel64 --download-mpich=/cluster/share/apps/mpich2-1.4.1p1.tar.gz --download-mumps=/cluster/share/apps/MUMPS_4.10.0-p3.tar.gz --download-scalapack=/cluster/share/apps/SCALAPACK-1.7.tar.gz --download-blacs=/cluster/share/apps/blacs-dev.tar.gz --download-parmetis=/cluster/share/apps/parmetis-4.0.2-p3.tar.gz --download-metis=/cluster/share/apps/metis-5.0.2-p3.tar.gz --download-cmake=/cluster/share/apps/cmake-2.8.7.tar.gz --with-scalar-type=complex --with-fortran-kernels=generic --with-large-file-io=1<br>
> [0]PETSC ERROR: ------------------------------------------------------------------------<br>
> [0]PETSC ERROR: MatPivotCheck_none() line 583 in /cluster/share/venkatesh/petsc-3.3-p7/include/petsc-private/matimpl.h<br>
> [0]PETSC ERROR: MatPivotCheck() line 602 in /cluster/share/venkatesh/petsc-3.3-p7/include/petsc-private/matimpl.h<br>
> [0]PETSC ERROR: MatLUFactorNumeric_SeqAIJ() line 570 in /cluster/share/venkatesh/petsc-3.3-p7/src/mat/impls/aij/seq/aijfact.c<br>
> [0]PETSC ERROR: MatLUFactorNumeric() line 2803 in /cluster/share/venkatesh/petsc-3.3-p7/src/mat/interface/matrix.c<br>
> [0]PETSC ERROR: PCSetUp_LU() line 160 in /cluster/share/venkatesh/petsc-3.3-p7/src/ksp/pc/impls/factor/lu/lu.c<br>
> [0]PETSC ERROR: PCSetUp() line 832 in /cluster/share/venkatesh/petsc-3.3-p7/src/ksp/pc/interface/precon.c<br>
> [0]PETSC ERROR: KSPSetUp() line 278 in /cluster/share/venkatesh/petsc-3.3-p7/src/ksp/ksp/interface/itfunc.c<br>
> [0]PETSC ERROR: PCSetUp_Redundant() line 176 in /cluster/share/venkatesh/petsc-3.3-p7/src/ksp/pc/impls/redundant/redundant.c<br>
> [0]PETSC ERROR: PCSetUp() line 832 in /cluster/share/venkatesh/petsc-3.3-p7/src/ksp/pc/interface/precon.c<br>
> [0]PETSC ERROR: KSPSetUp() line 278 in /cluster/share/venkatesh/petsc-3.3-p7/src/ksp/ksp/interface/itfunc.c<br>
> [0]PETSC ERROR: STSetUp_Sinvert() line 153 in /cluster/share/venkatesh/petsc-3.3-p7/slepc-3.3-p4/src/st/impls/sinvert/sinvert.c<br>
> [0]PETSC ERROR: STSetUp() line 280 in /cluster/share/venkatesh/petsc-3.3-p7/slepc-3.3-p4/src/st/interface/stsolve.c<br>
> [0]PETSC ERROR: EPSSetUp() line 204 in /cluster/share/venkatesh/petsc-3.3-p7/slepc-3.3-p4/src/eps/interface/setup.c<br>
> [0]PETSC ERROR: EPSSolve() line 109 in /cluster/share/venkatesh/petsc-3.3-p7/slepc-3.3-p4/src/eps/interface/solve.c<br>
> [0]PETSC ERROR: main() line 148 in ex7.c<br>
> application called MPI_Abort(MPI_COMM_WORLD, 71) - process 0<br>
> [unset]: aborting job:<br>
> application called MPI_Abort(MPI_COMM_WORLD, 71) - process 0<br>
><br>
> Venkatesh<br>
><br>
> On Fri, May 15, 2015 at 4:22 PM, Jose E. Roman <<a href="mailto:jroman@dsic.upv.es" target="_blank">jroman@dsic.upv.es</a>> wrote:<br>
><br>
> El 15/05/2015, a las 12:49, venkatesh g escribió:<br>
><br>
> > Hi,<br>
> > 'make test' was successful on the $PETSC_DIR.<br>
> > I configured MUMPS, metis, parmetis, scalapack through '--download-mumps=/path/MUMPS_4.10.0-p3.tar.gz' and etc..<br>
> ><br>
> > I will anyway download the new version and start installing it now., Is it due to some flag issue ? I have attached my makefile.<br>
> ><br>
> > Venkatesh<br>
> ><br>
><br>
> You must link with all libraries, otherwise you will get undefined symbols. The easiest way is to include PETSc and SLEPc makefile definitions. See for instance section 1.4.2 of SLEPc's users manual.<br>
><br>
> Jose<br>
><br>
><br>
<br>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>