<html><head><meta http-equiv="Content-Type" content="text/html; charset=UTF-8"></head><body><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 17px;">Thank you for your fast reply. However, I was already using 64 bit pointer.</span></font><br><br><br>-------- Original message --------<br>Subject: Re: [petsc-users] Errors from large matrices <br>From: Barry Smith <bsmith@mcs.anl.gov> <br>To: Joon hee Choi <choi240@purdue.edu> <br>CC: petsc-users@mcs.anl.gov <br><br><br><div style="word-break:break-all;"><br>On May 26, 2013, at 6:19 AM, Joon hee Choi <choi240@purdue.edu> wrote:<br><br>> Hello all,<br>> <br>> I need to multiply a large seqaij matrix(X1) and a maij(or baij) matrix(CC). I set up X1 (size:4273949x108965941330383, nonzeros:143599552) and C (size:25495389x10, nonzeros:254953890) and created a maij matrix CC from C. However, I got errors such as out of memory and Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range. Is this memory problem, and do I have to change seqaij into mpiaij and use multi processors? Or do I have another methods fixing it? If you know the method, then please let me know it. Thank you.<br><br>http://www.mcs.anl.gov/petsc/documentation/faq.html#with-64-bit-indices<br><br><br>> <br>> Joon<br>> <br>> <br>> Code:<br>> ...<br>> ierr = MatCreate(PETSC_COMM_SELF, &X1); CHKERRQ(ierr);<br>> ierr = MatSetSizes(X1, PETSC_DECIDE, PETSC_DECIDE, I, J*K); CHKERRQ(ierr);<br>> ierr = MatSetBlockSizes(X1, I, J); CHKERRQ(ierr);<br>> ierr = MatSetType(X1, MATSEQAIJ); CHKERRQ(ierr);<br>> ierr = MatSeqAIJSetPreallocation(X1, 0, nnz); CHKERRQ(ierr);<br>> <br>> for (int x=0; x<tups.size(); x++) {<br>> i = std::tr1::get<0>(tups[x]);<br>> j = std::tr1::get<2>(tups[x]) + std::tr1::get<1>(tups[x])*J;<br>> val = std::tr1::get<3>(tups[x]);<br>> ierr = MatSetValues(X1, 1, &i, 1, &j, &val, INSERT_VALUES); CHKERRQ(ierr);<br>> }<br>> ierr = MatAssemblyBegin(X1, MAT_FINAL_ASSEMBLY);<br>> ierr = MatAssemblyEnd(X1, MAT_FINAL_ASSEMBLY);<br>> ierr = PetscGetTime(&v1); CHKERRQ(ierr);<br>> ierr = PetscPrintf(PETSC_COMM_WORLD, "Setup Time: %2.1e \n", v1-v); CHKERRQ(ierr);<br>> <br>> // Create a matrix C (K x R) with all values 1<br>> ierr = MatCreateSeqAIJ(PETSC_COMM_SELF, K, R, R, NULL, &C); CHKERRQ(ierr);<br>> for (k=0; k<K; k++) {<br>> for (r=0; r<R; r++) {<br>> ierr = MatSetValues(C, 1, &k, 1, &r, &one, INSERT_VALUES); CHKERRQ(ierr);<br>> }<br>> }<br>> ierr = MatAssemblyBegin(C, MAT_FINAL_ASSEMBLY); CHKERRQ(ierr);<br>> ierr = MatAssemblyEnd(C, MAT_FINAL_ASSEMBLY); CHKERRQ(ierr);<br>> <br>> ierr = M...</div> </body>