256 250 MHZ IP27 Processors ------------------------------------------------------- Amber 6 SANDER Scripps/UCSF 1999 ------------------------------------------------------- | Sat Dec 9 19:35:25 2000 [-O]verwriting output File Assignments: |MDIN : /u/ac/jphillip/amber/pme.amber |MDOUT: /u/ac/jphillip/amber/bench_amber_252.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= 521791 | Real memory used by local nonb = 386115 | MAX NONBOND PAIRS = 8000000 | Memory Use Allocated Used | Real 4000000 1460123 | Hollerith 1200000 148532 | Integer 4000000 1187456 | 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 97 188 286 374 483 576 659 | 762 838 938 1031 1131 1223 1306 1401 | 1493 1586 1688 1769 1872 1957 2057 2143 | 2240 2341 2425 2513 2606 2699 2792 2885 | 2978 3071 3164 3257 3350 3443 3536 3629 | 3722 3815 3908 4001 4094 4187 4280 4373 | 4466 4559 4652 4745 4838 4931 5024 5117 | 5210 5303 5396 5489 5582 5675 5768 5861 | 5954 6047 6140 6233 6326 6419 6512 6605 | 6698 6791 6884 6977 7070 7163 7256 7349 | 7442 7535 7628 7721 7814 7907 8000 8093 | 8186 8279 8372 8465 8558 8651 8744 8837 | 8930 9023 9116 9209 9302 9395 9488 9581 | 9674 9767 9860 9953 10046 10139 10232 10325 | 10418 10511 10604 10697 10790 10883 10976 11069 | 11162 11255 11348 11441 11534 11627 11720 11813 | 11906 11999 12092 12185 12278 12371 12464 12557 | 12650 12743 12836 12929 13022 13115 13208 13301 | 13394 13487 13580 13673 13766 13859 13952 14045 | 14138 14231 14324 14417 14510 14603 14696 14789 | 14882 14975 15068 15161 15254 15347 15440 15533 | 15626 15719 15812 15905 15998 16091 16184 16277 | 16370 16463 16556 16649 16742 16835 16928 17021 | 17114 17207 17300 17393 17486 17579 17672 17765 | 17858 17951 18044 18137 18230 18323 18416 18509 | 18602 18695 18788 18881 18974 19067 19160 19253 | 19346 19439 19532 19625 19718 19811 19904 19997 | 20090 20183 20276 20369 20462 20555 20648 20741 | 20834 20927 21020 21113 21206 21299 21392 21485 | 21578 21671 21764 21857 21950 22043 22136 22229 | 22322 22415 22508 22601 22694 22787 22880 22973 | 23066 23159 23252 23345 23558 | Atom division among processors for gb: | 0 93 186 279 372 465 558 651 | 744 837 930 1023 1116 1209 1302 1395 | 1488 1581 1674 1767 1860 1953 2046 2139 | 2232 2325 2418 2511 2604 2697 2790 2883 | 2976 3069 3162 3255 3348 3441 3534 3627 | 3720 3813 3906 3999 4092 4185 4278 4371 | 4464 4557 4650 4743 4836 4929 5022 5115 | 5208 5301 5394 5487 5580 5673 5766 5859 | 5952 6045 6138 6231 6324 6417 6510 6603 | 6696 6789 6882 6975 7068 7161 7254 7347 | 7440 7533 7626 7719 7812 7905 7998 8091 | 8184 8277 8370 8463 8556 8649 8742 8835 | 8928 9021 9114 9207 9300 9393 9486 9579 | 9672 9765 9858 9951 10044 10137 10230 10323 | 10416 10509 10602 10695 10788 10881 10974 11067 | 11160 11253 11346 11439 11532 11625 11718 11811 | 11904 11997 12090 12183 12276 12369 12462 12555 | 12648 12741 12834 12927 13020 13113 13206 13299 | 13392 13485 13578 13671 13764 13857 13950 14043 | 14136 14229 14322 14415 14508 14601 14694 14787 | 14880 14973 15066 15159 15252 15345 15438 15531 | 15624 15717 15810 15903 15996 16089 16182 16275 | 16368 16461 16554 16647 16740 16833 16926 17019 | 17112 17205 17298 17391 17484 17577 17670 17763 | 17856 17949 18042 18135 18228 18321 18414 18507 | 18600 18693 18786 18879 18972 19065 19158 19251 | 19344 19437 19530 19623 19716 19809 19902 19995 | 20088 20181 20274 20367 20460 20553 20646 20739 | 20832 20925 21018 21111 21204 21297 21390 21483 | 21576 21669 21762 21855 21948 22041 22134 22227 | 22320 22413 22506 22599 22692 22785 22878 22971 | 23064 23157 23250 23343 23558 | Running AMBER/MPI version on 252 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.73142052 ------------------------------------------------------------------------------ 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 = 376.675 TOTAL TIME = 376.675 | Nonbond Ewald pairlist time components: | map, save_t = 0.53 0.14 | grid setup_t = 0.00 0.00 | grid ucell_t = 0.13 0.04 | grid image_t = 0.05 0.01 | build list_t = 0.47 0.12 | total list_t = 1.18 0.31 Nonbond Ewald force time components: | zero ene,force = 0.43 0.12 | map = 0.00 0.00 | map,adjust = 3.22 0.85 | self energy = 0.02 0.01 | 1-dim b-spline= 2.30 0.61 | grid charges = 1.58 0.42 | scalar sum = 1.20 0.32 | grad sum = 2.64 0.70 | FFT = 156.07 41.43 | Recip force = 163.81 43.49 | Direct force = 4.89 1.30 | adjust masked = 0.53 0.14 | accum force = 1.55 0.41 | finish virial = 4.55 1.21 | mfct4 = 0.51 0.14 | f_wait = 2.37 0.63 | stack_time = 0.46 0.12 | total force = 182.36 48.41 | | Routine Sec % | ---------------------------- | Nonbond 190.55 50.59 | Bond 0.52 0.14 | Angle 0.75 0.20 | Dihedral 1.61 0.43 | Shake 1.52 0.40 | GBrad 0.00 0.00 | Force 0.00 0.00 | GBraddist 0.00 0.00 | F,Xdist 153.31 40.70 | Other 16.03 4.26 | ---------------------------- | Total 376.67 0.10 Hours | Nonsetup 362.65 96.28% | Highest rstack allocated: 70674 | MAX_RSTACK = 10000000 | Highest istack allocated: 1963 | MAX_ISTACK = 200000 | Setup wallclock 22 seconds | Nonsetup wallclock 378 seconds