Commit 0b4e005e authored by Claus Kleinwort's avatar Claus Kleinwort
Browse files

fixed uninitialized values in QL decomposition of constraints matrix

git-svn-id: http://svnsrv.desy.de/public/MillepedeII/trunk@218 3547b9b0-65b8-46d3-b95d-921b3f43af62
parent 51a8651d
...@@ -322,7 +322,8 @@ SUBROUTINE qldecb(a,bpar,bcon,cons) ...@@ -322,7 +322,8 @@ SUBROUTINE qldecb(a,bpar,bcon,cons)
! number of used columns ! number of used columns
ncol=ilast-ipoff ncol=ilast-ipoff
! get column ! get column
vecN(kn)=0.0_mpd in=iblast+k1-iclast
vecN(in:kn)=0.0_mpd
vecN(ifirst:ilast)=matV(ioff1+1:ioff1+ncol) vecN(ifirst:ilast)=matV(ioff1+1:ioff1+ncol)
nrm = SQRT(dot_product(vecN(ifirst:ilast),vecN(ifirst:ilast))) nrm = SQRT(dot_product(vecN(ifirst:ilast),vecN(ifirst:ilast)))
IF (nrm == 0.0_mpd) CYCLE IF (nrm == 0.0_mpd) CYCLE
......
...@@ -51,7 +51,7 @@ ...@@ -51,7 +51,7 @@
!! 1. Download the software package from the DESY \c svn server to !! 1. Download the software package from the DESY \c svn server to
!! \a target directory, e.g.: !! \a target directory, e.g.:
!! !!
!! svn checkout http://svnsrv.desy.de/public/MillepedeII/tags/V04-09-04 target !! svn checkout http://svnsrv.desy.de/public/MillepedeII/tags/V04-09-05 target
!! !!
!! 2. Create **Pede** executable (in \a target directory): !! 2. Create **Pede** executable (in \a target directory):
!! !!
...@@ -1662,7 +1662,7 @@ SUBROUTINE feasma ...@@ -1662,7 +1662,7 @@ SUBROUTINE feasma
WRITE(lunlog,*) 'QL decomposition of constraints matrix' WRITE(lunlog,*) 'QL decomposition of constraints matrix'
CALL monini(lunlog,monpg1,monpg2) CALL monini(lunlog,monpg1,monpg2)
END IF END IF
IF(matsto > 0) THEN ! True unless LAPACK IF(matsto > 0) THEN ! True unless unpacked LAPACK
! QL decomposition ! QL decomposition
CALL qlini(nvgb,ncgb,npblck,mszcon,monpg1) CALL qlini(nvgb,ncgb,npblck,mszcon,monpg1)
! loop over parameter blocks ! loop over parameter blocks
...@@ -3392,8 +3392,8 @@ SUBROUTINE mupdat(i,j,add) ! ...@@ -3392,8 +3392,8 @@ SUBROUTINE mupdat(i,j,add) !
END IF END IF
! bad index? ! bad index?
IF(ij < 0.OR.ij > size(matPreCond)) THEN IF(ij < 0.OR.ij > size(matPreCond)) THEN
CALL peend(23,'Aborted, bad matrix index') CALL peend(23,'Aborted, bad matrix index')
STOP 'mupdat: bad index' STOP 'mupdat: bad index'
END IF END IF
! update? ! update?
IF(ij /= 0) matPreCond(ij)=matPreCond(ij)+add IF(ij /= 0) matPreCond(ij)=matPreCond(ij)+add
...@@ -4976,7 +4976,7 @@ SUBROUTINE avprds(n,l,x,is,ie,b) ...@@ -4976,7 +4976,7 @@ SUBROUTINE avprds(n,l,x,is,ie,b)
lj=1 ! index (in group region) lj=1 ! index (in group region)
DO i=ia,ib DO i=ia,ib
b(i)=b(i)+dot_product(REAL(globalMatF(indij+lj+ja2-ja:indij+lj+jb2-ja),mpd),x(ja:jb)) b(i)=b(i)+dot_product(REAL(globalMatF(indij+lj+ja2-ja:indij+lj+jb2-ja),mpd),x(ja:jb))
lj=lj+jn lj=lj+jn
END DO END DO
END IF END IF
IF (mextnd == 0.AND.ia2 <= ib2) THEN IF (mextnd == 0.AND.ia2 <= ib2) THEN
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment