<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><div id="yiv1195565870"><div class="yiv1195565870gmail_quote">On Thu, May 19, 2011 at 11:09, Bartłomiej W <span dir="ltr"><<a rel="nofollow" ymailto="mailto:bartlomiej.wach@yahoo.pl" target="_blank" href="/mc/compose?to=bartlomiej.wach@yahoo.pl">bartlomiej.wach@yahoo.pl</a>></span> wrote:<br><blockquote class="yiv1195565870gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<span style="font-style: italic;">[0]Total space allocated 2514,564,304 bytes<br>(...)<br>[0]PETSC ERROR: Memory allocated<span style="font-weight: bold;"> -<a rel="nofollow">2147483648</a> </span>Memory used by process 1848942592</span><span style="font-style: italic;"></span><br style="font-style: italic;">
<span style="font-style: italic;">[0]PETSC ERROR: Memory requested 1578,397,300!</span></blockquote><div><br></div><div><div>Please always send the whole error trace. I have updated petsc-dev so the sizes are not incorrectly truncated like you see here. Are you sure you have built a 64-bit executable (run "file your-executable")?</div>
</div><div> </div></div></div></blockquote>This might be the problem, it occurs it is 32 bit, tried to compile with -m64 though but got <br><br style="font-style: italic;"><span style="font-style: italic;">/usr/bin/ld: skipping incompatible .../petsc-3.1-p8/linux-gnu-c-debug/lib/libpetsc.a when searching for -lpetsc</span><br style="font-style: italic;"><br><div id="yiv1195565870"><div class="yiv1195565870gmail_quote"><br style="font-style: italic;"></div></div><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><div id="yiv1195565870"><div class="yiv1195565870gmail_quote"><blockquote class="yiv1195565870gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><br>Caused by the code<br><br><span style="font-style: italic;"> ierr = MatCreate(PETSC_COMM_WORLD,&Dinv);CHKERRQ(ierr);</span><br style="font-style: italic;">
<span style="font-style: italic;"> ierr = MatSetSizes(Dinv,PETSC_DECIDE,PETSC_DECIDE,n,n);CHKERRQ(ierr);</span><br style="font-style: italic;"><span style="font-style: italic;"> ierr = MatSetFromOptions(Dinv);CHKERRQ(ierr);</span><br style="font-style: italic;">
<span style="font-style: italic;"> ierr = MatAssemblyBegin(Dinv,MAT_FINAL_ASSEMBLY);CHKERRQ(ierr);</span></blockquote><div><br></div><div>Is this in parallel? How many processors? Does the name "Dinv" imply that this matrix is supposed to be the inverse of something? Why are you assembling without putting anything in the matrix?</div>
<div><br></div></div></div></blockquote><div id="yiv1195565870"><div class="yiv1195565870gmail_quote"><br></div></div><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><div id="yiv1195565870"><div class="yiv1195565870gmail_quote"><blockquote class="yiv1195565870gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><br>where n = 39459925</blockquote><div><br></div><div>How many nonzeros per row?<br><br></div></div></div></blockquote>I try to run in on single processor right now. Yes, Dinv is supposed to be an inverse of a diagonal matrix (1 nnz/row)<br>
<br>
Made it sparse and the app went on but stopped a bit further on next allocation, problem is the 32 bit build while it should be 64 but I cant get it right.<div id="yiv1195565870"><div class="yiv1195565870gmail_quote"><div><br></div></div></div><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><div id="yiv1195565870">
</div></blockquote></td></tr></table>