FA</span>
    r_psolve = @(state) incompTPFA_DFM(state, g_fine, S, fluid, <span class="string">'Wells'</span>, W,<span class="string">'src'</span>,src,<span class="string">'c2cTrans'</span>,S2,<span class="string">'MatrixOutput'</span>,true);


    <span class="comment">% MsFV using a few GMRES iterations to improve flux error</span>
    psolve   = @(state) solveMSFV_TPFA_Incomp_DFM(state, g_fine, CG, S, fluid, <span class="keyword">...</span>
        <span class="string">'Reconstruct'</span>, true, <span c