Changeset 38d77eb for partem_p.f
- Timestamp:
- 11/19/09 11:29:48 (14 years ago)
- Branches:
- master
- Children:
- 7137e5d
- Parents:
- 32289cd
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
partem_p.f
r32289cd r38d77eb 95 95 trackID = 1 96 96 odd = 1 97 write ( *,*) 'Starting parallel tempering.'98 write ( *,*) 'parameters, ',switch,newsta,nmes,nswp,nmes,97 write (logString, *) 'Starting parallel tempering.' 98 write (logString, *) 'parameters, ',switch,newsta,nmes,nswp,nmes, 99 99 & rep_id, num_rep, partem_comm, myrank 100 100 call flush(6) … … 167 167 do i=1,num_rep 168 168 read(13,*) j,inode(i),intem(i),yol(i),e_minp(i),h_maxp(i) 169 write ( *,*) "par_R.in:",i,j169 write (logString, *) "par_R.in:",i,j 170 170 end do 171 171 jold=(iold/nmes)*num_rep … … 178 178 & ,idum1, idum2, idum3, e_min 179 179 & ,dummy, dummy 180 write ( *,*) i180 write (logString, *) i 181 181 call flush(6) 182 182 end do … … 197 197 e_min = e_minp(rep_id+1) 198 198 h_max = h_maxp(rep_id+1) 199 write ( *,*) "E_min=",e_min," for ", rep_id + 1199 write (logString, *) "E_min=",e_min," for ", rep_id + 1 200 200 eol=energy() 201 201 if(.not.newsta.and.abs(yol(rep_id + 1) - eol).gt.0.1) then 202 write (*,*) rep_id, ' Warning: yol(rep_id).ne.eol:'203 write (*,*) rep_id, yol(rep_id + 1), eol202 write (logString, *) rep_id, ' Warning: yol(rep_id).ne.eol:' 203 write (logString, *) rep_id, yol(rep_id + 1), eol 204 204 endif 205 205 ! Start of simulation 206 write ( *,*) '[',rep_id, myrank, beta, partem_comm,206 write (logString, *) '[',rep_id, myrank, beta, partem_comm, 207 207 & '] Energy before equilibration:', eol 208 208 ! =====================Equilibration by canonical Metropolis … … 211 211 end do 212 212 CALL MPI_BARRIER(MPI_COMM_WORLD,IERR) 213 write (*,*) '[',rep_id,'] Energy after equilibration:', eol 213 write (logString, *) '[',rep_id,'] Energy after equilibration:', 214 & eol 214 215 call flush(6) 215 216 ! … … 222 223 eol0 = eol 223 224 if (myrank.eq.0.and.rep_id.eq.0) then 224 write ( *,*) "Finished sweep", nsw225 write (logString, *) "Finished sweep", nsw 225 226 call flush(6) 226 227 endif … … 494 495 & IERR) 495 496 if (myrank.ne.0) then 496 ! write ( *,*) '[', myrank,'] Missed', randomCount,497 ! write (logString, *) '[', myrank,'] Missed', randomCount, 497 498 ! & 'random numbers.' 498 499 do i = 1, randomCount 499 500 rd = grnd() 500 ! write ( *,*) '[', myrank,'] rd=', rd501 ! write (logString, *) '[', myrank,'] rd=', rd 501 502 enddo 502 503 endif … … 519 520 e_final=energy() 520 521 if (partem_comm.ne.MPI_COMM_NULL) then 521 write ( *,*) rep_id, ' E_final', e_final522 write (logString, *) rep_id, ' E_final', e_final 522 523 endif 523 524 eol0 = eol
Note:
See TracChangeset
for help on using the changeset viewer.