<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Exchange Server">
<!-- converted from rtf -->
<style><!-- .EmailQuote { margin-left: 1pt; padding-left: 4pt; border-left: #800000 2px solid; } --></style>
</head>
<body>
<font face="Calibri" size="2"><span style="font-size:11pt;">
<div>I&#8217;m running into time and memory problems creating an array to store multiple KSP objects.</div>
<div>&nbsp;</div>
<div>We are using a Fourier method for the temporal portion of our operator; currently, this is anywhere between 1024-2048 time samples equating to 1024-2048 frequencies to advance our solution. Each frequency has its own matrix, meaning there are 1024-2048
linear equations of the form Ax = b. So, while the matrices are different for each frequency, when we propagate the solution one step, the 1024-2048 matrices don&#8217;t change since the frequencies haven&#8217;t changed between steps.</div>
<div>&nbsp;</div>
<div>We would like to store each KSP solution (all 1024-2048 of them) in an array so that we can reuse them for the next step since this saves us time factoring the matrix. The problem I am seeing is that storing these KSP objects in any kind of array increases
the runtime and storage space exponentially &#8211; I am seeing times on the order of 30 mins to an hour and memory consumed hover around 20-30GB! Is there something inherent about the KSP object that I am missing and/or a better way to store these objects? I don&#8217;t
have any problems if I use one KSP object and set the operators each time to the different matrices, except that it takes longer than it normally would after the first step since it has to recalculate the preconditioner matrix.</div>
<div>&nbsp;</div>
<div>Note: We are using distributed matrices and vectors for our linear system.</div>
<div>&nbsp;</div>
<div>I have included some pseudocode to illustrate our linear system, which works:</div>
<div>&nbsp;</div>
<div>KSP k;</div>
<div>for steps = 1 to 100 do:</div>
<div style="text-indent:36pt;">for time = 1 to 1024 do:</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; create distributed matrix A and fill with local values (A[step=1,time=1] = A[step=2,time=1], etc)</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; create distributed vector b and fill with local values for given time</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; create distributed vector x for result</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; create KSP k and set matrices</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; KSPSolve(k, b, x) =&gt; takes a little extra time each step</div>
<div>&nbsp;</div>
<div>However, the following creates time/storage problems:</div>
<div>&nbsp;</div>
<div>KSP* kArray;</div>
<div>create kArray =&gt; takes inordinate amout of time and storage</div>
<div>for steps = 1 to 100 do:</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; for time=1 to 1024 do:</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; create distributed matrix A and fill with local values (A[step=1,time=1] = A[step=2,time=1], etc)</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; create distributed vector b and fill with local values for given time</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; create distributed vector x for result</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if time == 1 do:</div>
<div style="text-indent:36pt;padding-left:72pt;">set matrices on kArray[time] =&gt; also takes an inordinate amount of time and storage</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; KSPSolve(kArray[time], b, x)</div>
<div>&nbsp;</div>
<div>Any insight would be greatly appreciated!</div>
<div>Brenna</div>
<div>&nbsp;</div>
</span></font>
</body>
</html>