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