> On Feb 28, 2016, at 11:29 AM, Barry Smith <[email protected]> wrote: > > >> On Feb 28, 2016, at 11:18 AM, Guo, Hanqi <[email protected]> wrote: >> >> Yes. It works without MPIIO. > > We've never added support for saving sparse matrices with MPIIO; likely it > would not be much faster and it would be a lot of work to code. If you are > just saving matrices for testing it doesn't matter anyways; I'd never save > matrices in production runs.
Hi Barry, The matrices are the final products of my application. I’ll use serial IO instead. Thanks! -Hanqi > > Barry > >> >> -Hanqi >> >> --------------- >> Hanqi Guo >> Postdoctoral Appointee >> Mathematics and Computer Science Division >> Argonne National Laboratory >> 9700 S. Cass Ave., Bldg 240 >> Argonne, IL 60439 >> e-mail: [email protected] >> web: http://www.mcs.anl.gov/~hguo >> phone: 630-252-7225 >> >>> On Feb 28, 2016, at 11:05 AM, Barry Smith <[email protected]> wrote: >>> >>> >>> Does it work without the PetscViewerBinarySetUseMPIIO(viewer, PETSC_TRUE); >>> ? >>> >>> Barry >>> >>>> On Feb 27, 2016, at 9:57 PM, Guo, Hanqi <[email protected]> wrote: >>>> >>>> Hi, >>>> >>>> I am trying to write a sparse matrix to file with MPIIO but failed. I got >>>> the same error on vesta (gpfs), my Macbook (local disk), and a linux >>>> machine (nfs). The followings are my code snippet and error output. >>>> Could anyone help me check the code? Thanks! >>>> >>>> ... >>>> MatCreate(PETSC_COMM_WORLD, &mat); >>>> MatSetType(mat, MATMPIAIJ); >>>> MatSetSizes(mat, PETSC_DECIDE, PETSC_DECIDE, n, n); >>>> MatSetUp(mat); >>>> ... >>>> // write matrix >>>> ... >>>> MatAssemblyBegin(mat, MAT_FINAL_ASSEMBLY); >>>> MatAssemblyEnd(mat, MAT_FINAL_ASSEMBLY); >>>> >>>> PetscViewer viewer; >>>> PetscViewerCreate(PETSC_COMM_WORLD, &viewer); >>>> PetscViewerSetType(viewer, PETSCVIEWERBINARY); >>>> PetscViewerFileSetMode(viewer, FILE_MODE_WRITE); >>>> PetscViewerFileSetName(viewer, "sfm.mat"); >>>> PetscViewerBinarySetUseMPIIO(viewer, PETSC_TRUE); >>>> MatView(mat, viewer); >>>> PetscViewerDestroy(&viewer); >>>> ... >>>> >>>> >>>> [0]PETSC ERROR: --------------------- Error Message >>>> -------------------------------------------------------------- >>>> [0]PETSC ERROR: Write to file failed >>>> [0]PETSC ERROR: Error writing to file total size 4 err -1 wsize 16 >>>> [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html >>>> for trouble shooting. >>>> [0]PETSC ERROR: Petsc Release Version 3.6.3, Dec, 03, 2015 >>>> [0]PETSC ERROR: Unknown Name on a arch-linux2-c-debug named compute001 by >>>> hguo Sat Feb 27 21:40:04 2016 >>>> [0]PETSC ERROR: Configure options —prefix=(hidden) --with-fc=0 >>>> --download-f2cblaslapack --with-mpi-dir=(hidden) --with-precision=single >>>> [0]PETSC ERROR: #1 PetscBinaryWrite() line 405 in >>>> (hidden)/petsc-3.6.3/src/sys/fileio/sysio.c >>>> [0]PETSC ERROR: #2 MatView_MPIAIJ_Binary() line 1219 in >>>> (hidden)/petsc-3.6.3/src/mat/impls/aij/mpi/mpiaij.c >>>> [0]PETSC ERROR: #3 MatView_MPIAIJ_ASCIIorDraworSocket() line 1385 in >>>> (hidden)/petsc-3.6.3/src/mat/impls/aij/mpi/mpiaij.c >>>> [0]PETSC ERROR: #4 MatView_MPIAIJ() line 1470 in >>>> (hidden)/petsc-3.6.3/src/mat/impls/aij/mpi/mpiaij.c >>>> [0]PETSC ERROR: #5 MatView() line 886 in >>>> (hidden)/petsc-3.6.3/src/mat/interface/matrix.c >>>> >>>> >>>> --------------- >>>> Hanqi Guo >>>> Postdoctoral Appointee >>>> Mathematics and Computer Science Division >>>> Argonne National Laboratory >>>> 9700 S. Cass Ave., Bldg 240 >>>> Argonne, IL 60439 >>>> e-mail: [email protected] >>>> web: http://www.mcs.anl.gov/~hguo >>>> phone: 630-252-7225 >>>> >>> >> >
