1 2 3c$$$ subroutine get_a_not_hex(xc,anot) 4c$$$ 5c$$$ include "common.h" 6c$$$ 7c$$$ dimension xc(npro,nenl,nsd), anot(npro,nenl,nsd) 8c$$$ 9c$$$ 10c$$$ do i = 1, nsd 11c$$$ 12c$$$ anot(:,1,i) = pt125*(xc(:,1,i)+xc(:,2,i)+xc(:,3,i)+xc(:,4,i) 13c$$$ & +xc(:,5,i)+xc(:,6,i)+xc(:,7,i)+xc(:,8,i)) 14c$$$ 15c$$$ anot(:,2,i) = pt125*(-xc(:,1,i)+xc(:,2,i)+xc(:,3,i)-xc(:,4,i) 16c$$$ & -xc(:,5,i)+xc(:,6,i)+xc(:,7,i)-xc(:,8,i)) 17c$$$ 18c$$$ anot(:,3,i) = pt125*(-xc(:,1,i)-xc(:,2,i)+xc(:,3,i)+xc(:,4,i) 19c$$$ & -xc(:,5,i)-xc(:,6,i)+xc(:,7,i)+xc(:,8,i)) 20c$$$ 21c$$$ anot(:,4,i) = pt125*(-xc(:,1,i)-xc(:,2,i)-xc(:,3,i)-xc(:,4,i) 22c$$$ & +xc(:,5,i)+xc(:,6,i)+xc(:,7,i)+xc(:,8,i)) 23c$$$ 24c$$$ anot(:,5,i) = pt125*(xc(:,1,i)-xc(:,2,i)+xc(:,3,i)-xc(:,4,i) 25c$$$ & +xc(:,5,i)-xc(:,6,i)+xc(:,7,i)-xc(:,8,i)) 26c$$$ 27c$$$ anot(:,6,i) = pt125*(xc(:,1,i)+xc(:,2,i)-xc(:,3,i)-xc(:,4,i) 28c$$$ & -xc(:,5,i)-xc(:,6,i)+xc(:,7,i)+xc(:,8,i)) 29c$$$ 30c$$$ anot(:,7,i) = pt125*(xc(:,1,i)-xc(:,2,i)-xc(:,3,i)+xc(:,4,i) 31c$$$ & -xc(:,5,i)+xc(:,6,i)+xc(:,7,i)-xc(:,8,i)) 32c$$$ 33c$$$ anot(:,8,i) = pt125*(-xc(:,1,i)+xc(:,2,i)-xc(:,3,i)+xc(:,4,i) 34c$$$ & +xc(:,5,i)-xc(:,6,i)+xc(:,7,i)-xc(:,8,i)) 35c$$$ 36c$$$ enddo 37c$$$ 38c$$$ return 39c$$$ end 40 41 42 subroutine get_coeff_tet(xc,anot) 43 44 use spebc 45 include "common.h" 46 47 48 dimension xc(nelint,nenl,nsd), anot(nelint,nenl,nsd) 49 50 51 do i = 1, nsd 52 53 anot(:,1,i) = xc(:,4,i) 54 anot(:,2,i) = xc(:,1,i)-xc(:,4,i) 55 anot(:,3,i) = xc(:,2,i)-xc(:,4,i) 56 anot(:,4,i) = xc(:,3,i)-xc(:,4,i) 57 58 enddo 59 60 return 61 end 62 63