="hl codehilite"><pre><span></span><span class="w">    </span><span class="k">subroutine </span><span class="n">Incomplete_Cholesky_decomposition</span><span class="p">(</span><span class="n">A</span><span class="p">,</span><span class="w"> </span><span class="n">L</span><span class="p">,</span><span class="w"> </span><span class="n">level</span><span class="p">)</span>

<span class="w">        </span><span class="kt">real</span><span class="p">(</span><span class="n">dp</span><span class="p">),</span><span