scaling of a dense symmetric positive definite matrix a+e using the
factorization a+e=l*d*trans(l) obtained by the subroutine mxdpgf.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer, | intent(in) | :: | n |
order of the matrix a. |
||
| real(kind=wp), | intent(inout) | :: | a(*) |
|
||
| real(kind=wp), | intent(in) | :: | alf |
scaling factor. |
pure subroutine mxdpgs(n,a,alf) real(wp),intent(in) :: alf !! scaling factor. integer,intent(in) :: n !! order of the matrix a. real(wp),intent(inout) :: a(*) !! `a(n*(n+1)/2)` factorization `a+e=l*d*trans(l)` !! obtained by the subroutine [[mxdpgf]]. integer :: i , j j = 0 do i = 1 , n j = j + i a(j) = a(j)*alf end do end subroutine mxdpgs