aselimov/cacmb
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Working changes. Working on the efill code

Browse source
This commit is contained in:
Alex Selimov 2020-10-23 17:17:40 -04:00
parent b5629b1563
commit 22e250093b
2 changed files with 153 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

4
src/elements.f90
View file

@ -747,8 +747,8 @@ module elements
s = (1.0_dp*(j-1)-(size_ele(ie)-1)/2)/(1.0_dp*(size_ele(ie)-1)/2)
t = (1.0_dp*(k-1)-(size_ele(ie)-1)/2)/(1.0_dp*(size_ele(ie)-1)/2)
r(:) = 0
do ibasis = 1, bnum
do inod = 1, 8
do ibasis = 1, basisnum(lat_ele(ie))
do inod = 1, ng_node(lat_ele(ie))
r(:,ibasis) = r(:,ibasis) + a_shape(inod) * r_node(:,ibasis,inod,ie)
end do
end do