Dear Amber users and developers,
I received the following error after attempt to build all-to-all (2D) RMSD
plot
from a long trajectory:
CPPTRAJ: Trajectory Analysis. V18.01
___ ___ ___ ___
| \/ | \/ | \/ |
_|_/\_|_/\_|_/\_|_
| Date/time: 01/09/19 14:03:45
| Available memory: 15.367 GB
Loading previous history from log 'cpptraj.log'
> parm abk_oct_spc2_1264_HMR.prmtop
[parm abk_oct_spc2_1264_HMR.prmtop]
Reading 'abk_oct_spc2_1264_HMR.prmtop' as Amber Topology
Radius Set: modified Bondi radii (mbondi)
Warning: Amber topology flag 'LENNARD_JONES_CCOEF' is unrecognized and will
be skipped.
> trajin Abk2_NPT_Prodrun8us_GPU_HMR.mdcrd
[trajin Abk2_NPT_Prodrun8us_GPU_HMR.mdcrd]
Reading 'Abk2_NPT_Prodrun8us_GPU_HMR.mdcrd' as Amber NetCDF
> strip !:1-42
[strip !:1-42]
STRIP: Stripping atoms in mask [!:1-42]
> rms2d .CA out mass Abk2_NPT_Prodrun8us_GPU_HMR.2drms.gnu
[rms2d .CA out mass Abk2_NPT_Prodrun8us_GPU_HMR.2drms.gnu]
RMS2D: COORDS set [_DEFAULTCRD_], mask [.CA], using RMSD (fit),
mass-weighted
Output to 'mass'
> go
[go]
Warning: One or more analyses requested creation of default COORDS DataSet.
CREATECRD: Saving coordinates from Top abk_oct_spc2_1264_HMR.prmtop to
"_DEFAULTCRD_"
---------- RUN BEGIN -------------------------------------------------
PARAMETER FILES (1 total):
0: abk_oct_spc2_1264_HMR.prmtop, 14799 atoms, 4746 res, box: Trunc. Oct.,
4705 mol, 4698 solvent
INPUT TRAJECTORIES (1 total):
0: 'Abk2_NPT_Prodrun8us_GPU_HMR.mdcrd' is a NetCDF AMBER trajectory with
coordinates, time, box, Parm abk_oct_spc2_1264_HMR.prmtop (Trunc. Oct. box)
(reading 800000 of 800000)
Coordinate processing will occur on 800000 frames.
BEGIN TRAJECTORY PROCESSING:
.....................................................
ACTION SETUP FOR PARM 'abk_oct_spc2_1264_HMR.prmtop' (2 actions):
0: [strip !:1-42]
Stripping 14100 atoms.
Stripped topology: 699 atoms, 42 res, box: Trunc. Oct., 1 mol
1: [createcrd _DEFAULTCRD_]
Warning: COORDS data sets do not store times.
Estimated memory usage (800000 frames): 6.730 GB
----- Abk2_NPT_Prodrun8us_GPU_HMR.mdcrd (1-800000, 1) -----
0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% Complete.
Read 800000 frames and processed 800000 frames.
TIME: Avg. throughput= 924.7276 frames / second.
ACTION OUTPUT:
ANALYSIS: Performing 1 analyses:
0: [rms2d .CA out mass Abk2_NPT_Prodrun8us_GPU_HMR.2drms.gnu]
Mask [.CA] corresponds to 42 atoms.
terminate called after throwing an instance of 'std::bad_alloc'
what(): std::bad_alloc
Aborted (core dumped)
The amount of memory should be sufficient.
Could you please tell me what is the possible source of trouble here?
Kind regards,
Nick
Nikolay Kuzmich
Department of Drug Safety,
Research Institute of Influenza,
WHO National Influenza Centre of Russia,
15/17 Professor Popov St.,
Saint-Petersburg
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Jan 09 2019 - 04:00:02 PST
