Lines Matching refs:intp
701 do intp = 1, ngaussf
713 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
714 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
715 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
716 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
717 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
718 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
719 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
720 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
721 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
751 wght = Qwtf(intp)
756 rho = rho+shp(i,intp)*yl(:,i,1) !density at qpt
765 u1 = u1 + shp(i,intp)*yl(:,i,2)
766 u2 = u2 + shp(i,intp)*yl(:,i,3)
767 u3 = u3 + shp(i,intp)*yl(:,i,4)
773 fresli(:,i,1) = rho * u1 * shp(i,intp) !rho * u1 * WdetJ
774 fresli(:,i,2) = rho * u2 * shp(i,intp) !rho * u2 * WdetJ
775 fresli(:,i,3) = rho * u3 * shp(i,intp) !rho * u3 * WdetJ
776 fresli(:,i,4) = rho * shp(i,intp) !rho * WdetJ
777 fresli(:,i,5) = WdetJ*shp(i,intp) !Integral of filter kernel
779 fresli(:,i,6) = rho * u1 * u1 * shp(i,intp) !rho * u1 * u1 * WdetJ
780 fresli(:,i,7) = rho * u2 * u2 * shp(i,intp) !rho * u2 * u2 * WdetJ
781 fresli(:,i,8) = rho * u3 * u3 * shp(i,intp) !rho * u3 * u3 * WdetJ
782 fresli(:,i,9) = rho * u1 * u2 * shp(i,intp) !rho * u1 * u2 * WdetJ
783 fresli(:,i,10)= rho * u1 * u3 * shp(i,intp) !rho * u1 * u3 * WdetJ
784 fresli(:,i,11)= rho * u2 * u3 * shp(i,intp) !rho * u2 * u3 * WdetJ
856 do intp = 1, ngauss
869 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
870 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
871 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
872 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
873 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
874 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
875 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
876 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
877 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
907 wght=Qwt(lcsyst,intp) ! may be different now
914 fresli(:,22) = fresli(:,22)+shp(i,intp)*yl(:,i,1) ! density at qpt
915 fresli(:,25) = fresli(:,25)+shp(i,intp)*hfresl(:,i,1) !bar(rho u1)
916 fresli(:,26) = fresli(:,26)+shp(i,intp)*hfresl(:,i,2) !bar(rho u2)
917 fresli(:,27) = fresli(:,27)+shp(i,intp)*hfresl(:,i,3) !bar(rho u3)
945 shg(:,n,1) = (shgl(1,n,intp) * dxidx(:,1,1)
946 & + shgl(2,n,intp) * dxidx(:,2,1)
947 & + shgl(3,n,intp) * dxidx(:,3,1))
948 shg(:,n,2) = (shgl(1,n,intp) * dxidx(:,1,2)
949 & + shgl(2,n,intp) * dxidx(:,2,2)
950 & + shgl(3,n,intp) * dxidx(:,3,2))
951 shg(:,n,3) = (shgl(1,n,intp) * dxidx(:,1,3)
952 & + shgl(2,n,intp) * dxidx(:,2,3)
953 & + shgl(3,n,intp) * dxidx(:,3,3))
980 strnrm(:,intp) = fresli(:,22) * sqrt(
1003 u1 = u1 + shp(i,intp)*yl(:,i,2)
1004 u2 = u2 + shp(i,intp)*yl(:,i,3)
1005 u3 = u3 + shp(i,intp)*yl(:,i,4)
1019 fresli(:,16) = strnrm(:,intp) * fresli(:,10) ! rho *|Eps| *Eps11 *WdetJ
1020 fresli(:,17) = strnrm(:,intp) * fresli(:,11) ! rho *|Eps| *Eps22 *WdetJ
1021 fresli(:,18) = strnrm(:,intp) * fresli(:,12) ! rho *|Eps| *Eps33 *WdetJ
1022 fresli(:,19) = strnrm(:,intp) * fresli(:,13) ! rho *|Eps| *Eps12 *WdetJ
1023 fresli(:,20) = strnrm(:,intp) * fresli(:,14) ! rho *|Eps| *Eps13 *WdetJ
1024 fresli(:,21) = strnrm(:,intp) * fresli(:,15) ! rho *|Eps| *Eps23 *WdetJ
1088 do intp = 1, ngauss
1102 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
1103 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
1104 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
1105 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
1106 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
1107 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
1108 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
1109 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
1110 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
1140 wght=Qwt(lcsyst,intp)
1146 shg(:,n,1) = (shgl(1,n,intp) * dxidx(:,1,1)
1147 & + shgl(2,n,intp) * dxidx(:,2,1)
1148 & + shgl(3,n,intp) * dxidx(:,3,1))
1149 shg(:,n,2) = (shgl(1,n,intp) * dxidx(:,1,2)
1150 & + shgl(2,n,intp) * dxidx(:,2,2)
1151 & + shgl(3,n,intp) * dxidx(:,3,2))
1152 shg(:,n,3) = (shgl(1,n,intp) * dxidx(:,1,3)
1153 & + shgl(2,n,intp) * dxidx(:,2,3)
1154 & + shgl(3,n,intp) * dxidx(:,3,3))
1161 fresli(:,22) = fresli(:,22)+shp(i,intp) ! unit density at qpt
1190 strnrm(:,intp) = fresli(:,22) * (
1205 strnrm(:,intp) = strnrm(:,intp) * WdetJ ! ( |Eps|^2 )*WdetJ
1207 epsli(:,1) = -xmudmi(:,intp)*strnrm(:,intp)
1209 epsli(:,2) = rlsl(:,intp,1)*fresli(:,10) +
1210 & rlsl(:,intp,2)*fresli(:,11) +
1211 & rlsl(:,intp,3)*fresli(:,12) +
1212 & two*( rlsl(:,intp,4)*fresli(:,13)+
1213 & rlsl(:,intp,5)*fresli(:,14)+
1214 & rlsl(:,intp,6)*fresli(:,15) )