Lines Matching refs:fres
512 dimension fres(nshg,24), fwr(nshg),
545 fres = zero
578 call hfilterBB (yold, x, mien(iblk)%p, fres,
606 if(numpe>1) call commu (fres, ilwork, 24, 'in ')
613 fres(i,:) = fres(i,:) + fres(j,:)
620 fres(j,:) = zero
643 fres(isgbeg:isgend,:) = zero
660 call sparseCG (fres(:,i), pfres(:,i), lhsG,
669 fres(:,i) = pfres(:,i)
672 fres(:,22) = one
693 & two * (fres(:,10)**2 + fres(:,11)**2 + fres(:,12)**2)
694 & + four * ( fres(:,13)**2 + fres(:,14)**2 + fres(:,15)**2 ) )
696 fwr = fwr1 * fres(:,22) * strnrm
699 & * fres(:,10) + fres(:,16)
701 & * fres(:,11) + fres(:,17)
703 & * fres(:,12) + fres(:,18)
705 xmij(:,4) = -fwr * fres(:,13) + fres(:,19)
706 xmij(:,5) = -fwr * fres(:,14) + fres(:,20)
707 xmij(:,6) = -fwr * fres(:,15) + fres(:,21)
709 fres(:,22) = one / fres(:,22)
711 xlij(:,1) = fres(:,4) - fres(:,1) * fres(:,1) * fres(:,22)
712 xlij(:,2) = fres(:,5) - fres(:,2) * fres(:,2) * fres(:,22)
713 xlij(:,3) = fres(:,6) - fres(:,3) * fres(:,3) * fres(:,22)
714 xlij(:,4) = fres(:,7) - fres(:,1) * fres(:,2) * fres(:,22)
715 xlij(:,5) = fres(:,8) - fres(:,1) * fres(:,3) * fres(:,22)
716 xlij(:,6) = fres(:,9) - fres(:,2) * fres(:,3) * fres(:,22)
958 dimension fres(nshg,33), fwr(nshg),
975 fres = zero
1062 & fres, hfres, shgl(lcsyst,:,1:nshl,:),
1069 if(numpe>1) call commu (fres, ilwork, 33, 'in ')
1076 fres(i,:) = fres(i,:) + fres(j,:)
1083 fres(j,:) = fres(i,:)
1088 call commu (fres, ilwork, 33, 'out')
1091 fres(:,22) = one / fres(:,22)
1093 fres(:,j) = fres(:,j) * fres(:,22)
1096 fres(:,j) = fres(:,j) * fres(:,22)
1122 & two * (fres(:,10)**2 + fres(:,11)**2 + fres(:,12)**2)
1123 & + four * ( fres(:,13)**2 + fres(:,14)**2 + fres(:,15)**2 ) )
1128 & * fres(:,10) + fres(:,16)
1130 & * fres(:,11) + fres(:,17)
1132 & * fres(:,12) + fres(:,18)
1134 xmij(:,4) = -fwr * fres(:,13) + fres(:,19)
1135 xmij(:,5) = -fwr * fres(:,14) + fres(:,20)
1136 xmij(:,6) = -fwr * fres(:,15) + fres(:,21)
1139 xlij(:,1) = fres(:,4) - fres(:,1) * fres(:,1)
1140 xlij(:,2) = fres(:,5) - fres(:,2) * fres(:,2)
1141 xlij(:,3) = fres(:,6) - fres(:,3) * fres(:,3)
1142 xlij(:,4) = fres(:,7) - fres(:,1) * fres(:,2)
1143 xlij(:,5) = fres(:,8) - fres(:,1) * fres(:,3)
1144 xlij(:,6) = fres(:,9) - fres(:,2) * fres(:,3)
1642 & fres(nshg,24), pfres(nshg,24),
1657 fres = zero
1688 call resSij (yold, x, mien(iblk)%p, fres,
1717 if(numpe>1) call commu (fres, ilwork, 24, 'in ')
1724 fres(i,:) = fres(i,:) + fres(j,:)
1730 fres(j,:) = fres(i,:)
1734 if(numpe>1) call commu (fres, ilwork, 24, 'out')
1736 fres(:,22) = one / fres(:,22)
1738 fres(:,j) = fres(:,j) * fres(:,22)
1740 pfres = fres
2127 dimension fres(nshg,24), fwr(nshg),
2145 fres = zero
2175 call asithf (yold, x, strl(iel:inum,:), mien(iblk)%p, fres,
2182 if(numpe>1) call commu (fres, ilwork, 24, 'in ')
2189 fres(i,:) = fres(i,:) + fres(j,:)
2195 fres(j,:) = fres(i,:)
2199 if(numpe>1) call commu (fres, ilwork, 24, 'out')
2201 fres(:,23) = one / fres(:,23)
2203 fres(:,j) = fres(:,j) * fres(:,23)
2210 xlij(:,1) = fres(:,4) - fres(:,1)*fres(:,1)
2211 xlij(:,2) = fres(:,5) - fres(:,2)*fres(:,2)
2212 xlij(:,3) = fres(:,6) - fres(:,3)*fres(:,3)
2213 xlij(:,4) = fres(:,7) - fres(:,1)*fres(:,2)
2214 xlij(:,5) = fres(:,8) - fres(:,1)*fres(:,3)
2215 xlij(:,6) = fres(:,9) - fres(:,2)*fres(:,3)
2218 & two * (fres(:,10)**2 + fres(:,11)**2 + fres(:,12)**2)
2219 & + four * ( fres(:,13)**2 + fres(:,14)**2 + fres(:,15)**2 ) )
2221 xfac(:,1) = strnrm*strnrm*( fres(:,10)**2 + fres(:,11)**2 +
2222 & fres(:,12)**2
2223 & + two*( fres(:,13)**2 + fres(:,14)**2 + fres(:,15)**2 ) )
2225 xfac(:,2) = strnrm*( xlij(:,1)*fres(:,10) + xlij(:,2)*fres(:,11)
2226 & + xlij(:,3)*fres(:,12) +
2227 & two*(xlij(:,4)*fres(:,13) + xlij(:,5)*fres(:,14) +
2228 & xlij(:,6)*fres(:,15)) )
2230 xfac(:,3) = strnrm*( fres(:,10)*fres(:,16) + fres(:,11)*fres(:,17)
2231 & + fres(:,12)*fres(:,18) +
2232 & two*(fres(:,13)*fres(:,19) + fres(:,14)*fres(:,20) +
2233 & fres(:,15)*fres(:,21)) )
2235 xfac(:,4) = xlij(:,1)*fres(:,16) + xlij(:,2)*fres(:,17)
2236 & + xlij(:,3)*fres(:,18) +
2237 & two*(xlij(:,4)*fres(:,19) + xlij(:,5)*fres(:,20) +
2238 & xlij(:,6)*fres(:,21))
2240 xfac(:,5) = fres(:,16)*fres(:,16) + fres(:,17)*fres(:,17)
2241 & + fres(:,18)*fres(:,18) +
2242 & two*(fres(:,19)*fres(:,19) + fres(:,20)*fres(:,20) +
2243 & fres(:,21)*fres(:,21))
2478 dimension fres(nshg,24), fwr(nshg),
2528 fres = zero
2562 call asithf (yold, x, strl(iel:inum,:), mien(iblk)%p, fres,
2593 if(numpe>1) call commu (fres, ilwork, 24, 'in ')
2600 fres(i,:) = fres(i,:) + fres(j,:)
2606 fres(j,:) = fres(i,:)
2610 if(numpe>1) call commu (fres, ilwork, 24, 'out')
2612 fres(:,23) = one / fres(:,23)
2614 fres(:,j) = fres(:,j) * fres(:,23)
2623 & two * (fres(:,10)**2 + fres(:,11)**2 + fres(:,12)**2)
2624 & + four * ( fres(:,13)**2 + fres(:,14)**2 + fres(:,15)**2 ) )
2627 fwr = fwr3 * fres(:,22) * strnrm
2630 & * fres(:,10) + fres(:,16)
2632 & * fres(:,11) + fres(:,17)
2634 & * fres(:,12) + fres(:,18)
2636 xmij(:,4) = -fwr * fres(:,13) + fres(:,19)
2637 xmij(:,5) = -fwr * fres(:,14) + fres(:,20)
2638 xmij(:,6) = -fwr * fres(:,15) + fres(:,21)
2640 fres(:,22) = one / fres(:,22)
2642 xlij(:,1) = fres(:,4) - fres(:,1) * fres(:,1) * fres(:,22)
2643 xlij(:,2) = fres(:,5) - fres(:,2) * fres(:,2) * fres(:,22)
2644 xlij(:,3) = fres(:,6) - fres(:,3) * fres(:,3) * fres(:,22)
2645 xlij(:,4) = fres(:,7) - fres(:,1) * fres(:,2) * fres(:,22)
2646 xlij(:,5) = fres(:,8) - fres(:,1) * fres(:,3) * fres(:,22)
2647 xlij(:,6) = fres(:,9) - fres(:,2) * fres(:,3) * fres(:,22)
2661 xfac(:,1) = strnrm*strnrm*( fres(:,10)**2 + fres(:,11)**2 +
2662 & fres(:,12)**2
2663 & + two*( fres(:,13)**2 + fres(:,14)**2 + fres(:,15)**2 ) )
2665 xfac(:,2) = strnrm*( xlij(:,1)*fres(:,10) + xlij(:,2)*fres(:,11)
2666 & + xlij(:,3)*fres(:,12) +
2667 & two*(xlij(:,4)*fres(:,13) + xlij(:,5)*fres(:,14) +
2668 & xlij(:,6)*fres(:,15)) )
2670 xfac(:,3) = strnrm*( fres(:,10)*fres(:,16) + fres(:,11)*fres(:,17)
2671 & + fres(:,12)*fres(:,18) +
2672 & two*(fres(:,13)*fres(:,19) + fres(:,14)*fres(:,20) +
2673 & fres(:,15)*fres(:,21)) )
2675 xfac(:,4) = xlij(:,1)*fres(:,16) + xlij(:,2)*fres(:,17)
2676 & + xlij(:,3)*fres(:,18) +
2677 & two*(xlij(:,4)*fres(:,19) + xlij(:,5)*fres(:,20) +
2678 & xlij(:,6)*fres(:,21))
2680 xfac(:,5) = fres(:,16)*fres(:,16) + fres(:,17)*fres(:,17)
2681 & + fres(:,18)*fres(:,18) +
2682 & two*(fres(:,19)*fres(:,19) + fres(:,20)*fres(:,20) +
2683 & fres(:,21)*fres(:,21))
2694 fres(j,:) = zero
2721 fres(isgbeg:isgend,:) = zero
2777 xl1(ifath(i),1) = xl1(ifath(i),1) + fres(i,4)
2778 xl1(ifath(i),2) = xl1(ifath(i),2) + fres(i,5)
2779 xl1(ifath(i),3) = xl1(ifath(i),3) + fres(i,6)
2780 xl1(ifath(i),4) = xl1(ifath(i),4) + fres(i,7)
2781 xl1(ifath(i),5) = xl1(ifath(i),5) + fres(i,8)
2782 xl1(ifath(i),6) = xl1(ifath(i),6) + fres(i,9)
2784 xl2(ifath(i),1) = xl2(ifath(i),1) + fres(i,1)*fres(i,1)
2785 xl2(ifath(i),2) = xl2(ifath(i),2) + fres(i,2)*fres(i,2)
2786 xl2(ifath(i),3) = xl2(ifath(i),3) + fres(i,3)*fres(i,3)
2787 xl2(ifath(i),4) = xl2(ifath(i),4) + fres(i,1)*fres(i,2)
2788 xl2(ifath(i),5) = xl2(ifath(i),5) + fres(i,1)*fres(i,3)
2789 xl2(ifath(i),6) = xl2(ifath(i),6) + fres(i,2)*fres(i,3)
2791 ui(ifath(i),1) = ui(ifath(i),1) + fres(i,1)
2792 ui(ifath(i),2) = ui(ifath(i),2) + fres(i,2)
2793 ui(ifath(i),3) = ui(ifath(i),3) + fres(i,3)
3159 dimension fres(nshg,33), fwr(nshg),
3219 fres = zero
3309 & fres, hfres, shglf(lcsyst,:,1:nshl,:),
3316 if(numpe>1) call commu (fres, ilwork, 33, 'in ')
3323 fres(i,:) = fres(i,:) + fres(j,:)
3330 fres(j,:) = fres(i,:)
3335 call commu (fres, ilwork, 33, 'out')
3338 fres(:,22) = one / fres(:,22)
3340 fres(:,j) = fres(:,j) * fres(:,22)
3343 fres(:,j) = fres(:,j) * fres(:,22)
3369 & two * (fres(:,10)**2 + fres(:,11)**2 + fres(:,12)**2)
3370 & + four * ( fres(:,13)**2 + fres(:,14)**2 + fres(:,15)**2 ) )
3377 & * fres(:,10) + fres(:,16)
3379 & * fres(:,11) + fres(:,17)
3381 & * fres(:,12) + fres(:,18)
3383 xmij(:,4) = -fwr * fres(:,13) + fres(:,19)
3384 xmij(:,5) = -fwr * fres(:,14) + fres(:,20)
3385 xmij(:,6) = -fwr * fres(:,15) + fres(:,21)
3388 xlij(:,1) = fres(:,4) - fres(:,1) * fres(:,1)
3389 xlij(:,2) = fres(:,5) - fres(:,2) * fres(:,2)
3390 xlij(:,3) = fres(:,6) - fres(:,3) * fres(:,3)
3391 xlij(:,4) = fres(:,7) - fres(:,1) * fres(:,2)
3392 xlij(:,5) = fres(:,8) - fres(:,1) * fres(:,3)
3393 xlij(:,6) = fres(:,9) - fres(:,2) * fres(:,3)
3407 xfac(:,1) = strnrm*strnrm*( fres(:,10)**2 + fres(:,11)**2 +
3408 & fres(:,12)**2
3409 & + two*( fres(:,13)**2 + fres(:,14)**2 + fres(:,15)**2 ) )
3411 xfac(:,2) = strnrm*( xlij(:,1)*fres(:,10) + xlij(:,2)*fres(:,11)
3412 & + xlij(:,3)*fres(:,12) +
3413 & two*(xlij(:,4)*fres(:,13) + xlij(:,5)*fres(:,14) +
3414 & xlij(:,6)*fres(:,15)) )
3416 xfac(:,3) = strnrm*( fres(:,10)*fres(:,16) + fres(:,11)*fres(:,17)
3417 & + fres(:,12)*fres(:,18) +
3418 & two*(fres(:,13)*fres(:,19) + fres(:,14)*fres(:,20) +
3419 & fres(:,15)*fres(:,21)) )
3421 xfac(:,4) = xlij(:,1)*fres(:,16) + xlij(:,2)*fres(:,17)
3422 & + xlij(:,3)*fres(:,18) +
3423 & two*(xlij(:,4)*fres(:,19) + xlij(:,5)*fres(:,20) +
3424 & xlij(:,6)*fres(:,21))
3426 xfac(:,5) = fres(:,16)*fres(:,16) + fres(:,17)*fres(:,17)
3427 & + fres(:,18)*fres(:,18) +
3428 & two*(fres(:,19)*fres(:,19) + fres(:,20)*fres(:,20) +
3429 & fres(:,21)*fres(:,21))
3440 fres(j,:) = zero
3467 fres(isgbeg:isgend,:) = zero
3523 xl1(ifath(i),1) = xl1(ifath(i),1) + fres(i,4)
3524 xl1(ifath(i),2) = xl1(ifath(i),2) + fres(i,5)
3525 xl1(ifath(i),3) = xl1(ifath(i),3) + fres(i,6)
3526 xl1(ifath(i),4) = xl1(ifath(i),4) + fres(i,7)
3527 xl1(ifath(i),5) = xl1(ifath(i),5) + fres(i,8)
3528 xl1(ifath(i),6) = xl1(ifath(i),6) + fres(i,9)
3530 xl2(ifath(i),1) = xl2(ifath(i),1) + fres(i,1)*fres(i,1)
3531 xl2(ifath(i),2) = xl2(ifath(i),2) + fres(i,2)*fres(i,2)
3532 xl2(ifath(i),3) = xl2(ifath(i),3) + fres(i,3)*fres(i,3)
3533 xl2(ifath(i),4) = xl2(ifath(i),4) + fres(i,1)*fres(i,2)
3534 xl2(ifath(i),5) = xl2(ifath(i),5) + fres(i,1)*fres(i,3)
3535 xl2(ifath(i),6) = xl2(ifath(i),6) + fres(i,2)*fres(i,3)
3537 ui(ifath(i),1) = ui(ifath(i),1) + fres(i,1)
3538 ui(ifath(i),2) = ui(ifath(i),2) + fres(i,2)
3539 ui(ifath(i),3) = ui(ifath(i),3) + fres(i,3)