
<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;">Dear Jed,</span></font><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 17px;"><br></span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 17px;">I wanted to multiply those large matrices. I thought it was fast because I can get the results by just one multiplication. I think I was wrong and so I will try to get the result by the multiplication of a matrix and each vector. Thank you very much.</span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 17px;"><br></span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 17px;">Joon</span></font></div><br><br><br>-------- Original message --------<br>Subject: Re: [petsc-users] Errors from large matrices <br>From: Jed Brown <jedbrown@mcs.anl.gov> <br>To: Joon hee Choi <choi240@purdue.edu>,Barry Smith <bsmith@mcs.anl.gov> <br>CC: petsc-users@mcs.anl.gov <br><br><br><div style="word-break:break-all;">Joon hee Choi <choi240@purdue.edu> writes:<br><br>> Thank you very much. I was confused between 64 bit pointers and 64 bit<br>> indices. The module I am using supports 64 bit indices. I succeeded in<br>> setting up the matrix with 4273949 x 108965941330383 size (non-zeros:<br>> 143599552), but failed to set up the maij matrix with 108965941330383<br>> x 42739470 size and 4273947 block size using MatCreateMAIJ(Mat A,<br>> PetscInt dof, Mat *maij). The number of non-zeros are 25495389 x 10 x<br>> 4273947 = 1.1 x 10^15. Do I have the method to set up the matrix<br>> successfully?<br><br>What do you want to do with a matrix with this enormous number of columns?<br><br>PETSc matrix data structures are not intended to be used in this way.<br>You should reduce the dimension of the column space so that the<br>corresponding dense vector can be stored. </div> </body>