From 9c3296aad2b0598636115c82638c356f6df71d59 Mon Sep 17 00:00:00 2001
From: Yossi Eliaz <eliaz123@gmail.com>
Date: Mon, 12 Mar 2018 13:45:13 -0500
Subject: [PATCH] Tools/Matlab: Allow to read LAMPPS output field
Some output fields have attributes attached on the same
line. e.g. "ITEM: BOX BOUNDS pp pp pp". This patch replaced all
the strcmpi to strncmpi in order to limit the number of character
compared with LAMPPS outputs.
Signed-off-by: Yossi Eliaz <eliaz123@gmail.com>
---
tools/matlab/readdump_all.m | 6 +++---
tools/matlab/readdump_one.m | 8 ++++----
2 files changed, 7 insertions(+), 7 deletions(-)
diff --git a/tools/matlab/readdump_all.m b/tools/matlab/readdump_all.m
index 1187beeb73..5a6588507e 100644
--- a/tools/matlab/readdump_all.m
+++ b/tools/matlab/readdump_all.m
@@ -29,13 +29,13 @@ end
i=1;
while feof(dump) == 0
id = fgetl(dump);
- if (strcmpi(id,'ITEM: TIMESTEP'))
+ if (strncmpi(id,'ITEM: TIMESTEP',numel('ITEM: TIMESTEP')))
timestep(i) = str2num(fgetl(dump));
else
- if (strcmpi(id,'ITEM: NUMBER OF ATOMS'))
+ if (strncmpi(id,'ITEM: NUMBER OF ATOMS',numel('ITEM: NUMBER OF ATOMS')))
Natoms(i) = str2num(fgetl(dump));
else
- if (strcmpi(id,'ITEM: BOX BOUNDS'))
+ if (strncmpi(id,'ITEM: BOX BOUNDS',numel('ITEM: BOX BOUNDS')))
x_bound(i,:) = str2num(fgetl(dump));
y_bound(i,:) = str2num(fgetl(dump));
z_bound(i,:) = str2num(fgetl(dump));
diff --git a/tools/matlab/readdump_one.m b/tools/matlab/readdump_one.m
index fabefec909..762ca9bb20 100644
--- a/tools/matlab/readdump_one.m
+++ b/tools/matlab/readdump_one.m
@@ -46,21 +46,21 @@ else
end
while done == 0 & last_status == 0
id = fgetl(dump);
- if (strcmpi(id,'ITEM: TIMESTEP'))
+ if (strncmpi(id,'ITEM: TIMESTEP',numel('ITEM: TIMESTEP')))
if t == 0
timestep(i) = str2num(fgetl(dump));
t=1;
end
else
- if (strcmpi(id,'ITEM: NUMBER OF ATOMS'))
+ if (strcmpi(id,'ITEM: NUMBER OF ATOMS',numel('ITEM: NUMBER OF ATOMS')))
Natoms = str2num(fgetl(dump));
else
- if (strcmpi(id,'ITEM: BOX BOUNDS'))
+ if (strcmpi(id,'ITEM: BOX BOUNDS',numel('ITEM: BOX BOUNDS')))
x_bound(1,:) = str2num(fgetl(dump));
y_bound(1,:) = str2num(fgetl(dump));
z_bound(1,:) = str2num(fgetl(dump));
else
- if (strcmpi(id(1:11),'ITEM: ATOMS'))
+ if (strncmpi('ITEM: ATOMS',numel('ITEM: ATOMS')))
atom_data = zeros(Natoms,ncol);%Allocate memory for atom data
for j = 1 : 1: Natoms
atom_data(j,:) = str2num(fgetl(dump));
--
GitLab