<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">Mahir:</div><div class="gmail_quote">Direct solvers consume large amount of memory. Suggest to try followings:</div><div class="gmail_quote"><br></div><div class="gmail_quote">1. A sparse iterative solver if <span style="font-family:Calibri,sans-serif;font-size:14.6666669845581px"> [-omega^2M + K] is not too ill-conditioned. You may test it using the small matrix.</span></div><div class="gmail_quote"><font face="Calibri, sans-serif"><span style="font-size:14.6666669845581px"><br></span></font></div><div class="gmail_quote"><font face="Calibri, sans-serif"><span style="font-size:14.6666669845581px">2. Incrementally increase your matrix sizes. Try different matrix orderings.</span></font></div><div class="gmail_quote"><font face="Calibri, sans-serif"><span style="font-size:14.6666669845581px">Do you get memory crash in the 1st symbolic factorization? </span></font></div><div class="gmail_quote"><font face="Calibri, sans-serif"><span style="font-size:14.6666669845581px">In your case, matrix data structure stays same when omega changes, so you only need to do one matrix symbolic factorization and reuse it.</span></font></div><div class="gmail_quote"><font face="Calibri, sans-serif"><span style="font-size:14.6666669845581px"><br></span></font></div><div class="gmail_quote"><font face="Calibri, sans-serif"><span style="font-size:14.6666669845581px">3. Use a machine that gives larger memory.</span></font></div><div class="gmail_quote"><font face="Calibri, sans-serif"><span style="font-size:14.6666669845581px"><br></span></font></div><div class="gmail_quote"><font face="Calibri, sans-serif"><span style="font-size:14.6666669845581px">Hong</span></font></div><div class="gmail_quote"><font face="Calibri, sans-serif"><span style="font-size:14.6666669845581px"><br></span></font><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">





<div lang="SV" link="blue" vlink="purple">
<div>
<p class="MsoNormal">Dear Petsc-Users,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><span lang="EN-US">I am trying to use PETSc to solve a set of linear equations arising from Naviers equation (elastodynamics) in the frequency domain.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">The frequency dependency of the problem requires that the system<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">                             [-omega^2M + K]u = F<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">where M and K are constant, square, positive definite matrices (mass and stiffness respectively) is solved for each frequency omega of interest.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">K is a complex matrix, including material damping.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">I have written a PETSc program which solves this problem for a small (1000 degrees of freedom) test problem on one or several processors, but it keeps crashing when I try it on my full scale (in the order of 10^6 degrees
 of freedom) problem.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">The program crashes at KSPSetUp() and from what I can see in the error messages, it appears as if it consumes too much memory.
<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">I would guess that similar problems have occurred in this mail-list, so I am hoping that someone can push  me in the right direction…<span><font color="#888888"><u></u><u></u></font></span></span></p><span><font color="#888888">
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">Mahir<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
</font></span></div>
</div>

</blockquote></div><br></div></div>