minor improvements 10Nov2023

mualal · Nov 10, 2023 · 5e3b3c1 · 5e3b3c1
1 parent 63d8bde
commit 5e3b3c1
Show file tree

Hide file tree

Showing 2 changed files with 33,402 additions and 38,635 deletions.
diff --git a/matlab/lattice_lattice_interface.m b/matlab/lattice_lattice_interface.m
@@ -1,3 +1,4 @@
+%clc;
 clear;
 
 % time step to save disp
@@ -98,7 +99,6 @@
 
 % equations of motion integration
 for t=times
-    disp_copy = disp;
     %disp(int32(length(num_y)/a), int32(length(num_x)/a)) = sin(omega*t);
     %vel(int32(length(num_y)/a), int32(length(num_x)/a)) = omega*cos(omega*t);
     vel = vel + c./mass.*(circshift(disp,[-1 0])+circshift(disp,[1 0])+circshift(disp,[0 -1])+circshift(disp,[0 1])-4*disp).*dt;