 row and column indexes are the same, stiffness terms are added</span>
         <span class="n">endif</span>
         <span class="n">jr1</span> <span class="o">=</span> <span class="n">jr2</span>                                 <span class="c">! stores (i-1) and (j-1) for further comparison</span>
         <span class="n">ir1</span> <span class="o">=</span> <span class="n">ir2</span>
      <span class="n">enddo</span>
      <span class="n">nz</span> <span class="o">=</span> <span class="n">innz</span>
    