Commit d70f3d46 authored by Claus Kleinwort's avatar Claus Kleinwort
Browse files

protection for records without global derivatives

git-svn-id: http://svnsrv.desy.de/public/MillepedeII/trunk@179 3547b9b0-65b8-46d3-b95d-921b3f43af62
parent 57e2c28f
...@@ -5296,6 +5296,7 @@ SUBROUTINE loop2 ...@@ -5296,6 +5296,7 @@ SUBROUTINE loop2
INTEGER(mpi) :: nmatmo INTEGER(mpi) :: nmatmo
INTEGER(mpi) :: noff INTEGER(mpi) :: noff
INTEGER(mpi) :: nr INTEGER(mpi) :: nr
INTEGER(mpi) :: nrece
INTEGER(mpi) :: nrecf INTEGER(mpi) :: nrecf
INTEGER(mpi) :: nrecmm INTEGER(mpi) :: nrecmm
INTEGER(mpi) :: nst INTEGER(mpi) :: nst
...@@ -5398,6 +5399,7 @@ SUBROUTINE loop2 ...@@ -5398,6 +5399,7 @@ SUBROUTINE loop2
ENDIF ENDIF
! reading events===reading events===reading events===reading events= ! reading events===reading events===reading events===reading events=
nrece =0 ! empty records (no global derivatives)
nrecf =0 ! records with fixed global parameters nrecf =0 ! records with fixed global parameters
naeqng=0 ! count number of equations (with global der.) naeqng=0 ! count number of equations (with global der.)
naeqnf=0 ! count number of equations ( " , fixed) naeqnf=0 ! count number of equations ( " , fixed)
...@@ -5560,6 +5562,11 @@ SUBROUTINE loop2 ...@@ -5560,6 +5562,11 @@ SUBROUTINE loop2
dstat(1)=dstat(1)+REAL((nwrd+2)*2,mpd) ! record size dstat(1)=dstat(1)+REAL((nwrd+2)*2,mpd) ! record size
dstat(2)=dstat(2)+REAL(nagbn+2,mpd) ! indices, dstat(2)=dstat(2)+REAL(nagbn+2,mpd) ! indices,
dstat(3)=dstat(3)+REAL(nagbn*nagbn+nagbn,mpd) ! data for MUPDAT dstat(3)=dstat(3)+REAL(nagbn*nagbn+nagbn,mpd) ! data for MUPDAT
IF (nagbn == 0) THEN
nrece=nrece+1
CYCLE
ENDIF
CALL sort1k(globalIndexUsage,nagbn) ! sort global par. CALL sort1k(globalIndexUsage,nagbn) ! sort global par.
! overwrite read buffer with lists of global labels ! overwrite read buffer with lists of global labels
...@@ -5837,6 +5844,10 @@ SUBROUTINE loop2 ...@@ -5837,6 +5844,10 @@ SUBROUTINE loop2
WRITE(lu,101) 'NRECF',nrecf, & WRITE(lu,101) 'NRECF',nrecf, &
'number of records with fixed global derivatives' 'number of records with fixed global derivatives'
END IF END IF
IF (nrece > 0) THEN
WRITE(lu,101) 'NRECE',nrece, &
'number of records without global derivatives'
END IF
IF (ncache > 0) THEN IF (ncache > 0) THEN
WRITE(lu,101) 'NCACHE',ncache,'number of words for caching' WRITE(lu,101) 'NCACHE',ncache,'number of words for caching'
WRITE(lu,111) (fcache(k)*100.0,k=1,3) WRITE(lu,111) (fcache(k)*100.0,k=1,3)
......
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