<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">You might want to give a look at KSPSetDMActive() to use a DM description of the problem without forming your operator.<div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; ">Best,<br>Luc</span>
</div>
<br><div><div>On 29 janv. 2014, at 11:09, Konstantinos Kontzialis wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>
<div>   Dear all,</div>
<div> </div>
<div>  I am using a matrix free method with a user defined preconditioner. </div>
<div>  I run with the option -snes_mf_operator, but I get the following:</div>
<div> </div>
<div>  Must call DMShellSetMatrix() or DMShellSetCreateMatrix()!</div>
<div> </div>
<div>  Why is this happening? I use the TS library.</div>
<div> </div>
  Kostas</div>

</blockquote></div><br></div></body></html>