<div>
<div>Dear all,</div>
<div> </div>
<div>I am trying to use preconditioning in SNES within a matrix free conext. I use petsc 3.3 and whenever I use the option</div>
<div>-snes_mf_operator I get the following error:</div>
<div> </div>
<div>Must call DMShellSetMatrix() or DMShellSetCreateMatrix()<br /><br data-mce-bogus="1" /></div>
<div> </div>
<div>I code the following:</div>
<div> </div>
<div>      call TSCreate (petsc_comm_world, ts_mhd, ierpetsc)      <br />c           <br />      call TSSetProblemType (ts_mhd, TS_NONLINEAR, ierpetsc)     <br />c<br />      call TSSetIFunction ( ts_mhd, res_mhd, residual_mag,      <br />     @  PETSC_NULL_OBJECT, ierpetsc )      <br />c     <br />      call TSSetSolution( ts_mhd, Bmagnetic_pet, ierpetsc )<br />c<br />      call TSSetMaxSNESFailures ( ts_mhd, -1, ierpetsc)<br />c              <br />      call TSGetSNES (ts_mhd, snes_mhd, ierpetsc )</div>
<div> </div>
<div>      call MatCreateSNESMF ( snes_mhd, J_mf, ierpetsc )      <br />c      <br />      call SNESSetJacobian ( snes_mhd, J_mf, M_mhd,<br />    @ SNESDefaultComputeJacobianColor, fdcoloring,<br />    @ ierpetsc )</div>
<div> </div>
<div>Has anyone any suggestions on what to do?</div>
<div> </div>
<div>Kostas</div>
</div>