[petsc-users] multiplication of large matrices using MatMatMult

[Reza Madankan]

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>