From bd7b7b126a32993789d296a9be5ee368ad76b156 Mon Sep 17 00:00:00 2001 From: dadaoqiuzhi <56991198+dadaoqiuzhi@users.noreply.github.com> Date: Sun, 21 Jan 2024 10:08:19 +0800 Subject: [PATCH] debug --- chemi_mechanism/species_bonds_lammpstrj_framecheck.m | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/chemi_mechanism/species_bonds_lammpstrj_framecheck.m b/chemi_mechanism/species_bonds_lammpstrj_framecheck.m index e937ab1..bf911ff 100644 --- a/chemi_mechanism/species_bonds_lammpstrj_framecheck.m +++ b/chemi_mechanism/species_bonds_lammpstrj_framecheck.m @@ -176,7 +176,7 @@ species_frame_check2='n'; if strcmpi(species_frame_checkans,'y') for iijj=3:size(outputdatanew,1) - if mod(outputdatanew{iijj,1}-outputdatanew{iijj-1,1},trajper)~=0 + if mod(outputdatanew{iijj,1}-outputdatanew{iijj-1,1},trajper(1))~=0 fprintf('From the %d th frame, viz from %d timestep the arithmetical progression is damaged\n',iijj,outputdatanew{iijj,1}) fprintf('Frame No. of the %d th frame to %d th frame are:%d %d %d\n',iijj-1,iijj+1,outputdatanew{iijj-1,1},outputdatanew{iijj,1},outputdatanew{iijj+1,1}) msgbox('The arithmetical progression is damaged in species.* file. How to cope with the problem?') @@ -190,7 +190,7 @@ end if strcmpi(species_frame_check2,'y') for jjkk=iijj+1:size(outputdatanew,1) - if mod(outputdatanew{jjkk,1}-outputdatanew{jjkk-1,1},trajper)~=0 + if mod(outputdatanew{jjkk,1}-outputdatanew{jjkk-1,1},trajper(1))~=0 fprintf('For the destructive arithmetical progression, from %d th frame on, viz %d timestep, the arithmetical progression is damage again\n',jjkk,outputdatanew{jjkk,1}) msgbox('The arithmetical progression is damage again!') break @@ -204,11 +204,11 @@ if strcmpi(species_frame_check3,'y') for jjkk=ijk:size(outputdatanew,1) if tartrajectory==outputdatanew{jjkk,1} - tartrajectory=outputdatanew{jjkk-1,1}+trajper; - promptans3=outputdatanew{jjkk-1,1}+trajper; + tartrajectory=outputdatanew{jjkk-1,1}+trajper(1); + promptans3=outputdatanew{jjkk-1,1}+trajper(1); fprintf('The initial target frame No. is modified to %d once more\n',tartrajectory) end - outputdatanew{jjkk,1}=outputdatanew{jjkk-1,1}+trajper; + outputdatanew{jjkk,1}=outputdatanew{jjkk-1,1}+trajper(1); end fprintf('The arithmetical progression problem is repaired for species.* file (outputnew), but not always guaranteed for 100% sure\n') end