Lines Matching refs:fres
177 subroutine twohfilterB (y, x, strnrm, ien, fres, argument
182 dimension y(nshg,ndof), fres(nshg,33)
394 fres(ien(nel,j),:) = fres(ien(nel,j),:) + fresl(nel,:)
593 subroutine twohfilterBB (y, x, strnrm, ien, fres, argument
598 dimension y(nshg,ndof), fres(nshg,33)
771 fres(ien(nel,j),:) = fres(ien(nel,j),:) + fresl(nel,:)
1659 dimension fres(nshg,24), fwr(nshg),
1708 fres = zero
1743 call asithf (yold, x, strl(iel:inum,:), mien(iblk)%p, fres,
1779 if(numpe>1) call commu (fres, ilwork, 24, 'in ')
1786 fres(i,:) = fres(i,:) + fres(j,:)
1792 fres(j,:) = fres(i,:)
1796 if(numpe>1) call commu (fres, ilwork, 24, 'out')
1798 fres(:,23) = one / fres(:,23)
1800 fres(:,j) = fres(:,j) * fres(:,23)
1809 & two * (fres(:,10)**2 + fres(:,11)**2 + fres(:,12)**2)
1810 & + four * ( fres(:,13)**2 + fres(:,14)**2 + fres(:,15)**2 ) )
1812 fwr = fwr1 * fres(:,22) * strnrm
1815 & * fres(:,10) + fres(:,16)
1817 & * fres(:,11) + fres(:,17)
1819 & * fres(:,12) + fres(:,18)
1821 xmij(:,4) = -fwr * fres(:,13) + fres(:,19)
1822 xmij(:,5) = -fwr * fres(:,14) + fres(:,20)
1823 xmij(:,6) = -fwr * fres(:,15) + fres(:,21)
1825 fres(:,22) = one / fres(:,22)
1827 xlij(:,1) = fres(:,4) - fres(:,1) * fres(:,1) * fres(:,22)
1828 xlij(:,2) = fres(:,5) - fres(:,2) * fres(:,2) * fres(:,22)
1829 xlij(:,3) = fres(:,6) - fres(:,3) * fres(:,3) * fres(:,22)
1830 xlij(:,4) = fres(:,7) - fres(:,1) * fres(:,2) * fres(:,22)
1831 xlij(:,5) = fres(:,8) - fres(:,1) * fres(:,3) * fres(:,22)
1832 xlij(:,6) = fres(:,9) - fres(:,2) * fres(:,3) * fres(:,22)
1846 xfac(:,1) = strnrm*strnrm*( fres(:,10)**2 + fres(:,11)**2 +
1847 & fres(:,12)**2
1848 & + two*( fres(:,13)**2 + fres(:,14)**2 + fres(:,15)**2 ) )
1850 xfac(:,2) = strnrm*( xlij(:,1)*fres(:,10) + xlij(:,2)*fres(:,11)
1851 & + xlij(:,3)*fres(:,12) +
1852 & two*(xlij(:,4)*fres(:,13) + xlij(:,5)*fres(:,14) +
1853 & xlij(:,6)*fres(:,15)) )
1855 xfac(:,3) = strnrm*( fres(:,10)*fres(:,16) + fres(:,11)*fres(:,17)
1856 & + fres(:,12)*fres(:,18) +
1857 & two*(fres(:,13)*fres(:,19) + fres(:,14)*fres(:,20) +
1858 & fres(:,15)*fres(:,21)) )
1860 xfac(:,4) = xlij(:,1)*fres(:,16) + xlij(:,2)*fres(:,17)
1861 & + xlij(:,3)*fres(:,18) +
1862 & two*(xlij(:,4)*fres(:,19) + xlij(:,5)*fres(:,20) +
1863 & xlij(:,6)*fres(:,21))
1865 xfac(:,5) = fres(:,16)*fres(:,16) + fres(:,17)*fres(:,17)
1866 & + fres(:,18)*fres(:,18) +
1867 & two*(fres(:,19)*fres(:,19) + fres(:,20)*fres(:,20) +
1868 & fres(:,21)*fres(:,21))
1879 fres(j,:) = zero
1906 fres(isgbeg:isgend,:) = zero
1962 xl1(ifath(i),1) = xl1(ifath(i),1) + fres(i,4)
1963 xl1(ifath(i),2) = xl1(ifath(i),2) + fres(i,5)
1964 xl1(ifath(i),3) = xl1(ifath(i),3) + fres(i,6)
1965 xl1(ifath(i),4) = xl1(ifath(i),4) + fres(i,7)
1966 xl1(ifath(i),5) = xl1(ifath(i),5) + fres(i,8)
1967 xl1(ifath(i),6) = xl1(ifath(i),6) + fres(i,9)
1969 xl2(ifath(i),1) = xl2(ifath(i),1) + fres(i,1)*fres(i,1)
1970 xl2(ifath(i),2) = xl2(ifath(i),2) + fres(i,2)*fres(i,2)
1971 xl2(ifath(i),3) = xl2(ifath(i),3) + fres(i,3)*fres(i,3)
1972 xl2(ifath(i),4) = xl2(ifath(i),4) + fres(i,1)*fres(i,2)
1973 xl2(ifath(i),5) = xl2(ifath(i),5) + fres(i,1)*fres(i,3)
1974 xl2(ifath(i),6) = xl2(ifath(i),6) + fres(i,2)*fres(i,3)
1976 ui(ifath(i),1) = ui(ifath(i),1) + fres(i,1)
1977 ui(ifath(i),2) = ui(ifath(i),2) + fres(i,2)
1978 ui(ifath(i),3) = ui(ifath(i),3) + fres(i,3)
2344 dimension fres(nshg,33), fwr(nshg),
2401 fres = zero
2491 & fres, hfres, shglf(lcsyst,:,1:nshl,:),
2498 if(numpe>1) call commu (fres, ilwork, 33, 'in ')
2505 fres(i,:) = fres(i,:) + fres(j,:)
2512 fres(j,:) = fres(i,:)
2517 call commu (fres, ilwork, 33, 'out')
2520 fres(:,22) = one / fres(:,22)
2522 fres(:,j) = fres(:,j) * fres(:,22)
2525 fres(:,j) = fres(:,j) * fres(:,22)
2551 & two * (fres(:,10)**2 + fres(:,11)**2 + fres(:,12)**2)
2552 & + four * ( fres(:,13)**2 + fres(:,14)**2 + fres(:,15)**2 ) )
2557 & * fres(:,10) + fres(:,16)
2559 & * fres(:,11) + fres(:,17)
2561 & * fres(:,12) + fres(:,18)
2563 xmij(:,4) = -fwr * fres(:,13) + fres(:,19)
2564 xmij(:,5) = -fwr * fres(:,14) + fres(:,20)
2565 xmij(:,6) = -fwr * fres(:,15) + fres(:,21)
2568 xlij(:,1) = fres(:,4) - fres(:,1) * fres(:,1)
2569 xlij(:,2) = fres(:,5) - fres(:,2) * fres(:,2)
2570 xlij(:,3) = fres(:,6) - fres(:,3) * fres(:,3)
2571 xlij(:,4) = fres(:,7) - fres(:,1) * fres(:,2)
2572 xlij(:,5) = fres(:,8) - fres(:,1) * fres(:,3)
2573 xlij(:,6) = fres(:,9) - fres(:,2) * fres(:,3)
2587 xfac(:,1) = strnrm*strnrm*( fres(:,10)**2 + fres(:,11)**2 +
2588 & fres(:,12)**2
2589 & + two*( fres(:,13)**2 + fres(:,14)**2 + fres(:,15)**2 ) )
2591 xfac(:,2) = strnrm*( xlij(:,1)*fres(:,10) + xlij(:,2)*fres(:,11)
2592 & + xlij(:,3)*fres(:,12) +
2593 & two*(xlij(:,4)*fres(:,13) + xlij(:,5)*fres(:,14) +
2594 & xlij(:,6)*fres(:,15)) )
2596 xfac(:,3) = strnrm*( fres(:,10)*fres(:,16) + fres(:,11)*fres(:,17)
2597 & + fres(:,12)*fres(:,18) +
2598 & two*(fres(:,13)*fres(:,19) + fres(:,14)*fres(:,20) +
2599 & fres(:,15)*fres(:,21)) )
2601 xfac(:,4) = xlij(:,1)*fres(:,16) + xlij(:,2)*fres(:,17)
2602 & + xlij(:,3)*fres(:,18) +
2603 & two*(xlij(:,4)*fres(:,19) + xlij(:,5)*fres(:,20) +
2604 & xlij(:,6)*fres(:,21))
2606 xfac(:,5) = fres(:,16)*fres(:,16) + fres(:,17)*fres(:,17)
2607 & + fres(:,18)*fres(:,18) +
2608 & two*(fres(:,19)*fres(:,19) + fres(:,20)*fres(:,20) +
2609 & fres(:,21)*fres(:,21))
2620 fres(j,:) = zero
2647 fres(isgbeg:isgend,:) = zero
2703 xl1(ifath(i),1) = xl1(ifath(i),1) + fres(i,4)
2704 xl1(ifath(i),2) = xl1(ifath(i),2) + fres(i,5)
2705 xl1(ifath(i),3) = xl1(ifath(i),3) + fres(i,6)
2706 xl1(ifath(i),4) = xl1(ifath(i),4) + fres(i,7)
2707 xl1(ifath(i),5) = xl1(ifath(i),5) + fres(i,8)
2708 xl1(ifath(i),6) = xl1(ifath(i),6) + fres(i,9)
2710 xl2(ifath(i),1) = xl2(ifath(i),1) + fres(i,1)*fres(i,1)
2711 xl2(ifath(i),2) = xl2(ifath(i),2) + fres(i,2)*fres(i,2)
2712 xl2(ifath(i),3) = xl2(ifath(i),3) + fres(i,3)*fres(i,3)
2713 xl2(ifath(i),4) = xl2(ifath(i),4) + fres(i,1)*fres(i,2)
2714 xl2(ifath(i),5) = xl2(ifath(i),5) + fres(i,1)*fres(i,3)
2715 xl2(ifath(i),6) = xl2(ifath(i),6) + fres(i,2)*fres(i,3)
2717 ui(ifath(i),1) = ui(ifath(i),1) + fres(i,1)
2718 ui(ifath(i),2) = ui(ifath(i),2) + fres(i,2)
2719 ui(ifath(i),3) = ui(ifath(i),3) + fres(i,3)