Changeset 38d77eb for partem_p.f


Ignore:
Timestamp:
11/19/09 11:29:48 (14 years ago)
Author:
baerbaer <baerbaer@…>
Branches:
master
Children:
7137e5d
Parents:
32289cd
Message:

Redirected standard out to logString.

SMMP produced a lot of log messages. This became an issue when run in massively
parallel environments. I replaced all writes to standard out to a write to logString.
The next step is to pass this string to a function that writes the messages to a log
file according to their importance and the chosen log level.

git-svn-id: svn+ssh://svn.berlios.de/svnroot/repos/smmp/trunk@34 26dc1dd8-5c4e-0410-9ffe-d298b4865968

File:
1 edited

Legend:

Unmodified
Added
Removed
  • partem_p.f

    r32289cd r38d77eb  
    9595      trackID = 1
    9696      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,
    9999     &            rep_id, num_rep, partem_comm, myrank
    100100      call flush(6)
     
    167167            do i=1,num_rep
    168168               read(13,*) j,inode(i),intem(i),yol(i),e_minp(i),h_maxp(i)
    169                write (*,*) "par_R.in:",i,j
     169               write (logString, *) "par_R.in:",i,j
    170170            end do
    171171            jold=(iold/nmes)*num_rep
     
    178178     &              ,idum1, idum2, idum3, e_min
    179179     &              ,dummy, dummy
    180                write (*,*) i
     180               write (logString, *) i
    181181               call flush(6)
    182182            end do
     
    197197      e_min = e_minp(rep_id+1)
    198198      h_max = h_maxp(rep_id+1)
    199       write (*,*) "E_min=",e_min," for ", rep_id + 1
     199      write (logString, *) "E_min=",e_min," for ", rep_id + 1
    200200      eol=energy()
    201201      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), eol
     202         write (logString, *) rep_id, ' Warning: yol(rep_id).ne.eol:'
     203         write (logString, *) rep_id, yol(rep_id + 1), eol
    204204      endif
    205205!     Start of simulation
    206       write (*,*) '[',rep_id, myrank, beta, partem_comm,
     206      write (logString, *) '[',rep_id, myrank, beta, partem_comm,
    207207     &            '] Energy before equilibration:', eol
    208208!     =====================Equilibration by canonical Metropolis
     
    211211      end do
    212212      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
    214215      call flush(6)
    215216!     
     
    222223         eol0 = eol
    223224         if (myrank.eq.0.and.rep_id.eq.0) then
    224             write (*,*) "Finished sweep", nsw
     225            write (logString, *) "Finished sweep", nsw
    225226            call flush(6)
    226227         endif
     
    494495     &                IERR)
    495496               if (myrank.ne.0) then
    496 !                  write (*,*) '[', myrank,'] Missed', randomCount,
     497!                  write (logString, *) '[', myrank,'] Missed', randomCount,
    497498!     &                            'random numbers.'
    498499                  do i = 1, randomCount
    499500                     rd = grnd()
    500 !                     write (*,*) '[', myrank,'] rd=', rd
     501!                     write (logString, *) '[', myrank,'] rd=', rd
    501502                  enddo
    502503               endif
     
    519520      e_final=energy()
    520521      if (partem_comm.ne.MPI_COMM_NULL) then
    521          write (*,*) rep_id, ' E_final', e_final
     522         write (logString, *) rep_id, ' E_final', e_final
    522523      endif
    523524      eol0 = eol
Note: See TracChangeset for help on using the changeset viewer.