<div dir="ltr">So I was adding a new test in src/mat/examples/tests/ using the file ex125.c, this test required data from <a href="http://ftp.mcs.anl.gov/" style="font-family:Times;font-size:medium;background-color:rgb(213,234,255)">ftp.mcs.anl.gov</a>:/pub/petsc/matrices, I download the data, but I get the following error:<div>
<br></div><div><div>> [0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div><div>> [0]PETSC ERROR: Read from file failed</div><div>> [0]PETSC ERROR: Read past end of file</div>
<div>> [0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.</div><div>> [0]PETSC ERROR: Petsc Development GIT revision:   GIT Date: </div>
<div>> [0]PETSC ERROR: ./ex125 on a linux-complex named <a href="http://conte-fe00.rcac.purdue.edu">conte-fe00.rcac.purdue.edu</a> by jbermeol Thu Jun 26 21:36:20 2014</div><div>> [0]PETSC ERROR: Configure options CC=mpiicc CXX=mpiicpc FC=mpiifort COPTFLAGS=-O3 CXXOPTFLAGS=-O3 FOPTFLAGS=-O3 --with-scalar-type=complex --with-shared-libraries=0 --with-debugging=1 --with-pic=1 --with-clanguage=C++ --with-fortran=1 --with-fortran-kernels=0 --download-mumps --download-metis --download-parmetis --with-scalapack-lib="-L/apps/rhel6/intel/composer_xe_2013.3.163/mkl/lib/intel64 -lmkl_scalapack_lp64 -lmkl_blacs_intelmpi_lp64" --with-scalapack-include=/apps/rhel6/intel/composer_xe_2013.3.163/mkl/include --with-blas-lapack-dir=/apps/rhel6/intel/composer_xe_2013.3.163/mkl --with-mkl-pardiso-dir=/apps/rhel6/intel/composer_xe_2013.3.163/mkl</div>
<div>> [0]PETSC ERROR: #1 PetscBinaryRead() line 274 in /home/jbermeol/software/petsc_cplx/src/sys/fileio/sysio.c</div><div>> [0]PETSC ERROR: #2 MatLoad_SeqAIJ() line 4312 in /home/jbermeol/software/petsc_cplx/src/mat/impls/aij/seq/aij.c</div>
<div>> [0]PETSC ERROR: #3 MatLoad() line 997 in /home/jbermeol/software/petsc_cplx/src/mat/interface/matrix.c</div><div>> [0]PETSC ERROR: #4 main() line 36 in ex125.c</div><div>> [0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint@mcs.anl.gov----------</div>
<div>> application called MPI_Abort(MPI_COMM_WORLD, 66) - process 0</div></div><div><br></div><div>Not sure if the file is corrupted, could you please help me with that. Thanks </div></div><div class="gmail_extra"><br>
<br><div class="gmail_quote">On Thu, Jun 26, 2014 at 3:38 PM, Satish Balay <span dir="ltr"><<a href="mailto:balay@mcs.anl.gov" target="_blank">balay@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">
<div class="">On Thu, 26 Jun 2014, José David Bermeo López wrote:<br>
<br>
> Hi, I would like to add a test to the repository for mkl_pardiso, looking<br>
> for mumps test, I sow test ex125.c that uses MUMPS or other solvers<br>
> depending on the option, so i have some questions:<br>
><br>
> 1. Would be a good idea to add mkl_pardiso in that test(ex125.c)<br>
<br>
</div>I've added one test. But more tests would be good.<br>
[check runex2_mkl_pardiso in src/ksp/ksp/examples/tutorials/makefile]<br>
<div class=""><br>
> 2. For that test you load matrix A, is it a distributed matrix or a<br>
> sequential matrix.<br>
<br>
</div>Hm - this example appears to also test Superlu [which is sequential] -<br>
so adding mkl_pardiso similarly might work..<br>
<br>
[Hong might have better suggestions]<br>
<div class=""><br>
<br>
> 3. You had html pages, for documentation of each solver, is this<br>
> automatically generated base on comments, or do I had to added my self.<br>
<br>
</div>Ah - I didn't check that. You can look at the end of<br>
mat/impls/aij/mpi/superlu_dist/superlu_dist.c for its<br>
documentation. If you add something similar to mkl_pardiso.c - the<br>
docs for it will get autogenerated.<br>
<span class="HOEnZb"><font color="#888888"><br>
Satish</font></span></blockquote></div><br></div>