Hello everyone; I am trying to compose a covariance matrix which is of size 72576 x 72576, by multiplication of a vector and its transpose, i.e.
MatMatMult(Ypcq,YpcqT,MAT_INITIAL_MATRIX,PETSC_DEFAULT,&InnProd) where, Ypcq is a vector with 72576 elements and YpcqT is its transpose. Unfortunately Petsc returns out of memory message while running the code: [0]PETSC ERROR: --------------------- Error Message ------------------------------------ [0]PETSC ERROR: Out of memory. This could be due to allocating [0]PETSC ERROR: too large an object or bleeding by not properly [0]PETSC ERROR: destroying unneeded objects. [0]PETSC ERROR: Memory allocated 0 Memory used by process 15732801536 [0]PETSC ERROR: Try running with -malloc_dump or -malloc_log for info. [0]PETSC ERROR: Memory requested 18446744068420534272! [0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: Petsc Release Version 3.2.0, Patch 3, Fri Sep 30 10:28:33 CDT 2011 [0]PETSC ERROR: See docs/changes/index.html for recent updates. [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. [0]PETSC ERROR: See docs/index.html for manual pages. [0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: ./dec5pbst23 on a linux-imp named k07n14.ccr.buffalo.edu by rm93 Mon Jan 2 18:10:28 2012 [0]PETSC ERROR: Libraries linked from /util/petsc/petsc-3.2-p3/linux-impi-mkl/lib [0]PETSC ERROR: Configure run at Fri Oct 21 08:36:23 2011 [0]PETSC ERROR: Configure options --CC=/util/intel/impi/ 4.0.3.008/intel64/bin/mpiicc --FC=/util/intel/impi/ 4.0.3.008/intel64/bin/mpiifort --CXX=/util/intel/impi/ 4.0.3.008/intel64/bin/mpiicpc--with-blas-lapack-dir=/util/intel/composer_xe_2011_sp1/mkl/lib/intel64 --download-hypre=1 --with-debugging=0 -PETSC_ARCH=linux-impi-mkl --with-shared-libraries=1 [0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: PetscMallocAlign() line 49 in src/sys/memory/mal.c [0]PETSC ERROR: PetscFreeSpaceGet() line 13 in src/mat/utils/freespace.c [0]PETSC ERROR: MatMatMultSymbolic_SeqAIJ_SeqAIJ() line 76 in src/mat/impls/aij/seq/matmatmult.c [0]PETSC ERROR: MatMatMult_SeqAIJ_SeqAIJ() line 21 in src/mat/impls/aij/seq/matmatmult.c [0]PETSC ERROR: MatMatMult() line 8246 in src/mat/interface/matrix.c [0]PETSC ERROR: main() line 331 in "unknowndirectory/"dec5pbst23.c application called MPI_Abort(MPI_COMM_WORLD, 55) - process 0 Is there any way to get rid of this error? Thanks in advance, Reza -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120102/c4afe663/attachment.htm>
