<div><br></div><div><br><div class="gmail_quote"><div dir="ltr">On Tue, 9 Apr 2019 at 10:33, Mark Adams via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">PETSc solvers will cache any setup computations, like a matrix factorization, and will reuse them on subsequent solves. <div><br></div><div>If you call KSPSetOperators, PETSc will assume its cache is invalid and redo any required setup on subsequent solves.</div></div></blockquote><div dir="auto"><br></div><div dir="auto">Following from what Mark said, if you want an exact solve, use -pc_type lu. Note that this will not store the inverse, rather only the L, U factors. During setup, the L,U factors are created and subsequent solves will simply apply the forward and backward solves using these factors.</div><div dir="auto"><br></div><div dir="auto"><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div></div></div><div dir="ltr"><div><div><br></div><div>Mark</div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Apr 9, 2019 at 4:08 AM ztdepyahoo via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div style="font-family:微软雅黑,Verdana,"Microsoft Yahei",SimSun,sans-serif;font-size:14px;line-height:1.6">
<div></div><div>
<div>
<span>Dear group:
</span></div><div><span> I need to solve a Laplacian problem Ax=b every time step where A is a constant . I want to store the inverse of A at the first time step. </span></div><div><span>how to do with Petsc</span></div><div><span> Regards</span></div>
<div>
<span>
<br>
</span>
</div>
<div id="m_-4361030859792712552gmail-m_2754168156661115478ntes-pcmac-signature" style="font-family:微软雅黑">
<div style="font-size:14px;padding:0px;margin:0px;line-height:14px">
<div style="padding-bottom:6px;margin-bottom:10px;border-bottom:1px solid rgb(230,230,230);display:inline-block">
<a href="https://maas.mail.163.com/dashi-web-extend/html/proSignature.html?ftlId=1&name=ztdepyahoo&uid=ztdepyahoo%40gmail.com&iconUrl=https%3A%2F%2Fmail-online.nosdn.127.net%2Fqiyelogo%2FdefaultAvatar.png&items=%5B%22ztdepyahoo%40gmail.com%22%5D" style="display:block;background:rgb(255,255,255);max-width:400px;padding:15px 0px 10px;text-decoration:none;outline:none" target="_blank">
<table cellpadding="0" style="width:100%;max-width:100%;table-layout:fixed;border-collapse:collapse;color:rgb(155,158,161);font-size:14px;line-height:1.3">
<tbody style="word-break:break-all;font-family:"PingFang SC","Hiragino Sans GB","WenQuanYi Micro Hei","Microsoft Yahei",微软雅黑,verdana">
<tr>
<td width="38" style="padding:0px;box-sizing:border-box;width:38px">
<img width="38" height="38" style="vertical-align:middle;width:38px;height:38px;border-radius:50%" src="https://mail-online.nosdn.127.net/qiyelogo/defaultAvatar.png">
</td>
<td style="padding:0px 0px 0px 10px;color:rgb(49,53,59)">
<div style="font-size:16px;font-weight:bold;width:100%;white-space:nowrap;overflow:hidden;text-overflow:ellipsis">ztdepyahoo</div>
</td>
</tr>
<tr width="100%" style="width:100%;font-size:14px">
<td colspan="2" style="padding:10px 0px 0px;width:100%;font-size:14px">
<div style="width:100%;word-break:break-all;font-size:14px">ztdepyahoo@gmail.com</div>
</td>
</tr>
</tbody>
</table>
</a>
</div>
</div>
<div style="font-size:12px;color:rgb(181,185,189);line-height:18px">
<span>签名由</span>
<a style="text-decoration:none;color:rgb(65,150,255);padding:0px 5px" href="https://mail.163.com/dashi/dlpro.html?from=mail81" target="_blank">网易邮箱大师</a>
<span>定制</span>
</div>
</div>
</div>
</div>
</div>
</blockquote></div>
</blockquote></div></div>