128 250 MHZ IP27 Processors ------------------------------------------------------- Amber 6 SANDER Scripps/UCSF 1999 ------------------------------------------------------- | Tue Sep 26 17:58:54 2000 [-O]verwriting output File Assignments: |MDIN : /u/ac/jphillip/amber/pme.amber |MDOUT: /u/ac/jphillip/amber/bench_amber_4.log.tmp |INPCR: inpcrd |PARM : prmtop |RESTR: restrt |REFC : refc |MDVEL: mdvel |MDEN : mden |MDCRD: mdcrd |MDINF: mdinfo Here is the input file: short md, nve ensemble &cntrl ntx=1, irest=0, ntc=2, ntf=2, tol=0.0000001, nstlim=500, ntt=0, ntpr=100, ntwr=10000, dt=0.001, vlimit=10.0, tempi=300., temp0 = 300., cut=9., &end &ewald a=62.23, b=62.23, c=62.23, nfft1=64,nfft2=64,nfft3=64, skinnb=2., &end ------------------------------------------------------------------------------- short md, nve ensemble 1. RESOURCE USE: getting box info from bottom of parm getting new box info from bottom of inpcrd | peek_ewald_inpcrd: Box info found EWALD SPECIFIC INPUT: | Using the SGI specific (ZFFT3D) Fast Fourier Transform | Using the T3D specific (FFT3D0) Fast Fourier Transform Largest sphere to fit in unit cell has radius = 31.115 Calculating ew_coeff from dsum_tol,cutoff Box X = 62.230 Box Y = 62.230 Box Z = 62.230 Alpha = 90.000 Beta = 90.000 Gamma = 90.000 NFFT1 = 64 NFFT2 = 64 NFFT3 = 64 Cutoff= 9.000 Tol =0.100E-04 Ewald Coefficient = 0.30768 Interpolation order = 4 NATOM = 23558 NTYPES = 16 NBONH = 22290 MBONA = 1302 NTHETH = 2789 MTHETA = 1772 NPHIH = 5026 MPHIA = 3140 NHPARM = 0 NPARM = 0 NNB = 41859 NRES = 7182 NBONA = 1302 NTHETA = 1772 NPHIA = 3140 NUMBND = 45 NUMANG = 91 NPTRA = 33 NATYP = 31 NPHB = 1 IFBOX = 1 NMXRS = 25 IFCAP = 0 EWALD MEMORY USE: | Total heap storage needed = 1005 | Adjacent nonbond minimum mask = 41859 | Max number of pointers = 25 | List build maxmask = 83718 | Maximage = 34473 EWALD LOCMEM POINTER OFFSETS | Real memory needed by PME = 1005 | Size of EEDTABLE = 20768 | Real memory needed by EEDTABLE = 83072 | Integer memory needed by ADJ = 83718 | Integer memory used by local nonb= 588791 | Real memory used by local nonb = 386115 | MAX NONBOND PAIRS = 8000000 | Memory Use Allocated Used | Real 4000000 1460123 | Hollerith 1200000 148532 | Integer 4000000 1254456 | Max Nonbonded Pairs: 8000000 BOX TYPE: RECTILINEAR 2. CONTROL DATA FOR THE RUN TIMLIM= 999999. IREST = 0 IBELLY= 0 IMIN = 0 IPOL = 0 NTX = 1 NTXO = 1 IG = 71277 TEMPI = 300.00 HEAT = 0.000 NTB = 1 BOXX = 62.230 BOXY = 62.230 BOXZ = 62.230 NTT = 0 TEMP0 = 300.000 DTEMP = 0.000 TAUTP = 1.000 VLIMIT= 10.000 NTP = 0 PRES0 = 1.000 COMP = 44.600 TAUP = 0.200 NPSCAL= 0 NTCM = 0 NSCM = 9999999 NSTLIM= 500 NTU = 1 T = 0.000 DT = 0.00100 NTC = 2 TOL = 0.00000 JFASTW = 0 NTF = 2 NSNB = 25 CUT = 9.000 SCNB = 2.000 SCEE = 1.200 DIELC = 1.000 NTPR = 100 NTWR = 10000 NTWX = 0 NTWV = 0 NTWE = 0 IOUTFM= 0 NTWPRT= 0 NTWPR0= 0 NTAVE= 0 NTR = 0 NTRX = 1 TAUR = 0.00000 NMROPT= 0 PENCUT= 0.10000 IVCAP = 0 MATCAP= 0 FCAP = 1.500 OTHER DATA: IFCAP = 0 NATCAP= 0 CUTCAP= 0.000 XCAP = 0.000 YCAP = 0.000 ZCAP = 0.000 VRAND= 0 NATOM = 23558 NRES = 7182 Water definition for fast triangulated model: Resname = WAT ; Oxygen_name = O ; Hyd1_name = H1 ; Hyd2_name = H2 | PLEVEL = 1: runmd parallelization, no EKCMR 3. ATOMIC COORDINATES AND VELOCITIES begin time read from input coords = 0.000 ps Number of triangulated 3-point waters found: 7023 | Atom division among processors: | 0 5891 11780 17669 23558 | Atom division among processors for gb: | 0 5889 11778 17667 23558 | Running AMBER/MPI version on 4 nodes Sum of charges from parm topology file = -11.00000006 Forcing neutrality... --------------------------------------------------- APPROXIMATING switch and d/dx switch using CUBIC SPLINE INTERPOLATION using 5000.0 points per unit in tabled values TESTING RELATIVE ERROR over r ranging from 0.0 to cutoff | CHECK switch(x): max rel err = 0.2088E-14 at 2.598900 | CHECK d/dx switch(x): max rel err = 0.7671E-11 at 2.757160 --------------------------------------------------- Total number of mask terms = 34709 Total number of mask terms = 69418 | Total Ewald setup time = 0.19996834 ------------------------------------------------------------------------------ NSTEP = 0 TIME(PS) = 0.000 TEMP(K) = 436.89 PRESS = 0.00 Etot = -45424.6591 EKtot = 21002.9065 EPtot = -66427.5656 BOND = 98.5019 ANGLE = 410.5387 DIHED = 805.7590 1-4 NB = 614.1164 1-4 EEL = 7104.4424 VDWAALS = 4790.9648 EELEC = -80251.8887 EHBOND = 0.0000 CONSTRAINT = 0.0000 Ewald error estimate: 0.1425E-02 ------------------------------------------------------------------------------ NSTEP = 100 TIME(PS) = 0.100 TEMP(K) = 341.94 PRESS = 0.00 Etot = -52070.7160 EKtot = 16438.2467 EPtot = -68508.9627 BOND = 424.6502 ANGLE = 1145.5619 DIHED = 1030.9361 1-4 NB = 533.8425 1-4 EEL = 6477.1604 VDWAALS = 8251.7818 EELEC = -86372.8957 EHBOND = 0.0000 CONSTRAINT = 0.0000 Ewald error estimate: 0.2890E-03 ------------------------------------------------------------------------------ NSTEP = 200 TIME(PS) = 0.200 TEMP(K) = 342.89 PRESS = 0.00 Etot = -52070.7614 EKtot = 16484.2530 EPtot = -68555.0144 BOND = 361.6936 ANGLE = 1146.3981 DIHED = 1045.0411 1-4 NB = 566.2224 1-4 EEL = 6564.9541 VDWAALS = 8047.9505 EELEC = -86287.2742 EHBOND = 0.0000 CONSTRAINT = 0.0000 Ewald error estimate: 0.3991E-04 ------------------------------------------------------------------------------ NSTEP = 300 TIME(PS) = 0.300 TEMP(K) = 339.45 PRESS = 0.00 Etot = -52071.0780 EKtot = 16318.7018 EPtot = -68389.7798 BOND = 370.7418 ANGLE = 1132.3021 DIHED = 1063.8323 1-4 NB = 547.3142 1-4 EEL = 6549.3791 VDWAALS = 7698.6647 EELEC = -85752.0141 EHBOND = 0.0000 CONSTRAINT = 0.0000 Ewald error estimate: 0.8891E-04 ------------------------------------------------------------------------------ NSTEP = 400 TIME(PS) = 0.400 TEMP(K) = 339.83 PRESS = 0.00 Etot = -52071.6976 EKtot = 16337.0121 EPtot = -68408.7097 BOND = 388.1232 ANGLE = 1145.1296 DIHED = 1027.5420 1-4 NB = 560.6767 1-4 EEL = 6597.0788 VDWAALS = 7791.9710 EELEC = -85919.2310 EHBOND = 0.0000 CONSTRAINT = 0.0000 Ewald error estimate: 0.2532E-04 ------------------------------------------------------------------------------ NSTEP = 500 TIME(PS) = 0.500 TEMP(K) = 340.83 PRESS = 0.00 Etot = -52070.5632 EKtot = 16385.3016 EPtot = -68455.8647 BOND = 415.5202 ANGLE = 1186.8791 DIHED = 1014.8376 1-4 NB = 551.1868 1-4 EEL = 6547.1283 VDWAALS = 7569.0142 EELEC = -85740.4309 EHBOND = 0.0000 CONSTRAINT = 0.0000 Ewald error estimate: 0.4445E-04 ------------------------------------------------------------------------------ A V E R A G E S O V E R 500 S T E P S NSTEP = 500 TIME(PS) = 0.500 TEMP(K) = 339.13 PRESS = 0.00 Etot = -52057.9046 EKtot = 16303.3183 EPtot = -68361.2229 BOND = 378.8697 ANGLE = 1118.0702 DIHED = 1038.3148 1-4 NB = 556.9583 1-4 EEL = 6581.3192 VDWAALS = 7537.2669 EELEC = -85572.0220 EHBOND = 0.0000 CONSTRAINT = 0.0000 Ewald error estimate: 0.1715E-03 ------------------------------------------------------------------------------ R M S F L U C T U A T I O N S NSTEP = 500 TIME(PS) = 0.500 TEMP(K) = 8.99 PRESS = 0.00 Etot = 297.1036 EKtot = 432.0046 EPtot = 388.8519 BOND = 35.3622 ANGLE = 70.9069 DIHED = 28.4969 1-4 NB = 27.5367 1-4 EEL = 83.8194 VDWAALS = 890.8266 EELEC = 1253.1135 EHBOND = 0.0000 CONSTRAINT = 0.0000 Ewald error estimate: 0.2284E-03 ------------------------------------------------------------------------------ | | ELAPSED TIME = 536.992 TOTAL TIME = 536.992 | Nonbond Ewald pairlist time components: | map, save_t = 0.53 0.10 | grid setup_t = 0.00 0.00 | grid ucell_t = 0.13 0.02 | grid image_t = 0.12 0.02 | build list_t = 26.81 4.99 | total list_t = 27.59 5.14 Nonbond Ewald force time components: | zero ene,force = 0.46 0.08 | map = 0.00 0.00 | map,adjust = 6.53 1.22 | self energy = 0.01 0.00 | 1-dim b-spline= 6.62 1.23 | grid charges = 22.59 4.21 | scalar sum = 16.76 3.12 | grad sum = 34.32 6.39 | FFT = 62.20 11.58 | Recip force = 142.48 26.53 | Direct force = 271.97 50.65 | adjust masked = 3.05 0.57 | accum force = 0.17 0.03 | finish virial = 3.32 0.62 | mfct4 = 1.19 0.22 | f_wait = 21.89 4.08 | stack_time = 0.05 0.01 | total force = 451.10 84.01 | | Routine Sec % | ---------------------------- | Nonbond 484.42 90.21 | Bond 0.25 0.05 | Angle 1.69 0.31 | Dihedral 8.25 1.54 | Shake 3.35 0.62 | GBrad 0.00 0.00 | Force 0.00 0.00 | GBraddist 0.00 0.00 | F,Xdist 21.67 4.04 | Other 14.59 2.72 | ---------------------------- | Total 536.99 0.15 Hours | Nonsetup 524.31 97.64% | Highest rstack allocated: 595105 | MAX_RSTACK = 2400000 | Highest istack allocated: 9325 | MAX_ISTACK = 200000 | Setup wallclock 12 seconds | Nonsetup wallclock 526 seconds