;      <span class="keyword">subroutine </span>dch1up(n,R,ldr,u,w)</div>
<div class="line"><a name="l00022"></a><span class="lineno">   22</span>&#160;<span class="comment">c purpose:      given an upper triangular matrix R that is a Cholesky</span></div>
<div class="line"><a name="l00023"></a><span class="lineno">   23</span>&#160;<span class="comment">c               factor of a symmetric positive definite matrix A, i.e.</span></div>
<div class="line"><a name="l00024"></a><span class="lineno">   24</span>