From a528dc4c5245d4a6a4d1c033547b920dd1eed4f8 Mon Sep 17 00:00:00 2001 From: Alex Selimov Date: Wed, 28 Oct 2020 17:53:12 -0400 Subject: [PATCH] Final fix to get working opt_slip_plane --- src/opt_slip_plane.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/opt_slip_plane.f90 b/src/opt_slip_plane.f90 index e607ae8..80da930 100644 --- a/src/opt_slip_plane.f90 +++ b/src/opt_slip_plane.f90 @@ -86,7 +86,7 @@ module opt_slip_plane end do !Check to make sure all lattice points exist for the current element - if(any(lat_points(m:m+esize-1, n:n+esize-1, o:o+esize-1) == 0)) cycle latloop + if(any(.not.lat_points(m:m+esize-1, n:n+esize-1, o:o+esize-1))) cycle latloop !Check to see if the plane intersects this element if not then add it maxp = maxval(rfill(sdim,1:basisnum(lat_ele(ie)),1:ng_node(lat_ele(ie)))) minp = minval(rfill(sdim,1:basisnum(lat_ele(ie)),1:ng_node(lat_ele(ie))))