> On Feb 28, 2016, at 11:39 AM, Guo, Hanqi <[email protected]> wrote: > > >> 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! \
What do you use the matrices for? Just curious > > -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 >>>>> >>>> >>> >> >
