Changeset 6650a56
- Timestamp:
- 11/19/09 11:29:38 (14 years ago)
- Branches:
- master
- Children:
- 32289cd
- Parents:
- 30644a9
- Files:
-
- 15 edited
Legend:
- Unmodified
- Added
- Removed
-
EXAMPLES/Makefile
r30644a9 r6650a56 17 17 # Linker flags 18 18 # Default flags, e.g., for gfortran, g77 19 L_FLAGS=-O2 -g19 # L_FLAGS=-O2 -g 20 20 # Flags for Intel's ifort 21 #L_FLAGS=-O3 -axW -g21 L_FLAGS=-O3 -axW -g 22 22 # FLAGS for Cray 23 23 # L_FLAGS=-fastsse -Mipa … … 27 27 # Compiler flags 28 28 # Default flags, e.g., for gfortran, g77 29 F_FLAGS=-c -O2 -g29 # F_FLAGS=-c -O2 -g 30 30 # Flags for Intel's ifort 31 #F_FLAGS=-c -O3 -axW -g31 F_FLAGS=-c -O3 -axW -g 32 32 # Flags for debugging build 33 33 # F_FLAGS=-c -O0 -g … … 38 38 # F_FLAGS=-c -fPIC -O3 -funroll-loops -mmmx -msse2 -msse -march=i686 -malign-double -fomit-frame-pointer 39 39 # Select your compiler 40 F90=gfortran41 #F90=ifort40 # F90=gfortran 41 F90=ifort 42 42 #F90=pgf90 43 43 #F90=xlf -
EXAMPLES/multicanonical.f
r30644a9 r6650a56 31 31 ! Change the following directory path to where you want to put SMMP 32 32 ! libraries of residues. 33 libdir='../ SMMP/'33 libdir='../../SMMP/' 34 34 35 35 !! Choose energy type with the following switch … … 58 58 iabin = 1 ! =0: read from PDB-file 59 59 ! =1: ab Initio from sequence (& variables) 60 seqfile=' enkefa.seq'60 seqfile='1vp.seq' 61 61 ! varfile='EXAMPLES/enkefa.var' 62 62 varfile = ' ' … … 77 77 temp=1000.0 78 78 79 call mulcan_par(nsweeps, nup, temp, kmin, kmax, binWidth, l_iter)79 ! call mulcan_par(nsweeps, nup, temp, kmin, kmax, binWidth, l_iter) 80 80 ! Uncomment the following three lines to perform a multi canonical simulation with 81 81 ! the parameters stored in muca.d 82 !nmes = 100083 !call mulcan_sim(nmes, nsweeps, nmes, nup, kmin, kmax, binWidth,84 !& l_iter)82 nmes = 1000 83 call mulcan_sim(nmes, nsweeps, nmes, nup, kmin, kmax, binWidth, 84 & l_iter) 85 85 end program multicanonical_main 86 86 -
INCP.H
r30644a9 r6650a56 2 2 ! contents of a PDB file 3 3 ! ...................... 4 integer MXCHP, MXRSP, MXATP, nchp, nchrsp, nrsp, irsatp, nrsatp 5 integer natp, noatp 6 7 double precision xatp, yatp, zatp 4 8 5 9 parameter (MXCHP =100, ! max. no. of polypeptide chains -
anneal.f
r30644a9 r6650a56 32 32 ! external rand 33 33 external can_weight 34 35 double precision bmin, bmax, db, dv, grnd, vr, addang, eol, energy 36 double precision acz, ymin, vlvrm, rgy, ee, temp 37 38 integer nresi, i, iv, nsw, nemin, j 39 34 40 ! parameter(lrand=.true.) 35 41 ! parameter(nequi=100, nswp=100000,nmes=1000) -
canon.f
r30644a9 r6650a56 27 27 ! external rand 28 28 external can_weight 29 30 double precision dv, grnd, vr, addang, eol, energy, acz, rgy, ee 31 32 integer i, iv, nsw, nhel, mhel, nbet, mbet, mhb 33 29 34 30 35 logical lrand … … 132 137 ! CALLS: none 133 138 ! 139 implicit none 140 double precision beta, x 134 141 135 implicit real*8 (a-h,o-z)136 142 137 143 common/bet/beta -
doc/manual.lyx
r30644a9 r6650a56 1 #LyX 1.5. 0rc2created this file. For more info see http://www.lyx.org/1 #LyX 1.5.6 created this file. For more info see http://www.lyx.org/ 2 2 \lyxformat 276 3 3 \begin_document … … 39 39 \tracking_changes false 40 40 \output_changes false 41 \author "Jan H. Meinke" 41 \author "" 42 \author "" 42 43 \end_header 43 44 … … 2536 2537 original ECEPP/2/3 and FLEX potentials. 2537 2538 However, it should be noted that SMMP has its own numeration of the atom 2538 types and the types of torsions (see [1,2] for more details). 2539 types and the types of torsions (see 2540 \begin_inset LatexCommand cite 2541 key "Eisenmenger2001,Eisenmenger2006" 2542 2543 \end_inset 2544 2545 for more details). 2539 2546 \end_layout 2540 2547 -
helix.f
r30644a9 r6650a56 27 27 28 28 logical lhel,lbet 29 double precision hlim, philim, psilim, hlim2, philim2, psilim2 30 double precision xphi, xpsi 31 32 integer nhel, mhel, nbet, mbet, i, iv 29 33 30 34 parameter(hlim=30.0d0,philim=-70.0d0,psilim=-37.0d0) -
incl_lund.h
r30644a9 r6650a56 1 integer nhptyp, mxconr, mxcell, iN, iCa, iC, mlvr, iphi, ipsi 1 2 parameter (nhptyp=3,mxconr=40,mxcell=125000) 2 3 dimension iN(mxrs),iCa(mxrs),iC(mxrs),mlvr(mxvr) … … 13 14 common /bgs_i/ bgsnvar,iph 14 15 common /bgs_r/ abgs,bbgs,dph, bgsvar 16 integer ihpat, nhpat 15 17 dimension ihpat(mxrs,6),nhpat(mxrs) 16 18 double precision hpstrg … … 18 20 19 21 double precision exvk,exvcut,exvcut2 22 integer matcon 20 23 dimension matcon(-mxconr:mxconr,mxat) 21 24 … … 30 33 dimension sig2exv(mxtyat,mxtyat) 31 34 dimension asaexv(mxtyat,mxtyat),bsaexv(mxtyat,mxtyat) 32 35 integer lcp1, lcp2, ilpst, ilpnd 33 36 common /lundff/kbias, & 34 37 & epshb1,epshb2,powa,powb,sighb,cthb, & -
init_lund.f90
r30644a9 r6650a56 6 6 include 'INCL.H' 7 7 include 'incl_lund.h' 8 integer i, npprs, j, k 9 8 10 logical bgsposs 9 11 do i=1,mxrs -
minqsn.f
r30644a9 r6650a56 56 56 ! ............................................................. 57 57 58 implicit real*8 (a-h,o-z) 59 implicit integer*4 (i-n) 58 implicit none 59 60 double precision eps1, eps2, tiny, zero, dff, c, g, scal, h, fa, f 61 double precision xa, x, ga, d, w, dga, di, fmin, gmin, stmin 62 double precision stepub, steplb, acur, step, xb, fb, gb, gl1, gl2 63 double precision si, dgb, sig 64 65 integer izero, ione, mxn, nfun, itr, i, n, n1, i1, j, isfv, maxfun 66 integer ir 67 60 68 61 69 parameter ( eps1=0.1d0, … … 249 257 ! CALLS: none 250 258 ! 251 implicit real*8 (a-h,o-z) 252 implicit integer*4 (i-n) 259 implicit none 260 261 double precision sig, ti, w, z, a, v, eps, tim, al, r, b, gm, y 262 263 integer mxn, n, ir, np, ij, mk, i, j, mm 253 264 254 265 dimension a(mxn*(mxn+1)/2),z(mxn),w(mxn) … … 373 384 ! CALLS: none 374 385 ! 375 implicit real*8 (a-h,o-z) 376 implicit integer*4 (i-n) 386 implicit none 387 double precision w, z, v, a 388 389 integer mxn, ir, n, i, ij, i1, j, np, nip, ii, ip 377 390 378 391 dimension a(mxn*(mxn+1)/2),z(mxn),w(mxn) -
mklist.f
r30644a9 r6650a56 310 310 ! .................................................... 311 311 312 implicit integer*4 (i-n) 312 implicit none 313 integer mx, j, n, j1, l1, j2, l2, i1, i2, ja, k 313 314 314 315 dimension l1(mx),l2(mx) -
mulcan_par_mod.f90
r30644a9 r6650a56 44 44 real*8, intent(in) :: binWidth, temp 45 45 logical, intent(in) :: l_iter 46 46 47 double precision :: grnd, addang, eol, energy 48 47 49 integer :: i,j 48 50 integer :: iv,nsw,muold … … 214 216 real*8, allocatable :: xhist(:),g1(:),ent(:) 215 217 integer, allocatable :: ihist(:) 216 218 219 double precision :: eol_old, x, grnd, addang, eol, energy, dham 220 221 integer :: nresi, nswm, nhy, nhx 217 222 integer :: i,j 218 223 integer :: iv,nsw,muold -
opesol.f
r30644a9 r6650a56 32 32 33 33 include 'INCL.H' 34 34 double precision xfvr, yfvr, zfvr, xfrvr, yfrvr, zfrvr, esolan, dx 35 double precision dy, dz, xfat, yfat, zfat, xi, yi, zi, xfrat 36 double precision yfrat, zfrat, xb, yb, zb, ex, ey, ez, xfiv, yfiv 37 double precision zfiv, xfriv, yfriv, zfriv 38 39 integer ntlvr, nml, ix2, ifivr, ilavr, i, i1s, i1a, io, iv, it, ia 40 integer ib, i2s, ims, i1, i2, i2a, iad, lad, ivw, j, i14 41 35 42 dimension xfat(mxat),yfat(mxat),zfat(mxat), 36 43 & xfrat(mxat),yfrat(mxat),zfrat(mxat), … … 244 251 245 252 include 'INCL.H' 253 254 double precision del, vlvrx, ovr, eynw, esolan, gda, gdn 255 256 integer i, it, iv, nml 257 246 258 247 259 parameter (del=1.d-6) -
partem_p.f
r30644a9 r6650a56 38 38 ! nmes: number of Monte Carlo sweeps between measurments 39 39 ! newsta: .true. for new simulations, .false. for re-start 40 40 double precision temp, eavm, sph, geavm, gsph, dv, grnd, vr 41 double precision addang, dummy, eol, energy, acz, rmsv, rmsdfun 42 double precision rgy, ee, tmhb, dham, swp, wij, rd, e_final 43 44 integer ifrrm, nmes, nswp, num_rep, i, j, nresi, iold, inode 45 integer intem, iv, jold, idum1, idum2, idum3, mpi_integer 46 integer mpi_comm_world, ierr, mpi_double_precision, nsw, nequi 47 integer nml, nhel, mhel, nbet, mbet, mhb, imhb, nctot, ncnat 48 integer mpi_comm_null, k1, k, nu, no1, in, jn 49 41 50 dimension eavm(MAX_PROC),sph(MAX_PROC),intem(MAX_PROC), 42 51 & inode(MAX_PROC), geavm(MAX_PROC), gsph(MAX_PROC) -
redstr.f
r30644a9 r6650a56 27 27 ! .......................................................... 28 28 integer ib, ie, l 29 integer ln, is, ibegst, j, iendst 30 29 31 character spr,blnk,str*(*),strn*(*) 30 32 … … 91 93 ! ............................................................. 92 94 93 implicit integer*4 (i-n)94 95 implicit none 96 integer i 95 97 character blnk,str*(*) 96 98 data blnk/' '/ … … 117 119 ! .............................................................. 118 120 119 implicit integer*4 (i-n)120 121 implicit none 122 integer i 121 123 character blnk,str*(*) 122 124 data blnk/' '/ … … 146 148 ! .......................................... 147 149 148 implicit integer*4 (i-n)149 150 implicit none 151 integer mxd, ib, ibegst, ie, iendst, il, i0, i9, i, ii, in 150 152 parameter (mxd=9) ! max. # of digits 151 153 … … 186 188 ! .......................................... 187 189 188 implicit integer*4 (i-n)189 190 implicit none 191 integer mxd, mxap, mxip, ib, ibegst, ie, iendst, il, ip, ibp 190 192 parameter (mxd =17, ! max. # of digits 191 193 & mxap= 6, ! max. # of digits after period … … 292 294 293 295 integer lun 294 296 integer i, ibegst, j, k 295 297 logical exs 296 298 character*(*) filnam,stat,format
Note:
See TracChangeset
for help on using the changeset viewer.