xmd Code

55 lines (43 with data), 1.5 kB

2012-06-23
Add repot to particle structure to save reference total
potential energy; will be used by MC commands.

2012-10-17
ADDING MONTE CARLO

  Remove previously added 'repot' variable above.

  New Monte-Carlo comands
   (x)    MC START_RUN      -  Reset counters and energy sum, initializes energy, sets run flag
   (x)    MC START_TRIAL    -  Test run flag, save current particle coords and types
   (_)    SWAP args         -  Alters current particles.
   (x)    MC END_TRIAL temp -  Calc new energy, accepts or rejects move
   (x)    MC WRITE          -  Write MC statistics

  Usage:

   MC START_RUN
   REPEAT 1000
     MC START_TRIAL
     ..move some particles..
     MC END_TRIAL 1000
     MC WRITE
   END
   MC WRITE

2012-10-27
  New MC commands:  
    cdsubs.c :  read_swap(), read_mc()
    cdmc.c   :  (called by read_mc()  )  -> mc_start_run(), mc_start_trial(), mc_end_trial(), mc_write()
    sortsub.c:  (called by read_swap())  -> randchoice()
  Compile with out error.

  Next:  Create mcswap example

2012-10-30
  New MC command tested and working.
  (_)  Change version to 2.5.39

2012-11-06
  See test examples/mcswap/ilya/*_dir.
  xmd_dir reproduces simulation in prog_dir.
  See gnuplot graph.

2012-11-08
  (x) Send email to I. (be sure to include attachments).
  (x) Tighten up mcswap example
  (x) Remove test print statements from code
  (x) Make new package 
  (_) Update doc.

2012-11-20
  (x) Initialize indices to NULL in cdsubs.c, read_swap().