From aef9da49be23a10f9e9bbfcd8a4f1cbf7386966f Mon Sep 17 00:00:00 2001 From: Alex Selimov Date: Tue, 26 May 2020 17:43:12 -0400 Subject: [PATCH] Added greater range of element sizes --- src/mode_create.f90 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/mode_create.f90 b/src/mode_create.f90 index f33efdf..ca7331f 100644 --- a/src/mode_create.f90 +++ b/src/mode_create.f90 @@ -155,7 +155,7 @@ module mode_create call add_element(element_type, curr_esize, 1, 1, r_node_temp) exit end if - curr_esize=esize/(2**ei) - 1 + curr_esize=curr_esize-2 end do end do end if @@ -339,7 +339,7 @@ module mode_create esize_nums=0 do while (curr_esize >= 7) esize_nums=esize_nums+1 - curr_esize = esize/(2**esize_nums) -1 + curr_esize = curr_esize -2 end do else esize_nums=1 @@ -438,7 +438,7 @@ module mode_create !Now build the finite element region lat_ele_num = 0 lat_atom_num = 0 - curr_esize=esize/(2**(esize_nums-1)) + 1 + curr_esize=esize - 2*(esize_nums-1) allocate(r_lat(3,8,numlatpoints/curr_esize)) curr_esize=esize @@ -496,7 +496,7 @@ module mode_create end do end do esize_index(ei) = lat_ele_num - curr_esize=esize/(2**ei) - 1 + curr_esize=curr_esize-2 end do !Now figure out how many lattice points could not be contained in elements allocate(r_atom_lat(3,count(lat_points)))