aselimov/CACmb
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Added greater range of element sizes

Browse Source
This commit is contained in:
Alex Selimov 2020-05-26 17:43:12 -04:00
parent 080f19d644
commit aef9da49be
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
src/mode_create.f90
View File

@ -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)))